Я пытался запустить модель CAR (пакет CARBayes - функция S.CARbym в R) для прогнозирования риска заболеваний висцеральным лейшманиозом в штате Сан-Паулу, Бразилия. Ковариатами являются осадки, высота, температура и расстояние до автомагистралей. Когда я запускаю модель с всеми необработанными ковариатами , возникает ошибка: "Error in ar.yw.default(x, aic = aic, order.max = order.max, na.action = na.action, : zero-variance series"
Я пытался настроить значения ковариат с аналогичными диапазонами на основе
(наблюдаемые данные - среднее значение) / (стандартное отклонение наблюдаемых данных)
equation=LV~offset(log(Total_Pop))+precipitation+altitude+temperature+distance
model <- S.CARbym(formula=equation, family="zip", formula.omega=~1, data=spat.data, W=W2, burnin=100000, n.sample=300000, thin=20)
Тогда я мог бы запустить модель, но результаты значений DIC = NaN.
Мои данные (Y) содержат много "нулей", поэтому семейство zip кажется наиболее подходящим. Всего данных 643 функции. Я использовал население в качестве смещения.
Индивидуально, для каждого ковариата, это работает. Но я хотел бы запустить модель со всеми из них.
Как я могу это исправить?