Question: DESEq2 Principal Component Analysis with different genes biotypes (protein coding, miRNAs, lncRNAs, etc)
0
gravatar for JoannaF
17 months ago by
JoannaF0
JoannaF0 wrote:

Hi,

I would to generate different Principal Component Analysis using DESeq2, after rlog transformation. I would to generate one PCA for protein coding genes, one PCA for miRNAs, one PCA for lncRNAs, etc (for all gene_biotypes in my reference annotations). I have count files from 33 RNA-seq samples (11 conditions). Here is the code I used to generate a PCA from all my data (all gene_biotypes) :

directory <- "/mydirectory"

directory

sampleFiles <- grep("htseqcount",list.files(directory),value=TRUE)

sampleCondition=c("condition1","condition1","condition1","condition2","condition2","condition2","condition3","condition3","condition3",[...]"condition11","condition11","condition11")

sampleTable <- data.frame(sampleName = sampleFiles, fileName = sampleFiles, condition = sampleCondition)

library("DESeq2")

ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = directory, design= ~ condition)

dds <- DESeq(ddsHTSeq)

rld <- rlog(dds, blind=FALSE)

plotPCA(rld, intgroup=c("condition"))

How can I adapt this code to generate all PCAs I want to ? My issue is that I can extract all genes from one biotype (protein coding genes for example) and after run the code above and then do the same thing for others biotypes. But I think I must normalize data once with all biotypes and then generate PCAs. Have you some advices to generate theses PCA ?

Thanks a lot !

ADD COMMENTlink modified 17 months ago • written 17 months ago by JoannaF0
Answer: DESEq2 Principal Component Analysis with different genes biotypes (protein codin
1
gravatar for Michael Love
17 months ago by
Michael Love25k
United States
Michael Love25k wrote:

You should add an indexing step before calling plotPCA:

idx <- ... # identify the relevant gene sub-types

plotPCA(rld[idx,], ...)
ADD COMMENTlink written 17 months ago by Michael Love25k
Answer: DESEq2 Principal Component Analysis with different genes biotypes (protein codin
0
gravatar for JoannaF
17 months ago by
JoannaF0
JoannaF0 wrote:

Thanks a lot for your answer ! It helps me a lot !

ADD COMMENTlink written 17 months ago by JoannaF0
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 16.09
Traffic: 282 users visited in the last hour