我想做一个力量对效果大小的情节.我使用以下方法完成了功率与样本大小之一:
ptab<-cbind(NULL, NULL)
for (i in c(1:4588)){
pwrt<-power.t.test(d=.9,n=c(1:4588),sig.level=.05,type="one.sample",alternative="two.sided")
plot(pwrt$n,pwrt$power,type="b",xlab="sample size",ylab="power")
Run Code Online (Sandbox Code Playgroud)
但无法找到功率与效果大小的任何说明.有任何想法吗?
编辑:刚刚看了卡拉卡尔的回答(感谢队友).看起来很棒,但你的编码与我的尝试完全不同:
plot(pwrt$d,pwrt$power,type="b",xlab="effect size",ylab="power")
Run Code Online (Sandbox Code Playgroud)
虽然这只是产生一个直接在中间有一个点的图形.我试图将Caracal的编码转换为我的样本大小,尽管它有效,但图表却很疯狂!