最近涉及到画多边形的diversity的图
建立好的坐标系
set.seed(1234)
df <- data.frame(x = rnorm(20), y = rnorm(20), group = c(rep("Disease",12),rep("Healthy",8)))
数据举例:
# 计算凸点chull 以便于连线
find_hull <- function(df) df[chull(df$x, df$y),]
hulls <- ddply(df, "group", find_hull)
# 画图
geom_polygon 是ggplot2里面多边形连线的function
ggplot(df, aes(x=x, y=y, color=group)) + geom_polygon(data = hulls,alpha = 0.2) + geom_point(alpha=.7, size=2)
如下:
参考链接:
https://www.rdocumentation.org/packages/grDevices/versions/3.6.2/topics/chull