mateable (version 0.3.2)

dailyOMP: Calculate daily outcrossed mating potential

Description

dailyOMP generates an OMP object giving the daily outcrossed mating potential of individuals based on k nearest neighbors

Usage

dailyOMP(
  scene,
  k = 3,
  days = NULL,
  gamma = 1/13.3,
  nn.constant = FALSE,
  sum = FALSE,
  mean = FALSE
)

Value

a named matrix with a row for each id and a column for each day, and entries corresponding to ids' OMP each day

Arguments

scene

a matingScene object

k

integer, number of nearest neighbors to use in calculating OMP

days

the day or range of days to calculate OMP for, default is all days in a scene (see Details)

gamma

parameter of exponential decay to be used in calculating OMP (defaults to 1/13)

nn.constant

logical; indicates whether the nearest neighbors used in calculations should be the nearest on a given day (nn.constant = FALSE) or the nearest neighbors over an entire season (nn.constant = TRUE)

sum

logical; indicates if the return should be a sum of an individual's daily outcrossed mating potential over the range specified by days

mean

logical; indicates if the return should be the mean of an individual's daily outcrossed mating potential over the range of days that the individual was receptive to mating

Author

Amy Waananen

Details

Daily outcrossed mating potential is a weighted average of an individual's distance to their nearest neighbors on a given day (Wagenius et al. 2007). The days to calculate OMP for should be input as integers relative to the first day of flowering, as they are in the start and end columns of a matingScene object. If the number of ids receptive on a day is less than k, OMP will be calculated for the maximum number of neighbors.

References

Wagenius, S., E. Lonsdorf, and C. Neuhauser. 2007. Patch aging and the S-Allee effect: breeding system effects on the demographic response of plants to habitat fragmentation. American Naturalist 169:383-397.

See Also

makeScene, proximity, synchrony, receptivityByDay

Examples

pop <- simulateScene()
omp <- dailyOMP(pop)
omp.1 <- dailyOMP(pop, nn.constant = TRUE) # same nearest neighbors throughout the season
omp.2 <- dailyOMP(pop, nn.constant = FALSE) # nearest flowering neighbors