search for: bonhlp

..."))) v4 <- csimint(estpar=par4, df=n-6, # I'm not sure whether I found # the correct degrees of freedom covm=cov4, cmatrix=cm4, conf.level=0.95) sv4 <- summary(v4) # 2. Method: I found in Handbook of Statistics Vol 13, p.616, # same can be found in http://home.clara.net/sisa/bonhlp.htm # Bonferroni on correlated outcomes: raw.p <- sv4$p.value.raw co4 <- cor(df$y[df$v==0],df$y[df$v==1]) rho <- mean(c(1,co4,co4,1)) pai <- 1-(1-raw.p)^2^(1-rho) # The results of two methods are presented in the following lines: out <- cbind(raw.p, sv4$p.value.bon, sv4$p.value.ad...