mra (version 2.16.4)

F.update.df: 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" 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

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 )