# Analysis example to compute pooled power and so on 

if(!require(MBESS)){install.packages("MBESS")}
library(MBESS)
if(!require(pwr)){install.packages("pwr")}
library(pwr)

# # testset 1: file-drawer
# tVal = c(2.125, 2.801, 2.017, 2.398, 2.731)
# n1 = c(12, 23,25,14,13)
# n2 = n1
# gVal = c(0.837, 0.812, 0.562, 0.880, 1.037)

# # testset 3: valid
# tVal = c(4.286, 2.439, 2.253, 3.586, 3.791)
# n1 = c(22,21,22,30,24)
# n2 = n1
# gVal = c(1.269,0.739,0.667,0.914, 1.076)

# # testset 2: optional stopping
tVal = c(2.826,2.172,2.155,2.314,2.118)
n1 = c(10,17,12,10,10)
n2 = n1
gVal = c(1.210,0.727,0.850,0.991,0.907)

gVar = matrix(,ncol = 5)
gWeight = matrix(,ncol = 5)
pwrs = matrix(,ncol = 5)
pwrPool = matrix(,ncol = 5)

for  (ii in 1:5){
  smd = ci.smd(ncp = tVal[ii], n.1=n1[ii], n.2=n2[ii])
  dVal = smd$smd
  dVar = (n1[ii]+n2[ii])/(n1[ii]*n2[ii]) + dVal^2/2/(n1[ii]+n2[ii])
  gVar[ii] = (1-3/(4*(n1[ii]+n2[ii]-2)-1))^2*dVar
  gWeight[ii] = 1/gVar[ii]*gVal[ii]
  pp = pwr.t2n.test(n1=n1[ii], n2=n2[ii], d=dVal) 
  pwrs[ii] = pp$power
}
gPool = sum(gWeight)/sum(1/gVar) # pooled g*

prob = 1
for  (ii in 1:5){
  pp = pwr.t2n.test(n1=n1[ii], n2=n2[ii], d=gPool)
  pwrPool[ii] = pp$power
  prob = prob*pwrPool[ii]
}
prob # probability of rejecting all 5 hypothesis