Question: Help with edgeR
0
gravatar for Neu.S
6 months ago by
Neu.S10
Neu.S10 wrote:

Hi, i am using edgeR for differential expression study. I am providing raw read counts. In few cases i am getting logFC is negative even if the raw read counts are larger. Am I doing anything wrong. My commands are as below.

countdata<-read.delim ("raw_reads.txt",sep="\t",row.names=1)
group<-factor(c("clt","clt","clt","trt","trt","trt"))
dge = DGEList (counts=countdata,group=group)
dge <- estemateCommonDisp(dge)
dge<- estimateTagwiseDisp(dge)
et <- exacttest(dge)
etp<- toptags(et, n=100000)
etp$table$logFC = -etp$table$logFC
write.csv(etp$table, "cltvstrt.csv")
edger • 130 views
ADD COMMENTlink modified 6 months ago by Gordon Smyth39k • written 6 months ago by Neu.S10
Answer: Help with edgeR
0
gravatar for Gordon Smyth
6 months ago by
Gordon Smyth39k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth39k wrote:

You code seems fine except that:

  1. Several of the function names are misspelled or not capitalized correctly,
  2. You don't seem to have filtered low expression genes, and
  3. You are changing the sign of the logFC for no apparent reason.

As far as results are concerned, you don't show any examples of anything that you think might be wrong.

Are you perhaps simply saying that your own simple fold-change calculation doesn't agree with edgeR's more sophisticated glm approach for some genes? You can't expect to reproduce edgeR's calculation in a few lines of your code. edgeR's calculation, which weights values inversely according to their variances, and updates the variances iteratively, is generally much preferable to simple averaging.

ADD COMMENTlink modified 6 months ago • written 6 months ago by Gordon Smyth39k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 519 users visited in the last hour