mastif is used to estimate fecundity and seed dispersal based on mapped tree and seed trap data.
mastif operates on inventory plot data or simulated data. Inventory data for seeds of the genus Pinus are shown below for random plots and years. In these maps, large green symbols indicate seed traps that accumulated large numbers of seeds. Large brown symbols indicate large trees.
The simulator produces data sets with known parameter values and latent variables. Analysis provides insight on identifiability of seed production and the variables that control it.
Prediction validates the model and shows how well parameters and latent states are identified for specific data sets.
Maps of predicted seed rain and uncertainty from the fitted model. Predicted mean seed rain at left and predictive coefficient of variation at right. In-sample prediction includes mapped stands fitted with the model. Out-of-sample prediction is available for mapped stands that are not fitted to the data.
Year effects show ‘masting’, the tendency for individuals to synchronize reproduction. Year effects can be fitted by species and by species and region.
mastifuses Gibbs sampling with direct sampling, Metropolis, and Hamiltonian updating. It is implemented in R and C++ with
Rcpp and the
RcppArmadillo library for fast, efficient linear algebra operations. For this development version you will need to install them. Occasionally this is challenging on some platforms–if you have trouble send me an email. Someone in the lab has probably figured it out for a machine like yours. This will not be an issue once it is on CRAN.
Installation in RStudio:
> install.packages('mastif_1.0.tar.gz',repos=NULL, type='source') > library('mast')
> help('mastif') > browseVignettes('mastif')