- Methodology
- Open access
- Published:

# Modeling time-to-cure from severe acute malnutrition: application of various parametric frailty models

*Archives of Public Health*
**volume 73**, Article number: 6 (2015)

## Abstract

### Background

In developing countries about 3.5% of children aged 0–5 years are victims of severe acute malnutrition (SAM). Once the morbidity has developed the cure process takes variable period depending on various factors. Knowledge of time-to-cure from SAM will enable health care providers to plan resources and monitor the progress of cases with SAM. The current analysis presents modeling time-to-cure from SAM starting from the day of diagnosis in Wolisso St. Luke Catholic hospital, southwest Ethiopia.

### Methods

With the aim of coming up with appropriate survival (time-to-event) model that describes the SAM dataset, various parametric clustered time-to-event (frailty) models were compared. Frailty model, which is an extension of the proportional hazards Cox survival model, was used to analyze time-to-cure from SAM. Kebeles (villages) of the children were considered as the clustering variable in all the models. We used exponential, weibull and log-logistic as baseline hazard functions and the gamma as well as inverse Gaussian for the frailty distributions and then based on AIC criteria, all models were compared for their performance.

### Results

The median time-to-cure from SAM cases was 14 days with the maximum of 63 days of which about 83% were cured. The log-logistic model with inverse Gaussian frailty has the minimum AIC value among the models compared. The clustering effect was significant in modeling time-to-cure from SAM. The results showed that *age* of a child and *co-infection* were the determinant prognostic factors for SAM, but *sex* of the child and the *type of malnutrition* were not significant.

### Conclusions

The log-logistic with inverse Gaussian frailty model described the SAM dataset better than other distributions used in this study. There is heterogeneity between the kebeles in the time-to-cure from SAM, indicating that one needs to account for this clustering variable using appropriate clustered time-to-event frailty models.

## Background

Malnutrition is a common cause of morbidity and mortality in developing countries
[1]. It is a risk factor for over 50% of the 11 million annual childhood deaths
[2]. Severe acute malnutrition (SAM) is defined as a very low weight for height less than 70% NCHS median or less than 115 mm of MUAC
[3]. In developing countries, about 3.5% of children aged 0–5 years are suffer from SAM
[4], which is the most important nutritional disease because of its high prevalence and relationship with child mortality
[5]. Currently, there are a number of studies on time of recovery from acute malnutrition in various parts of the world. Survival data is a term used for describing data that measure the time to a given event of interest. The term *survival data* is often used for describing data that measure the time to the occurrence of a given event of interest. The event of interest can be seen as a transition from one state to another. In this study, the event of interest was the time-to-cure from SAM from the day of diagnosis. One of the major aim of this analysis was to model the time-to-cure from SAM among less than five year inpatient children in Wolisso St. Luke Catholic hospital and to compare the efficiency of various parametric frailty models using the same dataset. The classical model for this kind of data is the proportional hazards model popularized by Cox
[6]. However, correct inference based on Cox’s model needs identically and independently distributed samples. One of the reasons why this model is so popular is the ease with which technical difficulties such as censoring and truncation are handled. This is due to the appealing interpretation of the hazard as a risk that changes over time. Naturally, the concept allows for entering covariates to describe their influence and model different levels of risk for different subgroups. Nonetheless, subjects may be exposed to different risk levels, even after controlling for known risk factors; as some relevant covariates were often unavailable to the researcher or even unknown (univariate frailty case). The study population may also be divided into clusters so that subjects from the same cluster behave more cohesively than subjects from different clusters (multivariate frailty case). The frailty model, introduced in the statistical literature by Vaupel et al.
[7], and discussed in details
[8–10], accounts for heterogeneity in baseline. It is an extension of the proportional hazards of Cox’s model in which the hazard function depends upon an unobservable random quantity, the so-called frailty that acts multiplicatively on it. Study subjects (children) in this study came from clustered community and hence clustered child survival data may be correlated at the Kebele level. In this research, shared frailty models were explored assuming that children within the same cluster (kebele) share similar risk factors, which will take care of the frailty term at kebele level. This model is a conditional independence model where the frailty is common to all individuals in a cluster, and therefore responsible for creating dependence between event times. This is because ignoring the full dependence among observations may lead to standard errors that are understated and parameter estimates that are both biased and inconsistent
[11]. Estimation of the frailty model can be parametric or semi-parametric. In the first case, a parametric density is assumed for the event times, resulting in a parametric baseline hazard function. Estimation is then conducted by maximizing the marginal log-likelihood function. In the second case, the baseline hazard is left unspecified and more complex techniques are available to approach that situation
[12]. Even though semi-parametric estimation offers more flexibility, the parametric estimation will be more powerful if the form of the baseline hazard is somehow known in advance
[13]. In this study, parametric frailty models were used to investigate the relationship between different potential covariates (*sex*, *age*, *type of malnutrition* and *co-infection*) and time to cure from SAM for clustered survival data with random right censoring. The choice of distribution for the hazard is very important than the choice of frailty distribution
[14]. Hence, in this research exponential, weibull and log-logistic hazard functions were used and compared for their efficiencies. Regarding the frailty distribution, we assumed gamma and inverse Gaussian distributions. For comparison of different distributions, the AIC criteria were used, but for comparing nested models, likelihood ratio test were used.

