Doğrusal karışık efekt modelleri için kaldıraç ve Cook mesafeleri nasıl ayıklanır / hesaplanır

Herkes bir mersınıf nesnesi ( lme4paket yoluyla elde edilen ) kaldıraç ve Cook mesafeleri hesaplamak (veya çıkarmak) biliyor mu ? Bunları artık analizi için planlamak istiyorum.

— Roey Angel
kaynak

R paketine bir göz atmalısınız influence.ME. Tarafından oluşturulan karışık efekt modelleri için etkili veri ölçümlerini hesaplamanızı sağlar lme4.

Örnek bir model:

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

İşlev influence, tüm diğer adımların temelini oluşturur:

library(influence.ME)
infl <- influence(model, obs = TRUE)

Cook mesafesini hesapla:

cooks.distance(infl)

Cook'un uzaklığını çizin:

plot(infl, which = "cook")

resim açıklamasını buraya girin

— Sven Hohenstein
kaynak

Teşekkürler! Bu kesinlikle yardımcı olur. Bir Cook'un mesafe ve kaldıraç planının kaldıraçını hesaplamaya ne dersiniz?

— Roey Angel

@RoeyAngel Sanırım bu influence.MEpaketle mümkün değil . Ne yazık ki, bu görev için bir çözümüm yok.

— Sven Hohenstein

Olmamalı infl <- influence(model, group = "cyl"), çünkü rastgele etkiyi(1|cyl) ? Bilmiyorum, bunu hiç anlamadım, sadece etki yükledim ... ama ne zaman obs = TRUEve ne zaman kullanılacağını gerçekten bilmiyorum group...

— Meraklı

Aşağıdakileri eklemek istiyorum: Cook'un D mesafelerinin gerçekleştiği satır numarasını almak isterseniz - çizimde çizimde aynı sayı ortaya çıkıyorsa, Cooks 'D mesafeleri numaraları hakkında bir kesim ile aşağıdaki r formülünü kullanabilirsiniz. 0,1 değeri

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))]

— Elias Estatistics

Bu, burada önerilenhatvalues() işlevden daha iyi mi?

— Meraklı