赖江山的博客分享 http://blog.sciencenet.cn/u/laijiangshan 生态、统计与R语言

博文

很多人问twinspan在R怎么实现,我把代码放这里了。

已有 11006 次阅读 2016-3-17 11:33 |系统分类:科研笔记| TWINSPAN

install.packages ('devtools')

devtools::install_github("zdealveindy/twinspanR")

install.packages ('devtools')

devtools::install_github("zdealveindy/twinspanR")


library (twinspanR)

library (vegan)

data (danube)

res <- twinspan (danube$spe, modif = TRUE, clusters = 4)

k <- cut (res)

 

dca <- decorana (danube$spe)

par (mfrow = c(1,2))

ordiplot (dca, type = 'n', display = 'si', main = 'Modified TWINSPAN')

points (dca, col = k)

for (i in c(1,2,4)) ordihull (dca, groups = k, show.group = i, col = i,

 draw = 'polygon', label = TRUE)

ordiplot (dca, type = 'n', display = 'si', main = 'Original assignment\n (Ellenberg 1954)')

points (dca, col = danube$env$veg.type)

for (i in c(1:3)) ordihull (dca, groups = danube$env$veg.type,

 show.group = unique (danube$env$veg.type)[i], col = i,

 draw = 'polygon', label = TRUE)




https://blog.sciencenet.cn/blog-267448-963158.html

上一篇:画散点图光滑趋势线的R函数
下一篇:今天有人问到点格局分析在R里怎么实现,附上一个英文的参考资料
收藏 IP: 159.226.89.*| 热度|

1 王振亭

该博文允许注册用户评论 请点击登录 评论 (8 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-23 19:32

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部