Homogeneity pursuit and variable selection in regression models for multivariate abundance data.

When building regression models for multivariate abundance data in ecology, it is important to allow for the fact that the species are correlated with each other. Moreover, there is often evidence species exhibit some degree of homogeneity in their responses to each environmental predictor, and that most species are informed by only a subset of predictors. We propose a generalized estimating equation (GEE) approach for simultaneous homogeneity pursuit (ie, grouping species with similar coefficient values while allowing differing groups for different covariates) and variable selection in regression models for multivariate abundance data. Using GEEs allows us to straightforwardly account for between-response correlations through a (reduced-rank) working correlation matrix. We augment the GEE with both adaptive fused lasso- and adaptive lasso-type penalties, which aim to cluster the species-specific coefficients within each covariate and encourage differing levels of sparsity across the covariates, respectively. Numerical studies demonstrate the strong finite sample performance of the proposed method relative to several existing approaches for modeling multivariate abundance data. Applying the proposed method to presence-absence records collected along the Great Barrier Reef in Australia reveals both a substantial degree of homogeneity and sparsity in species-environmental relationships. We show this leads to a more parsimonious model for understanding the environmental drivers of seabed biodiversity, and results in stronger out-of-sample predictive performance relative to methods that do not accommodate such features.

Interval estimation of relative risks for combined unilateral and bilateral correlated data.

Measurements are generally collected as unilateral or bilateral data in clinical trials, epidemiology, or observational studies. For example, in ophthalmology studies, the primary outcome is often obtained from one eye or both eyes of an individual. In medical studies, the relative risk is usually the parameter of interest and is commonly used. In this article, we develop three confidence intervals for the relative risk for combined unilateral and bilateral correlated data under the equal dependence assumption. The proposed confidence intervals are based on maximum likelihood estimates of parameters derived using the Fisher scoring method. Simulation studies are conducted to evaluate the performance of proposed confidence intervals with respect to the empirical coverage probability, the mean interval width, and the ratio of mesial non-coverage probability to the distal non-coverage probability. We also compare the proposed methods with the confidence interval based on the method of variance estimates recovery and the confidence interval obtained from the modified Poisson regression model with correlated binary data. We recommend the score confidence interval for general applications because it best controls converge probabilities at the 95% level with reasonable mean interval width. We illustrate the methods with a real-world example.

A latent functional approach for modeling the effects of multidimensional exposures on disease risk.

Understanding the relationships between exposure and disease incidence is an important problem in environmental epidemiology. Typically, a large number of these exposures are measured, and it is found either that a few exposures transmit risk or that each exposure transmits a small amount of risk, but, taken together, these may pose a substantial disease risk. Further, these exposure effects can be nonlinear. We develop a latent functional approach, which assumes that the individual effect of each exposure can be characterized as one of a series of unobserved functions, where the number of latent functions is less than or equal to the number of exposures. We propose Bayesian methodology to fit models with a large number of exposures and show that existing Bayesian group LASSO approaches are a special case of the proposed model. An efficient Markov chain Monte Carlo sampling algorithm is developed for carrying out Bayesian inference. The deviance information criterion is used to choose an appropriate number of nonlinear latent functions. We demonstrate the good properties of the approach using simulation studies. Further, we show that complex exposure relationships can be represented with only a few latent functional curves. The proposed methodology is illustrated with an analysis of the effect of cumulative pesticide exposure on cancer risk in a large cohort of farmers.

Clinical trial design and analysis for comparing three treatments with intra-individual right- and left-hand data.

