## facility functions ####
Gmat <- function(nvec){
D <- length(nvec)
arealab <- rep(1:D, times = nvec)
unname(model.matrix(~as.factor(arealab)-1))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.