multistate - Stata development version history

Last updated on Aug 13, 2019

This page details the version history of the development version of multistate, hosted on this website, which can be installed by typing in Stata,

net install multistate, from(https://www.mjcrowther.co.uk/code/multistate)

or updated by typing adoupdate multistate.

bug fix: missed an edit in 2.3.1 fix preventing aj to work in predictms (merlin function names); now fixed
error check added for streg, d(ggamma) combined with aj - not supported

bug fix: edits in merlin function names (1.2.5) caused predictms to fail; now fixed
error check added for bhtime with strcs models - not currently supported

bug fix: stpm2 and strcs models failed when simulating, introduced in 2.2.0; now fixed
streg, dist(lognormal) now allowed with clock-forward models
bug fix: streg, dist(lognormal) with reset incorrectly assumed loglogistic; now fixed

bug fix: with novcv() when ci not specified introduced in 2.1.0; now fixed
bug in error check for log normal requiring reset; now fixed
Generalised gamma streg model now supported when used in the models() syntax
bug fix in msaj when number of transitions was not equal to dimension of transmatrix()

novcv() option added to use mean vector instead of draws in CI calculations, i.e. assumes a particular transition is estimated without uncertainty

bug fix: los was incorrectly scaled to 1, introduced in 1.2.0. Now fixed.
normal approximation for CIs is now the default. percentile option replaces normal.
bug fix: level() was ignored in CI calculations, always assumed 95%. Now fixed.
gen() removed
bug fix: in some cases mm_root_vec produced an error when solving the root. Now fixed.
all predictions now available in one call
infinite at#()s added (limited to 50 for error check reasons)
at() changed to at1()
difference option added
added atref(#) - default 1
standardise option added to calculate standardised (sample-averaged) predictions
bug fix: with Stata 15 and streg (apart from dist(exp)), with at2(), the ancilary parameter constant was ignored. Now fixed.
_time = 0 now allowed with normal ci’s
userfunction() added for user-defined prediction function, subroutines for probs and los
userlink() added for normal cis of userfunction() - default identity, can also be log or logit
cr option added to avoid specifying a transition matrix. For use with models() only.
bug fix: reset now synced with root-finding updates for stpm2 and strcs models
outsample added
stms added
enter() removed, now min(timevar or _t)
aj added for Aalen-Johansen estimator with Markov models
reversible transitions now allowed

some re-writes of source code for massive speed improvements when root-finding required
trans#() syntax removed, now parsing is done automatically and using covariates expanded using msset
model#() changed to models()
improvements to help files
covariates() added to msset to create transition-specific dummies, which must be used in model fitting (if required)
default n() now 100,000, or 10,000 with ci
survival added to calculate all predictions on a standard single event model