o_4p_p1p2: Model of tagloss based on Rivalan data

In phenology: Tools to Manage a Parametric Function that Describes Phenology and More

Model of tagloss based on Rivalan data

Description

Model of tagloss based on Rivalan data

Usage

o_4p_p1p2

Details

Model of tagloss based on Rivalan data

Author(s)

Marc Girondot

References

Rivalan, P., Godfrey, M.H., Prévot-Julliard, A.-C., Girondot, M., 2005. Maximum likelihood estimates of tag loss in leatherback sea turtles. Journal of Wildlife Management 69, 540-548.

See Also

Other Model of Tag-loss: Tagloss_LengthObs(), Tagloss_L(), Tagloss_cumul(), Tagloss_daymax(), Tagloss_fit(), Tagloss_format(), Tagloss_mcmc_p(), Tagloss_mcmc(), Tagloss_model(), Tagloss_simulate(), logLik.Tagloss(), plot.TaglossData(), plot.Tagloss()

Examples

## Not run: 
library(phenology)
# Example
data_f_21 <- Tagloss_format(outLR, model="21")

# Same data fitted with new model
par <- c(D1_1 = 100.15324837975547, A_1 = 5.9576927964120188, 
         B_1 = 8.769924225871069, B_2 = 8.2353860179664125)
pfixed <- c(D2D1_1 = 2568, D3D2_1 = 2568, D2D1_2 = 2568, D3D2_2 = 2568)
o_4p_p1p2 <- Tagloss_fit(data=data_f_21, fitted.parameters = par, 
                         fixed.parameters = pfixed, 
                         model_before = "par['C_1']=par['B_1'];
                         par['A_2']=par['A_1'];
                         par['C_2']=par['B_2'];
                         par['D1_2']=par['D1_1']", hessian=TRUE)
data(o_4p_p1p2)
plot(o_4p_p1p2, model="1", col="red")
plot(o_4p_p1p2, model="2", col="blue", add=TRUE)
legend("topright", legend=c("2->1", "1->0"), lty=1, col=c("blue", "red"))
plot(o_4p_p1p2, model="Cumul")

## End(Not run)

phenology documentation built on Oct. 16, 2023, 9:06 a.m.