#----------------------------------------------------
#永田（2003）第5章 1つの母平均の検定―母分散が未知の場合
#----------------------------------------------------

######################################
#1つの母平均の検定（検出力）
#(1)u=/ u_0　（精確値）
######################################

power.t <- function(sig.level,delta,n){
	pt(qt(p=sig.level/2,df=n-1,lower.tail=T),df=n-1,ncp=delta*sqrt(n),lower.tail=T)+
	pt(qt(p=sig.level/2,df=n-1,lower.tail=F),df=n-1,ncp=delta*sqrt(n),lower.tail=F)
	}

power.t(sig.level=.05,delta=.6,n=9)


######################################
#1つの母平均の検定（検出力）
#(2)u > u_0　（精確値）
######################################

power.t2 <- function(sig.level,delta,n){
	pt(qt(p=sig.level,df=n-1,lower.tail=F),df=n-1,ncp=delta*sqrt(n),lower.tail=F)
	}

power.t2(sig.level=.05,delta=.6,n=9)


######################################
#1つの母平均の検定（検出力）
#(3)u < u_0　（精確値）
######################################
power.t3 <- function(sig.level,delta,n){
	pt(qt(p=sig.level,df=n-1,lower.tail=T),df=n-1,ncp=delta*sqrt(n),lower.tail=T)
	}

power.t3(sig.level=.05,delta=.6,n=9)


######################################
#1つの母平均の検定（サンプルサイズ）
#(1)u =/ u_0　（精確値）
######################################
samplesize.t1 <- function(power,delta, sig.level){
	((qnorm(sig.level/2,lower.tail=F)-qnorm(power,lower.tail=F))/
	delta)^2+
	(qnorm(sig.level/2,lower.tail=F)^2)/2
}

samplesize.t1(power=.9,sig.level=.05,delta=1)


######################################
#1つの母平均の検定（サンプルサイズ）
#(2)u > u_0　（精確値）
######################################
samplesize.t2 <- function(power,delta, sig.level){
	((qnorm(sig.level,lower.tail=F)-qnorm(power,lower.tail=F))/
	delta)^2+
	(qnorm(sig.level,lower.tail=F)^2)/2
}
samplesize.t2(power=.95,sig.level=.05,delta=.5)


######################################
#1つの母平均の検定（サンプルサイズ）
#(3)u < u_0　（精確値）（2）と同一
######################################
samplesize.t3 <- function(power,delta, sig.level){
	((qnorm(sig.level,lower.tail=F)-qnorm(power,lower.tail=F))/
	delta)^2+
	(qnorm(sig.level,lower.tail=F)^2)/2
}
samplesize.t3(power=.8,sig.level=.05,delta=-.5)








##ex 5.2
delta <- seq(-1,1,by=.01)
plot(power.t(sig.level=.05,delta=delta,n=9)~delta,
	ylim=c(0,1),type="l")
lines(power.t(sig.level=.05,delta=delta,n=16)~delta)
lines(power.t(sig.level=.05,delta=delta,n=25)~delta)	

delta <- c(seq(-1,-.2,by=.2),seq(-.1,.1,by=.05),seq(.2,1,by=.2))
power.t(sig.level=.05,delta=delta,n=9)



##ex 5 4
delta <- seq(-1,1,by=.01)
plot(power.t2(sig.level=.05,delta=delta,n=9)~delta,
	ylim=c(0,1),type="l")
lines(power.t2(sig.level=.05,delta=delta,n=16)~delta)
lines(power.t2(sig.level=.05,delta=delta,n=25)~delta)	

delta <- c(seq(-1,-.2,by=.2),seq(-.1,.1,by=.05),seq(.2,1,by=.2))
power.t2(sig.level=.05,delta=delta,n=9)


##ex 5.6
delta <- seq(-1,1,by=.01)
plot(power.t3(sig.level=.05,delta=delta,n=9)~delta,
	ylim=c(0,1),type="l")
lines(power.t3(sig.level=.05,delta=delta,n=16)~delta)
lines(power.t3(sig.level=.05,delta=delta,n=25)~delta)	

delta <- c(seq(-1,-.2,by=.2),seq(-.1,.1,by=.05),seq(.2,1,by=.2))
power.t3(sig.level=.05,delta=delta,n=9)


##ex 5.7
samplesize.t1(power=.9,sig.level=.05,delta=1)
power.t(sig.level=.05,delta=1,n=13)
power.t(sig.level=.05,delta=1,n=12)



##ex 5.8
samplesize.t2(power=.95,sig.level=.05,delta=.5)
power.t2(sig.level=.05,delta=.5,n=45)
power.t2(sig.level=.05,delta=.5,n=44)


##ex 5.9
samplesize.t3(power=.8,sig.level=.05,delta=-.5)
power.t3(sig.level=.05,delta=-.5,n=26)