American Journal of Mathematics and Statistics

p-ISSN: 2162-948X    e-ISSN: 2162-8475

2019;  9(1): 17-22

doi:10.5923/j.ajms.20190901.03

 

Parametric Survival Modeling of Maternal Obstetric Risk; a Censored Study

Obubu Maxwell1, Babalola A. Mayowa2, Ikediuwa U. Chinedu1, Amadi E. Peace3

1Department of Statistics, Nnamdi Azikiwe University, Awka, Nigeria

2Department of Statistics, University of Ilorin, Ilorin, Nigeria

3Department of Statistics, Abia State Polytechnic, Aba, Nigeria

Correspondence to: Obubu Maxwell, Department of Statistics, Nnamdi Azikiwe University, Awka, Nigeria.

Email:

Copyright © 2019 The Author(s). Published by Scientific & Academic Publishing.

This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

Abstract

The world mortality rate has declined 45% since 1990, but still 800 women die every day from pregnancy or childbirth related causes. According to the United Nations Population Fund (UNFPA) this is equivalent to about one woman every two minutes and for every woman who dies, 20 or 30 encounter complications with serious or long-lasting consequences. Young mothers face higher risks of complications and death during pregnancy than older mothers, especially adolescents aged 15 years or younger. Adolescents have higher risks for postpartum hemorrhage, thus, an increased risk of death during pregnancy or childbirth compared with older women. The study is therefore focused on addressing the issue of good statistical estimators, is born out of the weaknesses of the estimators in use, and the apparent lack of research into the application of other methods. Survival analysis technique was employed, once the survival function has been developed, various tests and the modeling of Maternal Mortality, as well as the determination of the appropriate distributions that best describes maternal mortality was done, using both the parametric and non-parametric methods. These include the identification of prognostic factors through regression analysis and the determination of an appropriate distribution for maternal survival. Results of its application to data from Oyo State, Nigeria showed that while about 90% of pregnant women made it alive to delivery, only 86% of them survived to the end of the postpartum period. There were significant differentials by location, and Maternal Age: The Weibull distribution described maternal survival well.

Keywords: Maternal Mortality, Hazard, Survival Analysis, Weibull Distribution, Maternal Survival, Maternal Cox Regression

Cite this paper: Obubu Maxwell, Babalola A. Mayowa, Ikediuwa U. Chinedu, Amadi E. Peace, Parametric Survival Modeling of Maternal Obstetric Risk; a Censored Study, American Journal of Mathematics and Statistics, Vol. 9 No. 1, 2019, pp. 17-22. doi: 10.5923/j.ajms.20190901.03.

1. Introduction

Maternal death is defined as "The death of a woman while pregnant or within 42 days of termination of pregnancy, irrespective of the duration and the site of the pregnancy, from any cause related to or aggravated by the pregnancy or its management, but not from accidental or incidental causes [1]. The world mortality rate has declined 45% since 1990, but still 800 women die every day from pregnancy or childbirth related causes. According to the United Nations Population Fund (UNFPA) this is equivalent to "about one woman every two minutes and for every woman who dies, 20 or 30 encounter complications with serious or long-lasting consequences. Most of these deaths are injuries are entirely preventable" [2]. UNFPA estimated that 289000 women died of pregnancy or childbirth related causes in 2013 [3]. These causes rage from severe bleeding to obstructed labour, all of which have highly effective inventions. As women have gained access to family planning and skilled birth attendance with backup emergency obstetric care, the global maternal mortality ratio has fallen from 380 maternal death per 100000 live births in 1990 to 210 deaths per 100000 live births in 2013. This has resulted in many countries having their maternal death rates [4] while there has been a decline in world-wide mortality rates; much more has to be done. High rates still exist, particularly in impoverished communities with over 85% living in Africa [5]. The effect of a mother's death results in vulnerable families and their infants, if they survive, are more likely to die before reaching their second birthday [6]. Factors that increase maternal death can be direct or indirect. Generally, there is a distinction between a direct maternal death that is the result of a complication of the pregnancy, delivery, or management of the two, and an indirect maternal death [7] that is a pregnancy-related death in a patient with a preexisting or newly developed health problem unrelated to pregnancy. Fatalities during but unrelated to a pregnancy are termed accidental, incidental, or non-obstetrical maternal deaths. The most common cause of postpartum bleeding (15%), complications from unsafe abortion (15%), hypertensive disorders of pregnancy (10%), postpartum infections (8%), and obstructed labour (6%) [8]. Other causes of maternal death include blood clot (3%) and pre-existing conditions (28%) [9]. Indirect causes are malaria, anemia, [10] HIV/AIDS, and cardiovascular disease, all of which may complicate pregnancy or be aggravated by it. Socio-demographic factors such as age, access to resources and income level are significant indicators of maternal outcome. Young mothers face higher risks of complications and death during pregnancy than older mothers, [11] especially adolescents aged 15 years or younger. [12] Adolescents have higher risks for postpartum hemorrhage, puerperal endometritis, operative virginal delivery, episiotomy, low birth weight, preterm delivery, and small-for-gestational age infants, all of which can lead to maternal death [12]. Surgical support and family support influences maternal outcomes. Furthermore, social disadvantage and social isolation adversely affects maternal health which can lead to increase in maternal death [13]. Additionally, lack of access to skilled medical care during childbirth, the travel distance to the nearest clinic to receive proper care, number of prior births, barriers to accessing prenatal medical care and poor infrastructure all increases maternal deaths. Unsafe abortion is another major cause of maternal death. According to the World Health Organization, every eight minutes a woman dies from complications arising from unsafe abortions. Complications include hemorrhage, infection, sepsis and genital trauma [14]. Globally, preventable deaths from improperly performed procedure constitute 13% of maternal mortality, and 25% or more in some countries where maternal mortality from other causes is relatively low, making unsafe abortion the leading single cause of maternal mortality worldwide [15-20]. This study is therefore focused on addressing the issue of good statistical estimators, is born out of the weaknesses of the estimators in use, and the apparent lack of research into the application of other methods. As stated in [21], the challenge arising from having some censored data, at most one event occurring per subject and highly skewed observations, make the use of standard statistical techniques inappropriate: hence our proposal of the use of survival analysis technique provides more informative estimates with measurable precision and is a much better way of estimating obstetric risk with broader parametric and non-parametric prospects.

