The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: deseq2 analysis multiple factor and 4 different time points
gravatar for yusuf.zhc
4 months ago by
yusuf.zhc0 wrote:
I have multifactor data 
1. Participants are healthy population, having both males and females.
2. different excersice intensity for different legs. few participanrts have high intensity on left and low intensity in right, and few have high intensity of right and low in left.
3. samples were taken from both righ and left leg. 
4. samples were taken at 4 different time point
5. participants response to training, few participant show high muscle grwoth in high intensity some show high muscle growth in low intensity.

I am formulating this to normalize count table i got from htseq and rsem and trying to identify differentially expressed gene responsible for muscle growth and look if there is any biases based on leg, gender, excersice intensity. 

ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = ".", design= ~ condition)
dds <- DESeq(ddsHTSeq)
table_counts_normalized <- counts(dds, normalized=TRUE)
samples$Group <- factor(paste0(samples$participants,samples$response, samples$intensity, samples$time, samples$gender))
dds <- DESeqDataSetFromMatrix(countData = table_counts_normalized, colData=samples, design=formula(~Group))
dds <- DESeq(dds)
ADD COMMENTlink modified 4 months ago by swbarnes270 • written 4 months ago by yusuf.zhc0
Answer: deseq2 analysis multiple factor and 4 different time points
gravatar for swbarnes2
4 months ago by
swbarnes270 wrote:

You make dds from raw counts, not normalized.  I'm not even sure what happens if you try to make a dds object with non-integer counts as input.


That method of making a concatenated "group" column is useful for comparing subsets of samples to each other, while using the dispersion of the whole dataset, but I don't think that's what you want to do.  You probably do not want to compare one participant at one time with another participant at another time, but that's the kind of comparison you've set up there.  


You've also omitted the other key step in setting up DESeq; the call to results, which determines exactly what comparison you are making.





ADD COMMENTlink written 4 months ago by swbarnes270

thank you for your reply.

I understand how to give concatenated subsets but i am not able to add time course into that. how should i compare then. 

ADD REPLYlink written 4 months ago by yusuf.zhc0
Please log in to add an answer.


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