BACKGROUND: Chemotherapy-induced peripheral neuropathy can occur in the right and left hand. Studies on prevention treatments for chemotherapy-induced peripheral neuropathy have largely adopted either self-controlled designs or parallel designs to compare two preventive treatments. When three treatment options (two experimental treatments and a control treatment) are available, both designs can be extended. However, no clinical trials have adopted a self-controlled design to compare three prevention treatments for chemotherapy-induced peripheral neuropathy. The incomplete block crossover design for more than two treatments can be extended to compare three treatments in the self-controlled design. In simple extension, some of the participants receive two experimental treatments in both hands; however, it may be difficult to administer different experimental treatments in both hands for practical reasons, such as a concern for the different types of unexpected adverse events. This study proposes a design and analysis method appropriate for the situation where only one experimental treatment is provided to each participant. METHODS: We assume clinical trials to compare each of the two experimental treatments (E1 and E2) with the control treatment (C) and between two experimental treatments only when both experimental treatments are superior to the control treatment. We propose a self-controlled design, which equally randomizes to four arms to adjust for the dominant hand effect: Arm 1: E1 for right hand, C for left hand; Arm 2: C for right hand, E1 for left hand; Arm 3: E2 for right hand, C for left hand; and Arm 4: C for right hand, E2 for left hand. We compare operating characteristics of the proposed design with the three-arm parallel design in which the same treatment is performed in both hands by participants. We also assess three proposed analysis methods for comparisons between experimental treatments in the self-controlled design under several conditions of correlations between right and left hands using simulation studies. RESULTS: The simulation studies showed that the proposed design was more powerful than the three-arm parallel design when correlation was 0.3 or higher. For comparisons between experimental treatments, the methods based on the regression model, including the outcome of hands with C as a covariate, had the highest power under modest to high correlation among the analysis methods in the self-controlled design. CONCLUSION: The proposed design can improve the power for comparing between two experimental treatments and the control treatment. Our design is useful in situations where it is undesirable for participants to receive different experimental treatments in both hands for practical reasons.

A model selection criterion for clustered survival analysis with informative cluster size.

We propose a model selection criterion for correlated survival data when the cluster size is informative to the outcome. This approach, called Resampling Cluster Survival Information Criterion (RCSIC), uses the Cox proportional hazards model that is weighted with the inverse of the cluster size. The RCSIC based on the within-cluster resampling idea takes into account the possible variability of the within-cluster subsampling and the possible informativeness of cluster sizes. The RCSIC allows for easy execution for the within-cluster resampling idea without a large number of resamples of the data. In contrast with the traditional model selection method in survival analysis, the RCSIC has an additional penalization for the within-cluster subsampling variability. Our simulations show the satisfactory results where the RCSIC provides a more robust power for variable selection in terms of clustered survival analysis, regardless of whether informative cluster size exists or not. Applying the RCSIC method to a periodontal disease studies, we identify the tooth loss in patients associated with the risk factors, Age, Filled Tooth, Molar, Crown, Decayed Tooth, and Smoking Status, respectively.

lmerSeq: an R package for analyzing transformed RNA-Seq data with linear mixed effects models.

BACKGROUND: Studies that utilize RNA Sequencing (RNA-Seq) in conjunction with designs that introduce dependence between observations (e.g. longitudinal sampling) require specialized analysis tools to accommodate this additional complexity. This R package contains a set of utilities to fit linear mixed effects models to transformed RNA-Seq counts that properly account for this dependence when performing statistical analyses. RESULTS: In a simulation study comparing lmerSeq and two existing methodologies that also work with transformed RNA-Seq counts, we found that lmerSeq was comprehensively better in terms of nominal error rate control and statistical power. CONCLUSIONS: Existing R packages for analyzing transformed RNA-Seq data with linear mixed models are limited in the variance structures they allow and/or the transformation methods they support. The lmerSeq package offers more flexibility in both of these areas and gave substantially better results in our simulations.

Randomization-based confidence intervals for cluster randomized trials.

In a cluster randomized trial (CRT), groups of people are randomly assigned to different interventions. Existing parametric and semiparametric methods for CRTs rely on distributional assumptions or a large number of clusters to maintain nominal confidence interval (CI) coverage. Randomization-based inference is an alternative approach that is distribution-free and does not require a large number of clusters to be valid. Although it is well-known that a CI can be obtained by inverting a randomization test, this requires testing a non-zero null hypothesis, which is challenging with non-continuous and survival outcomes. In this article, we propose a general method for randomization-based CIs using individual-level data from a CRT. This approach accommodates various outcome types, can account for design features such as matching or stratification, and employs a computationally efficient algorithm. We evaluate this method's performance through simulations and apply it to the Botswana Combination Prevention Project, a large HIV prevention trial with an interval-censored time-to-event outcome.

A comparison of methods for multiple degree of freedom testing in repeated measures RNA-sequencing experiments.

