differential expression RNASeq DESeq2

I am using DESeq2 for differential expression analysis. I want to account for treatment and Sex. I have the following code. how and where in the code should I add sex?

dds <- DESeqDataSetFromMatrix(countData = rawdata,
                              colData = metadata,
                              design = ~ treatment)

dds <- DESeq(dds)
res <- results(dds)

rna-seq

• 1,762 views • link

written 4.5 years ago by Rob •

It is generally expected that you read the tools manual first. Even in the "quick start" section there is an example on adding covariates to the design: https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html

• 0 views • link

written 4.5 years ago by ATpoint

1 answer

Simply add it in the design formula:

dds <- DESeqDataSetFromMatrix(countData = rawdata,
                              colData = metadata,
                              design = ~ treatment + Sex )

• 0 views • link

written 4.5 years ago by nhaus •

Thank you nhaus Can also this be done for more than 2 covariates? such as: treatment, sex, stage of cancer?

• 0 views • link

written 4.5 years ago by Rob •

Yes, you can add it the same way.

• 0 views • link

written 4.5 years ago by nhaus •

Thanks for responding nhaus

• 0 views • link

written 4.5 years ago by Rob •

More posts like this

DESEQ2 Code
written 2.4 years ago by ananta.kapoor •

# Load the DESeq2 library library(DESeq2) # Read your count data counts <- read.csv("C:/Users/anant/OneDrive/Documents/RNASeq/count_data.csv", row.names = 1) # Read the metadata metadata <- read.csv("C:/Users/anant/OneDrive/Documents/RNASeq/Metadata.csv") # …
DESeq2 results function not working
written 5.1 years ago by mrk •

I'm trying to use DESeq2 to do differential expression analysis. I have a simple design where there are two groups. Running the function seems to …
Generating a heatmap in DEG analysis
written 5.4 years ago by Pranathi •

After performing DESEQ2 on my data, I could able to plot MA, PCA, and EnhancedVolcano plots including a (.xlsx) file consisting of log fold change …
Differential expression analysis with DESeq2
written 5.9 years ago by Rob •

Hi all I am using HT-Seq count data to analyze differential expression with DESeq2 in R. I have about 60000 genes(all coding+non-coding genes) and 159 …
DESeq2 for Differential expression Ensembl gene ID
written 5.9 years ago by Rob •

Hi friends I use DESeq2 for differential expression analysis for HTSeq read count data with about 20000 genes. my genes have ENSEMBL ID. I have …
Understanding DESeq2 padj and Volcano Plot
written 6.9 years ago by juan.crescente •

I'm using DESeq2 for a differential expression analysis of sRNA data. The R code looks like this: library("DESeq2") alpha <- 0.05 setwd(dirname(rstudioapi::getActiveDocumentContext()$path)); getwd() data_path <- …
Differential expression analysis using DESeq2
written 7.3 years ago by prameelaawale •

Hello everyone: I am new in using DEseq2. I am trying to analyze differential expression data using a publicly available data for my final project. …
lfcShrink function in DEseq2 is giving me an error
written 7.4 years ago by Fawzi Yassine •

Hi, lfcShrink function in DEseq2 is giving me the following error: Error in lfcShrink(dds = dds, res = res_deseq2) : res should be a DESeqResults …
DESeq2 log2 fold change
written 8.0 years ago by bioinfonerd •

Hi everyone, I did the DESeq2 on my Nanostring gene expression data using the below code: dds <- DESeqDataSetFromMatrix(countData = counts_set,colData = annotation,design = ~ …
heatmap using DESeq2
written 9.2 years ago by sumithrasank75

I am trying to generate a heat map of the top differentially expressed genes between two samples. I have used DESeq2 to identify these. My …