search for: pbh

...gin.table(ttab,1)/total) pb<-(margin.table(ttab,2)/total) ncate<-length(levels(as.factor(as.matrix(ratings)))) pp<-prop.table(ttab) sumagr<-sum(diag(ttab))/ns alpha<-(sumagr-sum(pa*pb))/(1-sum(pa*pb)) aickin <- function(theta,ratings) { alphanew<-theta[1] pah<-theta[2:3] pbh<-theta[4:5] ratings <- as.matrix(na.omit(ratings)) ns <- nrow(ratings) nr <- ncol(ratings) r1 <- ratings[, 1] r2 <- ratings[, 2] if (!is.factor(r1)) r1 <- factor(r1) if (!is.factor(r2)) r2 <- factor(r2) ifelse (length(levels(r1)) >= length(levels(r2)),lev <- c(level...

Dear List, I have been trying to use p.adjust() to do BH multiple test correction and have gotten some unexpected results. I thought that the equation for this was: pBH = p*n/i where p is the original p value, n is the number of tests and i is the rank of the p value. However when I try and recreate the corrected p from my most significant value it does not match up to the one computed by the method p.adjust: > setwd("C:/work/Methylation/IMA/GM/siteLists...