24小时热门版块排行榜    

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

wuchao924929

新虫 (初入文坛)

【答案】应助回帖

楼主请问下你有这个的原始数据吗?我想看下数据格式
5楼2019-08-29 21:01:09
已阅   回复此楼   关注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的回帖

changling

银虫 (小有名气)

看不懂 这个几个是啥意思,特请教~ subsdens + watrcont, ~ substrate + shrub + topo,
6楼2020-05-01 13:22:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 086000生物与医药319分求调剂 +4 Tolkien 2026-03-07 8/400 2026-03-10 21:34 by Tolkien
[考研] 313分生物学求调剂 +5 Yyt杨1 2026-03-09 7/350 2026-03-10 18:58 by 清风月
[考研] 考研材料与化工,求调剂 +6 戏精丹丹丹 2026-03-09 6/300 2026-03-10 17:02 by BruceLiu320
[考研] 材料专硕调剂 +7 慕辰123 2026-03-05 10/500 2026-03-10 14:22 by peike
[考研] 085701环境工程295求调剂,一志愿A区211 +8 等春来, 2026-03-06 9/450 2026-03-10 10:38 by 沙漠之狐994
[考研] 085600材料与化工 326 求调剂 +4 热爱生活ing 2026-03-09 4/200 2026-03-10 08:25 by Linda Hu
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
[考研] 307求调剂 +3 辛仁豆腐 2026-03-08 5/250 2026-03-09 14:09 by macy2011
[考研] 考研调剂,一志愿山东大学材料与化工,328分,政治51 +5 关你西红柿929 2026-03-08 6/300 2026-03-09 13:50 by 新篇章DFSS
[考研] 293一志愿华东理工 0817化学工程与技术 调剂 +5 fjj0912 2026-03-07 5/250 2026-03-09 09:13 by 30660438
[考研] 安徽农业大学材料与化学学院0856材料招收调剂 +3 akakk47 2026-03-05 9/450 2026-03-08 21:03 by deng.z.1314
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
[考研] 0703化学求调剂学校 +3 一志愿京区211 2026-03-04 3/150 2026-03-07 21:35 by yinhuanshun
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 085600,一志愿郑州大学,280分求调剂 +7 Wuqi725 2026-03-05 7/350 2026-03-06 08:12 by 杨杨杨紫
[考研] 274求调剂 +9 一个学习者 2026-03-04 9/450 2026-03-05 20:56 by 黑衣馒头人
[考研] 求调剂,学校研究所都可以,材料与化工267分 +6 wmx1 2026-03-05 6/300 2026-03-05 20:01 by houyaoxu
[考研] 271求调剂 +7 月色c 2026-03-05 8/400 2026-03-05 19:51 by wangjihu
[考研] 070300化学 280 一志愿太原理工 求调剂 +8 拾玖壹 2026-03-04 8/400 2026-03-05 08:52 by anny19840123
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见