2. Materials and Method

In applying Survival analysis techniques to Maternal Mortality, we make the assumption that the time of maternal deaths are a function of the random process and therefore the time of death is a random variable with a probability distribution and consequently, other related distributions.

2.1. Basic Concepts in Survival

Let T, denote the survival time of a woman from her time of conception. The distribution of T can be characterized by three equivalent functions [22].
Survival Function,
2.1.1. Probability Density Function
The probability density function of the survival time is defined as the limit of the probability that a woman dies in the short interval, t to , per unit width, ∆t, or simply the probability of dying in a small interval per unit time. It can be expressed as
Where denotes a woman who has conceived, dying after she has survived to time t. These functions are related by;
2.1.2. Estimating the Survival Functions
Survival functions can be estimated using a number of options; the Kaplan-Meier and the Life Table Method. In this study, the Life table method, a non-parametric method to estimate the survival functions was used. This method, which uses time intervals, estimates the survival functions up to the upper limit of each interval, and estimates the hazard and density functions at the midpoint of each interval [23]. For the interval, let be the end time and be the conditional probability of dying. Then;
Where,
is the mid-point of the interval.
is the number of women dying in the interval after their conception.
is the number of women exposed in the interval after conception.
is the conditional probability of a woman dying in the interval after conception
is the conditional probability of a woman dying in the interval after conception.
is the width of the interval.
The standard error of the survival function [18] is estimated by:
while that of the hazard function [19] is estimated by;
The probability density function [19] is estimated by;
2.1.3. Modeling and Other Tests Based on Survival Techniques
Once the survival function has been developed, various tests and the modeling of Maternal Mortality, as well as the determination of the appropriate distributions that best describes maternal mortality can be done, using both the parametric and non-parametric methods [24]. These include the identification of prognostic factors through regression analysis and the determination of an appropriate distribution for maternal survival.
2.1.4. The Log-Rank Test
The Log-rank test, [25] a non-parametric test of difference for survival functions, is the most widely used technique when data is censored. It measures the difference in the survival for the different groups at each of the given times. For a k factor group, this test the hypothesis that;
This is tested as a chi-square test which compares the observed numbers of failures to the expected number of failures under the hypothesis. Thus, given that and is the observed and expected number of deaths respectively for the group, the test statistic is given by;
Where
is the number of women still exposed to the risk of dying at time up to to t for the group.
is the total number of deaths for all groups at time t. Thus,
has approximately a chi-square distribution with degrees of freedom. A large chi-square value will lead to a rejection of the null hypothesis in favor of the alternative that the k groups do not have the same survival distribution.
2.1.5. Proportional Hazard Regression
The Cox proportional regression [26-28] was used to determine the effect of some socio economic and demographic factors on maternal survival based on the assumption of proportionality. In this model, the hazard for an individual is assumed to be related to the covariates through the equation;
Taking the logarithm of both sides, the model can also be written as
Where
The ratio of the hazard for the two individuals i and j (say rural women and urban women) is then given by;
Where measures the relative risk for the woman over the , with respect to the change in the covariate,
2.1.6. Determination of an Appropriate Distribution
Graphical and analytical method exists for the determination of an appropriate distribution for survival data. In this study, the hazard plotting technique [29], a graphical plotting technique designed to handle censored data was used. The shape of the hazard plot gives a good indication of the distribution that could best fit it: For instance, a constant hazard would suggest an exponential distribution; while a first increasing and then decreasing hazard would suggest a log normal or log-logistic model, a monotonically distribution could suggest a log-logistic model; The Weibull distribution could have an increasing, decreasing or constant hazard. The monotonically increasing shape of the hazard plot for maternal mortality, suggests that the assumption that maternal survival had a Weibull distribution. The parameters of the assumed Weibull distribution were then estimated using the hazard plotting technique. The adequacy of the fit was also assessed using R-squared value of the linear plot of against is the cumulative hazard at time t.