## Methods

### Study sample and setting

The data set used in this study were obtained from Wolisso St. Luke catholic hospital, Wolisso, south west Ethiopia. Children aged under five, having marasmus and/or kwashiorkor, diagnosed for severe acute malnutrition (SAM) according to the protocol for management of SAM [15] and admitted to inpatient, and started treatment during the period January 1, 2010 to January 31, 2012 were eligible for the analysis. A total of 929 children who came from 280 kebeles around Wolisso district were considered. Kebeles that contribute only a child were omitted since the shared frailty model should be done on at least two children per kebele. Therefore a total of 855 children with severe acute malnutrition from 206 kebeles were considered.

The response variable time-to-cure from SAM were obtained by calculating the difference (in day) from the start of treatment until the child were cured (recovered) or censored. Cured children according to the SAM treatment protocol are defined as children who have weight for height >85% and no bilateral edema. Children were considered to be cured and discharged, which is our event of interest, if they fulfilled the discharging criteria for SAM
[15]. However, the time-to-cure were censored for those children transferred to other hospital, dropped treatment, died, did not cure at January 31, 2012 (at the end of study). The following variables were considered for their influence on the time-to-cure from SAM; *sex*, *age*, *type of malnutrition* and *co-infection*. For *age* we used six categories; 0–5 months, 6–11 months, 12–23 months, 24–35 months, 36–47 months and 48–59 months. *Types of malnutrition* were categorised as Marasmus, Kwashiorkor and Marasmic-kwashiorkor. *Co-infection* was categorized based on whether the child has co-infection such as malaria, anemia, pneumonia, measles and giardiasis or not.

### Shared frailty model

Conditional on the random term, called the frailty denoted by *u*_{
i
}, the survival time in cluster *i* (1≤*i*≤*n*) are assumed to be independent, the proportional hazard model assumes

Whereas an alternative if the proportional hazards assumption does not hold is the accelerated failure time frailty model which assumes

Where *i* indicates the *i*^{th} cluster, *j* indicates the *j*^{th} individual in the *i*^{th} cluster, *h*_{0}(*t*) is the baseline hazard, *u*_{
i
} the random term for all subjects in cluster *i*, *X*_{
ij
} is the vector of covariates for subject *j* in cluster *i*, and ** β** the vector of regression coefficients.

We assumed that *Z* (where *Z* = *exp*(*u*_{
i
})) has the gamma or inverse Gaussian distribution so that the hazard function depends upon this frailty that acts multiplicatively on it. The main assumption of a shared frailty model is that all individuals in cluster *i* share the same value of frailty *Z*_{
i
}(*i* = 1,…,*n*), and that is why the model is called the shared frailty model. The survival time is assumed to be conditionally independent with respect to the shared (common survival times) frailty. This shared frailty is the cause of dependence between survival time within the clusters.

In order to investigate the effect of the candidate covariates on the time-to-cure from SAM, we first did a univariable analysis by fitting a separate model for each candidate covariates. Covariates that were found to be significant in the univariable analysis were included in the multivariable analysis.

