First 2 columns of Table 2:One Treatment per Cluster, Exchangeable working correlation
   K=20 clusters.

     d1= tau=0
     d4=tau=.5

Va= substitue Ui with Hi * Ui in sandwich estimator of variance where 
        H[i,,]<- diag( 1/sqrt( 1- pmin(.75,diag( omegaivm ))))



tau=0
> print.sim(out$pd1, out$cid1, out$erri.d1)
[1] "Number Converged= 991"
[1] "P-Values"
    Vm,Chi2   Vs,Chi2   Vs,F,K-p Vs,F,d.tilde.no.H Vs,F,d.hat.no.H    Va,Chi2 
 0.08980827 0.1029263 0.07870838        0.06559031      0.05247225 0.08779011
 Va,F,d.tilde.with.H Va,F,d.hat.with.H 
          0.05852674          0.04238143
[1] "Mean ci"
    Vm,Chi2    Vs,Chi2   Vs,F,K-p Vs,F,d.tilde.no.H Vs,F,d.hat.no.H 
 0.07946658 0.07855416 0.08420368        0.08730673      0.09736381
    Va,Chi2 Va,F,d.tilde.with.H Va,F,d.hat.with.H 
 0.08114371          0.09109155           0.1012976



tau=.5
> print.sim(out$pd4, out$cid4, out$erri.d4)
[1] "Number Converged= 1000"
[1] "P-Values"
 Vm,Chi2 Vs,Chi2 Vs,F,K-p Vs,F,d.tilde.no.H Vs,F,d.hat.no.H Va,Chi2 
   0.091   0.089    0.069             0.067            0.04   0.081
 Va,F,d.tilde.with.H Va,F,d.hat.with.H 
               0.064               0.039
[1] "Mean ci"
   Vm,Chi2   Vs,Chi2  Vs,F,K-p Vs,F,d.tilde.no.H Vs,F,d.hat.no.H   Va,Chi2 
 0.4226344 0.4208087 0.4510727          0.451229       0.5552596 0.4317172
 Va,F,d.tilde.with.H Va,F,d.hat.with.H 
           0.4629298           0.5695816