This function creates an inheritance cube to model a Killer-Rescue system. Killer-Rescue is a 2-locus system: one locus has a toxin and the other locus contains the antidote. The loci are assumed independent and are non-homing.
This drive has 3 alleles at locus 1 and 2 alleles and locus 2:
Locus 1
T: Wild-type allele
K: "Killer" toxin allele
R: Broken toxin allele
Locus 2
W: Wild-type allele
A: Antidote allele
cubeKillerRescue( eR =0, Keff =1, Aeff =1, eta =NULL, phi =NULL, omega =NULL, xiF =NULL, xiM =NULL, s =NULL)
Arguments
eR: Conversion of K allele to R allele, a basal mutation rate
Keff: Toxin efficacy
Aeff: Antidote efficacy
eta: Genotype-specific mating fitness
phi: Genotype-specific sex ratio at emergence
omega: Genotype-specific multiplicative modifier of adult mortality
xiF: Genotype-specific female pupatory success
xiM: Genotype-specific male pupatory success
s: Genotype-specific fractional reduction(increase) in fertility
Returns
Named list containing the inheritance cube, transition matrix, genotypes, wild-type allele, and all genotype-specific parameters.