Displaying 1 result from an estimated 1 matches for "rcnumber".
Did you mean:
cnumber
2011 Oct 05
4
SPlus to R
...,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(ne1)
+ }
> ###
> power1.f<-function(rc,nc,ne,pie,alpha=0.05){
+ #-------------------------------------
+ # rcnumber of response in historical control
+ # ncsample size in historical control
+ # ne sample size in experitment group
+ # pietrue response rate for experiment group
+ # alphasize of the test
+ #-------------------------------------
+
+ za<-qnorm(1-alpha)
+ re<-ne*pie
+ xe<-asin(sqrt((re+0....