Question: Consensus Of Hierarchical Clustering
2
gravatar for Diana
8.8 years ago by
Diana840
Germany
Diana840 wrote:

Hi all!

I am permuting (1000 times) the columns of my matrix containing gene expression data and then using hierarchical clustering to cluster the data in R. I don't understand how to make a consensus tree out of my results and plot it. Someone told me about consensus function in maanova but it works only for a specific type of object with fitted ANOVA model so I don't want to use that. Can anyone suggest anything else in R? Here's my code:

make.permuted.clust <- function(i)
{ 
  permuted <- data.matrix[,sample(ncol(data.matrix), 12, replace=TRUE)]
  d = dist(permuted, method = "euclidean", diag = FALSE, upper = FALSE, p = 2)
  clust = hclust(d, method = "complete", members=NULL)
  clust
}
all.clust <- lapply(1:1000, make.permuted.clust)

Thanks!!

clustering • 3.1k views
ADD COMMENTlink modified 8.8 years ago by Zev.Kronenberg11k • written 8.8 years ago by Diana840
4
gravatar for Zev.Kronenberg
8.8 years ago by
United States
Zev.Kronenberg11k wrote:

You should check out the package PVCLUST. It can constuct confidence intervals around clustered groups.

PVCLUST WEB

You could also try bootstrapping your original cluster using boot.phylo in the APE package. It's designed for phylogentics but you can use any distance function to generate a tree.

ADD COMMENTlink modified 8.8 years ago • written 8.8 years ago by Zev.Kronenberg11k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2302 users visited in the last hour
_