| 8 | 1/1 | 返回列表 |
| 查看: 12375 | 回復(fù): 7 | |||
[求助]
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 |
新蟲 (初入文壇)
新蟲 (小有名氣)
新蟲 (初入文壇)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0857 資源與環(huán)境 285分 +5 | 未名考生 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑,307分 +12 | 張泳銘1 2026-03-09 | 13/650 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +3 | KRSLSR 2026-03-10 | 3/150 |
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
[考研] 085602化工求調(diào)劑 +7 | 董boxing 2026-03-10 | 7/350 |
|
|
[考研] 求調(diào)劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 8/400 |
|
|
[考博] 找博導(dǎo) +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
[考博] 2026博士申請 +6 | 起泡酒 2026-03-08 | 6/300 |
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑,四六級(jí)已過,有比賽,有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
[考研] 中科大材料299求調(diào)劑 +10 | DAIjiayo 2026-03-05 | 16/800 |
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|