如何在包ngspatial中使用sparse.sglmm指定异构随机效应?

时间:2014-11-18 22:13:15

标签: r models spatial poisson

我使用包ngspatial(Hughes)创建了一个空间泊松模型。我的数据汇总到包含变量violentg5propconc_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)那样起作用。

谢谢!

0 个答案:

没有答案