vrsn <- "20120804"
Dirname <- "ExpMeanJ2"
FileName <- "RCTvalidateExpMeanJ2"
EPSILON <- 1e-10

RCTvalidateOutputDirname <- paste(RCTdirnameV,Dirname,sep="")
if(!file.exists(RCTvalidateOutputDirname))
	dir.create(RCTvalidateOutputDirname)
RCTvalidateOutputDirname <- paste(RCTvalidateOutputDirname,"/",sep="")
cat("#####", FileName, "-", vrsn, "\n")

####################
TestName <- "Expected Mean J=2 "
####################

rslt <- NULL
test <- 0

# Tests 1 - 3
z <- seqDesign(arms=1,J=2,sample.size=1:2,variance=1,P=1,EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 4 - 6
z <- seqDesign(arms=1,J=2,sample.size=c(1,3),variance=1,P=1,EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 7 - 9
z <- seqDesign(arms=1,J=2,sample.size=1:2,variance=2,P=1,EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 10 - 12
z <- seqDesign(arms=1,J=2,sample.size=1:2,variance=1,P=0.5,EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 13 - 15
z <- seqDesign(arms=1,J=2,sample.size=1:2,variance=1,P=c(0.5,1),EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 16 - 18
z <- seqDesign(arms=1,J=2,sample.size=1:2,variance=1,P=c(Inf,1),EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))

# Tests 19 - 21
z <- seqDesign(arms=1,J=2,sample.size=c(15,25),variance=3,P=c(.7,.8),EPSILON=EPSILON)
test <- test + 1
fn <- paste(RCTvalidateOutputDirname,"test",substring(format(10000+test),2),sep="")
dsnfn <- paste(fn,"dsn",sep="")
ocfn <- paste(fn,"oc",sep="")
if(!file.exists(dsnfn))
	dput(z,file=dsnfn)
check <- all.equal(z,dget(file=dsnfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
zOC <- seqOperatingChar (z,EPSILON=EPSILON)
if(!file.exists(ocfn))
	dput(zOC,file=ocfn)
check <- all.equal(zOC,dget(file=ocfn))
rslt <- c(rslt,ifelse1(is.logical(check),check,F))
rslt <- c(rslt,seqCheckExpMeanJ2(z,EPSILON))


if (all(rslt)) {cat("     ",TestName,": All",length(rslt),"tests  PASS\n")
} else cat("     ",TestName,": failing test numbers ", (1:length(rslt))[!rslt],"\n")


rm(vrsn,EPSILON,Dirname,FileName,TestName,RCTvalidateOutputDirname,rslt,check,z,zOC,test,fn,dsnfn,ocfn)