BACKGROUND: As the cost of RNA-sequencing decreases, complex study designs, including paired, longitudinal, and other correlated designs, become increasingly feasible. These studies often include multiple hypotheses and thus multiple degree of freedom tests, or tests that evaluate multiple hypotheses jointly, are often useful for filtering the gene list to a set of interesting features for further exploration while controlling the false discovery rate. Though there are several methods which have been proposed for analyzing correlated RNA-sequencing data, there has been little research evaluating and comparing the performance of multiple degree of freedom tests across methods. METHODS: We evaluated 11 different methods for modelling correlated RNA-sequencing data by performing a simulation study to compare the false discovery rate, power, and model convergence rate across several hypothesis tests and sample size scenarios. We also applied each method to a real longitudinal RNA-sequencing dataset. RESULTS: Linear mixed modelling using transformed data had the best false discovery rate control while maintaining relatively high power. However, this method had high model non-convergence, particularly at small sample sizes. No method had high power at the lowest sample size. We found a mix of conservative and anti-conservative behavior across the other methods, which was influenced by the sample size and the hypothesis being evaluated. The patterns observed in the simulation study were largely replicated in the analysis of a longitudinal study including data from intensive care unit patients experiencing cardiogenic or septic shock. CONCLUSIONS: Multiple degree of freedom testing is a valuable tool in longitudinal and other correlated RNA-sequencing experiments. Of the methods that we investigated, linear mixed modelling had the best overall combination of power and false discovery rate control. Other methods may also be appropriate in some scenarios.

Efficient gene-environment interaction tests for large biobank-scale sequencing studies.

Complex human diseases are affected by genetic and environmental risk factors and their interactions. Gene-environment interaction (GEI) tests for aggregate genetic variant sets have been developed in recent years. However, existing statistical methods become rate limiting for large biobank-scale sequencing studies with correlated samples. We propose efficient Mixed-model Association tests for GEne-Environment interactions (MAGEE), for testing GEI between an aggregate variant set and environmental exposures on quantitative and binary traits in large-scale sequencing studies with related individuals. Joint tests for the aggregate genetic main effects and GEI effects are also developed. A null generalized linear mixed model adjusting for covariates but without any genetic effects is fit only once in a whole genome GEI analysis, thereby vastly reducing the overall computational burden. Score tests for variant sets are performed as a combination of genetic burden and variance component tests by accounting for the genetic main effects using matrix projections. The computational complexity is dramatically reduced in a whole genome GEI analysis, which makes MAGEE scalable to hundreds of thousands of individuals. We applied MAGEE to the exome sequencing data of 41,144 related individuals from the UK Biobank, and the analysis of 18,970 protein coding genes finished within 10.4 CPU hours.

A permutation-based approach for heterogeneous meta-analyses of rare events.

The increasingly widespread use of meta-analysis has led to growing interest in meta-analytic methods for rare events and sparse data. Conventional approaches tend to perform very poorly in such settings. Recent work in this area has provided options for sparse data, but these are still often hampered when heterogeneity across the available studies differs based on treatment group. We propose a permutation-based approach based on conditional logistic regression that accommodates this common contingency, providing more reliable statistical tests when such patterns of heterogeneity are observed. We find that commonly used methods can yield highly inflated Type I error rates, low confidence interval coverage, and bias when events are rare and non-negligible heterogeneity is present. Our method often produces much lower Type I error rates and higher confidence interval coverage than traditional methods in these circumstances. We illustrate the utility of our method by comparing it to several other methods via a simulation study and analyzing an example data set, which assess the use of antibiotics to prevent acute rheumatic fever.

Frailty modeling under a selective sampling protocol: an application to type 1 diabetes related autoantibodies.

In studies following selective sampling protocols for secondary outcomes, conventional analyses regarding their appearance could provide misguided information. In the large type 1 diabetes prevention and prediction (DIPP) cohort study monitoring type 1 diabetes-associated autoantibodies, we propose to model their appearance via a multivariate frailty model, which incorporates a correlation component that is important for unbiased estimation of the baseline hazards under the selective sampling mechanism. As further advantages, the frailty model allows for systematic evaluation of the association and the differences in regression parameters among the autoantibodies. We demonstrate the properties of the model by a simulation study and the analysis of the autoantibodies and their association with background factors in the DIPP study, in which we found that high genetic risk is associated with the appearance of all the autoantibodies, whereas the association with sex and urban municipality was evident for IA-2A and IAA autoantibodies.

Fitting marginal models in small samples: A simulation study of marginalized multilevel models and generalized estimating equations.

