Description Usage Details Source Examples
A SpatialPointsDataFrame
containing the Rongelap data, in a UTM projection.
1 2 |
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 \Sexpr[results=rd]{c( '\\\code{rongelap} in the \\\CRANpkg{geoRglm} package', '\\\code{\\\\link[geoRglm]{rongelap}}' )[1+requireNamespace('geoRglm', quietly=TRUE)]}, or http://www.leg.ufpr.br/doku.php/pessoais:paulojus:mbgbook:datasets
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.