王胖的生信笔记第13期：随便搞点数据玩玩

这周有点累了，搞个简单的小花活儿给大家瞅瞅。

第9期的时候，搞了个08年NBA前五十的球员数据。

这次换个方式展示

library(ggplot2)
library(pheatmap)
library(data.table)
nba = read.csv("http://datasets.flowingdata.com/ppg2008.csv")
View(nba)
row <- nba[,1]
nba <- nba[,-1]
rownames(nba) <- row
nbanew <- t(nba)

这些和上次一样。

然后就用到了Y叔（如果搞生物的你没有用过他的clusterProfiler包，我觉得这是你的一大损失）的ggtree包。

library(ggtree)
mydata<-hclust(dist(nbanew))
plot(mydata)

这是平平无奇的一个树状分类图。

哎，让ggtree来施展一些魔法，搞一搞坐标变换。

ggtree(mydata, layout = "circular")

就变成这样了。

当然可以再来点花的

ggtree(mydata,layout = "circular",size = 0.8 )+ geom_tiplab2(offset=0.1, size=3)+ geom_highlight(node=24,fill="red",alpha=0.5)+ geom_highlight(node=26,fill="blue")+ geom_highlight(node=27,fill="orange")

再来点花的。

ggtree(mydata,layout = "circular",size = 0.8 )+ geom_tiplab2(offset=0.1, size=3)+ geom_highlight(node =24,fill="red",alpha=0.5)+ geom_highlight(node=26,fill="blue")+ geom_highlight(node=27,fill="orange") + geom_cladelabel(node=24,label="PT", offset=50,barsize =2,hjust=-1, vjust=-0.5,color="red")+ geom_cladelabel(node=26,label="Offensive",offset=50,barsize =2,hjust=1.5,vjust=-0.5,color="blue")+ geom_cladelabel(node=27,label="Defensive",offset=50,barsize=2,hjust=2,color="orange")