search for: lhood

...l1   the allowance between the expceted power and target power + #--------------------------- + pc<-(0:((1-d)*nc))/nc + ncl<-length(pc) + ne.old<-ne + ne.old1<-ne.old + ### sweeping forward + for(i in 1:ncl){ + cmin<-cc[1] + cmax<-cc[2] + ### fixed cci<-cmax bug + cci <-1 + lhood<-dbinom((i:ncl)-1,nc,pc[i]) + ne[i:ncl]<-(1+(cci-1)*(lhood/lhood[1])) * ne.old1[i:ncl] + Ep0 <-Epower(nc, d, ne, pc, alpha) + while(abs(Ep0[i]-power)>tol1){ + if(Ep0[i]<power) cmin<-cci + else cmax<-cci + cci<-(cmax+cmin)/2 + ne[i:ncl]<-(1+(cci-1)*(lhood/lhood[1])) * ne.o...