IPModel: Estimates the pattern of internesting intervals for a set of...

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

Estimates the pattern of internesting intervals for a set of parameters.

Description

This function fits a model of internesting period.
The parameters are:

• meanIP : The average number of days between two nesting processes

• DeltameanIP : The shift in days for IP at each new clutch.

• sdIP : The standard deviation of number of days between two nesting processes

• minIP : The minimum number of days between two nesting processes

• pAbort : The -logit of the probability to abort a nesting process

• meanAbort : The average of the number of days after the abortion of a nesting process

• sdAbort : The standard deviation of the number of days after the abortion of a nesting process

• pCapture : The -logit of the probability to capture a female on the beach

• meanECF : The average number of clutch a female will try to do being reprensented as ECF

• sdECF : The standard deviation of number of clutch a female will try to do

• N : The number of replicates to generate the distribution (default is 10000 if not indicated)

• ECF.x : The relative proportion of females nesting with ECF = x (ECF.1 being fixed to 1)

Usage

IPModel(
  par,
  parallel = TRUE,
  limits = list(meanIP = 40, meanECF = 4, minIP = 15, sdAbort = 1, sdIP = 1, sdECF = 1,
    DeltameanIP = 0.5, maxDays = 365)
)

Arguments

par	Set of parameters
parallel	If TRUE, will use parallel computing
limits	A list of limits for various parameters

Details

IPModel estimates the pattern of internesting intervals for a set of parameters.

Value

Return a list with two elements.

Author(s)

Marc Girondot [e-mail]

See Also

Other Model of Internesting Period: IPFit(), IPPredict(), plot.IP(), summary.IP()

Examples

## Not run: 
library(phenology)
# Example
par <- c(meanIP = 9.8, 
sdIP = 0.1, 
minIP = 7, 

pAbort = -logit(0.1), 
meanAbort = 2, 
sdAbort = 0.05, 

pCapture = -logit(0.8), 

meanECF = 4, 
sdECF = 0.1)

model <- IPModel(c(par, N=10000))

plot(model)


## End(Not run)

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