Question: Getting MDSplot from Cufflinks data [cummeRbund]
0
gravatar for anu014
2.3 years ago by
anu014180
India
anu014180 wrote:

Hello Biostars,

I am doing MDSplot like this :

library(cummeRbund)
cuff<-readCufflinks(gtfFile="merged.gtf",genome="grch38")
cuff
CuffSet instance with:
     2 samples
     56875 genes
     282953 isoforms
     131485 TSS
     81720 CDS
     56875 promoters
     131485 splicing
     18767 relCDS

try(genes.MDS<-MDSplot(genes(cuff)))

But it's throwing this error :

Error in cmdscale(d, eig = TRUE, k = 2) : 'k' must be in {1, 2, .. n - 1} Error in try(genes.MDS) : object 'genes.MDS' not found

How can I fix this error? Please help.. Thanks in advance :)

rna-seq next-gen R • 827 views
ADD COMMENTlink modified 2.3 years ago by Kevin Blighe60k • written 2.3 years ago by anu014180

I am not a huge expert on cummeRbund, but it looks like you first extract the genes from cuff object, then you use them to create MDS-plot. Usually, the input for MDS (or PCA) is some rectangular data, e.g. genes in rows, samples in columns. For me it looks like you are using just "genes". Am I right?

ADD REPLYlink written 2.3 years ago by e.rempel820
0
gravatar for Kevin Blighe
2.3 years ago by
Kevin Blighe60k
Kevin Blighe60k wrote:

Running MDS on just 2 samples will not tell you much. In any case, you have to specify the dimensions that you want to plot:

MDSplot(genes(cuff), "PC1", "PC2")
ADD COMMENTlink written 2.3 years ago by Kevin Blighe60k
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: 1005 users visited in the last hour