rRiskGEvec: Supplementary Code and Data Geoffrey R. Hosack, Adrien Ickowicz and Keith R. Hayes Commonwealth Scientific and Industrial Research Organisation Hobart, Tasmania, Australia Abstract The R package rRiskGEvec is provided as supporting electronic supplementary mate- rial for the study ‘Quantifying the risk of vector-borne disease transmission attributable to genetically modified vectors’, Royal Society Open Science, 8:201525. This R package provides the code and data to fully reproduce all analyses and figures. Keywords : R. S1. Introduction This vignette describes supplementary methods, code and data for Hosack, Ickowicz, and Hayes (2021). The results in this paper were obtained using R 3.6.2 with the rRiskGEvec 1.3 package. R itself and all packages used are available from the Comprehensive R Archive Network (CRAN) at https: //CRAN.R-project.org/. It is assumed that the user is able to install R and the rRiskGEvec package. The necessary R documentation is available on the R homepage at https://www.r-project.org/. To cite this R package, please refer to Hosack et al. (2021). S2. Estimation and Prediction a. Elicitation Session The following code chunk loads the package and R data objects needed for the statistical analyses. R> library(rRiskGEvec) R> data("elicitations") R> data("MRR") R> released <- MRR$released R> psc.dat <- MRR$psc.dat See ?elicitations for documentation of the elicited subjective probability distributions of the vec- torial capacity parameters included in this package and presented in the main text. See ?MRR for documentation of the summary statistics from the mark release recapture (MRR) experiment of fe- male mosquitoes in Burkina Faso conducted as part of larger MRR experiment by Epopa, Millogo, Collins, North, Tripet, Benedict, and Diabate (2017)). The model design matrices used for the elicitation session of all parameters except mortality is given by the following code chunk. R> X <- as.matrix(elicitations[[1]]$Z$design) # common to all elicitations R> X <- cbind(X, X[ , "GE"]*X[ , "G"], X[ , "GE"]*X[ , "G"]*X[ , "B"])