In correlated data settings, analysts typically choose between fitting conditional and marginal models, whose parameters come with distinct interpretations, and as such the choice between the two should be made on scientific grounds. For settings where interest lies in marginal-or population-averaged-parameters, the question of how best to estimate those parameters is a statistical one, and analysts have at their disposal two distinct modeling frameworks: generalized estimating equations (GEE) and marginalized multilevel models (MMMs). The two have been contrasted theoretically and in large sample settings, but asymptotic theory provides no guarantees in the small sample settings that are commonplace. In a comprehensive series of simulation studies, we shed light on the relative performance of GEE and MMMs in small-sample settings to help guide analysis decisions in practice. We find that both GEE and MMMs exhibit similar small-sample bias when the correct correlation structure is adopted (ie, when the random effects distribution is correctly specified or moderately misspecified)-but MMMs can be sensitive to misspecification of the correlation structure. When there are a small number of clusters, MMMs only slightly underestimate standard errors (SEs) for within-cluster associations but can severely underestimate SEs for between-cluster associations. By contrast, while GEE severely underestimates SEs, the Mancl and DeRouen correction provides approximately valid inference.

Fast QLB algorithm and hypothesis tests in logistic model for ophthalmologic bilateral correlated data.

In ophthalmologic or otolaryngologic studies, bilateral correlated data often arise when observations involving paired organs (e.g., eyes, ears) are measured from each subject. Based on Donner's model , in this paper, we focus on investigating the relationship between the disease probability and covariates (such as ages, weights, gender, and so on) via the logistic regression for the analysis of bilateral correlated data. We first propose a new minorization-maximization (MM) algorithm and a fast quadratic lower bound (QLB) algorithm to calculate the maximum likelihood estimates of the vector of regression coefficients, and then develop three large-sample tests (i.e., the likelihood ratio test, Wald test, and score test) to test if covariates have a significant impact on the disease probability. Simulation studies are conducted to evaluate the performance of the proposed fast QLB algorithm and three testing methods. A real ophthalmologic data set in Iran is used to illustrate the proposed methods.

Clinical risk prediction models and informative cluster size: Assessing the performance of a suicide risk prediction algorithm.

Clinical visit data are clustered within people, which complicates prediction modeling. Cluster size is often informative because people receiving more care are less healthy and at higher risk of poor outcomes. We used data from seven health systems on 1,518,968 outpatient mental health visits from January 1, 2012 to June 30, 2015 to predict suicide attempt within 90 days. We evaluated true performance of prediction models using a prospective validation set of 4,286,495 visits from October 1, 2015 to September 30, 2017. We examined dividing clustered data on the person or visit level for model training and cross-validation and considered a within cluster resampling approach for model estimation. We evaluated optimism by comparing estimated performance from a left-out testing dataset to performance in the prospective dataset. We used two prediction methods, logistic regression with least absolute shrinkage and selection operator (LASSO) and random forest. The random forest model using a visit-level split for model training and testing was optimistic; it overestimated discrimination (area under the curve, AUC = 0.95 in testing versus 0.84 in prospective validation) and classification accuracy (sensitivity = 0.48 in testing versus 0.19 in prospective validation, 95th percentile cut-off). Logistic regression and random forest models using a person-level split performed well, accurately estimating prospective discrimination and classification: estimated AUCs ranged from 0.85 to 0.87 in testing versus 0.85 in prospective validation, and sensitivity ranged from 0.15 to 0.20 in testing versus 0.17 to 0.19 in prospective validation. Within cluster resampling did not improve performance. We recommend dividing clustered data on the person level, rather than visit level, to ensure strong performance in prospective use and accurate estimation of future performance at the time of model development.

MCMSeq: Bayesian hierarchical modeling of clustered and repeated measures RNA sequencing experiments.

