| 5 | 1/1 | 返回列表 |
| 查看: 12361 | 回復: 7 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
VPA分析(Vegan包Varpart程序) 已有1人參與
|
||
|
這是varpart自帶example,求大神解讀:三個因子,為什么跑完程序圖中會出現(xiàn)四個圈? 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è),講解很耐心,小白也收獲很多。幫老師宣傳一下,有興趣的朋友們可以報名。http://mp.weixin.qq.com/s?__biz= ... 3&srcid=#rd |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 新疆大學地質(zhì)與礦業(yè)工程學院招生 +24 | another12 2026-03-04 | 32/1600 |
|
|---|---|---|---|---|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +4 | 每天散步 2026-03-09 | 5/250 |
|
|
[考研] 310 070300化學求調(diào)劑 +4 | 撲風鈴的貓 2026-03-08 | 5/250 |
|
|
[考研] 327分求調(diào)劑086 +3 | 西紅柿?小帥 2026-03-09 | 6/300 |
|
|
[考研] 083000環(huán)境科學與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 7/350 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 化學290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 085600材料與化工 298 調(diào)劑 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風了呢 2026-03-05 | 9/450 |
|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 274求調(diào)劑 +9 | 一個學習者 2026-03-04 | 9/450 |
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗 +5 | 錘某人 2026-03-04 | 8/400 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
[考研] 0703化學求調(diào)劑,總分320分,一志愿中南民族大學 +3 | 19286712440 2026-03-04 | 3/150 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|