Question: How to analyze DEGs other than makeContrasts method for time-point RNASeq experiment?
gravatar for anikng
9 months ago by
anikng0 wrote:

I am analyzing RNAseq datasets with a bit complex experimental setup. Specifically, I have 6 strains, each of which has 2 time point for control and drug treatment. ie, 642 biological replicates=48 samples. As I know, one way of analyzing the data is to make a time point based pairwise comparison by makeContrasts and I used (edgeR),

y<-DGEList(counts = raw_counts, group = group)

#Filterout low count genes



my.contrasts <- makeContrasts(DrugvsControl.2h.rk1=Drug.2h.rk1-Control.2h.rk1,  levels=design)

qlf<-glmQLFTest(fit,contrast = my.contrasts[,"DrugvsControl.2h.rk1"])

I would like to know if estimating the pairwise DEGs at each time points and intercepting with DEGs from other time point(as am planning to do) is a valid method to analyse the given experiment? Is there any dimension that I can apply to improve the analysis? Basically my aim includes identifying genes up or down regulated at each point across all genotypes, genes that show response at a specific time point across all genotypes and some contrasting response patterns.

edger rna-seq degs • 226 views
ADD COMMENTlink written 9 months ago by anikng0

Can anyone give suggestion on above situation? If this has obvious method, a link to the document would be helpful.

ADD REPLYlink written 9 months ago by anikng0
Please log in to add an answer.


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