Survival analysis is a fundamental tool in medical science for time-to-event data. However, its application to colony organisms like bees poses challenges due to their social nature. Traditional survival models may not accurately capture the interdependence among individuals within a colony. Frailty models, accounting for shared risks within groups, offer a promising alternative. This study evaluates the performance of semi-parametric shared frailty models (gamma, inverse normal, and positive stable-in comparison to the traditional Cox model using bees’ survival data). We examined the effect of misspecification of the frailty distribution on regression and heterogeneity parameters using simulation and concluded that the heterogeneity parameter was more sensitive to misspecification of the frailty distribution and choice of initial parameters (cluster size and true heterogeneity parameter) compared to the regression parameter. From the data, parameter estimates for covariates were close for the four models but slightly higher for the Cox model. The shared gamma frailty model provided a better fit to the data in comparison with the other models. Therefore, when focusing on regression parameters, the gamma frailty model is recommended. This research underscores the importance of tailored survival methodologies for accurately analyzing time-to-event data in social organisms.
References
[1]
Buchmann, A., Moore, K. and Fisher, D. (2010) EXPERIENCING FILM TOURISM: Authenticity Fellowship. Annals of Tourism Research, 37, 229-248. https://doi.org/10.1016/j.annals.2009.09.005
[2]
Brittain, C.A., Vighi, M., Bommarco, R., Settele, J. and Potts, S.G. (2010) Impacts of a Pesticide on Pollinator Species Richness at Different Spatial Scales. Basic and Applied Ecology, 11, 106-115. https://doi.org/10.1016/j.baae.2009.11.007
[3]
Wiest, L., Buleté, A., Giroud, B., Fratta, C., Amic, S., Lambert, O., Pouliquen, H. and Arnaudguilhem, C. (2011) Multi-Residue Analysis of 80 Environmental Contaminants in Honeys, Honeybees and Pollens by One Extraction Procedure Followed by Liquid and Gas Chromatography Coupled with Mass Spectrometric Detection. Journal of Chromatography A, 1218, 5743-5756. https://doi.org/10.1016/j.chroma.2011.06.079
[4]
Ndakidemi, B., Mtei, K. and Ndakidemi, P.A. (2016) Impacts of Synthetic and Botanical Pesticides on Beneficial Insects. Agricultural Sciences, 7, 364-372. https://doi.org/10.4236/as.2016.76038
[5]
Kumar, A., Joseph, S., Tsechansky, L., Privat, K., Schreiter, I.J., Schüth, C. and Graber, E.R. (2018) Biochar Aging in Contaminated Soil Promotes Zn Immobilization Due to Changes in Biochar Surface Structural and Chemical Properties. Science of the Total Environment, 626, 953-961. https://doi.org/10.1016/j.scitotenv.2018.01.157
[6]
Shah, P.A. and Pell, J.K. (2003) Entomopathogenic Fungi as Biological Control Agents. Applied Microbiology and Biotechnology, 61, 413-423. https://doi.org/10.1007/s00253-003-1240-8
[7]
Boff, S., Scheiner, R., Raizer, J. and Lupi, D. (2021) Survival Rate and Changes in Foraging Performances of Solitary Bees Exposed to a Novel Insecticide. Ecotoxicology and Environmental Safety, 211, Article 111869. https://doi.org/10.1016/j.ecoenv.2020.111869
[8]
Parish, J.B., Scott, E.S., Correll, R. and Hogendoorn, K. (2019) Survival and Probability of Transmission of Plant Pathogenic Fungi through the Digestive Tract of Honey Bee Workers. Apidologie, 50, 871-880. https://doi.org/10.1007/s13592-019-00697-6
[9]
Van Dooremalen, C., Gerritsen, L., Cornelissen, B., van der Steen, J.J.M., van Langevelde, F. and Blacquière, T. (2012) Winter Survival of Individual Honey Bees and Honey Bee Colonies Depends on Level of Varroa destructor Infestation. PLOS ONE, 7, e36285. https://doi.org/10.1371/journal.pone.0036285
[10]
Duchateau, L., Janssen, P., Lindsey, P., Legrand, C., Nguti, R. and Sylvester, R. (2002) The Shared Frailty Model and the Power for Heterogeneity Tests in Multicenter Trials. Computational Statistics Data Analysis, 40, 603-620. https://doi.org/10.1016/S0167-9473(02)00057-9
[11]
Clayton, D.G. (1978) A Model for Association in Bivariate Life Tables and Its Application in Epidemiological Studies of Familial Tendency in Chronic Disease Incidence. Biometrika, 65, 141-151. https://doi.org/10.1093/biomet/65.1.141
[12]
Vaupel, J.W., Manton, K.G. and Stallard, E. (1979) The Impact of Heterogeneity in Individual Frailty on the Dynamics of Mortality. Demography, 16, 439-454. https://doi.org/10.2307/2061224
[13]
Gorfine, M. and Zucker, D.M. (2023) Shared Frailty Methods for Complex Survival Data: A Review of Recent Advances. Annual Review of Statistics and Its Application, 10, 51-73. https://doi.org/10.1146/annurev-statistics-032921-021310
[14]
Hsu, L., Gorfine, M. and Malone, K. (2007) On Robustness of Marginal Regression Coefficient Estimates and Hazard Functions in Multivariate Survival Analysis of Family Data When the Frailty Distribution Is Mis-Specified. Statistics in Medicine, 26, 4657-4678. https://doi.org/10.1002/sim.2870
[15]
Duchateau, L. and Janssen, P. (2008) The Frailty Model. Springer Science & Business Media, New York.
[16]
Wienke, A. (2010) Frailty Models in Survival Analysis. Chapman and Hall/CRC, New York. https://www.taylorfrancis.com/books/9781420073911 https://doi.org/10.1201/9781420073911
[17]
Balan, T.A. and Putter, H. (2020) A Tutorial on Frailty Models. Statistical Methods in Medical Research, 29, 3424-3454. http://journals.sagepub.com/doi/10.1177/0962280220921889
[18]
Li, Y., Dick, A.W., Glance, L.G., Cai, X. and Mukamel, D.B. (2007) Misspecification Issues in Risk Adjustment and Construction of Outcome-Based Quality Indicators. Health Services and Outcomes Research Methodology, 7, 39-56. https://doi.org/10.1007/s10742-006-0014-z
[19]
Goethals, K., Ampe, B., Berkvens, D., et al. (2009) Modeling Interval-Censored, Clustered Cow Udder Quarter Infection Times through the Shared Gamma Frailty Model. Journal of Agricultural, Biological, and Environmental Statistics, 14, 1-14. https://doi.org/10.1198/jabes.2009.0001
[20]
Yavuz, A. and Lambert, P. (2016) Semi-Parametric Frailty Model for Clustered Interval-Censored Data. Statistical Modelling, 16, 360-391. https://doi.org/10.1177/1471082X16655631
[21]
Zike, D.T., Fenta, H., Workie, D. and Swain, P. (2018) Determinants of Under-Five Mortality in Ethiopia: An Application of Cox Proportional Hazard and Frailty Models. Biostatistics, 10, 123-136. https://doi.org/10.5336/biostatic.2018-60550
[22]
Cox, D.R. (1972) Regression Models and Life-Tables. Journal of the Royal Statistical Society: Series B (Methodological), 34, 187-202. https://doi.org/10.1111/j.2517-6161.1972.tb00899.x
[23]
Hougaard, P. (2000) Analysis of Multivariate Survival Data. Springer New York, New York. https://doi.org/10.1007/978-1-4612-1304-8
[24]
Hougaard, P. (1986) Survival Models for Heterogeneous Populations Derived from Stable Distributions. Biometrika, 73, 387-396. https://doi.org/10.1093/biomet/73.2.387
[25]
Fine, J.P., Glidden, D.V. and Lee, K.E. (2003) A Simple Estimator for a Shared Frailty Regression Model. Journal of the Royal Statistical Society: Series B (Methodological), 65, 317-329. https://doi.org/10.1111/1467-9868.00388
[26]
Oakes, D. (1989) Bivariate Survival Models Induced by Frailties. Journal of the American Statistical Association, 84, 487-493. https://doi.org/10.1080/01621459.1989.10478795
[27]
Tweedie, M.C.K. (1984) An Index Which Distinguishes between Some Important Exponential Families. In: Ghosh, J.K. and Roy, J., Eds., Statistics: Applications and New Directions. Proceedings of the Indian Statistical Institute Golden Jubilee International Conference, Indian Statistical Institute, Calcutta, 579-604.
[28]
Gorfine, M., Zucker, D.M. and Hsu, L. (2006) Prospective Survival Analysis with a General Semiparametric Shared Frailty Model: A Pseudo Full Likelihood Approach. Biometrika, 93, 735-741. https://doi.org/10.1093/biomet/93.3.735
[29]
Monaco, J.V., Gorfine, M. and Hsu, L. (2017) FrailtySurv: General Semi-Parametric Shared Frailty Model. R Package Version, 1.
[30]
Munda, M., Rotolo, F. and Legrand, C. (2012) Parfm: Parametric Frailty Models in R. Journal of Statistical Software, 51, 1-20. https://doi.org/10.18637/jss.v051.i11
[31]
Dempster, A.P., Laird, N.M. and Rubin, D.B. (1977) Maximum Likelihood from Incomplete Data via the EM Algorithm. Journal of the Royal Statistical Society: Series B (Methodological), 39, 1-22. https://doi.org/10.1111/j.2517-6161.1977.tb01600.x
[32]
Klein, J.P. (1992) Semiparametric Estimation of Random Effects Using the Cox Model Based on the EM Algorithm. Biometrics, 48, 795-806. https://doi.org/10.2307/2532345
[33]
Omuse, E.R., Niassy, S., Wagacha, J.M., Ong’Amo, G.O., Lattorff, H.M.G., Kiatoko, N., Mohamed, S.A., Subramanian, S., Akutse, K.S. and Dubois, T. (2022) Susceptibility of the Western Honey Bee Apismellifera and the African Stingless Bee Meliponulaferruginea (Hymenoptera: Apidae) to the Entomopathogenic Fungi Metarhiziumanisopliae and Beauveriabassiana. Journal of Economic Entomology, 115, 46-55. https://doi.org/10.1093/jee/toab211
[34]
Aalen, O.O., Borgan, Ø. and Gjessing, H.K. (2008) Survival and Event History Analysis. Springer, New York. https://doi.org/10.1007/978-0-387-68560-1
[35]
Balan, T.A. and Putter, H. (2019) FrailtyEM: An R Package for Estimating Semiparametric Shared Frailty Models. Journal of Statistical Software, 90, 1-29. https://doi.org/10.18637/jss.v090.i07
[36]
Alves, S.B., Marchini, L.C., Pereira, R.M. and Baumgratz, L.L. (1996) Effects of Some Insect Pathogens on the Africanized Honey Bee, ApisMellifera L. (Hym., Apidae). Journal of Applied Entomology, 120, 559-564. https://doi.org/10.1111/j.1439-0418.1996.tb01652.x
[37]
Seiedy, M., Saboori, A., Allahyari, H., Talaei-Hassanloui, R. and Tork, M. (2010) A Laboratory Investigation on the Virulence of Two Isolates of the Entomopathogenic Fungus Beauveriabassiana against the Twospotted Mite, Tetranychusurticae (Acari: Tetranychidae). International Journal of Acarology, 8, 527-532. https://doi.org/10.1080/01647954.2010.519718
[38]
Oreste, M., Bubici, G., Poliseno, M., Triggiani, O. and Tarasco, E. (2012) Pathogenicity of Beauveriabassiana (Bals.-Criv.) Vuill. and Metarhiziumanisopliae (metschn.) sorokin against Galleria mellonella L. and Tenebriomolitor L. in Laboratory Assays. Journal of Zoology, 95, 43-48.
[39]
Susan Gachau, W. (2014) Frailty Models with Applications in Medical Research: Observed and Simulated Data. Ph.D. Thesis, University of Nairobi, Nairobi. http://erepository.uonbi.ac.ke/handle/11295/73552
[40]
Adham, S.A. and AlAhmadi, A.A. (2016) Gamma and Inverse Gaussian Frailty Models: A Comparative Study. International Journal of Mathematics and Statistics Invention, 4, e4105.
[41]
Glidden, D.V. and Vittinghoff, E. (2004) Modelling Clustered Survival Data from Multicentre Clinical Trials. Statistics in Medicine, 23, 369-388. https://doi.org/10.1002/sim.1599
[42]
Su, Y.-R. and Wang, J.-L. (2016) Semiparametric Efficient Estimation for Shared-Frailty Models with Doubly-Censored Clustered Data. Annals of Statistics, 44, 1298-1331. https://doi.org/10.1214/15-AOS1406
