Как сделать мой график более плавным в R с помощью функции построения графика?

Question

Это мой код и мой окончательный результат на графике. Но я не знаю, как сделать линию более гладкой. Есть какие-нибудь предложения? Я пробовал лесс и сплайн, но не знаю, как применить их к моему коду.

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. Затем вы можете использовать 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))