#----------------------------------------------------
#永田（2003）第11章 その他の手法
#----------------------------------------------------

sig.level <- 0.05; p0 <- 0.1; p <- 0.05; n=200


#######################################
#母不良率の検定（検出力）
#######################################
p.power <- function(sig.level,p0,p,n){
A <- sqrt((p0*(1-p0))/(p*(1-p)))
B <- (p-p0)/sqrt((p*(1-p))/n)
pnorm(q=-qnorm(p=sig.level/2,lower.tail=F)*A-B,lower.tail=T)+
pnorm(q=qnorm(p=sig.level/2,lower.tail=F)*A-B, lower.tail=F)
}

##ex 11.1
p.power(sig.level=0.05,p0=0.1,p=0.05,n=200)


#######################################
#母不良率の検定（サンプルサイズ）
#######################################
省略