Due to non-normality and/or homocedasticity of the data (even after transformation), we used non-parametric tests. The Kruskal-Wallis test was employed to analyze the time to reach the next stage, with the Bonferroni post-hoc test applied when significant differences were observed (p < 0.05). Analysis of variance (ANOVA) was used to analyze the number of dead individuals during each stage, with a Tukey’s post-hoc test to identify significantly different groups. The Kruskal-Wallis non-parametric test was applied to compare the weight of the adults, and the Bonferroni post-hoc test was performed subsequently. To compare overall differences in the duration of each developmental stage for the four treatment groups, we calculated the ratio of the duration in days of the Pupa stage to the Larvae stage (ratioPL) and then applied an ANOVA with a Tukey post-hoc test.