| 查看: 12159 | 回复: 7 | |||
[求助]
VPA分析（Vegan包Varpart程序）已有1人参与
|
|
这是varpart自带example，求大神解读：三个因子，为什么跑完程序图中会出现四个圈？ data(mite) data(mite.env) data(mite.pcnm) # two explanatory data frames -- hellinger-transform y mod <- varpart(mite, mite.env, mite.pcnm, transfo="hel" ) mod ## use fill colours showvarparts(2, bg = c("hotpink","skyblue" )) plot(mod, bg = c("hotpink","skyblue" )) ## test fraction [a] using partial rda, '~ .' in formula tells to use ## all variables of data mite.env. afrac <- rda(decostand(mite, "hel" ) , mite.env, mite.pcnm) anova(afrac) ## rsquareadj gives the same result as component [a] of varpart rsquareadj(afrac) ## partition bray-curtis dissimilarities varpart(vegdist(mite), mite.env, mite.pcnm) ## three explanatory tables with formula interface mod <- varpart(mite, ~ subsdens + watrcont, ~ substrate + shrub + topo, mite.pcnm, data=mite.env, transfo="hel" ) mod showvarparts(3, bg=2:4) plot(mod, bg=2:4) ## use rda to test fraction [a] ## matrix can be an argument in formula rda.result <- rda(decostand(mite, "hell" ) ~ subsdens + watrcont + condition(substrate + shrub + topo) + condition(as.matrix(mite.pcnm)), data = mite.env) anova(rda.result) ## four explanatory tables mod <- varpart(mite, ~ subsdens + watrcont, ~substrate + shrub + topo, mite.pcnm[,1:11], mite.pcnm[,12:22], data=mite.env, transfo="hel" ) mod plot(mod, bg=2:5) ## show values for all partitions by putting 'cutoff' low enough: plot(mod, cutoff = -inf, cex = 0.7, bg=2:5) [ last edited by jjdg on 2019-5-7 at 09:11 ] [ Last edited by jjdg on 2019-5-7 at 09:11 ] |
回复此楼» 猜你喜欢
博士申请都是内定的吗？
已经有3人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
博士读完未来一定会好吗
已经有27人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
xp198766
铁杆木虫 (著名写手)
小木虫职业打酱油滴~~！
- 应助: 5 (幼儿园)
- 金币: 6525.1
- 散金: 262
- 红花: 25
- 帖子: 2193
- 在线: 554小时
- 虫号: 1046859
- 注册: 2010-06-24
- 性别: GG
- 专业: 环境微生物学
2楼2019-04-22 08:24:09
3楼2019-05-06 15:44:13
| 一直看赖老师的书和文章，听课的感觉完全不同。老师很专业，讲解很耐心，小白也收获很多。帮老师宣传一下，有兴趣的朋友们可以报名。http://mp.weixin.qq.com/s?__biz= ... 3&srcid=#rd |
4楼2019-05-27 11:14:53
wuchao924929
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 46
- 散金: 5
- 红花: 2
- 帖子: 39
- 在线: 41.9小时
- 虫号: 1944264
- 注册: 2012-08-18
- 专业: 环境工程
5楼2019-08-29 21:01:09
6楼2020-05-01 13:22:25
chenhezhou
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 803.3
- 散金: 88
- 红花: 2
- 帖子: 228
- 在线: 105小时
- 虫号: 3773922
- 注册: 2015-03-30
- 性别: GG
- 专业: 环境工程
7楼2020-07-08 15:58:07
8楼2021-03-28 20:22:20