Recent Releases of mvgam
mvgam - v1.1.5
This release brings better plotting for several plot() methods, along with more useful guidance in summary() for diagnostics procedures. Default priors have also changed for variance parameters to provide better out of the box regularisation
- R
Published by nicholasjclark 11 months ago
mvgam - v1.1.4
New functionalities
- Added the
how_to_cite.mvgam()function to generate a scaffold methods description of fitted models, which can hopefully make it easier for users to fully describe their programming environment - Improved various plotting functions by returning
ggplotobjects in place of base plots (thanks to @mhollanders #38) - Added the brier score (
score = 'brier') as an option inscore.mvgam_forecast()for scoring forecasts of binary variables when usingfamily = bernoulli()(#80) - Added
augment()function to add residuals and fitted values to an mvgam object's observed data (thanks to @swpease #83) - Added support for approximate
gp()effects with more than one covariate and with different kernel functions (#79) - Added function
jsdgam()to estimate Joint Species Distribution Models in which both the latent factors and the observation model components can include any of mvgam's complex linear predictor effects. Also added a functionresidual_cor()to compute residual correlation, covariance and precision matrices fromjsdgammodels. See?mvgam::jsdgamand?mvgam::residual_corfor details - Added a
stability.mvgam()method to compute stability metrics from models fit with Vector Autoregressive dynamics (#21 and #76) - Added functionality to estimate hierarchical error correlations when using multivariate latent process models and when the data are nested among levels of a relevant grouping factor (#75); see
?mvgam::ARfor an example - Added
ZMVN()error models for estimating Zero-Mean Multivariate Normal errors; convenient for working with non time-series data where latent residuals are expected to be correlated (such as when fitting Joint Species Distribution Models); see?mvgam::ZMVNfor examples - Added a
fevd.mvgam()method to compute forecast error variance decompositions from models fit with Vector Autoregressive dynamics (#21 and #76)
Deprecations
- Arguments
use_stan,jags_path,data_train,data_test,adapt_delta,max_treedepthanddrifthave been removed from primary functions to streamline documentation and reflect the package's mission to deprecate 'JAGS' as a suitable backend. Bothadapt_deltaandmax_treedepthshould now be supplied in a namedlist()to the new argumentcontrol
Bug fixes
- Bug fix to ensure
marginaleffects::comparisonsfunctions appropriately recognise internalrowidvariables - Updates to ensure
ensembleprovides appropriate weighting of forecast draws (#98) - Not necessarily a "bug fix", but this update removes several dependencies to lighten installation and improve efficiency of the workflow (#93)
- Fixed a minor bug in the way
trend_maprecognises levels of theseriesfactor - Bug fix to ensure
lfo_cvrecognises the actual times intime, just in case the user supplies data that doesn't start att = 1. Also updated documentation to better reflect this - Bug fix to ensure
update.mvgamcaptures anyknotsortrend_knotsarguments that were passed to the original model call
- R
Published by nicholasjclark about 1 year ago
mvgam - v1.1.3
mvgam 1.1.3
New functionalities
- Allow intercepts to be included in process models when
trend_formulais supplied. This breaks the assumption that the process has to be zero-centred, adding more modelling flexibility but also potentially inducing nonidentifiabilities with respect to any observation model intercepts. Thoughtful priors are a must for these models - Added
standata.mvgam_prefit,stancode.mvgamandstancode.mvgam_prefitmethods for better alignment with 'brms' workflows - Added 'gratia' to Enhancements to allow popular methods such as
draw()to be used for 'mvgam' models if 'gratia' is already installed - Added an
ensemble.mvgam_forecastmethod to generate evenly weighted combinations of probabilistic forecast distributions - Added an
irf.mvgammethod to compute Generalized and Orthogonalized Impulse Response Functions (IRFs) from models fit with Vector Autoregressive dynamics
Deprecations
- The
driftargument has been deprecated. It is now recommended for users to include parametric fixed effects of "time" in their respective GAM formulae to capture any expected drift effects
Bug fixes
- Added a new check to ensure that exception messages are only suppressed by the
silentargument if the user's version of 'cmdstanr' is adequate - Updated dependency for 'brms' to version >= '2.21.0' so that
read_csv_as_stanfitcan be imported, which should future-proof the conversion of 'cmdstanr' models tostanfitobjects (#70)
- R
Published by nicholasjclark over 1 year ago
mvgam - v1.1.2
This version brings several new features and efficiency improvements
* Added options for silencing some of the 'Stan' compiler and modeling messages using the silent argument in mvgam()
* Moved a number of packages from 'Depends' to 'Imports' for simpler package loading and fewer potential masking conflicts
* Improved efficiency of the model initialisation by tweaking parameters of the underlying 'mgcv' gam object's convergence criteria, resulting in much faster model setups
* Added an option to use trend_model = 'None' in State-Space models, increasing flexibility by ensuring the process error evolves as white noise (#51)
* Added an option to use the non-centred parameterisation for some autoregressive trend models,
which speeds up mixing most of the time
* Updated support for multithreading so that all observation families (apart from nmix()) can now be modeled with multiple threads
* Changed default priors on autoregressive coefficients (AR1, AR2, AR3) to enforce
stationarity, which is a much more sensible prior in the majority of contexts
* Fixed a small bug that prevented conditional_effects.mvgam() from handling effects with three-way interactions
- R
Published by nicholasjclark over 1 year ago
mvgam - v1.08
This release brings piecewise linear and logistic trends with automatic changepoint selection, similar to what is available in Facebook's popular prophet package. It also includes a few bug fixes to make gp() effects more stable and versatile
- R
Published by nicholasjclark about 2 years ago