Question: Why is my Differential Gene Expression Analysis returning fold change values for samples with null expressions?
0
gravatar for tovereem
9 days ago by
tovereem0
tovereem0 wrote:

Hi there,

I've been noticing some unexpected results, and it has me questioning how missing values are handled in lm (in R). Specifically, I have 200 samples, all but 5 of which have a null expression value for gene A. Of the 5 samples with expression values, 4 are in group A, 1 is in group C, and none are in group B.

I'm trying to find the log fold change between groups A and C, A and B, and B and C; specifically by running

fit <- lmFit(m, design)
contrasts_matrix <- makeContrasts(contrasts=contrast_pairs, levels=design)
fit2 <- contrasts.fit(fit, contrasts_matrix)
fit2 <- eBayes(fit2)

The weird part is that this returns non-null values for A-B fold change and B-C fold change. Can anyone help me figure out why this is? It seems to me that since there are no samples in groupB, the fold changes should be null as well.

Thanks!

rna-seq R software error • 58 views
ADD COMMENTlink written 9 days ago by tovereem0

Can you please define 'null'? - zero?; NA? Elaborating on your processing pipeline and experimental set-up (prior to limma) may help.

ADD REPLYlink modified 8 days ago • written 8 days ago by Kevin Blighe42k

Yeah sorry, by null i mean NA -- missing data (not 0). Not sure why, but certain genes just weren't sequences for some of these samples

ADD REPLYlink written 8 days ago by tovereem0
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 2.3.0
Traffic: 1531 users visited in the last hour