24小时热门版块排行榜    

查看: 12373  |  回复: 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的回帖

chenhezhou

新虫 (小有名气)

您好,请问您的分析数据是怎样的啊,能给我看看吗
7楼2020-07-08 15:58:07
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0857 资源与环境 285分 +3 未名考生 2026-03-09 3/150 2026-03-10 23:25 by rock_sdust
[考研] 313求调剂 +4 一朵人头 2026-03-04 4/200 2026-03-10 23:05 by Equinoxhua
[考研] 材料调剂,307分 +11 张泳铭1 2026-03-09 12/600 2026-03-10 22:52 by 星空星月
[考研] 070300化学 280 一志愿太原理工 求调剂 +9 拾玖壹 2026-03-04 9/450 2026-03-10 22:51 by peike
[基金申请] PDF在线压缩软件 10+3 haxia 2026-03-08 4/200 2026-03-10 13:44 by lfy8008
[考博] 求材料读博院校 +9 yanglei131 2026-03-08 9/450 2026-03-10 11:54 by 天道酬勤2026686
[考研] 物理学求调剂 +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
[考研] 0703化学调剂 +7 小奶苍 2026-03-07 7/350 2026-03-10 05:17 by anny19840123
[硕博家园] 2026级硕士研究生招生/调剂 +3 知足常乐的乐 2026-03-06 5/250 2026-03-09 20:56 by 木瓜膏
[考研] 310 070300化学求调剂 +4 扑风铃的猫 2026-03-08 5/250 2026-03-09 19:46 by 30660438
[考研] 一志愿南大化学339分求调剂,四六级已过,有比赛,有文章 +7 Gallantzhou 2026-03-07 7/350 2026-03-09 18:38 by 30660438
[考研] 334求调剂 +5 Trying] 2026-03-07 5/250 2026-03-09 14:47 by ruiyingmiao
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
[考研] 材料与化工,291分,求调剂 +13 咕噜咕噜123123 2026-03-04 15/750 2026-03-09 09:32 by 30660438
[考研] 303求调剂 +8 forgman95 2026-03-05 10/500 2026-03-08 12:41 by 蓝莓都给你
[考研] 哈尔滨理工大学2026年研究生调剂,材料科学与化学工程学院研究生调剂 +3 xinliu866 2026-03-06 3/150 2026-03-06 19:42 by 祁东原神
[考研] 不限学校专业的调剂同学看过来 +5 啊摆啊摆 2026-03-05 9/450 2026-03-06 12:06 by 啊摆啊摆
[考研] 纺织、生物、化学、材料等专业 +3 Eember. 2026-03-05 7/350 2026-03-05 16:35 by Eember.
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见