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.