#
# A file for testing the use of the forest gap model
#
initial <- init_landscape(states=c("+","0"),
cover=c(.5,.5))
parms <- forestgap$parms # Import defaults
parms$delta <- .15
output <- ca(initial, parms, model=forestgap, t_max=1000)
library(ggplot2)
output_fmt <- fortify(output)
ggplot(output_fmt) +
geom_raster(aes(x,y,fill=state)) +
facet_wrap( ~ time ) +
scale_fill_manual(values=c('#111111','#FFFFFF'))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.