[BioC] report a possible buy in edgeR
wang peter
wng.peter at gmail.com
Sun Jul 8 23:44:08 CEST 2012
dear all:
i read the source coding of edgeR_2.4.6
the function of .calcFactorWeighted has some problem
.calcFactorWeighted <- function(obs, ref, logratioTrim=.3,
sumTrim=0.05, doWeighting=TRUE, Acutoff=-1e10) {
if( all(obs==ref) )
return(1)
obs <- as.numeric(obs)
ref <- as.numeric(ref)
nO <- sum(obs) # this is should be the sum of matrix
nR <- sum(ref)
from the formula in the cited paper, it seems it should be column sum
--
shan gao
Room 231(Dr.Fei lab)
Boyce Thompson Institute
Cornell University
Tower Road, Ithaca, NY 14853-1801
Office phone: 1-607-254-1267(day)
Official email:sg839 at cornell.edu
Facebook:http://www.facebook.com/profile.php?id=100001986532253
More information about the Bioconductor
mailing list