rongelapUTM: Rongelap island data

In glmmBUGS: Generalised Linear Mixed Models with BUGS and JAGS

Description

A SpatialPointsDataFrame containing the Rongelap data, in a UTM projection.

Usage

data("rongelapUTM")
data("rongelapResult")

Details

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.

Source

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

Examples

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)

glmmBUGS documentation built on May 2, 2018, 1:03 a.m.