This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Obtain genes in PC1 and PC2 from bulk RNA-seq using DESeq2

Hi everyone,

I'm a beginner to RNA-seq, have done a few courses and trying to figure my way around somethings. I have created a PCA by the following means and would like to know what genes are in the principle components;

autoplot(pcDat, data = sampleinfo, colour = "calcified" , shape = "calcified", size = 5, x = 4, y = 3)

I tried following the previous thread How to get genes that make up principal components and got stuck at point one, wherein i cant get passed the following error:

rv <- rowVars(assay(pcDat)) Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘assay’ for signature ‘"prcomp", "missing"’

Any advice, or link to another thread for advice would be much appreciated! Many thanks, Cengiz

pc2 genes pc1
• 1,007 views • link
updated by Basti • Basti • written by cengiz • cengiz

First I would like to precise that a principal component do not "contain" specific genes. All genes matter in a principal component but their relative contribution will be different. If I understand well your goal is to find the genes that contribute the most to separate your samples along your first principal components. To understand your question, we would need to know which code you ran to obtain your pcDat object. Did you use prcomp ?

• 0 views • link
written by Basti • Basti

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • GEO data sample information error
    written by Jakpa • Jakpa
    Hi, I ran the following command in R to download GEO data: library("GEOquery") library(affy) library(oligo) library("org.Hs.eg.db") asserId = "GSE32894" gse = getGEO(asserId) pData(gse) ## sample …
  • Why I get error while performing DeSeq2 ?
    written by prithvi.mastermind • prithvi.mastermind
    While perfroming DEsEq2 analysis, I get the following error: **Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘NSBS’ …
  • help for genNullSeqs()
    written by koyee_d • koyee_d
    when I used `genNullSeqs()` genNullSeqs('test.bed',genomeVersion = 'hg19',length_match_tol=500) somthing wrong as following: Error in (function (classes, fdef, mtable) : unable to find an inherited method for …
  • I want to create count matrix from expanded vcf
    written by nikitanaik224 • nikitanaik224
    I want to create a mutation count matrix from an expanded vcf file. I tried the genCountMatrixFromVcf() from signeR package. I am getting the following …
  • Error in GDCprepare
    written by Kasthuri Kasthuri
    I am working with RNA seq data from TCGA using TCGAbiolinks. Everything has been working fine so far and suddenly I get this error: > …
  • CummeRbund (CS heatmap error message)
    written by saj98 saj98
    Hello I am using cummeRbund package to analyse my data. When I use the following code h<-csHeatmap(mySigMat,cluster='both') I get the following error Error in (function …
  • NCBI analyze dataset
    written by lvitale • lvitale
    I downloaded the following dataset: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE68928 and read the file using getGEO function. multicancer <- getGEO(filename='Data/multicancer.soft.gz') like here https://warwick.ac.uk/fac/sci/moac/people/students/peter_cock/r/geo/ A this point how I can …
  • How to edit count matrix of a DESeq object?
    written by kjkjindal • kjkjindal
    Hi, In the counts documentation (https://www.rdocumentation.org/packages/DESeq2/versions/1.12.3/topics/counts) there is a line which implies the count matrix of a deseq object can be edited: ## S4 replacement …
  • error in drawing PCA plot
    written by hassanlou110 • hassanlou110
    Hi, I can't draw PCA plot by this command. The following error massage appeared. I appreciated for any one help. mydata <- ReadAffy() eset <- …
  • ERROR in loading data through rma
    written by muhammadzubair777 • muhammadzubair777
    i am facing problem in loading data through rma in R. when i have run the command rmaData = rma(gse28403preset) it is giving following error! …