Ошибка в validObject(.Object): недопустимый объект класса «dsparseModelMatrix»: суперкласс «Mnumeric»
Я использую байесовскую модель с INLA. Несколько месяцев назад все работало нормально. Сегодня я снова запускаю эту ошибку: Ошибка в validObject(.Object): недопустимый объект класса «dsparseModelMatrix»: суперкласс «Mnumeric» не определен в среде класса объекта
Мой код:
formula = casesp ~
f(ID.area, model="bym",
graph=g,
adjust.for.con.comp=TRUE) +
f(ID.area1,
ID.time,
model="iid",
constr=TRUE) +
ID.time
res0 = inla(formula, family="poisson", data=datos1, E=NULL,
control.predictor=list(compute=TRUE),
control.compute=list(dic=TRUE),
quantiles=c(0.025,0.975), verbose = FALSE)
str данных:
'data.frame': 4024 obs. of 34 variables:
$ DPA_PARROQ : chr "020150" "020150" "020150" "020150" ...
$ populationa : num 13768 14382 13987 16095 7 ...
$ casesa : num 3 6 0 0 0 0 0 0 0 0 ...
$ populationp : num 4608 4140 3933 3968 4 ...
$ casesp : num 2 2 0 0 0 0 0 0 0 0 ...
$ year : chr "2017" "2018" "2019" "2020" ...
$ vaccinated_pigs : num 19558 20560 24259 19376 0 ...
$ vacinators : num 87 67 79 66 0 2 3 1 7 6 ...
$ pop_a_na : num 13768 14382 13987 16095 7 ...
$ pop_p_na : num 4608 4140 3933 3968 4 ...
$ pop_a_vac_cor : num 13768 14382 13987 16095 7 ...
$ percentage_vac_cor : num 0.71 0.71 0.87 0.6 0 0.48 0.5 0.5 0.46 0.82 ...
$ pop_a_s : num [1:4024, 1] 3.16 3.313 3.215 3.742 -0.285 ...
$ pop_p_s : num [1:4024, 1] 16.326 14.627 13.875 14.002 -0.394 ...
$ vac_pigs_s : num [1:4024, 1] 1.269 1.343 1.616 1.255 -0.177 ...
$ vacin_s : num [1:4024, 1] 12.236 9.262 11.046 9.113 -0.704 ...
$ percentage_vac_cor_s: num [1:4024, 1] 0.2675 0.2675 0.7422 -0.0589 -1.839 ...
$ exp_pop_p_na : num 1.76394 1.58479 1.50555 1.51895 0.00153 ...
$ SIR_p_na : num 1.13 1.26 0 0 0 ...
$ id : int 1 1 1 1 2 2 2 2 3 3 ...
$ DPA_DESPRO : Factor w/ 24 levels "AZUAY","BOLIVAR",..: 2 2 2 2 2 2 2 2 2 2 ...
$ DPA_DESPAR : Factor w/ 928 levels "10 DE AGOSTO",..: 321 321 321 321 275 275 275 275 371 371 ...
$ timeformatted : chr "2017" "2018" "2019" "2020" ...
$ time : Date, format: "2017-01-01" ...
$ idtime : int 1 2 3 4 1 2 3 4 1 2 ...
$ Area : num 519 519 519 519 162 ...
$ vac_sqr_km : num 37.7 39.6 46.7 37.3 0 ...
$ vac_sqr_km_s : num [1:4024, 1] 0.171 0.195 0.282 0.167 -0.291 ...
..- attr(*, "scaled:center")= num 23.7
..- attr(*, "scaled:scale")= num 81.6
$ id2 : num 0 0 0 0 1 1 1 1 2 2 ...
$ regiao : chr "1Highlands" "1Highlands" "1Highlands" "1Highlands" ...
$ ID.area : int 1 1 1 1 2 2 2 2 3 3 ...
$ ID.area1 : int 1 1 1 1 2 2 2 2 3 3 ...
$ ID.time : int 1 2 3 4 1 2 3 4 1 2 ...
$ ID.time1 : int 1 2 3 4 1 2 3 4 1 2 ...
Работает на x86_64-pc-linux-gnu,
arch x86_64
os linux-gnu
R версия 3.6.3
год 2020
'INLA' версия 21.02.23
Любая помощь!