我尝试了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位)
语言环境: [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
lme4
到最新的即将发布的版本(>1.0?)吗?这可能会破坏LMERConvenienceFunction
。我认为这将成为一个编程问题,而不是统计问题。请提供sessionInfo()
信息? - Ben Bolker