The multivariable survival analysis in the study was done by assuming the exponential, weibull and log-logistic distributions for the baseline hazard function; and the gamma and inverse Gaussian frailty distributions. It was performed using the three most significant covariates in the univariable analysis namely *age*, *type of malnutrition* and *co-infection*. However, we excluded *sex* which was not significant in univariable analysis.

## Results

Of all 855 malnourished patients, 711(83.16%) were cured and the median cure time from SAM was 14 days, while the minimum and the maximum cure times observed were 7 and 63 days, respectively (Table 1).

Using all the multivariable frailty models, the covariate *co-infection* was significant, indicating that it was the most important prognostic factor for the time-to-cure from SAM. *Age group* was significant in the three models namely, weibull-gamma, weibull-inverse Gaussian and log-logistic-inverse Gaussian frailty models. *Type of malnutrition* was not a significant factor for time-to-cure from SAM using all the models. The variance of the random effect (*θ*) was significant for the weibull and log-logistic baseline frailty models but not significant for the exponential baseline frailty models. It was highest when we assume the inverse Gaussian distribution (*θ* = 0.21) followed by the gamma distribution (*θ* = 0.172) with the log-logistic baseline hazard function. This term was again higher when we assume the inverse Gaussian frailty distribution (*θ* = 0.169) than the gamma distribution (*θ* = 0.163) for the weibull baseline hazard function. The Kendall’s tau (*τ*) was higher for the higher *θ* values. Accordingly the dependence within the clusters for the log-logistic-inverse Gaussian frailty model (*τ* = 0.081) was the maximum followed by the log-logistic-gamma frailty model (*τ* = 0.079). The AIC value of the log-logistic-inverse Gaussian model 4941.63, was the minimum among all the other AIC values of the models indicating that it was the most efficient model to describe the SAM dataset using various parametric frailty models (Table
2).

Analysis based on log-logistic-inverse Gaussian frailty model showed that *age group* of the children and presence of *co-infection* were significant. That is the confidence interval of acceleration factors(*ϕ*) for both *Age group* and *co-infection* did not include 1 (Table
3).

Children with age group between 12 and 23 months had significantly different curing time than the reference groups (aged 0–5months) with acceleration factor (*ϕ* = 1.18) and confidence interval (1.014, 1.374), did not includes 1, at 5% level of significance. An acceleration factor of greater than 1 indicates prolonging the time-to-cure from SAM. Therefore, children aged between 12 and 23 months had prolonged cure time by a factor of 1.18 than the reference groups. But the other age groups were not significantly different from the baseline age group at 5% level of significance. The confidence interval of the acceleration factor of *co-infection* was (1.093, 1.268), did not include 1, indicating *co-infection* is also significant prognostic factor for time-to-cure from SAM. It prolonged curing time by a factor of (*ϕ* = 1.177) at 5% level of significance.

The value of the shape parameter in the log-logistic-inverse Gaussian frailty model was (*ρ* = 3.56). This value greater than unity indicates that the shape of hazard function is unimodal, i.e., it increases up to some time and then decreases. The variability (heterogeneity) in the population of clusters (kebeles) estimated by our working model was *θ* = 0.21, and the dependence within clusters was about *τ* = 8.1*%*.

The predicted frailty values increases with range of 0.731 to 1.272 with increasing the median time of the cluster on log-logistic- inverse Gaussian frailty model (Figure 1). That is, these values are lower for lower values of event times and higher for higher value of event times. The median value of the frailty distribution is around 1.

The hazard functions given the 25^{th}, 50^{th} and 75^{th} quartiles of frailty distribution (*Z* = 0.837,*Z* = 1*and**Z* = 1.167, respectively) were plotted for the log-logistic-inverse Gaussian frailty model (Figure
2). The parameter estimates obtained by maximizing the likelihood of the log logistic-inverse Gaussian frailty model were *λ* = 1.306*e*-4, *ρ* = 3.56(*S*.*E* = 0.134) and *θ* = 0.21(*S*.*E* = 0.055).

