The Cry proteins (Cry1Ac and Cry2Ab) content in cotton pollen was compared among the treatments using one-way analysis of variance (ANOVA) followed by Tukey’s post-hoc test. The data from honey bees were analyzed with mixed models and used replicate (cage) as a random factor. The survival dynamics of honey bees were analyzed with Cox proportional hazards regression models, and the cumulative pollen consumption and THC results for honey bees were fitted to a log-linear model. The survival response of B. mori to different dietary treatments was analyzed using the Kaplan-Meier procedure and Logrank test. Nonparametric tests (K independent samples: Kruskal-Wallis H- tests; two independent samples: Mann-Whitney U tests) were performed on the developmental duration of B. mori larvae (from 1st instar to 5th instar), because the assumptions for parametric analyses were not fulfilled. The molter weight and THC results of B. mori larvae were compared using ANOVA, and means were compared by Tukey’s post-hoc test. All statistical tests were conducted using SAS Version 8.0 (SAS Institute Inc., Cary, NC, USA).