bfFixefLMER_t.fncorfitLMER.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 ビット)
ロケール: [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 グラフィックス grDevices utils データセット メソッド base
その他の添付パッケージ: [1] LMERConvenienceFunctions_2.0 lme4_0.99999911-8
[3] RcppEigen_0.3.1.2.1 Rcpp_0.10.4
[5] Matrix_1.0-12lattice_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