×

Location-scale mixed models and goodness-of-fit assessment applied to insect ecology. (English) Zbl 1521.62413

Summary: Survival models have been extensively used to analyse time-until-event data. There is a range of extended models that incorporate different aspects, such as overdispersion/frailty, mixtures, and flexible response functions through semi-parametric models. In this work, we show how a useful tool to assess goodness-of-fit, the half-normal plot of residuals with a simulated envelope, implemented in the hnp package in R, can be used on a location-scale modelling context. We fitted a range of survival models to time-until-event data, where the event was an insect predator attacking a larva in a biological control experiment. We started with the Weibull model and then fitted the exponentiated-Weibull location-scale model with regressors both for the location and scale parameters. We performed variable selection for each model and, by producing half-normal plots with simulated envelopes for the deviance residuals of the model fits, we found that the exponentiated-Weibull fitted the data better. We then included a random effect in the exponentiated-Weibull model to accommodate correlated observations. Finally, we discuss possible implications of the results found in the case study.

MSC:

62-XX Statistics

Software:

mpr; hnp; R
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Aitkin, M.; Francis, B.; Hinde, J.; Darnell, R., Statistical Modelling in R (2009), Oxford University Press: Oxford University Press, Oxford · Zbl 1211.62003
[2] Atkinson, A. C., Plots, Transormations, and Regression: An Introduction to Graphical Methods of Diagnostic Regression Analysis (1985), Clarendon: Clarendon, Oxford · Zbl 0582.62065
[3] Broyden, C. G.; Dennis, J. E.; Moré, J. J., On the local and superlinear convergence of quasi-Newton methods, J. Appl. Math., 12, 223-245 (1973) · Zbl 0282.65041
[4] Burke, K.; MacKenzie, G., Multi-parameter regression survival modelling: An alternative to proportional hazards, Biometrics, 73, 678-686 (2017) · Zbl 1372.62056 · doi:10.1111/biom.12625
[5] Cardinale, B. J.; Harvey, C. T.; Gross, K.; Ives, A. R., Biodiversity and biocontrol: Emergent impacts of a multi-enemy assemblage on pest suppression and crop yield in an agroecosystem, Ecol. Lett., 6, 857-865 (2003) · doi:10.1046/j.1461-0248.2003.00508.x
[6] Demétrio, C.G.B., Hinde, J. and Moral, R.A., Models for overdispersed data in entomology, in Ecological Modelling Applied to Entomology, C.P. Ferreira and W.A.C. Godoy, eds., Springer, Berlin, 2014, pp. 219-259.
[7] Kogan, M., Integrated pest management historical perspectives and contemporary developments, Annu. Rev. Entomol., 43, 243-270 (1998) · doi:10.1146/annurev.ento.43.1.243
[8] McNeish, D.; Stapleton, L. M., Modeling clustered data with very few clusters, Multivar. Behav. Res., 51, 495-518 (2016) · doi:10.1080/00273171.2016.1167008
[9] Moral, R. A.; Demétrio, C. G.B.; Hinde, J.; Godoy, W. A.C.; Fernandes, F. S., Parasitism-mediated prey selectivity in laboratory conditions and implications for biological control, Basic Appl. Ecol., 19, 67-75 (2017) · doi:10.1016/j.baae.2016.11.002
[10] Moral, R. A.; Hinde, J.; Demétrio, C. G.B., Half-normal plots and overdispersed models in R: The hnp package, J. Stat. Softw., 81, 1-23 (2017) · doi:10.18637/jss.v081.i10
[11] Nadarajah, S.; Cordeiro, G. M.; Ortega, E. M.M., The exponentiated Weibull distribution: A survey, Statist. Papers, 54, 839-877 (2013) · Zbl 1307.62033 · doi:10.1007/s00362-012-0466-x
[12] R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria; software available at: http://www.R-project.org.
[13] Stroud, A. H.; Secrest, D., Gaussian Quadrature Formulas (1966), Prentice-Hall · Zbl 0156.17002
[14] Therneau, T. M.; Grambsch, P. M.; Fleming, T. R., Martingale-based residuals for survival models, Biometrika, 77, 147-160 (1990) · Zbl 0692.62082 · doi:10.1093/biomet/77.1.147
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.