Thursday, January 26, 2012

edgeR tumor and matched normal

Software: edgeR bioconductor R package
Description: edgeR for paired significant test use generalized linear model and likelihood ratio test to calculate the p-value for each genes/transcript.
Error: when any zero value contain in the count matrix, it shows:
  Error in beta[k, ] <- betaj[decr, ] :
  NAs are not allowed in subscripted assignments
Solution: first filter out the all-zero rows in the matrix, then transform the matrix by matrx + 1 to create non-zero matrix.
Question: compared to the large size of the library, the different in library size effecting on 1 could be neglected?

No comments:

Post a Comment