This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Subsetting a set of genes of interest from DESeq2 res

Dear all,

This is probably a very basic R question but I have done RNA-seq on two groups of cell lines and I want to compare expressions of 50 genes of interest among these two cell lines. I have run DESeq2 and created the results file, and annotated the Ensembl ID with gene symbols. What is the best step to "intersect" my 50 genes of interest with the res file?

Help is very much appreciated!

deseq2 rnaseq gene of interest subsetting
• 2,842 views • link
updated by Haci Haci • written by alhafidzhamdan • alhafidzhamdan

1 answer

A head() of the result table that you want to subset would be helpful when answering your question, but anyway as far as I understand you have already added a gene_symbol column to your result table. If that is the case the following code would give you a table of "your genes of interest".

example_genes_of_interest <- c("EPCAM", "PTPRC", "COL1A1")
res[res$gene_symbol %in% example_genes_of_interest,]
• 0 views • link
written by Haci Haci

Log in to answer this question.

More posts like this