Displaying 1 result from an estimated 1 matches for "popdataadjmat".
2011 Jun 27
0
OpenBUGS and glmmBUGS package
...om the diseasemapping
package document:
----------------------------------------
library(diseasemapping)
data(popdata)
data(casedata)
model = getRates(casedata, popdata, ~age*sex, breaks=seq(0, 90, by=10) )
ontario = getSMR(popdata,model, casedata)
spplot(ontario, 'SMR')
library(spdep)
popDataAdjMat = poly2nb(ontario, ontario[["CSDUID"]])
library(glmmBUGS)
forBugs = glmmBUGS(formula=observed + logExpected ~ 1,
effects="CSDUID", family="poisson", spatial=popDataAdjMat,
data=ontario at data, modelFile="model.txt")
startingValues = forBugs$startingValues...