Study design and setting
This cross-sectional study utilized data collected in the Simiyu region from May 20th to June 3rd, 2021 by Amref Health Africa in Tanzania. The population of the Simiyu region according to the 2021 projection is approximately 2,418,495 with a growth rate of 5.0,(21). Agriculture is the most dominant economic activity in the region. Other economic activities include fishing, trade and commerce. The region has reported excessively worse situations for maternal and child health services coverage like postnatal check-ups (10%), facility delivery (40%), and newborns who received postnatal check-ups within two days (15%). Moreover among children under five years, 33.3% and 5% are reported to be stunted and wasted, respectively, and only 15% of children were exclusively breastfed in the region (10).
Sample Size And Sampling
A multi-stage cluster sampling was used. In the first stage, Enumeration Areas (EAs) in the Simiyu region from the 2012 Tanzania population and housing census were used as the sampling frame. A random number of clusters to be included were selected, and a total of 51 clusters were randomly selected for inclusion out of 67 clusters. From sampled clusters, enumeration of all households and their residents was carried out by a team of enumerators and mappers from the National Bureau of Statistics (NBS) to draw up a list of all the households in each cluster. In stage two, the number of households to be included in each cluster was selected, yielding a total 2035 of households. All women and men aged 15–49 residing within the selected households, who consented were interviewed; a total of 2020 women were interviewed. This study includes women of reproductive age (15–49 years) who gave birth two years preceding the survey. Out of 2020 women, 754 women gave birth less than two years before the survey. Among them 72 resorted to never breastfeeding their children and 13 had missing information on EIBF. As a result, 669 women were analysed (Fig. 1).
Data Collection Methods And Tools
An interviewer-administered questionnaire used for data collection was developed based on questions from the Tanzania 2015/16 Demographic and Health Survey (DHS). The questionnaires were translated into Swahili and installed on an android tablet using Open Data Kit (ODK) for data collection. The questionnaire captured information on pregnancy, deliveries, post-delivery care, family planning, sexual activities/behaviours, and decision-making on various issues such as healthcare utilization and management of household income. Information on EIBF was captured in the questionnaire. Data collection was done by a team of twenty research assistants (RAs). The head of the household provided permission to interview all eligible women available in the household.
Study Variables
The dependent variable was EIBF. EIBF refers to the initiation of breast milk feeding within one hour after delivery (11, 13, 14, 22) and was classified as “1” if breastfeeding initiation was within one hour and “0” if otherwise.
Independent variables included maternal and child demographic characteristics and reproductive and maternal health services. Maternal demographic characteristics included maternal age (15–19, 20–24, 25–29 and 30 + years), education level (non-formal education, primary, secondary,), marital status (single, married, cohabiting, divorced/ widowed), maternal employment (employed, not employed) and district of residence (Bariadi, Itimila, Busega, Maswa, Meatu). Child characteristics include child’s sex (male, female). Reproductive and maternal health services included; number of antenatal visits (< 4, 4 visits and more), place of delivery (health facility, home/others), mode of delivery (vaginal, cesarean section), placing of the baby on mother’s chest/abdomen after delivery (skin to skin contact) (yes, no), level of health facility of delivery (health centre, dispensary, hospital), birth attendant (health professionals, traditional birth, and friend/relative/ others).
Data Management And Analysis
For data cleaning and analysis, STATA software version 15 was used. Continuous variables were summarized using mean and standard deviation, while categorical variables were summarized into frequencies and percentages. To compare the proportion of EIBF between groups of independent variables, a cross-tabulation of EIBF against each independent factor along with a Chi-squared test was done. Since the study used data collected by a cluster sampling technique, the observations are dependent at several hierarchical levels, thus using a classical logistic regression may not be appropriate in analysing this data (23). Therefore, a multilevel logistic regression analysis model was used. To determine whether multilevel regression was required for the data set, a null model (a model without exposure variables) was first fitted. The intra-class correlation (ICC) a measure of the degree of variability between clusters (24), indicates that there was a significant cluster difference.
Moreover, two models were fitted: model 1 was a two-level random effects model with clusters (enumeration areas) as level 2. Model 2 was a three-level model with households as level 2 and clusters (enumeration areas) as level 3. In model 2 the ICC indicated that there was no significant difference between the households. As a result, the log odds of EIBF were best modelled using model 1 a two-level logistic regression with the following equation.
$$log\left[\frac{{\pi }_{ij}}{1-{\pi }_{ij}}\right]= {[\beta }_{o}+{\mu }_{j}]+ \sum _{k=1}^{n}{\beta }_{k}{x}_{ikj}+ {e}_{ij}$$
Where; \({{\beta }}_{{k}}\) Are the fixed coefficients, \({{\beta }}_{{o}}\)is the intercept-the odds of EIBF at the baseline levels of the explanatory variables, \({{\mu }}_{{j} }\)is the random effect for the jth cluster (enumeration areas, that is the effect of the enumeration areas on EIBF), \({{e}}_{{i}{j}}\)is random errors at the individual levels, i and j are the level 1 (women of reproductive age (15–49 years) and level 2 clusters (enumeration areas) units, respectively, n is the number of explanatory variables and\({{X}}_{{k}}\) is a set of explanatory variables both individual and community-level variables.
A bivariate analysis was performed to identify variables that are associated with EIBF. Variables with p-value < 0.10 analysis and those identified in the literature as potential confounders were considered for multivariable analysis. Akaike Information Criteria (AIC) together with the ROC analysis/curve were used for the final model selection. In multivariable analysis variables with a p-value of ≤ 0.05 were considered statistically significant.