F.update.df {mra}

F.update.df - Update degrees of freedom in a Cormack-Jolly-Seber fitted object

Description

Updates the degrees of freedom in a fitted object to either the value estimated from the rank of the variance-covariance matrix, the number of coefficients, or a user-specified value.

Usage

F.update.df(fit, df=NA)

Arguments

fit

An MRA fitted CJS model. Class must be c("cjs", "cr"). These are produced by F.cjs.estim.

df

The new value for degrees of freedom. If df == NA, the number of parameters estimated from the rank of the matrix of 2nd derivatives or Hessian, depending on cov.meth parameter. If df <= 0, the number of parameters will be set to NX+NY = the number of estimated coefficients. Otherwise, if df > 0, the supplied value is used. The penalty terms of AIC, QAIC, AICc, and QAICc are recomputed using this value.

Value

An object (list) of class c("cjs","cr") with degrees of freedom, AIC, QAIC, AICc, and QAICc updated.

Author(s)

Trent McDonald, WEST-INC, [e-mail]

See Also

F.cjs.estim

Examples

## Fit CJS model to dipper data, time-varying capture and survivals.
data(dipper.histories)
ct <- as.factor( paste("T",1:ncol(dipper.histories), sep=""))
attr(ct,"nan")<-nrow(dipper.histories)
dipper.cjs <- F.cjs.estim( ~tvar(ct,drop=c(1,2)), ~tvar(ct,drop=c(1,6,7)), dipper.histories )

## Update the degrees of freedom
dipper.cjs <- F.update.df( dipper.cjs, -1 )