similar to: Hclust tree to Figtree w/ branch lengths

Hi R helpers... I am having troubles because of the discrepancy between the dendrogram plotted from hclust and what is wrote in the hc2Newick file. I've got a matrix C: > hc <- hclust(dist(C)) > plot(hc) with the: > write(hc2Newick(hc),file='test.newick') both things draw completely different "trees"... I have also tried with the raw distance matrix D and

Dear R-helpers, I would like to export hierarchical clusters e.g. fit <- hclust(d, method="ward") in the newick format write(hc2Newick(fit),file='hclust.newick') Searching for a possibility I found this solution in the package ctc (http://www.bioconductor.org/packages/bioc/html/ctc.html). However, neither from this site nor via CRAN it is possible to install the package.

Hi, I have reconstructed ancestral character states on a phylogeny using MuSSE in the diversitree package and plotted the character state probabilities as pie charts on the nodes. I would, however, like to colour the character states of my extant species, i.e. the tip labels, the same colours as my pie charts, such that all species in state 1 are e.g. blue, species in state 2 red and species in

Hello, I've been trying to figure out how to access the individual elements from an object of class phylo. I am reading in 201 trees created by paup as below. > read.nexus("A_30knj_200t.txt", tree.names= NULL) -> anj30 > anj30[1] $tree1 $edge [,1] [,2] [1,] "-1" "-2" [2,] "-2" "-3" [3,] "-3" "1"

Hello, The following code does'nt work for me. The last command reports an error. I have created a consensus tree using the consensus comand from phylo but cannot manipulate the phylo object afterwards to create a dendogram , by transforming the phylo object into a hclust object and then into a dendogram ?? Thanks for any help library(ade4) library(cluster) library(stats) library(ape)

Hi, Use `check.names=FALSE` head(dd,2) #? Phylo.Tree Genesis.Tree #1????????? 1??????????? 2 #2????????? 2??????????? 3 ?dd <- data.frame("Phylo Tree"= c(1:10), "Genesis Tree"= c(2:11),check.names=FALSE) head(dd,2) #? Phylo Tree Genesis Tree #1????????? 1??????????? 2 #2????????? 2??????????? 3 write.csv(dd,"Crish.csv",row.names=FALSE) ?yy<-

Dear All, I get some problems using the 'phylo' and dissimilarity functions in R. I converted an output from 'hclust' into an order of phylo so as to be able to use the 'consensus' function on it. Each time I submit the consensus codes, my computer hangs. When I tried to see what the contents of the object converted into order phylo is, I get the message

Hi, I would be very grateful if anyone could give me a bit of help with this. I am trying to run Samba 3.0.5, as distributed with Mac OS X server 10.3.8, on an Apple Xserve. Things were working beautifully until I had problems with an IP address change and had to reinstall everything. At present, the situation seems to be that the smb protocol is operating, since I can make an smb

I am beginning to work with the 'ape' package in R, and have run into some trouble. I generated a UPGMA tree based on DNA sequence distance in Paup* and read it into R, where it became an object of class "phylo". However, I need it to be classified as a "dendrogram" for my purposes (to use it to order the layout of a heatmap). I get an error using as.hclust.phylo

Hi, I was recently splitting some massive phylo class objects with extract.clade() and noticed what appears to be a bug in how tip labels are copied from the full tree to the pruned tree. This possible bug was also mentioned here: http://www.mail-archive.com/r-sig-phylo at r-project.org/msg00537.html An example: library(ape) set.seed(5) x <- matrix(rnorm(100), ncol=10) p <-

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

I saw an example online of taking hclust dendrogram and plotting it using ggplot2 and thought I would give it a try to see what it would look like. I get an error when trying to use ggplot; Error: ggplot2 doesn't know how to deal with data of class phylo. Regular plot works fine but I can't get ggplot2 to work. see code below.... rows=100 columns=100 #create matrix

Let's say I would like to look inside the function corBrownian in library (ape). When I type in the function name I get the following, which is not nearly the detail that goes into this function. I am wondering how to begin cracking this function open (and others) so I can learn more about it and perhaps code my own corClass one day. Thanks. > corBrownian function (value = 1, phy, form

Hello all, I'm now working with hclust objects and was hoping to perform some basic editing on them like: - Joining = the merging of two hclust objects (so they will share one root) - Splicing = So to cut/extract a branch out of an hclust object - that by itself will be an hclust object. I noticed I could extract one element of an hclust object by turning it into a dendrogram,

Hi I try ? Iramuteq ? and when I?ll like ? methode reinert ? I have this message Erreur R Loading required package: Matrix There were 50 or more warnings (use warnings() to see the first 50) There were 50 or more warnings (use warnings() to see the first 50) Error in Ntip + Nnode : non-numeric argument to binary operator Calls: plot.dendropr -> plot.phylo Execution halted 1 None Loading

Dear all, The version 1.0 of ape (analysis of phylogenetics and evolution) is now on CRAN. The jump from version 0.2-1 to 1.0 is explained by the fact that the initial objectives of the project have been completed. The relevant part of the Changes file is shown below. All comments, suggestions, or bug reports are welcome. Emmanuel Paradis CHANGES IN APE VERSION 1.0 NEW FEATURES

Dear all, I'm trying to use phylo.pca function from phytools for the first time. I'm using an ultrametric tree with 167 tips and all branch lengths transformed to 1. I lunched the pPCA like this: pPCA<-phyl.pca(tree,data,method = "lambda") For some reason it takes for ever and never reaches the end of the process. I tried with method="BM" and the process runs

Hi, Using Ape, I have constructed an object of class "phylo", using the method 'nj' (lets call the object 'tree_ja'). I also have a given subset of 'tree_ja' in a vector (lets call the vector 'subspecies'). What I want to do, is construct a nj tree - plot(tree_ja) - but have the species in vector 'subspecies' shown as red at the tips of the

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the