如何使用 plot ZC1C425268E18385D14ZA974 在 R 中使我的 plot 更平滑？

Question

這是我的代碼和我的最終圖表結果。 但我不知道如何使線條更平滑。 有什么建議嗎？ 我嘗試了 loess 和 spline，但我不知道應該如何將它們應用到我的代碼中。

p2 <- seq(.1587, .18144, .00756)
np2 <- length(p2)
power <- seq(.4, .9, .1)
npower <- length(power)
samsize <- array(numeric(np2*npower), dim = c(np2,npower))
  for (i in 1:npower){
    for (j in 1:np2){
      result <- power.prop.test(n = NULL, p1 =.1511, p2 = p2[j], 
                                sig.level = .05, power = power[i],
                                alternative = c("two.sided","one.sided"),
                                strict = FALSE, tol = .Machine$double.eps^.25)
        samsize[j,i] <- ceiling(result$n)
      }
    }

    xrange <- range(p2)
    yrange <- round(range(samsize))
    colors <- rainbow(length(power))
    plot(xrange, yrange, type = "n",
         xlab = "Target Percentage (p2)",
         ylab = "Sample Size (n)")

    for (i in 1:npower){
      lines(p2, samsize[,i], lwd = 2, col=colors[i])
    }

    abline(h=seq(0,50000,5000), v=seq(xrange[1],xrange[2], .002),lty=20,col="gray89")
    title("Sample Size Estimation for Proportion Test")
    legend("topright", title= "Power", as.character(power),
           fill=colors)

謝謝！

Answer 1

您可以增加計算樣本量的p2的點數。

例如，如果您設置

p2 <- seq(.1587, .18144, length.out = 100)

然后你得到

Answer 2

您只需要將數組轉換為 ggplot 可以處理的 dataframe 即可。 然后你可以使用geom_smooth。

df <- data.frame(samsize)
colnames(df) <- as.character(power)

df <- tidyr::gather(df)
df$p2 <- p2

ggplot() + 
  geom_smooth(data=df, aes(x=p2, y=value, color=key)) +
  scale_color_discrete(name = 'Power') +
  labs(x = 'Target Percentage (p2)', y = 'Sample Size (n)', title="Sample Size Estimation for Proportion Test") +
  theme_bw() +
  theme(plot.title = element_text(hjust = 0.5))