Question: RNASeq biplot variables and samples
gravatar for ecg1g15
4 months ago by
ecg1g150 wrote:

I would like to generate a biplot

I have been following the Love et al (2019) version of RNASeq analysis using DESeq, and have been able to plot PCAs, however, they haven't designed a function of generating billets yet

DO you know how to transfer the PCA data from here to a biplot package?

I have checked frizzpcavar from the "facto extra" package but I am not able to input the variables, and is taking the gene IDs as variables (therefore a very messy graph).

DO you have other recommendation which might be easier to use?

My data looks like this (values are random)

Gene_name Sample1 Sample2 Sample 3 Sample4
Gene1 452435 23434 2342345 234235
Gene2 2435253 235245 56767 64365
gene3 4352345 56793 8437 46567
Gene4 9834757 98745 63404 92875

Environmental data looks like this

Samples Variable1 Variable2 Variable3 Variable4
Sample1 0.2 23 0.5 9 1
Sample2 0.2 23 0.5 9 2
Sample3 0.4 24 0.8 10 2
Sample4 0.3 24 0.7 10 1


rnaseq deseq2 • 126 views
ADD COMMENTlink modified 4 months ago by Michael Love26k • written 4 months ago by ecg1g150
Answer: RNASeq biplot variables and samples
gravatar for Michael Love
4 months ago by
Michael Love26k
United States
Michael Love26k wrote:

I'd recommend you take the transformed data itself and use another package to perform the dimension reduction, if you want a particular plot.

E.g., use:

mat <- t(assay(vsd))

As your matrix for dimension reduction. I've transposed it so that samples are now rows (typical for base R functions).

ADD COMMENTlink written 4 months ago by Michael Love26k
Please log in to add an answer.


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