Table2.5=function(m=15,n=15,sig2=2){
avec=c(.2,.4,.6,.8,1,1.2)
ax=-avec/sqrt(2*sig2)
p1=pnorm(ax)
if(!exists("pmnorm")) library(mnormt)
p2=NULL
for(i in 1:6){
p2[i]=pmnorm(c(ax[i],ax[i]),c(0,0),varcov=matrix(c(1,.5,.5,1),ncol=2))
}
EWXY=m*n*p1
varWXY=m*n*(p1*(1-p1)+(m+n-2)*(p2-p1^2))
Prob=2*pnorm(m*n*(.5-p1)/sqrt(varWXY))-1
alt.Prob=2*pnorm(sqrt(3*m*n/((m+n+1)*pi))*(avec/sqrt(sig2)))-1
out=round(rbind(p1,p2,EWXY,varWXY,Prob,alt.Prob),4)

out
}