R/ancJNmpcp.R

In musto101/wilcox_R: What the package does (short line)

ancJNmpcp <-
function(n1,n2,regfun=qreg,nreps=1000,alpha=.05,MC=FALSE,
SEED=TRUE,cov1=TRUE,DEEP=TRUE){
if(MC)library(parallel)
if(SEED)set.seed(2)
x=list()
n=max(c(n1,n2))
nmiss=n-min(c(n1,n2))
for(i in 1:nreps){
x[[i]]=rmul(n,p=6)
if(n1<n2){
if(nmiss>0)x[[i]][1:nmiss,1:3]=NA
}
if(n1>n2){
if(nmiss>0)x[[i]][1:nmiss,4:6]=NA
}
}
if(!MC)vals=lapply(x,ancJNmpcp.sub,regfun=regfun,cov1=cov1,DEEP=DEEP)
if(MC)vals=mclapply(x,ancJNmpcp.sub,regfun=regfun,cov1=cov1,DEEP=DEEP)
vals=as.vector(matl(vals))
pc.est=hd(vals,alpha)
list(pc.est=pc.est)
}