24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 12121  |  回复: 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的回帖

黑色幽灵

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by xp198766 at 2019-04-22 08:24:09
你最后都是  ## Four explanatory tables,有4项,肯定有4个圈啊

请问VPA分析结果解释度为负数是啥意思啊
8楼2021-03-28 20:22:20
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见