ECFOCF_f: Calculate a table of probabilities of ECF and OCF.

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/ECFOCF_f.R

Description

This function calculates a table of probabilities of ECF and OCF.
If p is lower or higher than 1E-100 or 1-1E-100, it is changed to 1E-100 and 1-(1E-100) respectively.
Names for p vector elements should be p, or px (with x=1:categories), or px.period.
If mu_season and sd_season are equal to NA, the model is not temporalized.
If mu_season and sd_season are not NA, the model is temporalized.

Usage

1
2
3
ECFOCF_f(mu, sd = NA, p, MaxNests = 15, mu_season = NA, sd_season = NA,
  MeanDaysBetween2Nests = 9.8,
  length_season = floor(365/MeanDaysBetween2Nests) + 1, parallel = TRUE)

Arguments

mu

The average of lognormal for clutch frequency.

sd

The sd parameter of lognormal for clutch frequency.

p

The capture probability for an individual nesting event.

MaxNests

Maximum number of nests by a female.

mu_season

The average of ordinal day for beginning of nesting season.

sd_season

The sd parameter of lognormal for ordinal day for beginning of nesting season.

MeanDaysBetween2Nests

Average number of days between two nests.

length_season

The total length of season based on groups of interclutch intervals.

parallel

If TRUE parallel computing is used.

Details

ECFOCF_f calculate a table of probabilities of ECF and OCF.

Value

Return a matrix of class TableECFOCF.

Author(s)

Marc Girondot

See Also

Other Model of Clutch Frequency: TableECFOCF, fitCF, lnLCF, logLik.ECFOCF, plot.ECFOCF, plot.TableECFOCF

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
## Not run: 
library(phenology)
# Example
modelECFOCF <- ECFOCF_f(mu=5.58013243236187, 
                    sd=1.225581130238, 
                    p=invlogit(1.3578137414575), 
                    MaxNests=15)
plot(modelECFOCF)
modelECFOCF <- ECFOCF_f(mu=5.58013243236187, 
                    sd=1.225581130238, 
                    mu_season=12, 
                    sd_season=2, 
                    p=c(p1=invlogit(1.3578137414575)), 
                    MaxNests=15, 
                    MeanDaysBetween2Nests=9.8, 
                    length_season=floor(365/9.8)+1
                    )
plot(modelECFOCF, period=2)

## End(Not run)

phenology documentation built on July 10, 2018, 5:07 p.m.