R's lmeまたはlmer関数(または の別のランダム効果ルーチン) でランダム効果の分散を一定に保つことができるかどうか、Rまたは少なくとも開始値を提供できるかどうか疑問に思っていました。
これは、 のparmsステートメントを使用して SAS で実行できるようPROC MIXEDです。Selyaらによる論文で。(2012)著者はこれを使用して、単純な固定効果構造を持つモデルの分散パラメーターを完全なモデルの分散パラメーターに設定します。
彼らが使用する PROC MIXED 内の特定の呼び出しは、parms/parmsdata = fullmodel.AB hold = ... 彼らの目標は、さまざまな固定効果構造を持つモデル間で分散推定値を一定に保つことです (ただし、これが SAS または R のいずれかで本当に可能かどうかは疑問です)。