Question: Getting MDSplot from Cufflinks data [cummeRbund]
0
gravatar for anu014
16 months ago by
anu014160
India
anu014160 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 • 553 views
ADD COMMENTlink modified 16 months ago by Kevin Blighe43k • written 16 months ago by anu014160

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 16 months ago by e.rempel770
0
gravatar for Kevin Blighe
16 months ago by
Kevin Blighe43k
Republic of Ireland
Kevin Blighe43k 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 16 months ago by Kevin Blighe43k
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: 1440 users visited in the last hour