R：繪製樣本大小對功率

Question

嗨我想繪製樣本大小對功率在R 我認爲這將需要一個for循環，但不知道如何正確實現這一點。任何幫助將不勝感激，謝謝。

for (i in c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1)){ 
    graph<-power.prop.test(p1=0.4,p2=0.24,power=c(0.1:1),sig.level=0.05)} 

    plot(graph$n,graph$power,type="b",xlab="sample size",ylab="power") 

Answer 1

R中的power.t.test和power.prop.test功能沒能很好的與量化的論據，你需要用電話的環狀構造你的建議行事。

下面我使用sapply函數進行迭代並提取樣本大小，然後用於繪圖。

powrange <- seq(0.4, .9, .1) 
n <- sapply(powrange, function(i) power.prop.test(p1=0.4,p2=0.24,power=i,sig.level=0.05)$n) 
plot(n, powrange, type="b",xlab="sample size",ylab="power") 

如果不提取n那麼你最終輸出結果的列表，你需要遍歷輸出要收到n。