The conditional hazard functions given the 75^{th} quartile frailty values was greater, followed by the conditional hazard functions given by the 50^{th} (median) and 25^{th} quartile frailty values of the clusters respectively. From this what we can observe is that the more frail groups (*Z* > 1) were more likely to get events earlier (cure earlier) and the less frail groups (*Z* < 1) had prolonged cure time. All the conditional hazard functions were almost equal at the beginning time (*t* = 0). But the gap widens through time, specifically at mid time. The pattern of all the hazard function is unimodal (increases up to some point and then decreases) as the shape parameter for the baseline hazard function is greater than unity (*ρ* = 3.56).

To check the adequacy of our baseline hazard, the exponential has been plotted by the cumulative hazard function with time-to-cure from SAM. Similarly, the weibull has been plotted by the logarithm cumulative hazard function with the logarithm of time-to-cure from SAM and log-logistic has been plotted by the logarithm of the failure odds with the logarithm of time-to-cure from SAM (Figure 3). The plot of log-logistic was more linear than the other plots, though only few observations were scattered at the beginning time. The patterns suggested that the log-logistic hazard function was appropriate in the model.

The Cox-Snell residuals together with their cumulative hazard function were obtained by fitting the exponential, weibull and log-logistic models to our dataset, via maximum likelihood estimation (Figure 4). The plots showed that the Cox-Snell residuals fitted to assess the log-logistic model for the dataset were nearest to the line through the origin as compared to the other models, again indicating that this model described the SAM dataset well.

A quantile-quantile or q-q plot was made to check if the accelerated failure time provided an adequate fit to the data using two different groups of population. We checked the adequacy of the accelerated failure-time model by comparing the significantly different age groups (children in the age group 0–5 months and 12–24 months); as well as the co-infected and non co-infected groups of patients (Figure
5). The figures appear to be approximately linear for both covariates *age group* and *co-infection* with slopes equivalent to the acceleration factors 1.180 and 1.177, respectively. Therefore the log-logistic baseline model used for time-to-cure from SAM was accelerated failure time model.

## Discussion

The main aim of the study was to model time-to-cure from SAM using appropriate survival model among various parametric frailty models. The comparison of distributions of the models was performed using the AIC criteria, where a model with minimum AIC is accepted to be the best [13]. Accordingly, the log-logistic-inverse Gaussian frailty model which has AIC value of 4941.63 was the most appropriate model to describe the SAM dataset.

This study also showed that there was a clustering (frailty) effect on modeling time-to-cure from SAM which might be due to the heterogeneity in kebele from which the child came, assuming children coming from the same kebele share similar risk factors related to SAM. Therefore, it was important considering the clustering effect in modeling the hazard function. Clusters with minimum median time have smaller frailties, so that these clusters are predicted to have a high hazard [9], more probable to cure in this case. Our data demonstrated that the nuisance (frailty) terms modified the hazard function, and therefore the the hazard function was evaluated conditionally on this effect. Kebeles that frail more were more likely to cure than the less frail kebeles (since the event is positive).

The inverse Gaussian frailty generates very strong dependence at mid time
[16]. According to the conditional hazard function given the 25^{th}, 50^{th} and 75^{th} quantile frailty values, the hazard function depends on these values especially at the mid time (Figure
1). This frailty distribution introduced as an alternative to the gamma distribution
[17] was better in this dataset compared to the gamma frailty distribution. The heterogeneity in the kebeles was estimated to be *θ* = 0.21, and the dependence within clusters is about *τ* = 8.1*%*. These values were the maximum among the variance of the random effects and the Kendall’s tau of all the models, which consolidates the idea that the better the model, the less unobserved heterogeneity will be
[18].

Nonetheless, the most acknowledged parametric model is the weibull, which allows the proportional hazards and accelerated life time model
[8]; the SAM data set was best described by the log-logistic baseline as compared to the exponential and weibull hazard functions. According to the diagnostic plots the *log failure odds* of log-logistic baseline with *log time* was more linear as compared to the plots of exponential (*cumulative hazard* versus *time*) and weibull (*log cumulative hazard* versus *log time*), showing the SAM dataset was best described by the log-logistic baseline. This result was also confirmed by the cumulative hazard plots for the Cox-Snell residuals of the exponential, weibull and the log-logistic models. The plot was more approached to the line in case of the log-logistic model, indicating that the log-logistic was best. A q-q plot was done to check if the accelerated failure time provided an adequate fit to the dataset and the log-logistic as baseline was accelerated failure time model. Hence, a survival model need not be chosen arbitrarily to fit event times, the baseline hazard function as well as the frailty distribution should be compared and the most appropriate model should be selected for appropriate inference.

