 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 12356 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
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 ] |
回復(fù)此樓» 猜你喜歡
0817找調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑的同學(xué)，走過路過，不要錯過
已經(jīng)有14人回復(fù)
-
288求調(diào)劑085600材料與化工
已經(jīng)有14人回復(fù)
-
【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工，服從調(diào)劑，抗壓能力強！
已經(jīng)有11人回復(fù)
-
0701-322 求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工求調(diào)劑
已經(jīng)有3人回復(fù)
-
298求調(diào)劑
已經(jīng)有7人回復(fù)
-
288求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料與化工290求調(diào)劑
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有13人回復(fù)
wuchao924929
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 46
- 散金: 5
- 紅花: 2
- 帖子: 40
- 在線: 42.3小時
- 蟲號: 1944264
- 注冊: 2012-08-18
- 專業(yè): 環(huán)境工程
5樓2019-08-29 21:01:09
xp198766
鐵桿木蟲 (著名寫手)
小木蟲職業(yè)打醬油滴~~！
- 應(yīng)助: 5 (幼兒園)
- 金幣: 6525.1
- 散金: 262
- 紅花: 25
- 帖子: 2193
- 在線: 554小時
- 蟲號: 1046859
- 注冊: 2010-06-24
- 性別: GG
- 專業(yè): 環(huán)境微生物學(xué)
2樓2019-04-22 08:24:09
| 一直看賴?yán)蠋煹臅臀恼拢犝n的感覺完全不同。老師很專業(yè)，講解很耐心，小白也收獲很多。幫老師宣傳一下，有興趣的朋友們可以報名。http://mp.weixin.qq.com/s?__biz= ... 3&srcid=#rd |
4樓2019-05-27 11:14:53
6樓2020-05-01 13:22:25
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑085600材料與化工 +12 | Daunrin 2026-03-07 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 5/250 |
|
|
[考研] 新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生 +15 | another12 2026-03-04 | 21/1050 |
|
|
|
[考研]
材料325求調(diào)劑
 30+6
|
mariusuki 2026-03-02 | 12/600 |
|
|
|
[考研] 歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +20 | lisimayy 2026-03-04 | 32/1600 |
|
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 7/350 |
|
|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
|
[論文投稿]
 主編終審卡者不動
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +6 | 鶴遨予卿 2026-03-04 | 8/400 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
|
[考研] 0856材料工程，初試313調(diào)劑 +7 | 賣個關(guān)子吧 2026-03-03 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|