問題タブ [multilevel-analysis]

現在、T 週間の I 店舗のパネル データを使用して階層ベイズモデルを実装しています。従属変数はブランドの売上です。使用している説明変数に店舗特性を関連付けることを目指していますが、2 つの別個の第 2 レベルが必要です。私のモデルは次のようになります。

y = 切片 + アルファ * X + ベータ * W + 誤差項、

ここで、2 番目のレベルを紹介します。

アルファ = ラムダ 1 * Z + エラー

ベータ = ラムダ 2 * Z + エラー、

ここで、Z には店舗特性が含まれます。ただし、私の知る限り、 in を使用する場合は 1 レベルしか使用できませrhierLinearModelんbayesm。

モデルを指定する方法やコードを調整して 2 つの第 2 レベルを取得する方法を知っている人はいますか?

事前にどうもありがとうございました！

2 つの反復測定の固定効果と変量効果のモデル化について質問があります。私のデータセットは次のようになります。

したがって、2 つの条件を持つタスクがあり、すべての参加者が両方の条件 (反復測定) を実行し、すべての条件が 3 つの試行 (反復測定) で存在します。そして連続DV。

(1) まず、変数 subject id、condition、trial が入れ子になっているのか交差しているのかが気になります。すべての科目は両方のレベルの条件を実行しているため、科目と条件が交差していると思いました。あれは正しいですか？
条件 1 と条件 2 は両方とも 3 つのレベルの試行を行っていますが、もちろん時間的には同じ試行ではありません (参加者はタスクの条件を同時にではなく、互いに連続してプレイするため)。これは、条件と試行がネストされていることを意味しますか? それとも、どちらも 3 つのレベルの試練があるため、交差していますか?
さらに、参加者全員が 6 段階の試練を経験しているため、科目と試練が交差していると思います。これは正しいです？

(2) 次に、マルチレベル モデルを実行したいのですが、固定効果と変量効果を指定する方法がわかりません。最低レベルは試行/時間になります。これは、両方の条件で時間の経過に伴う DV の曲線/成長軌道 (試行) をモデル化するためです。実際には、条件ごとに 20 回の試行があります。すべての参加者のランダム インターセプトをモデル化する次のコードが作成されました。

M2 は、M1 よりもデータへの適合性が大幅に向上しています。しかしこれは、参加者ごとに 6 回の試行すべてでランダムな切片がモデル化されているということでしょうか? 私は実際に条件ごとにこれを知りたいです。では、条件ごとに複数レベルのモデルを作成する必要があります (両方の条件のすべての分析を分離します)。または、コードに条件を含めることができますか?コードは次のようにする必要があります:

M3 が M1 よりもデータへの適合性が大幅に優れている場合、両方の条件にランダムな切片があることを意味しますか?

(3) 条件ごとに試行の効果があるかどうかを知りたい場合、これをどのようにモデル化すればよいですか? このような？

試験の有意な効果は、両方の条件内で試験の線形効果があることを意味しますか?

(4) 最後に、DV とトライアルの関係にランダムな勾配があるかどうかを条件ごとに知りたいです (つまり、6 つのトライアル全体ではありません)。6回の試行全体についてこれを知りたくありませんが、条件ごとに（実際には条件ごとに20回の試行）。モデルは次のようになります。

全体として、モデル内で繰り返される 2 つの変数 (状態と試行) をどのように扱うか混乱しています。誰かが助けてくれることを願っています！

ランダム インターセプト モデルを推定するコードを R から複製しようとしています。Rコードは次のとおりです。

lme4 パッケージの lmer コマンドを使用して、異なるグループ (groupid で定義) の観測値の変数 resid のランダムな切片を推定しています。「固定効果」部分がないため、(1|groupid) の前に変数はありません。さらに、各グループの切片を取得するために定数を推定する必要はありません。

Python で同様の推定を行う方法がわかりません。私は次のようなものを試しました：

A は前の例の resid ですが、groupid は同じです。

1) mdf.random_effects が探しているランダムな切片であるかどうかわかりません

2) 変数 B を削除できません。これは固定効果部分であると理解しています。私が試してみると：

「配列を空にすることはできません」というエラーが表示されます。

Rコードとまったく同じモデルを推定しようとしています。アドバイスをいただければ幸いです。