R/gof.Pearson.R

In TrialSize: R Functions for Chapter 3,4,6,7,9,10,11,12,14,15 of Sample Size Calculation in Clinical Research

gof.Pearson <-
function(alpha,beta,pk,pk0,r){
b=0;noncen=0;

for (i in 1:1000){
b[i+1]<-pchisq(qchisq(alpha,df=r-1,ncp=0,lower.tail=FALSE),df=r-1,ncp=i/100)
noncen=rbind(noncen,c(i,b[i],b[i+1]))
}

delta=noncen[which(noncen[,2]>beta & noncen[,3]<beta)]/100

n<-delta*(sum((pk-pk0)^2/pk0))^-1
return(n)
}