Displaying 1 result from an estimated 1 matches for "ncstar".
Did you mean:
cstar
2011 Oct 05
4
SPlus to R
...-d+rc/nc
+ ne1<-nef(rc,nc,pe*ne,ne,alpha,power)
+ ## if(is.na(ne1)) print(paste('rc=',rc,',nc=',nc,',pe=',pe,',ne=',ne))
+ }
+ if (ne1>100000) return(NA)
+ else return(ne=ne1)
+ }
+ ### for method 3
+ if (method==3) {
+ if (tol1 > tol2/10) tol1<-tol2/10
+ ncstar<-(1-d)*nc
+ pc<-(0:ncstar)/nc
+ ne<-rep(NA,ncstar + 1)
+ for (i in (0:ncstar))
+ { ne[i+1]<-ss.rand(i,nc,d,alpha=.05,power=.8,tol=.01)
+ }
+ plot(pc,ne,type='l',ylim=c(0,max(ne)*1.5))
+ ans<-c.searchd(nc, d, ne, alpha, power, cc, tol1)
+ ### check overall absolute deviance
+...