2 つの IV、時間 (被験者内の 3 つのレベル、t1、t2、t3)、および修正タイプ (被験者間の 3 つのレベル) の実験があります。DV = 態度 (連続)。完全な lmer モデルは次のようになります。
agreement ~ correction * time + (1 + time|subject) + (1 + correction + time|item)
モデル全体から単純な効果を得るにはどうすればよいですか?
時点 3 で 3 つのレベルの修正すべてを比較したいと思いますが、時点 1 での個々の参加者のスコア (ベースライン メジャー) を考慮に入れます。これは、時点 1 からのベースライン スコアを継続的な尺度として使用した、時点 3 の ANCOVA のようなものです。
私がすることができます：
library(lsmeans)
ref_levs <- emmeans(mod2, "correction", by = "time", at = list(time = "t3"))
pairs(ref_levs)
しかし、t1 での個々の参加者のスコアを制御する比較の出力を得ることができません。
このようなモデルに合わせて、必要なものを得るためにデータセットをサブセット化する必要がありますか?
agreement_t3 ~ correction + time_t1 + (1|subject) + (1 + correction|item)
または、データセットをサブセット化することなくそれを取得する方法が emmeans にありますか?