Блог

Здравствуйте, я пытаюсь выполнить анализ соответствия для модели занятости за один сезон, созданной с использованием пакета без маркировки и оцененной с помощью пакета AICcmodavg. Вы можете загрузить исходную модель в виде RDS здесь

 library(unmarked)
library(AICcmodavg)
BestMylu <- readRDS("best2.My.Lu2.rds")
obs.boot <- mb.gof.test(BestMylu, nsim = 5000)
  

Я получаю следующую ошибку:

Ошибка в data.frame (det.hist, preds.psi) : аргументы подразумевают разное количество строк: 123, 111 Кроме того: Предупреждающие сообщения: 1: Некоторые наблюдения были отброшены, поскольку отсутствовали соответствующие ковариаты. 2: 12 сайты были удалены из-за отсутствия данных. 3: Некоторые наблюдения были отброшены, поскольку отсутствовали соответствующие ковариаты. 4: 12 сайты были удалены из-за отсутствия данных.

Я знаю, что я получил это из-за данных NA в моем исходном data.frame, и я могу вернуться и удалить эти строки, но мне пришлось бы переделать 25 разных моделей, и я бы предпочел преодолеть эту ошибку.

Есть ли какой-либо способ устранить ошибку в этой функции, или можно использовать другую функцию, чтобы получить хорошее соответствие?

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