[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