Reads coefficients, standard errors, log-likelihoods, maximum gradients, correlation and variance-covariance matrices from AD Model Builder output files
1 2 3 4 5 6
(character) Base name of AD Model Builder
(logical) drop negative-phase (fixed) parameters from results?
(character) file name for covariance matrix information
warn if report file is in nonstandard format?
(character) Names of variables
Given the output from an ADMB run on FOO.tpl,
read_pars reads the
files FOO.par (parameters, log-likelihood, max gradient); FOO.std (standard
deviations); FOO.cor (correlations); FOO.rep (report variables);
admodel.hes for hessian; and admodel.cov for
read_psv reads the output of MCMC runs.
read_rep (called by
read_admb to read the
.rep file) first
checks if the report file is in a standard format: first line starts with a
comment character (#); thereafter, each block starts with a single commented
line containing the name of the parameter (possibly ending with a colon),
followed by a block of all-numeric lines, which are read as a single vector.
If the report file is in a standard format, the values are added to the end
of the coefficients list. Otherwise, the numeric values from the report file
are included in the results as a single, concantenated numeric vector.
List containing the following elements
coeflistparameter estimates in list format, with proper shape (vectors, matrices, etc.)
seestimated standard errors of coefficients
maxgradmaximum gradient of log-likelihood surface
nparnumber of parameters
heshessian matrix (only if no vcov matrix)
reportvalues from report file (if non-standard report file)
coeflist component is untested for data
structures more complicated than scalars, vectors or matrices (i.e. higher-dimensional or ragged arrays)
Because ADMB hard-codes the file name for covariance matrix information (
admodel.cov), care is necessary when running different models in the same directory; users may want to rename this file by hand and use the
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.