也就是说,使用spdep包中nb类的region.id,而不是像spdep::nb2INLA那样忽略它?我一直试图将包含这些区域的数据中的一列作为一个因素链接到INLA模型,并提供描述它们的空间排列的图表。#something like this model = "bym2",
graph = inla_graphs$gb_regio
我正在使用R-INLA运行以下模型(Treatment,Animal.1和Animal.2是因子,Encounter.Length是连续的): formula <- Encounter.LengthTreatment +f(Animal.1, model = "iid", n = n.animal) +
m.1 <- inla(formula, data = inla