search for: pair_kfwer

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...