similar to: estimate the number of clusters

Hello all, This might have been fixed in later versions (I am using R1.7.0), r-help archive contains messages reporting similar problems but no reports of codes fixes. I have encountered a couple of problems using the silhouette function. one occurs when the clustering contains clusters composed of 1 element (Martin Maechler posted code few months ago that fixes a similar problem that occurs

Hi list, When I was using 'silhouette' from the 'cluster' package to calculate clustering performances, R crashed. I traced the problem to the fact that my clustering labels only have 2's and 3's. when I replaced them with 1's and 2's, the problem was solved. Is the function purposely written in this way so when I have clustering labels, "2" and

Dear all, I am trying (without much success) to use the fuction 'silhouette'. Would anyone encountered that before (or would know where I am wrong ?) Please find below the R ouput. Thanks in advance, L. > s <- silhouette(ct, as.dist(metric)) Error in "[<-"(*tmp*, iC, "sil_width", value = s.i) : number of items to replace is not a multiple of replacement

After ordering the table of membership degrees , i must get the difference between the first and second coloumns , between the first and second largest membership degree of object i. This for K=2,K=3,....to K.max=6. This difference is multiplyed by the Crisp silhouette index vector (si). Too it dependending on K=2,...,K.max=6; the result divided by the sum of these differences I need a final

Dear R users, I am mailing you about the graphical output of silhouette (cluster package) From the example of silhouette in help(silhouette): > ar <- agnes(ruspini) > si3 <- silhouette(cutree(ar, k = 5), # k = 4 gave the same as pam() above + daisy(ruspini)) > plot(si3, nmax = 80, cex.names = 0.5) from which one may conclude that group 1 is composed by

I'm fairly new to the silhouette functionality in the cluster package, so apologize if I'm asking something naive. If I run the 'agnes(ruspini)' example from the silhouette section of the cluster package vignette, and assign colours to clusters, two clusters have what appear to be incorrect colours in the silhouette plot. library(cluster) data(ruspini) ar<- agnes(ruspini)

Dear R users, Just wondering how to extract clustered data in a file. More precisely after running any clustering algorithms, you get information about the clusterd vectors. I want to save that data after clustering is done in the same order as the clusters are produced. How can I do that After runining algorithms you get information which vectors are clustered. for example, after running pam you

I've written the function below to simulate the mean 1st through nth nearest neighbor distances for a random spatial pattern using the functions nndist() and runifpoint() from spatsat. It works, but runs relatively slowly - would appreciate suggestions on how to speed up this function. Thanks. --Dale library(spatstat) sim.nth.mdist <- function(nth,nsim) { D <- matrix(ncol=nth,

I don't know if there is some fault in compiling or a bug of the new R-1.7.0 version: cl.pam.2 <- pam(as.dist(1-cor(mel.data)),2) plot(cl.pam.2) perform a right partitioning and silhouette plot on the old R-1.6.2 instead "Error in clusplot.default(x$diss,...... ; x is not numeric" is the output on the new R-1.7.0. Same platform: RH8.0 i686. Some suggestions? A.S.

I would like to take the values of observations and map them to a new index. I am not sure how to accomplish this. The result would look like so: x[1,2,3,4,5,6,7,8,9,10] becomes y[2,4,6,8,10,12,14,16,18,20] The "newindex" would not necessarily be this sequence, but a sequence I have stored in a vector, so it could be all kinds of values. here is what happens: > x <- rnorm(10)

Hello! # I have a data frame: mydf<-data.frame(c1=rep(NA,5),c2=rep(NA,5),c3=rep(NA,5)) # I have an index whose length is always the same as nrow(mydf): myindex<-c(1,2,3,2,1) # I need c1 to have 1s in rows 1 and 5 (based on the information in myindex) # I need c2 to have 1s in rows 2 and 4 (also based on myindex) # I need c3 to have 1 in row 3 # In other words, I am trying to achieve this

Hello, after splitting a data-frame I want to access the results. Maybe the problem is, that the factor/index is a string... ...or do I miss knowing details of the index-uasge? Please look and help: ======================================= > weblog <- read_weblog("web.log") > > > str(weblog) 'data.frame': 2247 obs. of 18 variables: $ host : Factor w/ 77

On Tuesday 28 June 2005 15:30, r-help-request at stat.math.ethz.ch wrote: Re : 37. Re: A. Mani : colours in Silhouette (Mulholland, Tom) > > Message: 37 > Date: Tue, 28 Jun 2005 09:08:24 +0800 > From: "Mulholland, Tom" <Tom.Mulholland at dpi.wa.gov.au> > Subject: Re: [R] A. Mani : colours in Silhouette > To: <a_mani_sc_gs at vsnl.net>, <r-help at

Hi there, I am working on clustering genomic data. I am using HOPACH method for clustering. I have a query here. How can i get SS (split silhouette) value for individual clusters. What I am getting as MSS value, is Mean of all SS value of individual cluster. Hence I am getting single value for a gene set. But I would be very helpful for me if I get silhouette value for individual clusters.

Hi, I am using fanny and I have estrange results. I am wondering if someone out there can help me understand why this happens. First of all in most of my tries, it gives me a result in which each object has equal membership in all clusters. I have read that that means "the clustering is entirely fuzzy". Looking at the graphics it is really difficult to understand how objects with so

Hi dear all It's the first time for me to install a developmental version of R, I came across following errors, my system is DISTRIB_ID=Ubuntu DISTRIB_RELEASE=10.04 DISTRIB_CODENAME=lucid DISTRIB_DESCRIPTION="Ubuntu 10.04.1 LTS" I downloaded the dev version of R in cran R source, and downloaded the recommended packages by using wget as described in R manual, and run

Sorry for the long message. I'm doing my best to try to explain myself. I have fitted a spline to my data, I have fitted a spline, filled in the missing data by replicating the spline coefficients associated to the last node. I obtained a number of dendograms by different combination of distance and link-method by calling DIST and AGNES. The agglomerative coefficient is very high (~ 0.99) for

Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by calculating the distance "as the crow flies", but will probably eventually need to calculate the

Hi, I am using code below to get a plot that will show me on the X axis the number of clusters and on the Y axis the cluster average widths. However I am getting this error: Error in summary(silhouette(cutree(d, x), dist(iris[, -5])))$si.summary : $ operator is invalid for atomic vectors the code I am using is: avgs<-sapply(1:20,function(x) summary(silhouette(cutree(d,x),

Hallo Sebastian, >>>>> "SP" == Sebastian P?lsterl <sebp at k-d-w.org> >>>>> on Sun, 14 Jun 2009 14:04:52 +0200 writes: SP> Hello Martin, SP> I plotting the silhouette of a clustering and storing it as png. When I SP> try to store the image as png the bars are missing. The bars are plotted SP> when I use x11 or