Displaying 1 result from an estimated 1 matches for "pair_kfwer".
2010 Jun 23
2
question about a program
...vnorm((k*(k-1))/(m*(m-1))*alpha, tail="upper",
sigma=var)$quantile} else
{cc_z[i] <- qmvnorm((k*(k-1))/((m-i+k)*(m-i+k-1))*alpha,
tail="upper", sigma=var)$quantile}
}
cc <- 1-pnorm(cc_z)
return(cc)
}
## second step
pair_kFWER=function(m,crit_cons,pvaluessort){
k=0
while((k <m)&&(crit_cons[m-k] < pvaluessort[m-k])){
k=k+1
}
return(m-k)
}
###########################################
##############Program with one function ##############
pair_kFWER=function(m,alpha,rho,k,pvaluessort){
library(mvtnorm)
cc_z &l...