Question: Why different tSNE plot with Seurat and Monocle if they'r using the same R package?
0
gravatar for jgarces
11 months ago by
jgarces10
Spain
jgarces10 wrote:

Hi all,

I'm doing my first serious analysis with single-cell data and I've a (maybe naïve, sorry) question.

I've analysed the data with Seurat and Monocle R packages... when I draw the tSNE plot for my sample, they are totally different and I don't understand it because internally both are using Rtsne package with defaults parameters. Any idea of this differences? tSNEs from Seurat and Monocle

seurat tsne monocle R single-cell • 1.6k views
ADD COMMENTlink modified 11 months ago by dppb0570 • written 11 months ago by jgarces10

Each time you run tSNE you'll get different results. I think you should be able to see this if you do another plot with Seurat and Monocle, you should see that it would change (unless they keep the same seed and random numbers). Also are you sure it is the same data? The first one is painted with two colors and the second one with 3.

ADD REPLYlink written 11 months ago by Lluís R.830
1
gravatar for dppb05
11 months ago by
dppb0570
dppb0570 wrote:

Despite both Seurat and monocle using Rtsne there are a few reasons the plots you got are different:

  • Assuming you have used their respective standard pipelines on your data: they have different pipelines which will alter the data considerably, especially in the QC part.

  • They are using different seeds: in Seurat::RunTSNE you can change the seed from the default value of 1, whereas in monocle::reduceDimension the seed used is 2016. You could try changing Seurat's seed to 2016 and see if you can get something closer to that of monocle, although, due to the previous point, it may not happen.

Lastly, keep in mind that in Seurat::RunTSNE you can change the parameter tsne.method in order to use a different function other than the default Rtsne.

ADD COMMENTlink written 11 months ago by dppb0570

Thanks both, guys! Nothing simpler that set.seed()...

ADD REPLYlink written 11 months ago by jgarces10
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: 1111 users visited in the last hour