PROSPER (version 0.3.3)

sel_resist: Calculating resistance values

Description

Calculates the phenotypic resistance value (Renton et al. 2011) for each genotype in dfgenotype.

Usage

sel_resist(Rmx, epis, dom)

Arguments

Rmx

maximum resistance value, if all gene loci under consideration are homozygous resistant. \(\code{numeric}, \ge 1\).

epis

epistasis value, describing the interaction between resistance alleles. numeric. epis = 0: no interaction i.e. additive effects of resistance alleles, epis < 0: effect of resistance alleles is smaller than additive, epis > 0: effect of resistance alleles is higher than additive.

dom

dominance of resistance alleles. numeric vector with length of n_loci (number of genes) and elements in [0,1].

Value

numeric vector with resistance values

Details

This function is used in struc_preparation2 to calculate the resistance value using the following term: $$ 1 + (Rmx - 1) * (sum(dom)/n\_loci)^{2^{epis}}). $$ Resistance values range from 1 to Rmx. The higher the resistance value is, the less the plant is susceptible to the herbicide. sel_resist is only used by struc_preparation2.

References

Renton, M.; Diggle, A.; Manalil, S. & Powles, S. (2011): Does cutting herbicide rates threaten the sustainability of weed management in cropping systems? Journal of Theoretical Biology, 283, 14-27.

See Also

sel_herb intern_herbicide

Examples

# NOT RUN {
#dfgenotype is usually generated by the function 'struc_preparation2'.
#Here, a simple example is done by hand.
var1 <- c("00","01","02","10","11","12","20","21","21")
var2 <- c(0,0,0,1,1,1,2,2,2)
var3 <- c(0,1,2,0,1,2,0,1,2)
dfgenotype <- data.frame(genotype=var1,l1=var2,l2=var3, stringsAsFactors = TRUE)
sel_resist(Rmx=10, epis=0, dom=1)
# }