R/Qmcp.R

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

Qmcp <-
function(x,q=.5,con=0,SEED=TRUE,nboot=NA,alpha=.05,HOCH=FALSE){
#
#  Multiple comparisons among independent groups
#  based on the quantile indicated by the argument
#  q
#
#   Familywise error is controlled with the Bonferroni method.
#   For n small, this is better than Hochberg or Rom methods.
#   For n sufficiently large, presumably Hochberg or Rom
#   gives better results, but uncertain how large n must be.
#
# The Harrell--Davis estimator is used in order to deal with tied values
#
res=pbmcp(x,est=hd,q=q,nboot=nboot,SEED=SEED,con=con)
if(!HOCH){
res$output[,4]=alpha/nrow(res$output)
res=res[1:2]
res$sum.sig=sum(res$output[,3]<=res$output[,4])   
}
res
}