Statistical analysis

The dose-mortality response was assessed by the R-script BioRssay [18]. This computes the doses of insecticides killing 50% and 95% of the tested colony or strain. After developing a protocol for data exploration [19], a Generalized Linear Mixed Model (GLMM) was applied with negative binomial distribution or zero-inflated negative binomial (as the data were over-dispersed) using the glmm ADMB package [20]. The response variables were Ae. albopictus female and egg abundance and parity rates. Explanatory variables are Control/Treatment and Pre/Post treatment and the interaction of both, while random variable was trap. Significant values were corrected by False Discovery Rate. Statistical analysis was performed in R version 2.14.2 [18]. The interpretation of the variables used in our statistical analysis is summarized in Additional file 5: Table S1.