similar to: clustering methods in R

Dear All, I am attempting to use the FANNY fuzzy clustering function in R (Kaufman & Rousseeuw, 1990), found in the "cluster" package. I have run into a variety of difficulties; the two most crucial difficulties are enumerated below. 1. Where is the 'm' parameter in FANNY? In _Finding Groups in Data: An Introduction to Cluster Analysis_ (1990) by Kaufman & Rousseeuw,

Hi, I'm trying to use a similarity matrix (triangular) as input for pam() or fanny() clustering algorithms. The problem is that this algorithms can only accept a dissimilarity matrix, normally generated by daisy(). However, daisy only accept 'data matrix or dataframe. Dissimilarities will be computed between the rows of x'. Is there any way to say to that your data are already a

Hello, I would like to use my own dissimilarity matrix in a PAM clustering with method "pam" (cluster package) instead of a dissimilarity matrix created by daisy. I read data from a file containing the dissimilarity values using "read.csv". This creates a matrix (alternatively: an array or vector) which is not accepted by "pam": A call

Hi, I'm using the library cluster to cluster a set of figures (method CLARA). Somebody that it work with clustering would know informs what I make to evaluate the clustering? Tks VM, Francisco. ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Francisco Júnior, Computer Science - UFPE-Brazil "One life has more value that the world whole" ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

I use hierclust (hierarchical clustering) in multiv package. In the documentation it is said that plclust (plotting a dendrogram) is available in S-plus. Can I find it anywhere (I have searched through http://lib.stat.cmu.edu/S and found only quotations of plclust in multiv) or is it only part of the S package (which I don't have)? Thanks --------------- Charles RAUX, Laboratoire

Dear People, Consider the following fragment of R code choose <- function(n,r) { if( is.integer(n) && is.integer(r) && n > 0 && r >= 0 ) { .C("choose",as.double(n),as.double(r),comb = double(1))$comb } else stop("n must be a positive integer and r a non-negative integer.") } This is a practice function (n choose r), which I

Hi, Is there a way to check whether a file is empty in R. I did the customary searches, but did not find anything. Please cc me on any reply. Thanks, Faheem.

Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem is trying to figure out how the (model) identifier (e.g, EII, VII, VVI, etc.) relates to the covariance matrix. The parameterization of the covariance matrix makes use of the method of decomposition in Banfield and Rraftery (1993) and Fraley and Raftery (2002) where

Dear People, After considerable effort, I managed to install gcc-3.0 on my Debian system. However, Debian keeps gcc-2.95 around as the default compiler, and gcc, g++ point to the 2.95 version. The 3.0 versions need to be explicitly invoked as gcc-3.0, g++-3.0 etc. How can I get R CMD SHLIB to invoke gcc-3.0 instead of gcc and g++-3.0 instead of g++? I would like this to be done by default but

Dear People, In a qqplot I am doing, I get lines/points that are very thick. I've tried setting the lwd variable to 0.1, but it doesn't seem to have any effect. Also, I have set the value of lty to dashed, but I still get dots. The command looks like qqplot(cdf.inv(seq(0,1,length=size),theta,pos,len),empmargdistvec(len,theta,pos,size), xlim=c(-theta,theta), ylim=c(-theta,theta),

Hello, As suggested in "De'ath, 2002. Multivariate regression trees: A new technique for modelling species-environment relationships. Ecology, 83 (4):1105-1117" (for those interested), I am trying to compare the performance of a multivariate regression tree to a cluster analysis. A simple partitioning with k clusters (as done by `pam`) seemed straightforward and appropriate

A small comment: The code of Agnes is written in Fortran. Following book give more details: Kaufman, L. and Rousseeuw, P.J. (1990). _Finding Groups in Data: An Introduction to Cluster Analysis_. Wiley, New York. The 'hclust' function is based an Algorithm contributed to STATLIB by F. Murtagh and the code is written in C. Probably the differences between agnes and hclust causes

I've been asked (privately) >>>>> "CarlosJ" == jaramilloc <jaramilloc at si.edu> >>>>> on Wed, 14 Sep 2005 09:40:22 -0400 writes: .......... CarlosJ> In Kaufman & Rousseeuw 2000 book on Cluster Analysis, it says that CarlosJ> Daisy can compute Pearson correlation between variables and then CarlosJ> transform

This code, added to drivers/genericups.h, supports the CyberPower SL series UPSes. I've tested with my 725SL - the online/on battery transition, and low battery status, and finally the shutdown signal all work as expected. *** genericups.h.orig Sun Nov 2 16:54:24 2003 --- genericups.h Mon Oct 3 13:49:00 2005 *************** *** 220,225 **** --- 220,234 ---- TIOCM_ST

Hi, Consider the following > x = c(1,2) > y = c(3,4) > d = data.frame(cbind(x,y)) > d$x [1] 1 2 > d$"x" [1] 1 2 > > foo = function(val) + { + return(d$val) + } > > bar = function() + { + return(d$"x") + } > > foo("x") NULL > bar() [1] 1 2 I'm a little surprised that R accepts both the form d$x and d$"x", but

Hi. I'm checking into robust methods for principal components analysis. There seem to be several floating around. I'm currently focusing my attention on a method of Hubert, Rousseeuw, and Vanden Branden (http://wis.kuleuven.be/stat/Papers/robpca.pdf) mainly because I'm familiar with other work by Rousseeuw and Hubert in robust methodologies. Of course, I'd like to obtain

Dear People, I am trying to use R's Math library as standalone, as documented in /src/nmath/standalone. I am using C++ in Debian testing, and the versions are as follows: ii g++-3.0 3.0.3-1 The GNU C++ compiler. ii r-mathlib 1.4.0-1 `GNU S' - Standalone R math library I have a file (rand.cc) as follows. I don't think that lattice.hh or mh.hh are very

Dear rsync people, I have just started using rsync for backups. I have had a couple of issues. Note I'm trying to use rsync as user using ssh between two machines both running Debian GNU/Linux potato (2.2r3). The local machine is currently running 2.4.6-1 and the remote 2.3.2-1.2. 1) When I run rsync with the vv option, stuff scrolls of my screen faster than I can read it. I was wondering if

Dear People, I am trying to include a barplot (see code below), in a latex document. However, the plot appears in totally the wrong position. I know this is not an R question as such, but the eps file is produced by R and I thought that perhaps other people had had similar problems and could give me pointers. Am I doing anything obviously wrong? Sincerely, Faheem Mitha. PS.

Dear People, I have been trying to perform concatenation operations on lists without success. Consider the following example. > x <- runif(2) > y <- runif(2) > z <- runif(2) > xy <- list(x,y) > xyz <- list(x,y,z) I want to have a function foo such that foo(xy,z) is the same as xyz. This appears not to be as easy as I thought, though there must be a way. For example