> require(NSM3)
> a=c(1,2,3,4)
> b=c(5,6,7,8)
> c=c(9,10,11,12)
> dat=c(a,b,c)
> grp=rep(1:3,c(4,4,4))
> pSDCFlig(dat,grp,method="Exact")