R/power.z.test.R

In asbio: A Collection of Statistical Tools for Biologists

power.z.test<- function(sigma=1,n=NULL,power=NULL,alpha=0.05,effect=NULL,test=c("two.tail", "one.tail"),strict=FALSE){
if(is.null(n)&test=="one.tail"){
    n<-((qnorm(1-alpha)+qnorm(power))^2*sigma^2)/effect^2}

if(is.null(n)&test=="two.tail"){
    n<-((qnorm(1-(alpha/2))+qnorm(power))^2*sigma^2)/effect^2}

if(is.null(power)&test=="one.tail"){
    power<-pnorm(qnorm(1-alpha)*(sigma/sqrt(n)), mean=effect,sd=sigma/sqrt(n),lower.tail=FALSE)}

if(is.null(power)&test=="two.tail"){
    power<-pnorm(qnorm(1-(alpha/2))*(sigma/sqrt(n)), mean=effect,sd=sigma/sqrt(n),lower.tail=FALSE)
    if(strict==TRUE)power<-power+pnorm(-1*(qnorm(1-(alpha/2))*(sigma/sqrt(n))), mean=effect,sd=sigma/sqrt(n),lower.tail=TRUE)}
    
if(is.null(effect)&test=="one.tail"){
    effect<-sqrt(((qnorm(1-alpha)+qnorm(power))^2*sigma^2)/n)}

if(is.null(effect)&test=="two.tail"){
    effect<-sqrt(((qnorm(1-(alpha/2))+qnorm(power))^2*sigma^2)/n)}

res<-list(sigma=sigma,n=n,power=power,alpha=alpha,effect=effect,test=test)
res
}