3. Results and Discussion

Table 1. Maternal Survival Estimates
     
Table 1 summarizes the survival estimates with their measures of precision. It shows that, in Delta State Nigeria, approximately 90 percent of pregnant women make it through to delivery while 86 percent make it to the end of the postpartum period: This means a maternal mortality of 14 percent. Approximately 10 percent of women, who conceived, lost their lives during the first month after delivery: representing over 50 percent of the total maternal deaths. This makes the first month after delivery the most risky period for a pregnant woman. Thus, a woman, once pregnant, has approximately a 1 in 10 chance of dying before delivery and a 2 in 10 chance of dying by the end of the ten month after conception (up to one month after delivery). Accordingly, if a pregnant Woman delivered safely, she still had approximately a 1 in 10 chance of dying within one month after delivery. A focus therefore on a mother's health immediately after she was delivered and up to one month thereafter, could reduce maternal mortality by up to half.
Figure 1. Hazard plot for Maternal Survival against time
From the Hazard plot in figure 1 above, we clearly see that maternal mortality increased at an increasing rate. This plot also suggests a possible two part mortality trend: the period before delivery and the period after delivery, with the latter being more perilous.
Table 2. Log-Rank Test of Differences in Maternal Survival
     
The results of the Log-rank tests shown in Table 2, show that, Type of town, ever given birth, marital status and Age of Mother at birth are all associated with maternal mortality.
Table 3. Maternal Cox Regression
     
From Table 3, the results of the Cox regression shows that a mother who has ever given birth has a 67 percent less risk of dying due to pregnancy related issues than a woman who has never given birth. Generally, as a woman grows older, the risk of dying decreases by 4.5 percent for every additional year of age: These results are statistically significant. Area differences were also significant with rural folks having a 67% less risk than urban folks: Thus, the more urban the locality, the higher the risk of maternal deaths. This could be a reflection of interventions in rural communities to make motherhood safer or of the data inadequacies or unreliability.
Table 4. Summary of fitted distribution for Maternal Survival
     
The estimated parameters of the fitted Weibull distribution, and its appropriateness, are displayed on Table 4, the results shows that the Weibull distribution describes Maternal Survival very well.

4. Conclusions and Recommendations

This study has clearly shows that Survival Techniques can be applied to the study of Maternal Survival. Mortality is adequately described by a Weibull model. The Weibull distribution which described maternal survival, had shape and scale parameters , indicating an increasing risk of mortality with time and thus a hazard that increases at an increasing rate: The prognosis factors that influences the hazard, and hence survival, are the determinants of the Weibull distribution that described maternal survival for this data was determined to be and respectively.

References

