Y0<-1; N0<-49; Y1<-5; N1<-50
 YC<-47; NC<-1176; S2<-0.00065
 

  Wc.stat(Y0,N0,Y1,N1,YC,NC,S2)
       ac      bc     muC   phiC     Tc  se.Tc     Wc   pvalue
 1 2.4417 58.6537 0.03997 0.0164 2.3635 1.0158 2.3267 0.009991