绘制功效分析图形

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

以下程序绘制了一幅功效分析图形：
> library(pwr)
> r<-seq(.1,.5,.01)  
> nr<-length(r)     
> p<-seq(.4,.9,.1)
> np<-length(p)
> samsize<-array(numeric(nr*np),dim=c(nr,np))   
> for(i in 1:np){      
+ for(j in 1:nr){
+ result<-pwr.r.test(n=NULL,r=r[j],sig.level=.05,power=p,
+ alternative='two.sided')
+ samsize[j,i]<-ceiling(result$n)
+ }}
> xrange<-range(r)
> yrange<-range(range(samsize))
> colors<-rainbow(length(p))
> plot(xrange,yrange,type='n')
> for(i in 1:np){
+ lines(r,samsize[,i],type='l',lwd=2,col=colors)
+ }
> abline(v=0,h=seq(0,yrange[2],50),lty=2,col='grey89')
> abline(h=0,v=seq(xrange[1],xrange[2],.02),lty=2,col='grey89')
> legend('topright',title='power',as.character(p),fill=colors)

请问：其中的seq(0,yrange[2],50)中的yrange[2]中的[2]是什么意思，我没想明白，望高手指点，谢谢。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：分析图 Character Library abline Result library legend result title power