我使用包ngspatial(Hughes)创建了一个空间泊松模型。我的数据汇总到包含变量violent
,g5prop
和conc_dis
的人口普查区域。我想适应异质性随机效应(即每个人口普查区的随机效应)。这是如何完成的?
这是我对我的模型所拥有的,它的工作非常好:
library(ngspatial)
Z = stl$violent
X = cbind(1, stl$g5prop, stl$conc_dis)
load("A.Rda") # adjacency matrix
set.seed(123456)
fit = sparse.sglmm(Z ~ X - 1 + offset(log(stl$total_pop)),
random = 1|stl$id, family=poisson, A=A,
attractive=40,tune=list(sigma.s = 0.02),
verbose = TRUE)
summary(fit)
现在,如何指定随机效果?它不像lme4
(例如(1|subject)
)或其他公式(random = 1~subject
)那样起作用。
谢谢!