Displaying 20 results from an estimated 22 matches for "cophenetic".
2016 Apr 21
2
"cophenetic" function for objects of class "dendrogram"
Hello,
I have been using the "cophenetic" function for objects of class "dendrogram" and I have realised that it gives different results when it is used with objects of class "hclust". For instance, running the first example in the help file of the "cophenetic" function,
d1 <- dist(USArrests)
hc &l...
2016 Apr 21
1
"cophenetic" function for objects of class "dendrogram"
Note that cophenetic.default (which works on the output of hclust(dist(X)))
uses the
row names of X as labels. as.dendrogram.hclust does not retain those row
names
so cophenetic.dendrogram cannot use them (so it orders them based on the
topology of the dendrogram).
Bill Dunlap
TIBCO Software
wdunlap tibco.com
On Thu...
2016 Apr 21
0
"cophenetic" function for objects of class "dendrogram"
...o, drop = FALSE]
}
identical(canonicalize.dist(d2), canonicalize.dist(d3))
[1] TRUE
Bill Dunlap
TIBCO Software
wdunlap tibco.com
On Thu, Apr 21, 2016 at 2:37 AM, Alberto Fern?ndez Sabater <
alberto.fernandez at urv.cat> wrote:
> Hello,
>
>
> I have been using the "cophenetic" function for objects of class
> "dendrogram" and I have realised that it gives different results when it is
> used with objects of class "hclust". For instance, running the first
> example in the help file of the "cophenetic" function,
>
>
> d1...
2007 Feb 17
1
help with cluster stopping rules
Is there a function available in R that implements Mojena's Upper Tail Rule or that draws a Mojena plot?
I would also like to find a function that implements Duda and Hart's stopping rule.
Finally with function cophenet how can I achieve a straightforward Cophenet correlation coefficient (i.e. one number for example 0.876) that gives the agreement between the structure of the
2006 May 09
1
problem accessing trees after read.nexus from ape package
...0.00913 0.01506
$tip.label
[1] "H" "G" "F" "E" "D" "C" "B"
attr(,"class")
[1] "phylo"
The only problem is i can't access or save the individual trees after
importing. What i want to do is calculate the cophenetic distance for
each tree and take the correlation between pairs of trees. I keep
getting errors like this:
x<- anj30[1]
> class(x)
[1] "list"
> as.phylo(x)
Error in as.phylo(x) : no applicable method for "as.phylo"
> cophenetic.phylo(x)
Error in cophenetic.phy...
2001 Jun 12
1
cophenetic matrix
Hello,
I analyse some free-sorting data so I use hierarchical
clustering.
I want to compare my proximity matrix with the tree
representation to evalute the fitting. (stress, cophenetic correlation
(pearson's correlation)...)
"The cophenetic similarity of two objects a and b is defined as the
similarity level at wich objects a and b become members of the same
cluster during the course of clustering" Legendre, P and Legendre, L Numeri
Ecology, 1998.
I've made a...
2010 Feb 21
1
How to: Compare Two dendrograms (Hierarchical Clusterings) ?
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
2009 Nov 10
2
All possible combinations of functions within a function
Dear All,
I wrote a function for cluster analysis to compute cophenetic correlations
between dissimilarity matrices (using the VEGAN library) and cluster
analyses of every possible clustering algorithm (SEE ATTACHED)
http://old.nabble.com/file/p26288610/cor.coef.R cor.coef.R . As it is now,
it is extremely long, and for the future I was hoping to find a more
efficient...
2009 Jul 30
1
stepwise variable selection method wanted
Hi List,
I am looking for a variable selection procedure with a forward-backward selection method.
Firstly, it is meant to work with the cophenetic
correlation coefficient (CPCC) and intended to find the variable combination with the
highest cophenetic correlation. Secondly, it is aimed at Gower metric with
wards method (though this could be easily extended) aimed at categorical data.
What I have so far is a function for backward selection th...
2007 May 22
1
Goodness of fit for hclust?
I'd like to get a measure of goodness of fit for a heirarchical
clustering result from hclust. Something that would indicate the
extent to which the dendrogram accurately represents the original
dissimilarity matrix. Is there an easy way to do this?
Or, does anyone have code for computing distances between nodes given
an hclust structure? So far, my searches have come up dry.
-- David
2011 Sep 13
2
help with hclust
Hello,
how can I get the similarity value (i.e., the inner cluster similarity) that was used to cut a hierarchical tree at a specific height?
I would appreciate your help!
Best regards,
Madeleine
2013 Feb 07
1
Feature selection for kmeans
I know that within sum of squares, DB, sillhouette and cophenetic are
indicators of clustering quality, but what indicators I need to observe when
I choose attributes for kmeans?
--
View this message in context: http://r.789695.n4.nabble.com/Feature-selection-for-kmeans-tp4657830.html
Sent from the R help mailing list archive at Nabble.com.
2010 Jun 02
1
compute the associate vector of distances between leaves in a binary non-rooted tree
Hello.
I'd like to compute the associate vector of distances between leaves in a binary non-rooted tree. The definition of a distance between two leaves in a binary non-rooted tree is the number of edges in the path joining the two leaves.
I've tried the ape package but I'm unable to find this vector.
For example, using rtree(5,rooted=F) I've obtained the following tree:
$edge
2003 Aug 04
1
hclust() and agnes() method="average" divergence (PR#3648)
...library(cluster)
library(mva)
daisemol<-daisy(emol,type=list(asymm=c(1:21)))
layout(matrix(c(1,1,2,2),2,2,byrow=TRUE))
hclustmol<-hclust(daisemol,method="average")
plot(hclustmol,cex=0.6)
agnesmol<-agnes(daisemol,method="average")
pltree(agnesmol,cex=0.6)
coagnes<-cophenetic(agnesmol)
cor(coagnes,daisemol)
cohclust<-cophenetic(hclustmol)
cor(cohclust,daisemol)
I'm using R 1.7.1 (2003-06-16) on Windows XP. Library "cluster" version
1.7.4 and library "mva" version 1.7.1.
Best wishes,
Mikkel
Mikkel Grum
International Plant Genetic Resource...
2006 Aug 05
0
cor of two matrices whose columns got shuffled
...322347 1.289899 1.553109 1.872578 1.452274 1.456126
1.335005 1.184658 1.018210
The problem is that these matrices are the endpoint of an extensive
simulation involving phylogenetic trees and the order of the columns
and rows shifts around according to how they were ordered in the
tree. I used cophenetic.phylo to get the distance matrix and now I
need to find a way to reorder or sort the columns and rows such that
the two matrices match and cor is taking the correlation of the
corresponding column and row and not just the one in the same
position. I don't want to sort the values in the col...
2008 Mar 09
0
Matrix Row name Compatibility
...ng the partial mantel
method. However, i am having trouble creating compatible matrices. I can
create matrices from my numerical data using the 'distance(x, method)'
function from ecodist, but the row/column names are out put as numbers.
While when i create a phyo-distance matrix using 'cophenetic(x)' from APE,
the row/column names are the taxa names. Different data sets incorporate
different taxa from the tree, which requires me to prune it for each data
set, meaning I cannot just convert the taxa names to corresponding numbers
in the data matrix (because the numbers will change for eac...
2009 Jan 07
1
troubles performing Moran.I test
dear R users,
I have troubles performing Moran.I test as suggested on
http://www.ats.ucla.edu/stat/r/faq/morans_i.htm
my spatial data are longitude and lattitide of communities. The
calculation of the inverse distance matrix according to the homepage
(using my data)
datAL <- read.csv2("C:\\Konvergenz AL.csv", header=T)
ALdist <- as.matrix(dist(cbind(datAL$L?nge,
2008 Nov 28
1
Unable to plot or analyse phylogeny (PR#13345)
...2)
*** error: can't allocate region
*** set a breakpoint in malloc_error_break to debug
The same problem comes up when I try to run an autocorrelative model with some
dummy continuous data. (I've tried many different functions, all produce similar
error responses):
> Moran.I(dummy$dist, cophenetic(sa))
Error: cannot allocate vector of size 1.4 Gb
R(539,0xa0810fa0) malloc: *** mmap(size=1488809984) failed (error code=12)
*** error: can't allocate region
*** set a breakpoint in malloc_error_break to debug
R(539,0xa0810fa0) malloc: *** mmap(size=1488809984) failed (error code=12)
*** error:...
2002 Oct 01
6
R-1.6.0 is released
...`priority'.
o termplot() uses factor levels rather than 1,2,3... for x-axis.
o Workaround for optimization bugs on gcc 3.1/2 on 32-bit Solaris.
o The trace() function has been robustified and a new function
tracingState() added to turn tracing temporarily on and off.
o New cophenetic() in "mva" as utility for hierarchical clustering.
o p.adjust() has two new methods, 'Hommel' and 'FDR', contributed
by Gordon Smyth <smyth at wehi.edu.au>.
o stars() now has add and plot arguments.
DEPRECATED & DEFUNCT
o The assignment operator...
2002 Oct 01
6
R-1.6.0 is released
...`priority'.
o termplot() uses factor levels rather than 1,2,3... for x-axis.
o Workaround for optimization bugs on gcc 3.1/2 on 32-bit Solaris.
o The trace() function has been robustified and a new function
tracingState() added to turn tracing temporarily on and off.
o New cophenetic() in "mva" as utility for hierarchical clustering.
o p.adjust() has two new methods, 'Hommel' and 'FDR', contributed
by Gordon Smyth <smyth at wehi.edu.au>.
o stars() now has add and plot arguments.
DEPRECATED & DEFUNCT
o The assignment operator...