if (FALSE) {
library(pvclust)
data(lung) # 916 genes for 73 subjects
set.seed(13134)
result <- pvclust(lung[, 1:20], method.dist = "cor", method.hclust = "average", nboot = 100)
dend <- as.dendrogram(result)
result %>%
as.dendrogram() %>%
hang.dendrogram() %>%
plot(main = "Cluster dendrogram with AU/BP values (%)")
result %>% text()
result %>% pvrect(alpha = 0.95)
dend %>%
pvclust_show_signif(result) %>%
plot()
dend %>%
pvclust_show_signif(result, show_type = "lwd") %>%
plot()
result %>% text()
result %>% pvrect(alpha = 0.95)
dend %>%
pvclust_show_signif_gradient(result) %>%
plot()
dend %>%
pvclust_show_signif_gradient(result) %>%
pvclust_show_signif(result) %>%
plot(main = "Cluster dendrogram with AU/BP values (%)\n bp values are highlighted by signif")
result %>% text()
result %>% pvrect(alpha = 0.95)
}
Run the code above in your browser using DataLab