predictEN: Expectation of the Number of Ticks on a Roe Deer as a...

In ClementCalenge/tickTF: Infection of Roe Deer Fawns by Ticks from 1992 to 2018 in the Trois-Fontaines Forest

Description

This function calculates the expected number of ticks on an average-sized roe deer fawn of a given age in average environmental conditions, at various ages, as predicted by an infection model.

Usage

predictEN(mcmcList, TicksData, TicksConst)

Arguments

mcmcList	An object of class mcmc.list containing the results of the fit of a model.
TicksData	The list containing the data passed to nimbleMCMC. WARNING: the function does not check that this list actually corresponds to the list passed to nimbleMCMC and trusts the user.
TicksConst	The list containing the constants passed to nimbleMCMC. WARNING: the function does not check that this list actually corresponds to the list passed to nimbleMCMC and trusts the user.

Value

A list containing the following elements:

age	a numeric vector giving the ages, from 1 to 8 days, at which the expected number is calculated
EN	a matrix with R rows (corresponding to the R MCMC iterations) and N columns (corresponding to the N values of age for which the expected number is calculated) giving the expected number of ticks.

Author(s)

Clement Calenge [e-mail]

References

Bariod L., Said S., Calenge C., Chabot S., Badeaud V. \& Bourgoin G. in prep. Infection of Roe Deer Fawns by Ticks from 1992 to 2018 in the Trois-Fontaines Forest (France).

Examples

## Not run: 
    en <- predictEN(resultsModel$samples, TicksData, TicksConst)

    lam <- en$EN
    lr <- range(unlist(lam))
    age <- en$age
    plot(age, lam[1,], ty="n", ylim=lr, xlab="Age (days)", ylab="Mean number of ticks")
    tmp <- lapply(1:nrow(lam), function(i) lines(age, lam[i,], col=rgb(0.1,0.1,0.1, 0.02)))
    mo <- apply(lam,2,mean)
    lines(age, mo, lwd=5)
    lines(age, mo, lwd=3, col="yellow")


## End(Not run)

ClementCalenge/tickTF documentation built on Dec. 17, 2021, 2:06 p.m.