INLAutils
INLAutils copied to clipboard
Nice way to plot mean and stdv of spatial projection
http://people.bath.ac.uk/fl353/isba/isbaspde.R
levelplot(row.values=proj$x,
column.values=proj$y,
x=inla.mesh.project(proj,
inla.result$summary.linear.predictor$mean[index[mesh.index$field.repl==1]]),
xlim=c(0,1),
ylim=c(0,1),
col.regions=cp, at=at+5,
aspect="iso",
contour=FALSE, labels=FALSE, pretty=TRUE,
xlab=NULL,ylab=NULL,scales=list(draw=FALSE))