問題タブ [lmertest]

いくつかの選択基準を使用して 2 つのモデルを比較したいと思います。ただし、BIC および CAIC の結果と比較して、AICC では異なる結果が得られます。Gurka (2006) によると、これらの 2 つのタイプの選択基準は、モデルの定式化に応じて異なる反応を示します。どうやら、Gurka (2006) によると、モデルの定式化に対数尤度関数の定数が含まれているかどうかによって、それらの適切性は異なります。SAS procmixed では、定数は含まれません。ただ、モデル式がよくわからないので、lmerwithのモデル計算に含まれているのかわかりません。lmerTest誰かが答えを知っていれば、私はそれを非常に感謝しています:)

文献: Gurka, MJ (2006): REML での最適な線形混合モデルの選択。中: アメリカの統計学者 60 (1)、S. 19–26。

以外の統計検定で対になった観測をどのように説明するのt-testですか? 以下では、混合効果アプローチでそうしようとして失敗した 2 つの例について説明します。

t.test(..., paired=T)例 1:での再現lm()方法

観測値がペアになっている場合、それを説明できますが、 （可能であれば）t.test()どのように行うのですか？lm()混合効果モデルのアプローチを使用しようとしましたが、次のようになります。

エラーが発生します：

その間：

実行されますが、固定効果パラメーターの推定値と関連する標準。エラーと t 値は、 によって生成されるものと同じですlm()。

例 2: 対応のある観測による線形回帰

私が作成したデータセットの観察が、30 日間隔で測定された被験者からのものであると想像してみましょう。つまり、100 人の被験者のそれぞれが 0 日目に測定され、次に 30 日目に測定されました。時間の経過に伴う変化率を推定したかったのです。

データは次のようになります (黒の線形回帰、赤い線で結ばれた対のデータ):

lm1観測のペアの性質を説明していません。データの各ペアの切片と勾配が異なることを可能にする混合効果モデルを実行することを考えましたが、R は、あまりにも多くのパラメーターを推定しようとしていることに再び抗議します。

データ ペアの切片が異なるが勾配が異なることを可能にする単純なモデル。つまり、次のようになります。

次のように訴えます。

しかし、 によって生成されるものと同一の固定効果推定値が実行され、生成されlm()ます。

[アップデート]

@Limey は、対応のある t 検定は、2 つのグループ間の対の差がゼロではないかどうかを評価する t 検定に他ならないことを思い出させてくれました。このような一対の差は、検定のほかに任意の対応のある統計的検定を実行するために使用できます。これを確認するために、3 つの異なる「応答」変数を作成しました。これらは、さまざまな方法で組み合わせx1てx2順序付けされています (それぞれ、元のランダムな順序x1、昇順とx2降順、両方とも昇順)。

対応する違いを計算しました：

そして、それらを使用して線形モデルを実行しました。

勾配の推定値が異なると予想していましたが、すべて同じでした。

それらの傾きの推定値は、対応する「対応のない」線形モデル ( lm(response~time)、lm(response2~time)、およびlm(response3~time)) から推定されたものと同じです。私は何が欠けていますか？