Bugzilla – Bug 15364
Last modified: 2013-07-19 14:45:45 UTC
Created attachment 1459 [details]
data stored in pca.pred[,1]
The kmeans() function of the stats package stalls/hangs under the following circumstances:
tmp <- kmeans(pca.pred[,1], centers=2, nstart=10)
I have attached a copy of the data stored in pca.pred[,1]. It appears that the lack of variance in the samples might be causing the calculation to stall for whatever reason.
The same data and function appears to work OK under R.2.15...
I've just installed R v2.15.3 on a different machine (same OS, Linux 12.04 LTS), and it stalls. So the situation seems more complicated than first thought.
The machine on which this works is actually running R v2.15.2 not 2.15.3 as previously thought. I've since tried R v2.15.2 on the other machine and it still stalls.
This all seems reminiscent of an issue that was previously reported:
Thank you; your example is (strictly speaking not completely) reproducible
for me on my Linux desktop....
so it should be easy to fix.
I'm having a go at it.
(In reply to comment #4)
> Thank you; your example is (strictly speaking not completely) reproducible
> for me on my Linux desktop....
> so it should be easy to fix.
> I'm having a go at it.
It's a rounding / precision problem that I can reproduce .. only on 64-bit BTW.
I've decided to *not* add numerical fuzz .. which would potentially also change cases that have been working previously.
Rather, we now catch too many 'steps' in the QTRAns routine.
Committed to R-devel and R '3.0.1 patched"