The prognostic factors considered were the *age group* of the child, *type of malnutrition* and presence or absence of *co-infection*, which were significant covariates using univariable analysis. Analysis using the best model, log-logistic-inverse Gaussian frailty model showed that the *age group* of children and presence of *co-infection(s)* were the determinant factors for the time-to-cure from SAM. Children aged between 12–23 months had prolonged cure time as compared to older age groups. This *age group* is known to be the time when the prevalence of SAM is the highest
[19–21]. This may be related to the fact that they start sub-optimal complementary feeding and compromise breast feeding practice which is important in preventing malnutrition among children
[22]. Literatures like
[23–25] identified infection as a prognostic indicators, likewise, co-infection(s) prolonged the time-to-cure from SAM in this study. However, our findings showed that the time of curing did not depend on *type of malnutrition*. Similarly, Efrem et al.
[26] showed that there was no significant difference in average length of stay among patients with severe wasting and edematous malnutrition. Therefore, a special attention should be given to children with SAM between one to two years old and/or co-infected by other disease(s) in order to control excess numbers of children in the program at any one time which might increase the cost of the program. A program is well functioning and considered *acceptable* if the length of stay in a hospital is less than 4 weeks with acceptable cure proportion of greater than 75%
[15]. This study revealed that the median time-to-cure from SAM for patients in Wolliso St. Luke Catholic hospital was 14 days with maximum cure time of 63 days of which (83.16%) were cured. This implies that the program was *acceptable* as per the above standard.

We classified SAM cases as marasmus, kwashiorkor and marasmic-kwashiorkor, a naming might imply that the cause is protein and calorie deficiency perse. However, the intent was to see if there are differential in recovery time between edematous and non-edematous cases of SAM. We acknowledge the limitation of our study for not being able to isolate SAM cases with micronutrient deficiency.

## Conclusion

The most appropriate statistical model for our dataset among various parametric frailty models, which well described the time-to-cure from SAM of the patients who were diagnosed in Wolisso St. Luke Catholic hospital is the log-logistic-inverse Gaussian frailty model. There is a frailty (clustering) effect on time-to-cure from SAM that arises due to heterogeneity between the kebeles of the children. The median curing time of the children is about 14 days with maximum cure time of 63 days of which 83.16% were cured. These values show *acceptable* functioning of the program in the hospital. Children who are aged between a year and two years and/or co-infected by other disease(s) prolonged their curing time as compared to the other groups of the patients among children under the age of five years.

## References

Black R, Morris S, Bryce J:

**Where and why are 10 million children dying every year?***Lancet*2003,**361:**937.Black R, Allen L, Bhutta Z, Caulfield L, Onis M, Ezzati M:

**Maternal and child under nutrition: global and regional exposures and health consequences.***Lancet*2008,**371:**243–60. 10.1016/S0140-6736(07)61690-0Asayehegn T, Mekitie W, Girma A, Kebede D:

**Cost effectiveness of community-based and in-patient therapeutic feeding programs to treat severe acute malnutrition in Ethiopia.***Biomed Central*2012. Available from (accessed on: July 11, 2013) http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3323427/Kerac M:

**Improving the treatment of severe acute malnutrition in childhood. Centre for international health & development, institute of child health, UCL.**2009.http://discovery.ucl.ac.uk/1306755/1/1306755.pdfZelalem T, Tsinuel G, Ayalew T:

**Treatment outcome of severe acute malnutrition in children with and without hiv infection: a historical cohort study in South-West Ethiopia.***Ethiop J Pediatr Child Health*2010,**6:**14–26.Cox D:

**Regression models and lifetables (with discussions).***J Roy Stat Soc*1972,**34:**187–220.Vaupel J, Manton K, Stallard E:

**The impact of heterogeneity in individual frailty on the dynamics of mortality.***Demography*1979,**16:**439–454. 10.2307/2061224Hougaard P:

