Return an x and y 3-D array for estimation of a traditional time-variant Cormack-Jolly-Seber capture-recapture model.
Number of individuals/animals.
Number of trap/mark occasions.
\$x are designed to useful
for fitting classical CJS models with time-variant, but individual-invariant effects.
To fit a CJS model using this function, the commands would be
tmp<-F.cjs.covars(nan,ns);F.cjs.estim(capture=~tmp\$x[,,2]+tmp\$x[,,3]+ ..., survival=
~tmp\$x[,,1]+tmp\$x[,,2]+ ..., histories=my.histories)
A list containing a single component,
\$x, that can be used to
estimate a classical CJS model when included in a subsequent call to
F.cjs.estim. The returned component,
is a 3-D array containing 0's everywhere, except for 1's in certain columns.
\$x has dimension
Element [i,j,k] of
\$x is 1 if j == k, and 0 otherwise. I.e., the k-th "page" of
the 3-D array has 1's in the k-th column, 0's elsewhere.
Trent McDonald, WEST Inc., [email protected]
1 2 3 4 5 6
## Fit CJS model to dipper data, time-varying capture and survivals. data(dipper.histories) xy <- F.cjs.covars( nrow(dipper.histories), ncol(dipper.histories) ) dipper.cjs <- F.cjs.estim( capture=~xy$x[,,2]+xy$x[,,3]+xy$x[,,4]+xy$x[,,5]+xy$x[,,6], survival=~xy$x[,,1]+xy$x[,,2]+xy$x[,,3]+xy$x[,,4]+xy$x[,,5], dipper.histories ) print(dipper.cjs)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.