Re: [R] Need help to plot clara results [clustering]

2013-03-25 Thread Martin Maechler 
> capricy gao 
> on Thu, 21 Mar 2013 06:07:12 -0700 writes:

> I am going to use clara for� gene expression analysis,
> so tried to play around with the examples from R document:


> http://127.0.0.1:10699/library/cluster/html/clara.html

> Everything looked fine until I tried to plot the results:

> it says:  waiting to confirm page change...

> I waited for more than 10 min and NO plot came out...

:-) :-)  

I really had a great chuckle!!

*Who* do you think is waiting when R says
  " waiting to confirm page change..."
??

Of course, R is waiting for *you* to confirm the page change...


> Should I wait longer? Anything wrong like this?

:-) ;-)  

... I can hardly stop chuckling away... sorry ...

> Thanks for any input:)

You're welcome.
Martin

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Need help to plot clara results [clustering]

2013-03-21 Thread capricy gao 

I am going to use clara for  gene expression analysis, so tried to play around 
with the examples from R document:


http://127.0.0.1:10699/library/cluster/html/clara.html

Everything looked fine until I tried to plot the results: 

it says: waiting to confirm page change...

I waited for more than 10 min and NO plot came out...

Should I wait longer? Anything wrong like this?

Thanks for any input:)

And the code and results:
=
> x <- rbind(cbind(rnorm(200,0,8), rnorm(200,0,8)),
+    cbind(rnorm(300,50,8), rnorm(300,50,8)))
> clarax <- clara(x, 2, samples=50)
> clarax
Call:    clara(x = x, k = 2, samples = 50) 
Medoids:
   [,1]   [,2]
[1,] -0.9695444  0.3985362
[2,] 49.4829294 49.9229633
Objective function:  9.786942
Clustering vector:   int [1:500] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ...
Cluster sizes:   200 300 
Best sample:
 [1]   2  20  21  89  91 104 106 107 129 151 163 188 189 193 217 222 231 235 244
[20] 265 266 288 304 305 313 323 324 339 352 371 387 406 407 412 423 436 441 443
[39] 448 459 461 462 490 494

Available components:
 [1] "sample" "medoids"    "i.med"  "clustering" "objective" 
 [6] "clusinfo"   "diss"   "call"   "silinfo"    "data"  
> clarax$clusinfo
 size max_diss  av_diss isolation
[1,]  200 24.39542 9.871896 0.3450682
[2,]  300 27.28630 9.730307 0.3859591
> all.equal(clarax[-8],
+   clara(x, 2, samples=50, pamLike = TRUE)[-8])
[1] TRUE
> plot(clarax)
Waiting to confirm page change...
[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.