Ошибка при усреднении моделей GLM с использованием `model.avg ()` из пакета `MuMIn` - PullRequest
Новая
       37

Ошибка при усреднении моделей GLM с использованием `model.avg ()` из пакета `MuMIn`

0 голосов
/ 30 апреля 2019

Я анализирую некоторые данные, используя glmer(). Пожалуйста, потерпите меня, потому что я изучаю и R, и моделирование.

Мои данные состоят из двух переменных-предикторов Area и Day_Night, которые я использую для объяснения возникновения данного поведения среди Standing, Feeding или Foraging

Поскольку у меня есть данные по нескольким лицам, я бы хотел усреднить лучшую модель, выбранную для каждого человека, которая оказалась с таким же аддитивным эффектом, как я уже рассчитал AICc значения. Модель: 

best_model_individual1<-glmer(cbind(Standing,Feeding_Foraging) ~ Day_Night + Area+(1|ID) , data=GLM_df , family=binomial)

Однако при усреднении наилучшей модели для каждого человека я получаю ошибку: 

> models<- list(best_model_individual1,best_model_individual2,best_model_individual3)
> averages<-model.avg(models, beta = c("none"),
+           rank = NULL, rank.args = NULL, revised.var = TRUE,
+           dispersion = NULL, ct.args = NULL)
Error in model.avg.default(models, beta = c("none"), rank = NULL, rank.args = NULL,  : 
  models are not all fitted to the same data

Я не уверен, связано ли это с тем, что я мог допустить ошибку в аргументах функции, или с тем, что такое среднее значение не может быть получено статистически.

Я надеюсь, что кто-нибудь может дать мне толчок! Любой вклад приветствуется.

Я предоставляю dput() образец данных, используемых в модели для индивидуума1 ниже: 

> dput(GLM_df)
structure(list(V1 = 1:144, hour = c(0L, 1L, 2L, 3L, 4L, 5L, 6L, 
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 
20L, 21L, 22L, 23L, 0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 
11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 
0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 
15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 0L, 1L, 2L, 3L, 
4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 
18L, 19L, 20L, 21L, 22L, 23L, 0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 
8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 
21L, 22L, 23L, 0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L), 
    Feeding = c(4L, 3L, 2L, 7L, 8L, 7L, 7L, 7L, 6L, 7L, 8L, 7L, 
    7L, 7L, 9L, 13L, 5L, 4L, 3L, 5L, 5L, 7L, 5L, 5L, 7L, 5L, 
    5L, 11L, 14L, 10L, 9L, 9L, 9L, 9L, 10L, 9L, 10L, 9L, 10L, 
    9L, 8L, 5L, 4L, 8L, 11L, 11L, 11L, 9L, 6L, 3L, 3L, 8L, 11L, 
    11L, 9L, 9L, 9L, 9L, 10L, 9L, 10L, 9L, 7L, 9L, 8L, 3L, 4L, 
    6L, 7L, 10L, 7L, 7L, 4L, 4L, 2L, 8L, 11L, 11L, 9L, 10L, 9L, 
    8L, 8L, 7L, 6L, 6L, 5L, 2L, 4L, 3L, 3L, 4L, 4L, 6L, 7L, 6L, 
    4L, 4L, 3L, 6L, 7L, 6L, 5L, 6L, 5L, 4L, 4L, 5L, 4L, 3L, 3L, 
    2L, 4L, 3L, 3L, 4L, 4L, 4L, 5L, 4L, 7L, 3L, 4L, 6L, 5L, 6L, 
    7L, 9L, 8L, 8L, 10L, 10L, 10L, 8L, 6L, 5L, 8L, 5L, 3L, 6L, 
    6L, 7L, 6L, 6L), Foraging = c(23L, 24L, 24L, 19L, 18L, 20L, 
    18L, 19L, 19L, 17L, 16L, 17L, 17L, 17L, 14L, 10L, 20L, 24L, 
    24L, 22L, 21L, 21L, 21L, 22L, 18L, 20L, 21L, 9L, 5L, 9L, 
    10L, 12L, 12L, 10L, 10L, 10L, 9L, 9L, 7L, 5L, 14L, 20L, 20L, 
    16L, 11L, 12L, 12L, 15L, 15L, 19L, 18L, 7L, 5L, 5L, 9L, 9L, 
    9L, 8L, 6L, 8L, 7L, 8L, 6L, 5L, 10L, 17L, 17L, 15L, 13L, 
    12L, 13L, 13L, 22L, 24L, 26L, 10L, 3L, 6L, 8L, 9L, 9L, 6L, 
    7L, 8L, 8L, 6L, 3L, 2L, 8L, 24L, 23L, 18L, 15L, 14L, 15L, 
    19L, 19L, 21L, 24L, 10L, 4L, 6L, 8L, 8L, 6L, 6L, 7L, 6L, 
    6L, 6L, 3L, 3L, 11L, 23L, 24L, 21L, 18L, 17L, 11L, 16L, 19L, 
    23L, 22L, 6L, 5L, 5L, 7L, 7L, 8L, 8L, 7L, 9L, 6L, 4L, 3L, 
    2L, 7L, 22L, 20L, 18L, 15L, 12L, 12L, 15L), Standing = c(1L, 
    1L, 2L, 4L, 4L, 2L, 5L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 
    2L, 1L, 1L, 2L, 1L, 2L, 3L, 2L, 2L, 2L, 1L, 6L, 7L, 7L, 6L, 
    5L, 5L, 6L, 4L, 6L, 6L, 8L, 9L, 12L, 3L, 2L, 2L, 3L, 4L, 
    5L, 3L, 4L, 4L, 4L, 4L, 8L, 9L, 8L, 5L, 6L, 6L, 5L, 8L, 7L, 
    7L, 6L, 10L, 12L, 6L, 3L, 3L, 4L, 5L, 5L, 5L, 5L, 6L, 5L, 
    3L, 12L, 16L, 14L, 12L, 11L, 13L, 15L, 15L, 15L, 16L, 18L, 
    23L, 28L, 17L, 5L, 5L, 9L, 12L, 12L, 13L, 8L, 9L, 5L, 4L, 
    13L, 19L, 19L, 17L, 17L, 19L, 20L, 20L, 19L, 20L, 22L, 24L, 
    28L, 15L, 5L, 4L, 7L, 12L, 10L, 14L, 11L, 7L, 8L, 5L, 19L, 
    23L, 20L, 16L, 13L, 14L, 14L, 13L, 13L, 13L, 18L, 21L, 24L, 
    16L, 6L, 8L, 9L, 12L, 14L, 13L, 11L), ID = c("41361", "41361", 
    "41361", "41361", "41361", "41361", "41361", "41361", "41361", 
    "41361", "41361", "41361", "41361", "41361", "41361", "41361", 
    "41361", "41361", "41361", "41361", "41361", "41361", "41361", 
    "41361", "41365", "41365", "41365", "41365", "41365", "41365", 
    "41365", "41365", "41365", "41365", "41365", "41365", "41365", 
    "41365", "41365", "41365", "41365", "41365", "41365", "41365", 
    "41365", "41365", "41365", "41365", "41366", "41366", "41366", 
    "41366", "41366", "41366", "41366", "41366", "41366", "41366", 
    "41366", "41366", "41366", "41366", "41366", "41366", "41366", 
    "41366", "41366", "41366", "41366", "41366", "41366", "41366", 
    "41366bis", "41366bis", "41366bis", "41366bis", "41366bis", 
    "41366bis", "41366bis", "41366bis", "41366bis", "41366bis", 
    "41366bis", "41366bis", "41366bis", "41366bis", "41366bis", 
    "41366bis", "41366bis", "41366bis", "41366bis", "41366bis", 
    "41366bis", "41366bis", "41366bis", "41366bis", "41367", 
    "41367", "41367", "41367", "41367", "41367", "41367", "41367", 
    "41367", "41367", "41367", "41367", "41367", "41367", "41367", 
    "41367", "41367", "41367", "41367", "41367", "41367", "41367", 
    "41367", "41367", "41368", "41368", "41368", "41368", "41368", 
    "41368", "41368", "41368", "41368", "41368", "41368", "41368", 
    "41368", "41368", "41368", "41368", "41368", "41368", "41368", 
    "41368", "41368", "41368", "41368", "41368"), Area = structure(c(2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Loliondo", "Seronera"
    ), class = "factor"), Feeding_Foraging = c(27L, 27L, 27L, 
    26L, 25L, 27L, 25L, 26L, 25L, 24L, 24L, 24L, 24L, 24L, 23L, 
    23L, 25L, 28L, 27L, 27L, 27L, 27L, 27L, 27L, 25L, 25L, 26L, 
    19L, 19L, 19L, 20L, 21L, 20L, 20L, 20L, 19L, 19L, 18L, 16L, 
    14L, 22L, 25L, 24L, 24L, 22L, 22L, 23L, 23L, 21L, 22L, 21L, 
    16L, 15L, 16L, 18L, 18L, 17L, 17L, 16L, 17L, 17L, 17L, 13L, 
    14L, 18L, 21L, 22L, 21L, 20L, 22L, 20L, 20L, 10L, 9L, 6L, 
    20L, 27L, 25L, 22L, 21L, 21L, 23L, 23L, 22L, 22L, 24L, 28L, 
    30L, 21L, 8L, 8L, 13L, 17L, 18L, 19L, 15L, 22L, 26L, 27L, 
    16L, 11L, 11L, 13L, 13L, 11L, 10L, 11L, 11L, 10L, 9L, 6L, 
    5L, 15L, 26L, 27L, 25L, 22L, 21L, 16L, 20L, 26L, 26L, 26L, 
    11L, 10L, 12L, 15L, 16L, 16L, 16L, 17L, 18L, 17L, 12L, 10L, 
    8L, 15L, 27L, 24L, 24L, 21L, 19L, 18L, 21L), Feeding_Standing = c(4L, 
    4L, 4L, 11L, 11L, 10L, 12L, 10L, 10L, 10L, 11L, 10L, 10L, 
    11L, 13L, 17L, 7L, 5L, 4L, 7L, 6L, 9L, 8L, 7L, 9L, 7L, 6L, 
    16L, 21L, 16L, 15L, 14L, 14L, 15L, 15L, 15L, 16L, 17L, 19L, 
    21L, 11L, 7L, 6L, 11L, 16L, 15L, 14L, 13L, 10L, 7L, 7L, 16L, 
    19L, 18L, 14L, 15L, 14L, 14L, 18L, 16L, 17L, 15L, 17L, 21L, 
    14L, 6L, 8L, 10L, 13L, 15L, 12L, 12L, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, 12L, 10L, 7L, 19L, 26L, 25L, 22L, 23L, 24L, 
    25L, 24L, 24L, 25L, 25L, 27L, 29L, 19L, 7L, 7L, 11L, 15L, 
    14L, 19L, 15L, 14L, 11L, 10L, 25L, 28L, 26L, 23L, 22L, 22L, 
    22L, 24L, 23L, 24L, 26L, 28L, 29L, 24L, 10L, 12L, 15L, 17L, 
    20L, 18L, 18L), Standing_Foraging = c(24L, 25L, 26L, 23L, 
    21L, 22L, 22L, 22L, 22L, 20L, 18L, 20L, 21L, 21L, 18L, 14L, 
    22L, 25L, 25L, 24L, 23L, 23L, 24L, 24L, 20L, 22L, 22L, 14L, 
    12L, 16L, 16L, 17L, 17L, 16L, 15L, 16L, 15L, 17L, 16L, 17L, 
    17L, 22L, 22L, 18L, 15L, 16L, 15L, 19L, 20L, 23L, 22L, 15L, 
    13L, 13L, 14L, 15L, 14L, 13L, 14L, 15L, 14L, 13L, 16L, 17L, 
    16L, 20L, 21L, 19L, 18L, 16L, 19L, 18L, 28L, 29L, 30L, 22L, 
    19L, 20L, 21L, 20L, 21L, 21L, 22L, 22L, 24L, 24L, 26L, 30L, 
    26L, 29L, 28L, 27L, 28L, 26L, 27L, 27L, 28L, 27L, 28L, 23L, 
    23L, 25L, 25L, 25L, 26L, 26L, 26L, 25L, 26L, 27L, 28L, 30L, 
    26L, 28L, 28L, 28L, 30L, 26L, 25L, 27L, 27L, 31L, 28L, 25L, 
    28L, 25L, 23L, 21L, 22L, 23L, 20L, 21L, 20L, 22L, 24L, 27L, 
    23L, 28L, 28L, 27L, 26L, 26L, 24L, 26L), Day_Night = structure(c(2L, 
    2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L), .Label = c("Day", "Night"
    ), class = "factor")), row.names = c(NA, -144L), class = c("data.table", 
"data.frame"), .internal.selfref = <pointer: 0x00000000025e1ef0>)
...