変量効果ロジスティック回帰のICCの計算

次の形式でロジスティック回帰モデルを実行しています。

lmer(response~1+(1|site), family=binomial, REML = FALSE)

通常、切片と残差からICCを計算しますが、モデルの要約には残差は含まれません。これをどのように計算しますか？

logistic lme4-nlme glmm

— ミーガン
ソース

なぜICCを計算しているのですか？

— AdamO

GLMMを使用する必要があるという証拠として、通常のロジスティック回帰がこれらのデータに対して有効でないという仮定をテストするため。ICClogit =切片の分散^ 2 /（切片の分散^ 2 + pi ^ 2/3）という方程式を見つけました。これは合理的ですか？

— ミーガン

完全な最尤法を使用しています。固定効果モデルに対して1自由度で尤度比検定を行うことはできませんか？

— AdamO

@ミーガン：あなたは正しい。実際には、Zeger等。（1988）が示唆している

作品よりも良い 2は非常に接近しているものの、ロジスティック回帰モデルの残差分散として。SL Zeger、KY Liang、およびPS Albertを参照してください。縦断データのモデル：一般化された推定方程式アプローチ。生体認証、44：1049-1060 1988.

(15 / 16)^{2} π^{2} / 3

$(15/16)^2\pi^2/3$

π^{2} / 3

$\pi^2/3$

— ランデル

@Megan：それはintercept_variance / (intercept_variance + pi^2/3)-分散を二乗しないでください。

— ヴォルフガング

ヘルプファイル?sjstats::iccには、バイナリ応答を含む混合モデルの数式への参照があります。

Wu S、Crespi CM、Wong WK。2012.がん予防クラスター無作為化試験におけるバイナリ応答のクラス内相関係数を推定する方法の比較。現代臨床試験33：869-880（doi：10.1016 / j.cct.2012.05.004）

ロジスティック回帰の残差はに固定され(pi ^ 2) / 3ます。

— ダニエル
ソース

この式の参照はありますか？

— ジーニー

私のことですか？OPの投稿で最初にコメントしていませんでしたか？

— ダニエル

@ Jeanine- ICCの引用：Moineddin、R.、Matheson、FI、＆Glazier、RH（2007）。マルチレベルロジスティック回帰モデルのサンプルサイズのシミュレーション研究。BMC Medical Research Methodology、7、34。doi.org

— 10.1186/