rongelapUTM {glmmBUGS} | R Documentation |
A SpatialPointsDataFrame
containing the Rongelap data, in a UTM projection.
data("rongelapUTM") data("rongelapResult")
These coordinates were obtained by translating and rotating the original Rongelap data until all the coordinates fit into the Rongelap border given by www.gadm.org. So they are not exact.
See the help file for
rongelap
in the geoRglm package, or
http://www.leg.ufpr.br/doku.php/pessoais:paulojus:mbgbook:datasets
data("rongelapUTM") if(require("sp", quietly=TRUE)){ plot(rongelapUTM) } ## Not run: rongelapBorderLL = raster::getData("GADM", country="MHL",level=0) library("rgdal") rongelapBorderUTM = spTransform(rongelapBorderLL, CRS(proj4string(rongelapUTM))) plot(rongelapBorderUTM, add=TRUE) ## End(Not run) rongelapUTM$logOffset = log(rongelapUTM$time) rongelapUTM$site = seq(1, length(rongelapUTM$time)) forBugs = glmmBUGS( formula=count + logOffset ~ 1, family="poisson", data=rongelapUTM@data, effects="site", spatial=rongelapUTM, priors=list(phisite="dgamma(100,1)") ) startingValues = forBugs$startingValues startingValues$phi=list(site = 100) source("getInits.R") ## Not run: rongelapResult = bugs(forBugs$ragged, getInits, parameters.to.save = names(getInits()), model.file="model.bug", n.chain=2, n.iter=20, n.burnin=4, n.thin=2, program="winbugs", debug=TRUE, working.directory=getwd()) rongelapParams = restoreParams(rongelapResult, forBugs$ragged) ## End(Not run) data("rongelapResult") rongelapParams = restoreParams(rongelapResult) checkChain(rongelapParams) rongelapSummary=summaryChain(rongelapParams)