是否有办法使用包裹素食主义者的函数ordistep对部分规范排序(冗余分析或对应分析)执行变量减少?我检查了Borcard等人。 (2011)数字生态与R和我找不到答案。 我尝试了以下语法,它返回了下一个错误:
Error in formula.default(prda1) : invalid formula
prda1 <- rda(spp.h, env.partial, avu.rda)
rda.ordistep <- ordistep(rda(spp ~ 1, data = env),
scope = formula(partial.rda),
direction= "forward",
pstep=10000,
trace = FALSE)
答案 0 :(得分:1)
如果您在ordistep
中进行部分排序,则必须在较低范围内使用部分字词。以下示例适用于素食主义者
data(dune, dune.env)
require(vegan)
data(dune, dune.env)
m0 <- rda(dune ~ Condition(Management), dune.env)
m1 <- rda(dune ~ Condition(Management) + ., dune.env)
m <- ordistep(m0, scope = list(lower=m0, upper=m1))
如果您在较低范围内没有Condition(Management)
,则会将其视为任何其他变量并考虑删除,或者如果起始模型没有Condition(Management)
,则将其视为任何其他变量候选人,但永远不能添加。