similar to: clustering indexes/indices

Hi, I'm trying to debug a function I wrote to calculate the cindex for a hierarchical tree. For this it is useful to compare my calculations with those in output from the clustindex function, in the cclust library. There's no way, however, to have the cindex value for a given output of the cclust function, as a NA value is always returned. This happens almost surely because the cindex in

Hi- I'd like to apply some of the cluter indices, using 'clustindex' in the cclust package, but using output from hierarchical clustering using 'hclust' in the mva package. The two packages use a different format for storing the results of clustering; has anyone coerced one into the other? Thanks! -John Barnett

Hi all, I'd like to add the Bhattacharyya method to the dist function. What is the best way to do this? I'm using R 1.9.0 and I was looking for the code that defines the already existing distances, but I didn't manage. As far as I understand, dist is defined in mva that is part of the stats package now, but where is the code? Thank you very much, Giampiero

Hi, I am trying to produce incremental plots directly using the pdf device. I mean that I want to produce a plot, save it in a pdf file, then add details and save the new plot in a new file, and so on. Before I used to do this with x11, by just adding graphics to a plot and then using dev.copy2eps at the right times for each eps figure. Now I want to do this in batch mode and the x11 device is not

Hi all, once more about the old subj :-) My data has too much various distribution families and for every particular experiment I need just to decide whether the data is "quite homogeneous" or it has two or more clusters. I've revisited the following libraries: amap, clust, cclust, mclust, multiv, normix, survey. And I didn't find any ready-to-use general

Hi, is there a way to reference to a data object without copying it? For example I have a huge matrix called dist and I want two objects obj1 and obj2 to have a memeber dist that points to the matrix, but I don't want, for memory reasons, to copy the matrix twice. As far as I understand the following code will generate three copies of my data: dist <- some_code_that_generates_data

Hi, I was checking the source code to the function combn that "generates all combinations of the elements of 'x' taken 'm' at a time.", because I wished to modify it. I have a doubt about a statement. This is the main loop. ._1 <- 1:1 nmmp1 <- n - m + ._1 while (a[1] != nmmp1) { if (e < n - h) { h <- ._1 e <-

Hello, here is my problem, After looking at the mail archives, I found a description of the error I get when I use this package. At first I even tought that they were showing how to solve it. But the thing is that by saying "the programmer forgot drop=FALSE" doesn't show me how I should get rid of the problem I have looked inside the package very quickly and I found three

I have to solve a clustering problem. My first step is to determinate the number of clusters, that's why I 'm using the Calinski index ( [tr(b)/(k-1)]/[tr(w)/(k-1)] ) which i try to maximize to have the best number of clusters. A function is already implemented in R to calculate this index : clustIndex(cl,x, index="calinski") where cl is the result of a clustering method ,

> Well, that is in all good texts on R, together with the > solution: drop=FALSE. See ?"[" for the on-line details. OK. Thank you a lot. Now patched cclust and clustIndex work fine for 1D case. BTW, why not to apply the "drop=F" to these functions? I guess other users need 1D case as well. kind regards, Valery A.Khamenya

Hi all, Why subscripting a one column matrix drops one dimension? > x<- matrix(rnorm(100), ncol=1) > str(x) num [1:100, 1] -0.413 -0.845 -1.625 -1.393 0.507 ... > str(x[20:30,]) num [1:11] -0.315 -0.693 -0.771 0.448 0.204 ... > str(x[20:30]) num [1:11] -0.315 -0.693 -0.771 0.448 0.204 ... This breaks: > cov(x) [,1] [1,] 0.9600812 >

Dear R-friends, the dataset I am using (data.it) is organized as follows partner stp btp reg hk 0.64 1 s ger 0.27 1 d tur 0.27 1 s rom 0.24 1 s-f por 0.24 1 s spa 0.23 1 s gre 0.22 1 d-f aus 0.17 1 d uk 0.16 1 s be 0.16 1 d arg 0.15 1 s usa 0.13 1 d-f fra 0.13 1 s neth

Hi all, I created a data frame with three factors, plus the response that looks like this: x1 x2 x3 y a 1 1 0.3 a 2 1 0.1 b 1 1 0.4 c 4 3 0.1 ... I would like to analise the effect of two of them, keeping the third fixed (I already know the effect of the last). The reason why I don't create several data frames for each value of the thirs factor is simply convenience (I'd like to be able

Hi all, I would like to add phonetic symbols in my figures. Usually I typeset my documents in LaTeX and I use the tipa [1] package to get International Phonetic Alphabet (IPA) fonts. So, my problem would be solved if I could insert LaTeX commands in the text() function (I guess at least). I would like to avoid using psfrag (that is, inserting a string in the eps figure and then substituting it

Hi, I am trying to create a plot and save it as pdf to the file file1.pdf, then add some details and save the new plot as file2.pdf. I would like to avoid repeating all the instructions needed to create file1.pdf when I create file2.pdf. This is what I have tried: pdf(file="file1.pdf") plot(1:10) dev.copy(pdf,file="file2.pdf") points(2,3) dev.off() dev.off() When I run

I would like to see an example of a data matrix for cclust and how to import it to cclust. In fact, i don't know how to give my imput for cclust program! i test this file 1 0.23 1.52 2 0.52 1.25 3 0.13 1.89 4 0.78 1.11 i do >library(cclust) >x<-scan("test.matrice.phyl") >cclust(x,2,method="kmeans") i have this error message: Error in sample(length(x),

Dear R users, When I run the following code, R crashes: require(cclust) x <- matrix(c(0,0,0,1.5,1,-1), ncol=2, byrow=TRUE) cclust(x, centers=x[2:3,], dist="manhattan", method="kmeans") While this works: cclust(x, centers=x[2:3,], dist="euclidean", method="kmeans") I'm posting this here because I am not sure if it is a bug. I've been searching

I have checked that section already. Sorry, I should have mentioned that. Memory limit increase does not work. Installtion of msvcrt.dll does not work either. Thank you. -----Original Message----- From: ripley at stats.ox.ac.uk [mailto:ripley at stats.ox.ac.uk] Sent: Wednesday, March 05, 2003 2:44 PM To: Igor Oleinik Cc: r-help at stat.math.ethz.ch Subject: Re: [R] problem with cclust[er]

There is a new version of the 'CCLUST' package ,where i removed the extra command for the kmeans algorithm in the .R programm and also the comments about it in the .Rd help page. Now in the cclust library the kmeans algorithm can be applied only by using the cclust function. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

There is a new version of the 'CCLUST' package ,where i removed the extra command for the kmeans algorithm in the .R programm and also the comments about it in the .Rd help page. Now in the cclust library the kmeans algorithm can be applied only by using the cclust function. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read