Пространственное моделирование результата бинарного ответа - PullRequest
Новая
       112

Пространственное моделирование результата бинарного ответа

0 голосов
/ 06 августа 2020

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

Yik следует распределению Бернулли для дочернего k, находящегося в состоянии i = 1, ..., 37 с ответом либо 0, либо 1. Мы добавлены структурированные и неструктурированные случайные эффекты для учета ненаблюдаемых влияющих факторов, которые варьируются в зависимости от штата. Я запустил свой код; 

Formula<- yik~x1+...+xn+f(NAME_1a,model="besag",graph="Map.adj")+
+   f(NAME_1b,model="iid")
> Analysis<-inla(Formula
+                ,family="binomial"
+                ,data=dat
+                ,control.compute=list(dic=TRUE, cpo=TRUE)),

, но получил это сообщение; 

Error in inla(Formula, family = "binomial", data = dat, control.compute = list(dic = TRUE,  : 
In f(NAME_1a): 'covariate' must match 'values',  and both must either be 'numeric', or 'factor'/'character'.

Чтобы исправить ошибку, я также запускаю следующий код: 

m = get("inla.models", INLA:::inla.get.inlaEnv())
m$latent$rw2$min.diff <- NULL
assign("inla.models", m, INLA:::inla.get.inlaEnv())

Но вот этот, m$latent$rw2$min.diff <- NULL, выдал эту ошибку; 

"Error in `*tmp*`$latent :`  object of type 'closure' is not subsettable"
...