24小时热门版块排行榜    

查看: 12366  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhenlisha

银虫 (小有名气)

[求助] 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

changling

银虫 (小有名气)

看不懂 这个几个是啥意思,特请教~ subsdens + watrcont, ~ substrate + shrub + topo,
6楼2020-05-01 13:22:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!


jjdg: 金币+1, 感谢参与回帖 2019-05-07 09:11:52
你最后都是  ## Four explanatory tables,有4项,肯定有4个圈啊
2楼2019-04-22 08:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenlisha

银虫 (小有名气)

一直看赖老师的书和文章,听课的感觉完全不同。老师很专业,讲解很耐心,小白也收获很多。帮老师宣传一下,有兴趣的朋友们可以报名。http://mp.weixin.qq.com/s?__biz= ... 3&srcid=#rd
4楼2019-05-27 11:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchao924929

新虫 (初入文坛)

【答案】应助回帖

楼主请问下你有这个的原始数据吗?我想看下数据格式
5楼2019-08-29 21:01:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 读博申请 +4 感dd 2026-03-10 6/300 2026-03-10 17:01 by nxgogo
[考研] 311求调剂 +3 牛乳糖的卡卡 2026-03-10 3/150 2026-03-10 16:19 by 球场大飞机
[考研] 物理学求调剂 +4 sx宋 2026-03-05 4/200 2026-03-10 11:11 by changlu999
[考研] 考研一志愿长安大学材料与化工309分请求调剂 +6 dtdxzxx 2026-03-06 8/400 2026-03-10 08:31 by Linda Hu
[考研] 334求调剂 +5 Trying] 2026-03-07 5/250 2026-03-09 14:47 by ruiyingmiao
[考研] 【求调剂】293分环境工程求调剂材料/化工,服从调剂,抗压能力强! +13 xiiiia 2026-03-04 14/700 2026-03-09 14:06 by macy2011
[考研] 求调剂,数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 一志愿苏大材料工程专硕293求调剂 +5 钢铁大炮 2026-03-04 5/250 2026-03-09 09:49 by EBSD
[考研] 0817找调剂 +6 kk扛 2026-03-08 6/300 2026-03-09 06:38 by houyaoxu
[考研] 083000环境科学与工程调剂 +5 加油呀fxy 2026-03-07 6/300 2026-03-08 21:07 by sdgdsdgd
[考研] 招调剂 +20 帆船哥 2026-03-04 23/1150 2026-03-08 16:25 by taycat
[考研] 0817化学工程与技术312分求调剂 +7 T123 tt 2026-03-04 7/350 2026-03-08 11:18 by 我的船我的海
[考研] 材料科学(0805)338 求调剂 +7 xiaokang3286 2026-03-07 7/350 2026-03-08 03:42 by lfhuang
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 求调剂 +3 snw石 2026-03-07 3/150 2026-03-07 20:01 by houyaoxu
[考研] 第一志愿上海大学,专业化学工程与技术,总分288,求调剂 +3 1829197082 2026-03-07 3/150 2026-03-07 19:14 by houyaoxu
[考研] 085600材料与化工 292分 求调剂 +6 程晴之 2026-03-06 6/300 2026-03-07 09:22 by 斩魂滴兔子!
[考研] 材料调剂 +10 ounce. 2026-03-04 12/600 2026-03-07 09:14 by Mornach1988
[考研] 085600材料调剂 总分330 +6 池池丶 2026-03-03 6/300 2026-03-06 15:50 by @飒飒飒飒
[考研] 复试调剂 +5 呼呼?~+123456 2026-03-05 5/250 2026-03-05 19:43 by caszguilin
[考研] 材料调剂 +4 L9370 2026-03-05 4/200 2026-03-05 19:36 by caszguilin
信息提示
请填处理意见