similar to: Robust PCA?

Dear R people, I need to do some work with clustering, but know next to nothing about it at present. R has (at least) three clustering packages, cluster, mclust, cclust. I was wondering if someone can direct me to some good books where I could find documentation and background on the functions in these packages. The html help in these packages lists the following as references. Can people

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,

Has anyone ported the Bagplot function by Rousseeuw, Ruts, and Tukey from S to R? The S function comprises a script and a FORTRAN function. I assume porting is relatively uncomplicated, but since I have not done any porting before I would not want to invest the effort if a port is readily available.

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

Hi. I have been looking at the PCAproj function in package pcaPP (R 2.4.1) for robust principal components, and I'm trying to interpret the results. I started with a data matrix of dimensions RxC (R is the number of rows / observations, C the number of columns / variables). PCAproj returns a list of class princomp, similar to the output of the function princomp. In a case where I can

Hello to you all, Here''s the situation: While doing a case replacement in my home storage server I accidentally removed the post-it with the disk number from my three 1TB disks before connecting them back to the corresponding SATA connector. The issue now is that I don''t know in which order they should be connected. Do any of you know how can I _safely_ bring the zpool

Hi. I am aware of three different R packages for linear programming: glpk, linprog, lpSolve. From what I can tell, if there are N variables and M constraints, all these solvers require the full NxM constraint matrix. Some linear solvers I know of (not in R) have a sparse matrix input format. Are there any linear solvers in R that have a sparse matrix input format? (including the

Hi. I recently tried the following in R 2.5.1 on Windows XP: >ev2<-c(0.8,-0.6) >ev1<-c(0.6,0.8) >ev1%*%ev2 [,1] [1,] -2.664427e-17 >sum(ev1*ev2) [1] 0 > (I got the same result with R 2.4.1 on a different Windows XP machine.) I expect this issue is very familiar and probably has been discussed in this forum before. Can someone please point me to some

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 everybody, I want to use the cluster package (Cluster Analysis Extended Rousseeuw et al.). I downloaded it from the CRAN and installed it on my linux system (fedora core 4). All seemed to be allright. But when trying to launch examples, I obtained the following message : > library(cluster) > data(votes.repub) > agn1 <- agnes(votes.repub, metric = "manhattan",

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Date: Wed, 19 Aug 1998 23:30:58 +0200 (EEST) From: Marcin Bohosiewicz <marcus@venus.wis.pk.edu.pl> To: redhat-announce-list@redhat.com Subject: screen-3.7.4 (security update) Message-ID: <Pine.LNX.3.96.980819232924.15924A-100000@venus.wis.pk.edu.pl> I updated my package after BUGTRAQ fix for tmp-races in screen package. I built

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL

Hi. I have a procedure that reads a directory, loops through a set of particular .RData files, loading each one, and feeding its object(s) into a function, as follows: cvListFiles<-list.files(fnDir); for(i in grep(paste("^",pfnStub,".*\\.RData$",sep=""),cvListFiles)){ load(paste(fnDir,cvListFiles[i],sep="/")); myFunction(rliObject); rm(rliObject); };

Dear R-helpers, I hope someone can help me with the following problem: I derived a variable from many others and produced data.frame: toktempo When I look at the variable name I get: names(toktempo) [1] "otok.V5" But I want the variable name to be Tempo so I googled around to find info about renaming variables. I found a webpage which seemed to indicate to me that there was a

Hi. Suppose I have a vector that I partition into disjoint, contiguous subvectors. For example, let v = c(1,4,2,6,7,5), partition it into three subvectors, v1 = v[1:3], v2 = v[4], v3 = v[5:6]. I want to find the maximum element of each subvector. In this example, max(v1) is 4, max(v2) is 6, max(v3) is 7. If I knew that the successive subvector maxima would never decrease, as in the example,

Hi. I have the following code in a loop. It splits a vector into subvectors of equal size. But if the size of the original vector is not an exact multiple of the desired subvector size, then the first few subvectors have one more element than the last few. I know that the cut function could be used to determine where to break up the vector, but it doesn't seem to provide control over

Hi. I am using R 2.5.1 on a Windows XP machine. Here is an example of a piece of code I was running in older versions of R on the same machine. I am looking for underscores and replacing them with periods. This result is from R 2.4.1: >gsub ( "\\_+","\.","AAA_I") [1] "AAA.I" > Here is what I get in R 2.5.1: >gsub (

>---------- Forwarded message ---------- >Date: Tue, 13 May 1997 21:18:33 +0200 >From: Wojciech Swieboda <wojtek@ajax.umcs.lublin.pl> >To: BUGTRAQ@NETSPACE.ORG > >Hello, > I''ve lately found an overflow vulnerability in Elm (Elm is setgid >mail on linux, and perhaps on some other platforms aswell). I''ve tested >this bug on versions 2.3 and

Hi. If I have a vector, v_1, and another vector of positive integers, i_1, the same length as v_1, then rep(v_1,i_1) will repeat v_i[j] exactly i_1[j] times, like so: >rep(c(1,2,3),c(3,2,1)) [1] 1 1 1 2 2 3 > I'd like to do the same sort of thing where I replace v_1 with a matrix, and the jth row of the matrix is repeated i_1 times. Obviously, I could do this with for loops, like

Hi. I am using R 2.3.1 on WIndows XP, and I am having trouble with the rank function in the presence of numerical NA data. I want the NA's all to get the same rank, but they don't. Here is an example from my session: >ct_align_rets_f2$liq[6851:6859] [1] 115396 NA 362595 NA 242986 340805 NA 692905 251533