| 5 | 1/1 | 返回列表 |
| 查看: 12359 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
VPA分析(Vegan包Varpart程序) 已有1人參與
|
||
|
這是varpart自帶example,求大神解讀:三個(gè)因子,為什么跑完程序圖中會(huì)出現(xiàn)四個(gè)圈? 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 ] |
新蟲 (初入文壇)
鐵桿木蟲 (著名寫手)
小木蟲職業(yè)打醬油滴~~!
| 一直看賴?yán)蠋煹臅臀恼,聽課的感覺完全不同。老師很專業(yè),講解很耐心,小白也收獲很多。幫老師宣傳一下,有興趣的朋友們可以報(bào)名。http://mp.weixin.qq.com/s?__biz= ... 3&srcid=#rd |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑,307分 +4 | 張泳銘1 2026-03-09 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +3 | wyx79618 2026-03-06 | 6/300 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 7/350 |
|
|
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3 | 0703楊悅305分 2026-03-05 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|