"Reaction ~ Days + (Days | Subject)" のような式です。
データフレームです (列のオブジェクトです)。
Optionaloptions: LmerOptions
REML、および contrasts と levels は lm と同じです (LmerOptions で受け付けるキーを参照)。
// R: lme4::lmer(y ~ 1 + (1 | g), data)
const fit = lmer("y ~ 1 + (1 | g)", {
y: [9.8, 10.4, 10.1, 12.3, 12.9, 12.0, 8.1, 7.6, 8.4],
g: ["a", "a", "a", "b", "b", "b", "c", "c", "c"],
});
fit.fixef; // [10.177778]
fit.sigma; // 0.39299
fit.varCorr[0].sd; // [2.1726]
fit.ranef[0].levels; // ["a", "b", "c"]
fit.ranef[0].values; // [[-0.076939], [2.1982], [-2.1213]]
fit.REMLcrit; // 19.015
REML (既定) または最尤法で線形混合モデルを当てはめます。R の
lme4::lmer()と同等です。メソッドに対応する量として、次を持ちます:fixef()、coef(summary())(p 値なし)、vcov()、sigma()、VarCorr()、ranef()、fitted()、residuals()、logLik()、AIC()、BIC()、deviance()、REMLcrit()、isSingular()。式の変量効果の項は、
(lhs | group)(groupの各水準の中でのlhsの列の、相関のある効果) または(lhs || group)(相関なし) です。groupはa:b(組み合わせ)、a/b(入れ子:b:aとa) などでもよく、lme4 と同じです。固定効果の部分は lm と同じように符号化されます。縮退した列は (メッセージとともに) 除かれます。式の変数に欠測値がある行は除かれます (使った行はrowsに示されます)。共分散パラメータ θ は、項の相対共分散因子の下三角部分です。推定値は、プロファイルされた基準の厳密な最小化点です (lme4 の最適化はその近くで止まるため、lme4 の値は通常、5 から 9 桁目の有効数字が異なります)。境界上の要素は厳密に 0 になり、そのとき当てはめは特異になります。計算は疎です (lme4 と同じく、疎なコレスキー分解を使います)。そのため、多くのグループも速く扱えます。
lmerSummary、lmerAnova、lmerContest など、当てはめを受け取る関数には、
lmerが返したオブジェクトそのものが必要です。そのコピー (structuredClone、{ ...fit }、JSON、postMessage) は、それらが使う内部の値を持たず、TypeError になります。コピーのフィールドは、値を保ちます。