BACKGROUND: As the barriers to incorporating RNA sequencing (RNA-Seq) into biomedical studies continue to decrease, the complexity and size of RNA-Seq experiments are rapidly growing. Paired, longitudinal, and other correlated designs are becoming commonplace, and these studies offer immense potential for understanding how transcriptional changes within an individual over time differ depending on treatment or environmental conditions. While several methods have been proposed for dealing with repeated measures within RNA-Seq analyses, they are either restricted to handling only paired measurements, can only test for differences between two groups, and/or have issues with maintaining nominal false positive and false discovery rates. In this work, we propose a Bayesian hierarchical negative binomial generalized linear mixed model framework that can flexibly model RNA-Seq counts from studies with arbitrarily many repeated observations, can include covariates, and also maintains nominal false positive and false discovery rates in its posterior inference. RESULTS: In simulation studies, we showed that our proposed method (MCMSeq) best combines high statistical power (i.e. sensitivity or recall) with maintenance of nominal false positive and false discovery rates compared the other available strategies, especially at the smaller sample sizes investigated. This behavior was then replicated in an application to real RNA-Seq data where MCMSeq was able to find previously reported genes associated with tuberculosis infection in a cohort with longitudinal measurements. CONCLUSIONS: Failing to account for repeated measurements when analyzing RNA-Seq experiments can result in significantly inflated false positive and false discovery rates. Of the methods we investigated, whether they model RNA-Seq counts directly or worked on transformed values, the Bayesian hierarchical model implemented in the mcmseq R package (available at https://github.com/stop-pre16/mcmseq ) best combined sensitivity and nominal error rate control.

Critique of the pairwise method for estimating qPCR amplification efficiency: beware of correlated data!

BACKGROUND: A recently proposed method for estimating qPCR amplification efficiency E analyzes fluorescence intensity ratios from pairs of points deemed to lie in the exponential growth region on the amplification curves for all reactions in a dilution series. This method suffers from a serious problem: The resulting ratios are highly correlated, as they involve multiple use of the raw data, for example, yielding ~ 250 E estimates from ~ 25 intensity readings. The resulting statistics for such estimates are falsely optimistic in their assessment of the estimation precision. RESULTS: Monte Carlo simulations confirm that the correlated pairs method yields precision estimates that are better than actual by a factor of two or more. This result is further supported by estimating E by both pairwise and Cq calibration methods for the 16 replicate datasets from the critiqued work, and then comparing the ensemble statistics for these methods. CONCLUSION: Contrary to assertions in the proposing work, the pairwise method does not yield E estimates a factor of 2 more precise than estimates from Cq calibration fitting (the standard curve method). On the other hand, the statistically correct direct fit of the data to the model behind the pairwise method can yield E estimates of comparable precision. Ways in which the approach might be improved are discussed briefly.

Methods for generalized change-point models: with applications to human immunodeficiency virus surveillance and diabetes data.

In many epidemiological and biomedical studies, the association between a response variable and some covariates of interest may change at one or several thresholds of the covariates. Change-point models are suitable for investigating the relationship between the response and covariates in such situations. We present change-point models, with at least one unknown change-point occurring with respect to some covariates of a generalized linear model for independent or correlated data. We develop methods for the estimation of the model parameters and investigate their finite-sample performances in simulations. We apply the proposed methods to examine the trends in the reported estimates of the annual percentage of new human immunodeficiency virus (HIV) diagnoses linked to HIV-related medical care within 3 months after diagnosis using HIV surveillance data from the HIV prevention trial network 065 study. We also apply our methods to a dataset from the Pima Indian diabetes study to examine the effects of age and body mass index on the risk of being diagnosed with type 2 diabetes.

Use of partitioned GMM marginal regression model with time-dependent covariates: analysis of Chinese Longitudinal Healthy Longevity Study.

BACKGROUND: Elderly population's health is a major concern for most industrial nations. National health surveys provide a measure of the state of elderly health. One such survey is the Chinese Longitudinal Healthy Longevity Survey. It collects data on risk factors and outcomes on the elderly. We examine these longitudinal survey data to determine the changes in health and to identify risk factors as they impact health outcomes including the elderly's ability to do a physical check. METHODS: We use a Partitioned GMM logistic regression model to identify risk factors. The model also accounts for the correlation between lagged time-dependent covariates and the outcomes. It addresses present and past measures of time-dependent covariates on simultaneous outcomes. The relation produces additional regression coefficients as byproduct of the Partitioned model, identifying the immediate, delayed effects (lag - 1), further delayed (lag-2), etc. Therefore, the model presents the opportunity for decision makers to monitor the covariate over time. This technique is particularly useful in healthcare and health related research. We use the Chinese Longitudinal Health Longevity Survey data to identify those risk factors and to display the utility of the model. RESULTS: We found that one's ability to make own decisions, frequently consuming vegetables, exercise frequently, one's ability to transfer without assistance, having visual difficulties and being able to pick book from floor while standing had varying effects of significance on one's health and ability to complete physical checks as they get older. CONCLUSIONS: The partitioning of the covariates as immediate effect, delayed effect or further delayed effect are important measures in a declining population.

Group penalized generalized estimating equation for correlated event-related potentials and biomarker selection.

BACKGROUND: Event-related potentials (ERP) data are widely used in brain studies that measure brain responses to specific stimuli using electroencephalogram (EEG) with multiple electrodes. Previous ERP data analyses haven't accounted for the structured correlation among observations in ERP data from multiple electrodes, and therefore ignored the electrode-specific information and variation among the electrodes on the scalp. Our objective was to evaluate the impact of early adversity on brain connectivity by identifying risk factors and early-stage biomarkers associated with the ERP responses while properly accounting for structured correlation. METHODS: In this study, we extend a penalized generalized estimating equation (PGEE) method to accommodate structured correlation of ERPs that accounts for electrode-specific data and to enable group selection, such that grouped covariates can be evaluated together for their association with brain development in a birth cohort of urban-dwelling Bangladeshi children. The primary ERP responses of interest in our study are N290 amplitude and the difference in N290 amplitude. RESULTS: The selected early-stage biomarkers associated with the N290 responses are representatives of enteric inflammation (days of diarrhea, MIP1b, retinol binding protein (RBP), Zinc, myeloperoxidase (MPO), calprotectin, and neopterin), systemic inflammation (IL-5, IL-10, ferritin, C Reactive Protein (CRP)), socioeconomic status (household expenditure), maternal health (mother height) and sanitation (water treatment). CONCLUSIONS: Our proposed group penalized GEE estimator with structured correlation matrix can properly model the complex ERP data and simultaneously identify informative biomarkers associated with such brain connectivity. The selected early-stage biomarkers offer a potential explanation for the adversity of neurocognitive development in low-income countries and facilitate early identification of infants at risk, as well as potential pathways for intervention. TRIAL REGISTRATION: The related clinical study was retrospectively registered with https://doi.org/ClinicalTrials.gov , identifier NCT01375647, on June 3, 2011.

Bayesian frailty modeling of correlated survival data with application to under-five mortality.

BACKGROUND: There is high rate of under-five mortality in West Africa with little effort made to study determinants that significantly increase or decrease its risk across the West African sub-region. This is important since it will help in the design of effective intervention programs for each country or the entire region. The overall objective of this research evaluates the determinants of under-five mortality prior to the end of the 2015 Millennium Development Goals, to guide West African countries implement strategies that will aid them achieve the Sustainable Development Goal 3 by 2030. METHOD: This study used the Demographic and Health Survey (DHS) data from twelve (12) out of the eighteen West African countries; Ghana, Benin, Cote d' Ivoire, Guinea, Liberia, Mali, Niger, Nigeria, Sierra Leone, Burkina Faso, Gambia and Togo. Data were extracted from the children and women of reproductive age files as provided in the DHS report. The response or outcome variable of interest is under-five mortality rate. A Bayesian exponential, Weibull and Gompertz regression models via a gamma shared frailty model were used for the analysis. The deviance information criteria and Bayes factors were used to discriminate between models. These analyses were carried out using Stata version 15 software. RESULTS: The study recorded 101 (95% CI: 98.6-103.5) deaths per 1000 live births occurring among the twelve countries. Burkina Faso (124.4), Cote D'lvoire (110.1), Guinea (116.4), Nigeria (120.6) and Niger (118.3) recorded the highest child under-5 mortality rate. Gambia (48.1), Ghana (60.1) and Benin (70.4) recorded the least unde-5 mortality rate per 1000 livebirths. Multiple birth children were about two times more likely to die compared to singleton birth, in all except Gambia, Nigeria and Sierra Leone. We observed significantly higher hazard rates for male compared to female children in the combined data analysis (HR: 1.14, 95% CI: [1.10-1.18]). The country specific analysis in Benin, Cote D'lvoire, Guinea, Liberia, Mali and Nigeria showed higher under-5 mortality hazard rates among male children compared to female children whilst Niger was the only country to report significantly lower hazard rate of males compared to females. CONCLUSION: There is still quite a substantial amount of work to be done in order to meet the Sustainable Development Goal 3 in 2030 in West Africa. There exist variant differences among some of the countries with respect to mortality rates and determinants which require different interventions and policy decisions.