Two dimensional map from a Bayesian hierarchical model

This function defines the model for the spatial behaviour of the GEV parameter.

1 |

`data` |
A matrix representing the data. Each column corresponds to one location. |

`formula` |
A R formula. See details for further details. |

**need to be documented**

Mathieu Ribatet

1 2 3 4 5 6 7 8 9 10 11 12 | ```
## 1- A design matrix from a classical linear model
n.site <- 5
coord <- matrix(rnorm(2*n.site, sd = sqrt(.2)), ncol = 2)
colnames(coord) <- c("lon", "lat")
loc.form <- loc ~ lat + I(lon^2)
modeldef(coord, loc.form)
## 2- A design and penalization matrix from a penalized smoothin spline
x <- sort(runif(10, -2, 10))
n.knots <- 3
knots <- quantile(x, prob = 1:n.knots / (n.knots + 2))
modeldef(x, y ~ rb(x, knots = knots, degree = 3, penalty = 1))
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.