similar to: Fwd: RE: pvclust function

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 and perform hierarchical clustering with

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: ------------------------------------------------------------------------------------- s <-

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,

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 result2 = as.dendrogram(result) however I get the error message

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

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

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

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

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

Hello all, I posted a question to this list last week and received no response. I am unsure if this means no-one knows the answer or if I posed the question badly. I'm going to assume I posed the question badly and try again. I am new to R so it is quite likely it's a very naive question, however if there is something blindingly obvious that I am missing or if there is another resource I

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) In addition: Warning message: NAs

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

Dear R users, I'm still a beginner and I'm wondering whether S3 or S4 methods really differ for my use. I understand more or less the distinction between the 2 classes from the documentation I've read but the big question is: _*does it make a difference in practice**?*_ Up to now, I've worked without noticing anything, but it might be important to differentiate and to know

Dear users, I think my questions are pretty simple, but I got lost in the hundreds of par() and plot() arguments and plot functions, so I don't know in which direction I should go. Here is my sample dataset: test <- structure(list(DIET = structure(c(1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L), .Label = c("G", "GG", "L",

Dear users, When I boxplot(), the lines of the whiskers are dashed. However, when I save in an svg file, the dashed lines of the whiskers are not dashed anymore. How can I have the dashed lines in the svg file? I don't have this problem with a ps file, but I cannot edit such file as easily as an svg file. That's why I'd like to stick to the svg format. Thanks in advance, Ivan df

Dear Jim, Here is again exactly what I did and with the output of Rprof (with this reduced dataset and with a simpler function, it is here much faster than in real life). Thanks you again for your help! ## CODE ## mydata1<- structure(list(species = structure(1:8, .Label = c("alsen","gogor", "loalb", "mafas", "pacyn", "patro",

Dear users, I've noticed some time ago that some of the packages I use (especially doBy) could not run (yet) on R 2.11. So I kept on using R 2.10. But I think it should be fine by now. So my question is how do I check whether a package is compatible with the last version of R? How can I check whether it works perfectly or whether there are still some bugs to be fixed? I could ask the

Dear users, I'm trying to get a dot plot but always end up with a boxplot. Can someone please tell me what I am doing wrong? df <- structure(list(FACETTE = structure(c(1L, 1L, 1L, 1L, 2L, 2L, + 2L, 2L), .Label = c("base", "tip"), class = "factor"), Sq = c(274836, + 0.74182, 0.709205, 0.984552, 279869, 255712, 26566, 301464)), .Names = c("FACETTE",

Dear users, I have looked on different sources and found different functions to prompt the user to provide input. However, I couldn't find one that does exactly what I'm looking for. select.list() and menu() are nice because a graphic window appears to prompt the user. However, the user can only choose from a predefined list of choices. readline() and scan() are more free in the

Dear users, I would like to process all the lists from all *.rda files that I have in one folder. Up to now, I can load all the *.rda files without any problem. The problem is when I want to access the list saved within each *.rda file (only one list per rda file). Here is my code: fpath <- "D:/R" listnames <- list.files(path=fpath, pattern=glob2rx("*.rda"),