similar to: give PAM my own medoids

I am using the PAM algorithm in the CLUSTER library. When I allow PAM to seed the medoids using the default __build__ algorithm things work well: > pam(stats.table, metric="euclidean", stand=TRUE, k=5) But I have some clusters from a Hierarchical analysis that I would like to use as seeds for the PAM algorithm. I can't figure what the mediod argument wants. When I put in the

Hello I am using k medoids in R to generate sets of clusters for datasets through time. I can plot the individual clusters OK but what I cannot find is a way of pulling out the co-ordinates of the individual points in the cluster diagrams - none of the kmed$... info sets seems to be this. Beneath is an example of a k medoid prog using the built in US arrests dataset - this is not the data I am

Greetings, I have had good success using the clara() function to perform a simple cluster analysis on a large dataset (1 million+ records with 9 variables). Since the clara function is a wrapper to pam(), which will accept known medoid data - I am wondering if this too is possible with clara() ... The documentation does not suggest that this is possible. Essentially I am trying to

Hello R-community: Much of the time I want to use loops to look at graphs, etc. For example, I have 25 plots, for which the names are m.1$medoids, m.2$medoids, ..., m.25$medoids. I want to index the object number (1:25) as below (just to show concept). for (i in 1:25){ plot(m.i$medoids) } I've tried the following, with negative results for ...

Hello, I have a customized pairs () fonction as follows that displays correctely my data. ------------------------------------------------------------------------ panel.cor1 <- function (x, y, digits=2, prefix="") { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- cor(x, y,use="pairwise.complete.obs",

Hello all. I wish to run k-means with "manhattan" distance. Since this is not supported by the function "kmeans", I turned to the "pam" function in the "fpc" package. Yet, when I tried to have the algorithm run with different starting points, I found that pam ignores and keep on starting the algorithm from the same starting-points (medoids). For my

Hello Martin Maechler and All, A simple question (I hope): How can I compute the "sum of the dissimilarities" that appears in the pam command (from the cluster package) ? Is it the "manhattan" distance (such as the one implemented by "dist") ? I am asking since I am running clustering on a dataset. I found 7 medoids with the pam command, and from it I have the

Thanks very much to all who answered my question yesterday on loops. Just a further quick question: If I have this: for(i in 1:n) { subset <- data[which(data$id==i),] which returns the rows where I have the variable ID=1, how can I tell the same statement to return just column Y where ID=i? I have tried this for(i in 1:n) { subset <- data[which(data$id==i),data$Y] but it

Hello list! I am new in * but i want to learn about its possibilities. I want somebody to tell me if what I want to do is possible with *. I have a teleconference tool which uses SIP and now I am using Asterisk as POTS gateway. When I dial certain number from a telephone I connect with asterisk which asks me for an extension. When I dial certain extension I connect with my SIP application

Hi, I recently installed an IMAP server with dovecot and the dovecot antispam plugin with dspam backend. If I configure dspam to not include the signature in the X-DSPAM-Signature header will I be able to move mail into SPAM folder? Regards, Isabel

Hi all Perhaps someone could help me? I am running a loop from i=1 to n In my data.frame I have a list of i.d's (i=1 to n) and each i.d has 2 rows of data (representing 2 time points) e.g. id Y X1 X2 1 1 0 0 1 0 0 0 2 1 1 0 2 1 1 0 I am wanting to put into my 'for' loop, a line which can pick out the rows where id=i and call this a

Hello. The laboratory where I work has the possibility to get a Sun machine with 12 processors. My question is: can R be used to implememt parallel algorithms which take into account the several processors? Is there a way to run *threads* in R? If this is possible, would somebody share his/her experiences? Thanks. Isabel -- Isabel Cañette Iguá 4225 Centro de Matemática 11400

Hola a todos, es la primera vez que escribo a la lista, espero explicarme correctamente. Trabajo en cronobiología de mamíferos y quiero comenzar a analizar mis datos con R (haciendo análisis de varianza ya llevo unos años en el entorno). Me interesa saber si cada una de mis variables podría ajustarse a una función sinusoidal; tengo medidas cada dos horas para cada una. Querría extraer el mesor y

HI R users With clara function I get a data frame (maybe this is not the exact word, I'm new to R) with the following variables: > names(myclara) [1] "sample" "medoids" "i.med" "clustering" "objective" [6] "clusinfo" "diss" "call" "silinfo" "data" I want to

Hello everyone, I need to do k-medoids clustering for data which consists of 50,000 observations. I have computed distances between the observations separately and tried to use those with pam(). I got the "cannot allocate vector of length" error and I realize this job is too memory intensive. I am at a bit of a loss on what to do at this point. I can't use clara(), because I

Hello all, I am using the clustering functions in R in order to work with large masses of binary time series data, however the clustering functions do not seem able to fit this size of practical problem. Library 'hclust' is good (though it may be sub par for this size of problem, thus doubly poor for this application) in that I do not want to make assumptions about the number of

Hello, I am trying to do some parallel programming with R. I programmed with C and MPI before. I heard that there is a package called Rmpi and one called snow. What is the difference? I know the administrator installed snow in our system so I wonder if this mean there is Rmpi in it. I believe the Rmpi will have similar functions than MPI that is why I am specially interested in it. Thanks

Hello, I am writting a program in C++ and I need to use the function polygamma. I wonder if I can call this function that is included in the R library from my program in C++. Thanks a million, Isabel

Dear R user, Can you please help me. How do I convert part of a cluster analysis output under the heading “Clustering vector” as shown below, showing the clusters to which each respondent belongs to:      [1] 1 1 2 2 1 2 1 2 1 1 2 2 1 2 2 2 2 1 1 1 1 2 2 1 2 2 1 2 2 2 2 2 2 2 2 1 2   [38] 2 1 1 2 2 2 2 2 1 2 1 2 2 2 2 1 2 1 2 2 1 2 2 2 2 2 2 1 2 1 2 2 2 1 1 2 2   [75] 2 1 2 2 2 2 2 2 2 1 1 2

Yo! > -----Urspr?ngliche Nachricht----- > Von: Mar?a Isabel L?pez S?nchez-Huete [mailto:marisa@ugr.es] > Gesendet: Freitag, 14. M?rz 2003 10:00 > An: samba > Betreff: Re: [Samba] subdirectory of home > > > Thank you for your rapid answer, Tiago, but I suppose I don't well > explain my doubt. > > What I want is to share a subdirectory of home instead the