Question: Issue with DESeq: when using DESeqDataSet object what to put under argument design when there are 3 variables
gravatar for catglen012
4 months ago by
catglen0120 wrote:


My ultimate goal is for each locality, see how males and females are differentially expressed during each separate life stage.

Example: In Chicago (locality), how is my plant at the protonema stage, differentially expressed in males compared to females?

When I construct my object dds, I was wondering how can I make my dataset include only 1 locality but also include Sex, and life stages?

I am using the following so far:

dds <- DESeqDataSetFromMatrix(countData = countData, colData = colData, design = ~ Sex + Stage + Locality)

I have 3 replicates for each different locality, as well for each different sex, and life stage.

rna-seq deseq2 R • 255 views
ADD COMMENTlink modified 9 weeks ago • written 4 months ago by catglen0120
gravatar for catglen012
4 months ago by
catglen0120 wrote:

Would this be a solution to this:

dds$group <- factor(paste0(dds$celltype, dds$condition))
design(dds) <- ~ group
dds <- DESeq(dds)
# e.g. for condition KO cell type 2 vs cell type 1
results(dds, contrast=c("group","2KO","1KO")) 
# e.g. for cell type 1 KO vs WT
results(dds, contrast=c("group","1KO","1WT"))

except instead of 2 variables it is 3? Example: dds$group <- factor(paste0(dds$celltype, dds$condition, dds$place))

ADD COMMENTlink written 4 months ago by catglen0120
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: 752 users visited in the last hour