find_mate: Finds a mate for a female mosquitoes of breeding age.

In beeysian/cairnsmozzie: Agent-based model of Wolbachia release in a nice tidy package

Description

We draw a box of 'radius' k around the agent. k is an ABC parameter and usually between 11m to 100m. Note that here we mainly deal with ID number, an actual data.table variable, and not the agent's index in the main data.table.

Usage

find_mate(femID, k, max_daily_mates, mozzie.dt)

Arguments

femID	ID of female looking for mate.
k	Distance within which female mosquito 'searches' for mate. ABC parameter.
max_daily_mates	Number of times a male can mate in a day.
mozzie.dt	Data.table of adult mosquitoes.

Value

mateID of agent, or -1 if agent does not find a mate.

Algorithm

Algorithm for finding a mate:

• If there are no males within box, return -1 for no mate.

• If there is 1 male, that is the new mate.

• If there is >1 male, randomly choose one.

Note that there is a max. number of agents that a male can mate with in a day, given by max_daily_mates. It is important in testing to ensure that the gonoCycle variable for each male mate is properly updated on a global scope. We assume that females only mate once.

Who can mate?

Recall that in order to mate, a female mosquito must have:

• enzyme between 1.2 and 1.8 inclusive

• gonoCycle 1, 2 or 3

beeysian/cairnsmozzie documentation built on Feb. 15, 2021, 12:12 a.m.