search for: pvclust

Hi, I am new to clustering and was wondering why pvclust using "maximum" as distance measure nearly always results in p-values above 95%. I wrote an example programme which demonstrates this effect. I uploaded a PDF showing the results Here is the code which produces the PDF file: ------------------------------------------------------------...

Hi, if you look at the first image (Image1) you see that there are 2 main clusters 7 and 8 I wanted to use pvclust to calculate a p-value whether these clusters are due to chance or statistically significant. Unfortunately pvclust does not provide a p-value for the first brunch (7 and 8). So I added a row to my matrix which is very different to the rest of the data to create an additional brunch. See image her...

Hi all, I'm new to R and I'm struggling to decipher an error message. Briefly, I am trying to use the pvclust package to do hierarchical clustering of some CGH data. The data is from the Progenetix CGH database. It is arranged as a table where each column is a single case and each row is a single chromosome band. The value in each cell is either 0, 1, 2, or -1. Corresponding to no change, gain, high level...

Hi all I am writing to you with a question regarding the pvclust package. And yes, before the usual people produce their usual contact-the-package-maintainers line, ye, I tried that but the emails one can find on the web either bounce or are not responded to. Also, yes, this error has already been reported as a bug but been shot down as not reproducible (<htt...

Is there a way to assign color to nodes as with hclust/as.dendrogram/dendrapply when using pvclust? The problem is that as.dendrogram isn't working on the pvclust objects. library(pvclust) pvc <- pvclust(matrix, nboot=1000) plot(pvc) -- View this message in context: http://www.nabble.com/Node-colors-in-pvclust-tp24405329p24405329.html Sent from the R help mailing list archive at Nabbl...

Dear All I'm using PVClust to perform hierarchical clustering, for the output plot I can control most of the graphical I need, however the plot is large and I would like to split it vertically into two panels one above the other. Is there a way to plot only part of a PVClust plot, I tried to convert it to a dendrogram with...

You should reply to the list, not just me, and even more because I cannot really help you!! My guess (but I don't know this package, and even less this function) is that pvclust() is expecting a matrix as the first argument. However, "cluster" is no data, it is a function. Why, I don't know. Take a closer look at ?pvclust, especially the "usage", "arguments" and "examples" parts. HTH, Ivan -------- Message original -------...

hello, i'd be happy if someone could provide help with the following problem: i have a dist.matrix that comes from vegdist() function of the vegan package. the used method = "horn" is not accepted as argument in hvcluster(...,dist.method="..."). is there a way to incorporate the method "horn" in hvcluster()? thanks in advance! yours, kay -- View this

...another resource I should consult that I haven't seen would someone be kind enough to point it out because it isn't obvious to me. Although my data is from biological experiments I think my problem is with R rather than the nature of the data, but I may be wrong. I am attempting to use the pvclust package to do some hierarchical clustering on some CGH data I have downloaded from the Progenetix database (http://www.progenetix.de/~pgscripts/progenetix/Aboutprogenetix.html). The data is in tab delimited format, each column is a single sample each row is a chromosome band some example dummy data...

Hi all I am trying to perform the follwing: fit<-pvclust(wq, method.hclust="ward", method.dist="euclidean") but get a strange error message that I just cant figure out. Has anyone come across this? Any help would be most appricieated Error in hclust(distance, method = method.hclust) : NA/NaN/Inf in foreign function call (arg 11) I...

I realize questions about packages should go to the package maintainer, but perhaps I have an old email address (suzuki3 at is.titech.ac.jp) Also I have both a general, and a specific, question. 1) General question: i've used pvclust before to assess significance of clusters and got reasonable results. However, on a new data set (see below) the results seem odd. I wonder if pvclust is a generally used package to assess cluster signficance, or if another package/approach is considered standard? The "approximately unbiased&q...

Hello, I am attempting to assign significance levels to a UPGMA cluster analysis as part my doctoral research. The pvclust function works well but doesn't include the similarity index I need (morisita's) as an option for computing a distance matrix. Morisita's is available in vegdist in the VEGAN library but I am having a hard time getting the vegdist function to "direct drop" into pvclust as th...

Hello there, Here's a question regarding p-values on clusters produced by hierarchical cluster analysis. A web search led me to the program pvclust to tackle this problem. But when I run the problem I get strange results. The 'AU' (approximately unbiased) p-values are very different from the 'BP' values (ordinary boot-strap) p-values. The AUs commonly are in the 80-100 range where the BPs are in the 0-10 range! One clue as t...

Dear List, I have another question to bother you about how to do clustering. My data consists of 49 columns (49 variables) and 238804 rows. I would like to do hierarchical clustering (unsupervised clustering and PCA). So far I tried pvclust (www.is.titech.ac.jp/~shimo/prog/*pvclust* /) but I always had the problem like for R like "cannot allocate the memory". I am curious about what else packages can perform the clustering analysis while memory efficient. Meanwhile, is there any way that I can extract the features of each...

...a dataset composed by 364 subjects and 13 binary variables (0,1 = absence,presence). I am testing possible association (co-presence) of my variables. To do this, I was trying with cluster analysis. My main interest is to check for the significance of the obtained clusters. First, I tried with the pvclust() function, by using method.hclust="ward" and method.dist="binary". Altoghether it works (clusters and significance obtained). However, I'm not convinced by the distance matrix. Association between variables are indeed different from results obtained in PAST by using Ward on...

Hi there, I have been trying to use the "pvclust" package but have been having some difficulties. This is the first time I have used R so I am sure the mistake I am making is a basic one. The data I have is a distance matrix and I have been using the command; fit <- pvclust(cluster, nboot=1000, method.dist="euclidean") to try an...

I would like to get an hierarchical clustering tree with bootstrap values indicated on the nodes, as in pvclust. The problem is that I have only distance matrix instead of the raw data, required for pvclust. Is there a way to get it? fit1 <- hclust(dist) # an object of class '"dist" plot(fit1) # dendogram without p values library(pvclust) fit2 <- pvclust(raw.data, method.hclust="wa...

...ust(). But the clustering on the rows failed due to the size: > y<-hclust(dist(data),method='average') Error: cannot allocate vector of size 1.9 Gb I tried to increase the memory using memory.limit(size=3000), still got the same error. I also tried agnes() from cluster package and pvclust() from pvclust package without success. My computer has 2G memory. Is there a more memory efficient clustering packages available? Thanks John > sessionInfo() R version 2.11.1 (2010-05-31) i386-pc-mingw32 locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United State...

Hi R users,   I am doing multiscale bootstrapping for clustering through pvclust package. I have large data set (observations 182 and variables 5546). When i tried to make bootstrapping, then i got message as "Bootstrap (r = 0.5)............Error: cannot allocate vector of size 117.3 Mb". I am new R user and could not understand what is the problem and also don't...

...de esa forma fijas los resultados, pero no estoy seguro si otra agrupación funcione mejor. Es decir me interesa un método de agrupación que genere la "mejor" agrupación y como los resultados de kmeans cambian, no sé cual agrupación elegir. Utilicé otros métodos de agrupación como mclust y pvclust ( http://www.statmethods.net/advstats/cluster.html) que entiendo generan la "mejor" agrupación y sobretodo no varían en los resultados. A diferencia de kmeans y pvclust, para el paquete mclust no es necesario definir de antemano el número de grupos. Sobre el control del tamaño de los grup...