線形混合効果モデルのレバレッジとクックの距離を抽出/計算する方法
回答:
Rパッケージを見てくださいinfluence.ME。によって生成された混合効果モデルの影響力のあるデータの測定値を計算できますlme4。
モデルの例：
library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)この機能influenceは、以降のすべてのステップの基礎となります。
library(influence.ME)
infl <- influence(model, obs = TRUE)クックの距離を計算します。
cooks.distance(infl)プロットクックの距離：
plot(infl, which = "cook")
ありがとう！これは確かに役立ちます。クックの距離とレバレッジのプロットのレバレッジを計算してみませんか？
—
Roey Angel 2013
@RoeyAngelこれは
—
Sven Hohenstein 2013
influence.MEパッケージでは不可能だと思います。残念ながら、このタスクの解決策はありません。
infl <- influence(model, group = "cyl")ランダムな効果を次のように指定したので、そうではありません(1|cyl)か？わからない、まったくわからない、インフルエンスをインストールしただけだ...しかしobs = TRUE、いつ使用するか、いつ使用するかわからないgroup...
私は以下を追加したいと思います：クックのD距離が発生する行番号を取得したい場合-プロットせずにプロットで同じ番号が発生する場合、カット付きのクックのD距離番号について次のr式を使用できますたとえば0.1のオフ値
—
エリアス統計
cooksD_data<-as.data.frame(cooks.distance(ft1)) cooksD_data_select<-cooksd[cooksD_data>0.1,drop=FALSE,] cooksD_oultiers<-as.numeric(rownames(cooksD_data_select))]