Question: DESEq2 Principal Component Analysis : generate two PCA, one with condition1 samples and an other with condition2
0
gravatar for JoannaF
15 months ago by
JoannaF0
JoannaF0 wrote:

Hi,

I have 33 RNA-seq samples from two conditions which can be separated in other sub-conditions. I ran DESEq2 on all the dataset up to

rld <- rlog(dds, blind=FALSE)

Then, I would to run plotPCA into a subset of samples of each condition. Is this code right in order to do that ?

For the first condition :
plotPCA(rld[protein_coding_genes_id_vect,c(1,2,3,10,11,12,19,20,21,22,23,24,25,26,27,31,32,33)], intgroup=c("condition"))

For the second condition :

plotPCA(rld[protein_coding_genes_id_vect,c(4,5,6,7,8,9,13,14,15,16,17,18,28,29,30)], intgroup=c("condition"))

Thanks a lot for your answers !

 

deseq2 pca • 266 views
ADD COMMENTlink modified 15 months ago by Michael Love25k • written 15 months ago by JoannaF0
Answer: DESEq2 Principal Component Analysis : generate two PCA, one with condition1 samp
2
gravatar for Michael Love
15 months ago by
Michael Love25k
United States
Michael Love25k wrote:

Yes, you can subset the DESeqTransform object directly to perform PCA on subsets.

ADD COMMENTlink written 15 months ago by Michael Love25k

Thanks lot for your answer !

ADD REPLYlink written 15 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: 98 users visited in the last hour