This is an application motivated by this nice walkthrough.

library("treelapse")
library("knitr")
data(iris)
hc_iris <- hclust(dist(iris[, 1:4]), method = "complete")
edges <- ape::as.phylo(hc_iris)$edge
colnames(edges) <- c("parent", "child")
class(edges) <- "character"

X <- t(iris[, 1:4])
colnames(X) <- 1:150

values <- tree_fun_multi(
  edges,
  X,
  tree_mean
)

values$time <- rownames(X)[values$row]
timebox_tree(values, edges)
treebox(values, edges)


krisrs1128/treelapse documentation built on Jan. 3, 2020, 6:29 p.m.