##################################################################################
# Comments to sarisbro at uottawa dot ca - Thanks.                               #
##################################################################################

attach(zzzM8)
EBinter.no <- inter[ EB==0 ]
EBinter.yes <- inter[ EB==1 ]
EBcplx.no <- cplx[ EB==0 ]
EBcplx.yes <- cplx[ EB==1 ]
# summary stats
summary(EBinter.no)
summary(EBinter.yes)

# BS, fixed sizes for inter
n.no <- length(EBinter.no)
n.yes <- length(EBinter.yes)
diff.fixed.bs <- c()
for(b in 1:100000) {
	avg.no  <- mean(sample(EBinter.no,n.no,replace=T))
	avg.yes <- mean(sample(EBinter.yes,n.yes,replace=T))
	diff.fixed.bs <- c(diff.fixed.bs,avg.no-avg.yes)
}
quantile(diff.fixed.bs,c(.025,.975))

count <- 0
for(i in 1:100000){
	if(diff.fixed.bs[i]<0) count <- (count +1)
}
print( count/100000 )
rm(count)

quantile(diff.fixed.bs,probs=.01)

hist(diff.fixed.bs,nclass=80,xlab="# of interactions: avg.wo/PS-avg.w/PS",col=1)
export.graph("D:/data/Pandit62/intercations/BS_inter_fixed_size.wmf", Name="GSD2", ExportType = "WMF")
rm(diff.fixed.bs,n.no,n.yes,avg.no,avg.yes)

# BS, fixed sizes for cplx
n.no <- length(EBcplx.no)
n.yes <- length(EBcplx.yes)
diff.fixed.bs <- c()
for(b in 1:100000) {
	avg.no  <- mean(sample(EBcplx.no,n.no,replace=T))
	avg.yes <- mean(sample(EBcplx.yes,n.yes,replace=T))
	diff.fixed.bs <- c(diff.fixed.bs,avg.no-avg.yes)
}
quantile(diff.fixed.bs,c(.01,.025,.975))

count <- 0
for(i in 1:100000){
	if(diff.fixed.bs[i]<0) count <- (count +1)
}
print( count/100000 )
rm(count)

hist(diff.fixed.bs,nclass=80,xlab="# of complexes formed: avg.wo/PS-avg.w/PS",col=1)
export.graph("D:/data/Pandit62/intercations/BS_cplx_fixed_size.wmf", Name="GSD2", ExportType = "WMF")
rm(diff.fixed.bs,n.no,n.yes,avg.no,avg.yes)

rm(EBinter.no,EBinter.yes,EBcplx.no,EBcplx.yes)
detach()