*Analysis of Multivariate Survival Data*. New York: Springer-Verlag; 2000.Duchateau L, Janssen P:

*The Frailty Model*. New York: Springer-Verlag; 2008.Wienke A:

*Frailty Models in Survival Analysis*. Chapman and Hall/CRC biostatistics series; 2011. Available from: (accessed on: June may 2 2013) http://www.amazon.com/Frailty-Survival-Analysis-Chapman-Biostatistics/dp/1420073885#reader_1420073885Sastry N:

**A nested frailty model for survival data, with an application to the study of child Survivalin Northeast Brazil.***J Am Stat Assoc*1997,**92:**426–435. 10.1080/01621459.1997.10473994Abrahantes J, Legrand C, Burzykowski T, Janssen P, Ducrocq V, Duchateau L:

**Comparison of different estimation procedures for proportional hazards model with random effects.***Comput Stat Data Anal*2007,**51:**3919–3930.Munda M:

**Parametric frailty models in R.***Am Stat Assoc*2012,**55:**1–21.Do I, Lee Y:

**Estimating frailty models via Poisson hierarchical generalized linear models.***J Comput Graph Stat*2003,**12:**663–681. 10.1198/1061860032256Chamois S:

**Protocol for the management of severe acute malnutrition.***Ethiopia-Federal Minist Health*2007. Available from: [accessed (accessed on: June 10, 2013) http://motherchildnutrition.org/resources/pdf/mcn-protocol-for-the-management-of-severe-acute-malnutrition.pdfHougaard P:

**Frailty models for survival data.***Lifetime Data Anal*1995,**1:**255–273. 10.1007/BF00985760Hougaard P:

**Life table methods for heterogeneous populations: Distributions describing heterogeneity.***Biometrika*1984,**71:**75–83. 10.1093/biomet/71.1.75Zorn CJW:

**Unobserved heterogeneity and frailty models.***Polit Anal*2000,**2:**79–86.Thailand Burma BorderConsortium:

**Nutrition survey report to CCSDPT health agencies.***TBBC*2012. Available from: (accessed on: July 05, 2013) http://reliefweb.int/sites/reliefweb.int/files/resources/nutrition-2012-report-ccsdpt-health-agencies.pdfPlay Therapy Africa:

**Emotional stimulation in the context of emergency food interventions.***Biometrica*2009,**15:**5–20.Center for disease control:

**Malnutrition and micronutrient deficiencies among bhutanese refugee children Nepal.***MMWR*2007,**57:**370–373.Eko Z:

**The influence of weaning pattern on nutritional statusamong 12–24 months old children, in purworejo, central java, Indonesia.**1999.http://www.chnrl.net/publikasi/pdf/TES-EF.pdfSchofield C, Ashworth A:

**Why have mortality rates for severe malnutrition remained so high?***Bull World Health Organ*1996,**74:**223–229.Bitwe R, Dramaix M, Hennart P:

**Simplified prognostic model of overall intrahospital mortality of children in central Africa.***Trop Med Int Health*2006,**11:**73–80. 10.1111/j.1365-3156.2005.01540.xSunguya B, Koola J, Atkinson S:

**Infections associated with severe malnutrition among hospitalised children in East Africa.***Tanzan Health Res Bull*2006,**8:**189–92.Efrem T, Meskele L, Sahle S, Zerihun B, Daniel G, Mohammed A:

**Treatment outcome of children with severe acute malnutrition admitted to therapeutic feeding centers in Southern Region of Ethiopia.***Ethiop J Health Dev*2010,**24:**234–138.

## Acknowledgments

The authors gratefully acknowledge Wolisso St.Luke’s Catholic hospital and College of Nursing for allowing to use the data. This work was financially supported by the College of Natural Science, Jimma University.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

AB conceived the idea, AB DS TB BB YG contributed in the design analyses and interpretation, AB the corresponding author drafted the manuscript. All authors read and approved the final manuscript.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## About this article

### Cite this article

Banbeta, A., Seyoum, D., Belachew, T. *et al.* Modeling time-to-cure from severe acute malnutrition: application of various parametric frailty models.
*Arch Public Health* **73**, 6 (2015). https://doi.org/10.1186/2049-3258-73-6

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/2049-3258-73-6