rongelapUTM {glmmBUGS}R Documentation

Rongelap island data

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 rongelap in the geoRglm package, 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)

[Package glmmBUGS version 2.4.2 Index]