Question: Subsetting a set of genes of interest from DESeq2 res
0
gravatar for alhafidzhamdan
15 days ago by
alhafidzhamdan0 wrote:

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!

ADD COMMENTlink modified 15 days ago by Haci120 • written 15 days ago by alhafidzhamdan0
3
gravatar for Haci
15 days ago by
Haci120
Haci120 wrote:

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,]
ADD COMMENTlink modified 15 days ago • written 15 days ago by Haci120

Thank you! That solved my problem!

ADD REPLYlink written 15 days ago by alhafidzhamdan0
Please log in to add an answer.

Help
Access

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