我尝试使用bfFixefLMER_t.fnc
或fitLMER.fnc
从LMERConvenienceFunctions
包中。在这两种情况下,我都会收到“输入模型不是 mer 对象”的错误。
我尝试了http://artax.karlin.mff.cuni.cz/r-help/library/LMERConvenienceFunctions/html/00Index.html中的示例。我得到同样的错误。
例如,当我从示例中运行时
fitLMER.fnc(mB, backfit.on = "t", item = FALSE,
ran.effects = c("(FreqB | Subject)",
"(LengthB | Subject)", "(WMC | Item)"))
这是我得到的结果。
Warning in fitLMER.fnc(mB, backfit.on = "t", item = FALSE, ran.effects = c("(FreqB | Subject)", :resetting argument "method" to "t"
**backfitting fixed effects**
Warning in bfFixefLMER_t.fnc(model = model, item = item, method = method, :factor variable with more than two levels in model terms, backfitting on t-values is not appropriate, please use function "bfFixefLMER_F.fnc" instead.
Error in bfFixefLMER_t.fnc(model = model, item = item, method = method, : the input model is not a mer object
有人对这些功能有过这种经验吗?
有一些函数可以反向拟合固定效应和正向拟合随机效应。有没有办法对模型的固定效应进行前向拟合glmer
?还是这在统计上毫无意义?我正在做生态建模,所以我对高级统计的理解并不多,所以,如果有人可以用外行的话更好地解释一下
sessionInfo()
R 版本 3.0.1 (2013-05-16) 平台:x86_64-pc-linux-gnu (64-bit)
语言环境:[1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_GB.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_GB.UTF-8
[ 7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
附加的基础包:[1] stats graphics grDevices utils datasets methods base
其他附加包:[1] LMERConvenienceFunctions_2.0 lme4_0.99999911-8
[3] RcppEigen_0.3.1.2.1 Rcpp_0.10.4
[5] Matrix_1.0-12 lattice_0.20-23
[7] LCFdata_1.0
通过命名空间加载(未附加):[1] grid_3.0.1 MASS_7.3-28 minqa_1.2.1 nlme_3.1-111 rpart_4.1-2
[6] splines_3.0.1 tools_3.0.1