[1]  http://www.maternalmortalitydata.org/Definitions.html.
[2]  http://www.unfpa.org/maternal-health.
[3]  http://www.unfpa.org/maternal-health.
[4]  http://www.unfpa.org/maternal-health.
[5]  http://www.unfpa.org/maternal-health.
[6]  http://www.unfpa.org/maternal-health.
[7]  Amere Deribew. “Incidence, prevalence and mortality rates of malaria in Ethiopia from 1990 to 2015: analysis of the global burden of diseases 2015 https://malariajournal.biomedcentral.com/articles/10.1186/s12936-017-1919-4.
[8]  Obubu, M. Nwokolo P.C (2016) “Prevalence of Breast Cancer in Delta State, Nigeria.” World Journal of Probability and Statistics Reearch. Vol. 2, No. 2, Pp 1-9.
[9]  Osuji G.A., Obubu, M., Obiora-Ilouno H.O (2016) “Uterine Fibroid on Women's Fertility and Pregnancy Outcome in Delta State, Nigeria.” Journal of Natural Sciences Research, Vol. 6, No 2, pp. 27-33.
[10]  Osuji G.A., Obubu, M., Obiora-Ilouno H.O (2016) “An investigation on the causes of Low birth weight in Delta State, Nigeria” European Journal of Statistics and Probability Vol. 4, No 1, pp. 1-6.
[11]  Osuji G.A., Obubu, M., Nwosu, C.A (2016) “Stock Investment Decision in Nigeria; A PC Approach” World Journal of Multidisciplinary and Contemporary Research, Vol. 2, No 1, pp. 1-11. http://wjmcr.com.
[12]  Osuji G.A., Obubu, M., Nwosu, C.A (2016) “Preconception sex selection using proper ovulation timing” World Journal of Probability and Statistics Research U.S.A, Vol. 2, No 1, pp. 1-12. http://wjmcr.com.
[13]  Osuji, G.A., Okoro, C.N., Obubu, M., Obiora-Ilouno H.O. (2016) “Effect of Akaike Information Criterion on Model Selection in Analyzing Auto-Crash Variables.” International Journal of Sciences: Basic and Applied Research (IJSBAR). Vol. 26, No 1, pp. 98-109.
[14]  Osuji G.A., Obubu, M., Obiora-Ilouno H.O., Okoro, C.N (2015) “Post-Partum Hemorrhage in Delta State, Nigeria; A Logistic Approach.” International Journal of Sciences: Basic and Applied Research (IJSBAR),Vol. 24, No 6, pp. 45-53.
[15]  Osuji G.A., Obubu, M., Obiora-Ilouno H.O., Nwosu, D.F (2015) “Perinatal Mortality and Associated Obstetric Risk Factors in Urban Delta State, Nigeria; Rural-Urban Differences.” International Journal of Mathematics and Statistics Studies Vol. 3, No. 5, PP 32-46.
[16]  Osuji G.A., Obubu, M., Obiora-Ilouno H.O (2015) “An Investigation on Crime Rate in Southeastern Nigeria.” European Journal of Statistics and Probability. Vol. 3, No 4, pp. 1-9.
[17]  Obubu, M., Okoye Valentine, Omoruyi Frederick, Ngonadi Lilian Oluebube (2017) “Infant Mortality; a continuing social problem in Northern Nigeria: Cox Regression Approach. American Journal of Innovative Research and Applied Sciences. 2017; 5(5):1-5.
[18]  Conde-Agudelo A, Belizan JM, Lammers C. Maternal-perinatal morbidity and maternal mortality associated with adolescent pregnancy in Latin America: Cross-sectional study. American Journal of Obstetrics and Gynecology, 2004. 192:342-349 http://www.sciencedirect.com/science/article/pii/S000293780401779X.
[19]  Morgan, K. J., & Eastwood, J. G. (2014). Social determinants of maternal self-rated health in South Western Sydney, Australia. BMC Research Notes, 7(1), 1-12. doi:10.1186/1756-0500-7-51]. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3899616/.
[20]  Haddad, L. B., & Nour, N.M. (2009). Unsafe abortion: unnecessary maternal mortality. Reviews in obstetrics and gynecology, 2(2), 122. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2709326/.
[21]  Fertility Regulation and Reproductive Health in the Millennium Development Goals: The search for a perfect indicator.
[22]  Lee, E. T. & Wang, J. W., 2003. Statistical Methods for Survival Data Analysis. 3rd ed. New Jersey: John Wiley & Sons. http://eu.wiley.com/WileyCDA/WileyTitle/productCd-1118095022.html.
[23]  Gehan, E. A., 1969. Estimating Survival function from the life table. Journal of chronic diseases, 21: 629-644. http://www.sciepub.com/reference/10550.
[24]  Greenwood, M. (1926). “The natural duration of cancer”. Reports on Public Health and Medical subjects, (London: Her Majesty’s Stationary Office), 33: 1-26.
[25]  Gehan, E. A., 1972. Estimating Survival function from the life table. Journal of chronic diseases, 21: 629-644. http://www.sciepub.com/reference/10550.
[26]  Lee, E. T. & Wang, J. W., 2003. Statistical Methods for Survival Data Analysis. 3rd ed. New Jersey: John Wiley & Sons. http://eu.wiley.com/WileyCDA/WileyTitle/productCd-1118095022.html.
[27]  Peto, R. & Peto, J., 1972. Asymptotically Efficent rank Invariant Procedures. Journal of the Royal Statistical Society, 135(Series A): 185-207. https://www.jstor.org/stable/2344317.
[28]  Cox, D. R., 1972. Regression Models and Life Tables. Journal of the Royal Statistical Society, 34(SeriesB): 187-220. http://www.stat.cmu.edu/~ryantibs/journalclub/cox_1972.pdf.
[29]  Nelson, W. (1972). Theory and application of Hazard plotting for censored failure data. Technometrics, 14: 945-966. https://www.jstor.org/stable/1267144.