MARSS()
outputs marssMLE
objects. coef(object)
, where object
is the output from a MARSS()
call, will print out the estimated parameters. The default output is a list with values for each parameter, however the output can be altered using the type
argument to output a vector of all the estimated values (type="vector"
) or a list with the full parameter matrix with the estimated and fixed elements (type="matrix"
). For a summary of the parameter estimates with CIs from the estimated Hessian, use try tidy(object)
.
# S3 method for marssMLE
coef(object, ..., type = "list", form = NULL, what = "par")
A list of the estimated parameters for each model matrix.
A marssMLE
object.
Other arguments. Not used.
What to output. Default is "list". Options are
A list of only the estimated values in each matrix. Each model matrix has it's own list element.
A vector of all the estimated values in each matrix.
A list of the parameter matrices each parameter with fixed values at their fixed values and the estimated values at their estimated values. Time-varying parameters, including d and c in a marxss form model, are returned as an array with time in the 3rd dimension.
Returns the parameter matrix for that parameter with fixed values at their fixed values and the estimated values at their estimated values. Note, time-varying parameters, including d and c in a marxss form model, are returned as an array with time in the 3rd dimension.
This argument can be ignored. By default, the model form specified in the call to MARSS()
is used to determine how to display the coefficients. This information is in attr(object$model,"form")
. The default form is "marxss"
; see MARSS.marxss()
. However, the internal functions convert this to form "marss"
; see MARSS.marss()
. The marss form of the model is stored (in object$marss
). You can look at the coefficients in marss form by passing in form="marss"
.
By default, coef()
shows the parameter estimates. Other options are "par.se", "par.lowCI", "par.upCI", "par.bias", and "start".
Eli Holmes, NOAA, Seattle, USA.
tidy()
, print()
dat <- t(harborSeal)
dat <- dat[c(2, 11), ]
fit <- MARSS(dat)
coef(fit)
coef(fit, type = "vector")
coef(fit, type = "matrix")
# to retrieve just the Q matrix
coef(fit, type = "matrix")$Q
Run the code above in your browser using DataLab