Question: Determine Condition using LogFC values
0
gravatar for ApoorvaB
4.1 years ago by
ApoorvaB230
United States
ApoorvaB230 wrote:

Hi Everyone,

I am analyzing RNA-Seq data having 22 samples from 3 batches for differential expression of genes. My condition is test between defective and normal phenotype. For one of the samples, the condition is indeterminate. From the slides, it appears to be slightly defective but might not necessarily be so. Is there anyway I can determine what to label it as? Would it help to see changes in logFC values when i first label it as defective and then as normal?

 Thanks

fold changes rna-seq • 1.0k views
ADD COMMENTlink modified 4.1 years ago by dariober10k • written 4.1 years ago by ApoorvaB230

If you're unsure what the sample is, it is best to exclude it (not only for your interests but for others). On the other hand, if you want to look at how similar the replicates are then you can use a simple Pearson correlation values/plots to make the decision.

ADD REPLYlink modified 5 weeks ago by RamRS25k • written 4.1 years ago by arnstrm1.7k
1
gravatar for dariober
4.1 years ago by
dariober10k
WCIP | Glasgow | UK
dariober10k wrote:

Principal components analysis (PCA) is sometimes applied to expression levels from RNA-Seq data to spot outliers or otherwise unexpected sample behaviours. You could apply it to your case and see first, if samples cluster neatly by condition. Then see which cluster your undetermined sample best belongs to. Having said that, I would be careful with labelling this sample as one or the other group just on the bases of PCA.

If you are using edgeR look at function plotMDS, I think DESeq has some similar function.

ADD COMMENTlink modified 5 weeks ago by RamRS25k • written 4.1 years ago by dariober10k
PCs on all genes will classify the unknown sample correctly, but I bet you $100 the anomalous sample lies somewhere between the two and removed in the third PC.
ADD REPLYlink written 4.1 years ago by karl.stamm3.6k
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: 2011 users visited in the last hour