similar to: plotting dendrograms

Dear R user, I used from the library MASS the dataset cpus and the skript plot.tree.sequence.R, which are included in the tree package: >library(MASS) > data(cpus) > cpus.ltr <- tree(log10(perf) ~ syct + mmin + mmax + cach >+ + chmin + chmax, data=cpus) At home it works but at work I get this message: >Error in .C("BDRgrow1", as.double(X),

Dear R fans, I have got a difficult sounding problem. For fitting a linear model using continuous response and then for re-fitting the model after excluding every single variable, the following functions can be used. library(MASS) model = lm(perf ~ syct + mmin + mmax + cach + chmin + chmax, data = cpus) dropterm(model, test = "F") But I am not sure whether any similar functions is

Hi All, I am trying to use a neural network for my work, but I am not sure about my approach to select a parsimonious model. In R with nnet, the IAC has not been defined for a feed-forward neural network with a single hidden layer. Is this because it does not make sens mathematically in this case? For example, is this pseudo code sensible? Thanks in advance for your help. I am sorry if this

Hello all, I am using libraries (mva,cluster) to produce dendrograms. With 1000 examples the dendrogram gets too crowded, and i am wondering whether there is an option (which i cannot find) to set the number of leaf nodes, like in matlab, and return the plot and the assignment map examples -> leaf nodes. Any suggestion is appreciated. Thanks Edo

I'd like to be able to cut dendrograms at a height I specify and then plot the resulting subtrees. I wanted to use the dendrogram object for this purpose because there doesn't seem to be a canned way to cut a hclust object and get a list of hclust objects, but there is a function (cut) that does that for dendrograms. The problem I'm having is that when I plot a dendrogram, I

Hello fellow R users, I am currently learning to use R, so please forgive me if there is an obvious explanation for the following problem. My goal is to perform WGCNA on a dataset of 19776 genes, so I opted to follow the block-wise network construction (Section 2c) in the WGCNA R Tutorial by Peter Langfelder and Steve Horvath. I ended up with 10 blocks and am able to plot the dendrograms and

Good afternoon, I ran heatmap and hclust on the same matrix x (strictly, I ran heatmap(x), and hclust(dist(t(x))), and realized that the two dendrograms were slightly different, in that the left-right arrangement of one pair of subclusters (columns) was reversed in the two functions (but all individual columns were grouped correctly). Looking through the code for heatmap as a most definite

I am using hclust and plot to produce dendrograms. Using my input data I am able to complete an analysis and obtain a vertical plot. I want to be able to plot the dendrogram horizontally.I am using version 2.6 of R and have updated my packages recently. Using the sample script for dendrograms I can produce a horizontal plot using the instruction horiz = TRUE in plot(). When I use the same

Hi, I'm trying to build dendrograms to pass to heatmap(). The dendrograms I build plot properly, but when I pass them to heatmap() I get the error message "row dendrogram ordering gave index of wrong length" (see output log below). I looked in the code of heatmap() and saw that the error was due to a NULL return value from order.dendrogram(), which in turn got a NULL return value

I would like to be able to exert certain types of control over the plotting of dendrograms (representing hierarchical clusterings) that I think is best achieved by modifying the dendrogram object prior to plotting. I am using the "dendrogram" class and associated methods. Define the cluster number of each cluster formed as the corresponding row of the merge object. So, if you are

Hi, I have produced a dendrogram of categorical data in R using the hclust function, although the input was a dissimilarity matrix produced in SAS, as I have defined my own distances. The dendrogram is fine and I can view and use this. However, I was wondering if there is a method by which I can find out the optimal place to place groups, rather than relying on my visual analysis? I don't

Using Mathematica, I've performed a hierarchical clustering of 3,107 U. S. counties based on 1995-2000 intercounty migration, using an algorithm of my own devising (see http://arxiv.org/abs/0807.1550). I can also generate the associated 3,107 x 3,107 matrix of ultrametric distances. Unfortunately, as far as I have been able to determine, the Mathematica hierarchical clustering package only

I'm running libvirt under Fedora 28. I would like to attach a USB device to a VM, but when I select "Redirect USB Device" from the "Virtual Machine" menu in virt-manager and then select the device, I get the error: USB redirection error spice-client-error-quark: Could not redirect [device name] at 1-11: Error setting USB device node ACL: 'Not authorized' (0)

Ladies and Gentlemen, As Johan Lindberg points out, the documentation for handling dendrograms is sparse....Does anyone know who is responsible for or working on development of tree methods and objects? I've written a couple of scripts for my own use to translate between parenthetical (A(B(CD))) or binary A00 B10 C11 D11 tree formats and cluster objects in R, but as an inexperienced

Dear list, I am trying to plot a radial dendrogram using the ape package, which requires my data to be of class 'phylo'. Currently I have my dendrogram stored as an object of class 'dendrogram' which was produced from an outside bit of C code, but was made into an object of class 'igraph.eigenc' and converted to a dendrogram using 'as.dendrogram()' from the igraph

Hi, I am trying to generate a group of graphics with an iteration. Some thing like this... x=1 y=1 max=10 myfiles <- paste("foo", x:max, ".png", sep="") while (x =< max) { png(file=myfiles, pointsize = 20, width = 600, height = 600, units = "px", bg="#eaedd5") plot(x,y) dev.off() x=x+1 y=y+1

Dear List, I'd like to run a discriminate analysis on a data set, but have no idea how to go about this in R. I have attempted to locate info in the manuals, but may not be consulting the right sections or documents. Can anyone point me to appropriate documentation if such exists. Many thanks, David S. David White sdavidwhite at bigfoot.com Columbus, Ohio

Hello, I would like to use by to create a series of plots, but I'm not sure how to design the function. Here's what I've got: 1. A data frame of observations of measured value 2 measured vlaues and a pch code. V1 V2 pchCode .0045 123 1 .0034 145 2 .0045 123 1 .0046 167 3 ... So what I want to do is create a single scatterplot of the multiple conditions (specified by

Hello all, I wish to compare two dendrograms (representing Hierarchical Clusterings). My problems are several: 1) how do I manually create a dendrogram object ? That is, how can I reconstruct it as an "hclust" object that creates such a dendrogram, when all I have is the dendrogram image (but don't have the underlaying distance matrix that produced it) ? I see that there is a

Hi, I'm using the latticeExtra package's dendrogramGrob function to create a grob with both a dendrogram and a covariate bar (via the 'add' parameter). When the dendrogram has very short branches at its leaves, it appears as a line and becomes impossible to see when placed next to the covariate bar. For example: > x <- outer(-5:5, -5:5, "*") + matrix(nrow = 11,