Rのパッケージを使用してlsmeans/emmeans、treatA (バイナリ/因子変数) のレベル間の応答におけるペアワイズ比較のプロットを作成しています。lsmeans (コントラスト) を使用して差の推定値を取得できますが、信頼限界ではなく、推定値の SE のみが提供されます。信頼限界は、対比ではなく、個々の効果に対してのみ提供されます。平均差 (コントラスト) の CL を生成するのを手伝ってくれる人はいますか?
model = glmer.nb(response ~ treatA * treatB + (1|random), data, family=nbinom1)
lsm <- lsmeans(model, pairwise ~ treatA*treatB)
contrast(lsm, method = "pairwise")
最終的にはこのようなプロットを作成したいと考えていますが、SE ではなく CL を使用します。
