Почему cook.distance () не работает для результата impact_mixed () в glmmTMB - PullRequest
Новая
       16

Почему cook.distance () не работает для результата impact_mixed () в glmmTMB

0 голосов
/ 27 мая 2020

Это то, что я копирую model_evaluation.pdf с материалами glmmTMB. Я не знаю, почему приведенный ниже код работал правильно, а сегодня он не работает. Кто-нибудь подскажет, зачем и доработал? Спасибо за 


L <- load(system.file("vignette_data","model_evaluation.rda",
                      package="glmmTMB"))

source(system.file("other_methods","influence_mixed.R", package="glmmTMB"))



owls_nb1 <- glmmTMB(SiblingNegotiation ~ FoodTreatment*SexParent +
                      (1|Nest)+offset(log(BroodSize)),
                    contrasts=list(FoodTreatment="contr.sum",
                                   SexParent="contr.sum"),
                    family = nbinom1,
                    zi = ~1, data=Owls)



owls_nb1_influence_time <- system.time(
  owls_nb1_influence <- influence_mixed(owls_nb1, groups="Nest")
)

inf <- as.data.frame(owls_nb1_influence[["fixed.effects[-Nest]"]])
inf <- transform(inf,
                 nest=rownames(inf),
                 cooks=cooks.distance(owls_nb1_influence))

Каждый раз, когда я выполнял cook.distance(owls_nb1_influence), он показывал ошибку. Error in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), : 'data' must be of a vector type, was 'NULL'

Что делать, чтобы решить проблему? Спасибо заранее.

...