Objective This study examined the prevalence and predictors of maternal and newborn skin-to-skin contact at birth in Papua New Guinea. Design Data for the study was extracted from the 2016-18 Papua New Guinea Demographic and Health Survey. We included 6,044 women with birth history before the survey in the analysis. Percentages were used to summarise the prevalence of maternal and newborn skin-to-skin contact. A multivariable multilevel binary logistic regression was adopted to examine the predictors of maternal and newborn skin-to-skin contact. The results were presented using adjusted ORs (aORs), with their respective 95% confidence intervals (CIs). Statistical significance was set at p<0.05. Setting The study was conducted in Papua New Guinea. Participant Mothers with children under 5 years. Outcome measures Mother and newborn skin-to-skin contact. Results The prevalence of mother and newborn skin-to-skin contact was 45.2% (95% CI=42.4 to 48.0). The odds of mother and newborn skin-to-skin contact was higher among women with primary education (aOR=1.38; 95% CI=1.03 to 1.83), women with four or more antenatal care attendance (aOR=1.27; 95% CI=1.01 to 1.61), those who delivered at the health facility (aOR=1.27; 95% CI=1.01 to 1.61), and women from communities with high socioeconomic status (aOR=1.45; 95% CI=1.11 to 1.90). Conclusion The study has demonstrated that the prevalence of mother and newborn skin-to-skin contact in Papua New Guinea is low. Factors shown to be associated with mother and newborn skin-to-skin contact were maternal level of education, antenatal care attendance, health facility delivery, and community socioeconomic status. A concerted effort should be placed in improving maternal health service utilisation such as antenatal care attendance and skilled birth delivery, which subsequently lead to the practice of skin-to-skin contact. Also, women should be empowered through education as it has positive impact on their socioeconomic status and health service utilisation.
We performed a secondary analysis of data from the 2016–2018 Papua New Guinea Demographic and Health Survey (DHS). The data for the study was extracted from the Kid’s recode file (KR File) of the DHS. The DHS is a nationally representative survey conducted in over 85 countries worldwide.15 The survey captures data on men, women, and child indicators including SSC.15 The DHS employed a cross-sectional design. Standardised structured interviewer-administered questionnaires were used to collect the data from the respondents. A stratified two-stage cluster sampling design was used to recruit the samples for the survey. In the first stage, clusters were selected using a probability proportional to size sampling technique. In the second stage, a predetermined number of households (usually 28–30) were selected using a systematic sampling technique. The detailed study methodology can be found in the DHS report.16 We included 6,044 women with birth history before the survey who had complete data on all variables of interest in the study. The dataset can be assessed freely at https://dhsprogram.com/data/dataset/Papua-New-Guinea_Standard-DHS_2017.cfm?flag=1.17 Mother and newborn SSC was the outcome variable in the present study. This variable was assessed using the question ‘Was child put on mother’s chest and bare skin after birth?’. With this question, the response options were 0=no; 1=put on chest, touching bare skin; 2=put on chest, no touching of bare skin; 3=put on chest, do not know/missing on touching on bare skin and 8=do not know. For this study’s purpose and with reference to literature,18–20 the response options were further recoded into ‘1=practiced SSC’ for women who response category was “put on chest, touching bare skin” whilst the remaining response options were categorised as ‘0=not practiced SSC’. We included 20 explanatory variables in the study. We selected the variables based on their availability in the DHS dataset as well as their significant association with mother and newborn SSC from literature.18–22 The variables were grouped into individual level and community level. The individual level variables consisted of sex of child, birth order, birth weight, caesarean delivery, type of birth, mother’s age, maternal educational level, marital status, current working status, number of antenatal care visits, place of delivery, health insurance coverage, exposure to watching television, exposure to listening to radio, exposure to reading newspaper or magazine and wealth index. We maintained the existing coding as found in the DHS for sex of child, type of birth, mothers age, caesarean delivery, health insurance coverage, and wealth status. The remaining individual-level variables were coded as birth order (first, second, third, fourth, and fifth or more); birth weight (normal and low birth weight); maternal educational level (no education, primary, and secondary or higher), marital status (never married, married, cohabiting, and previously married); number of antenatal care visits (below four visits and four or more visits); place of delivery (home, health facility, and other); exposure to watching television (no and yes); exposure to listening to radio (no and yes) and exposure to reading newspaper or magazine (no and yes). The DHS devised a wealth index as a proxy measure of socioeconomic position. It was calculated using component rankings derived from principal component analysis on family asset ownership, such as access to drinking water, kind of toilet, type of cooking fuel and possession of a television and refrigerator. The community level variables consisted of place of residence (urban and rural), region (Southern, Highlands, Momase, and Islands), community literacy level (low, medium, and high) and community socioeconomic status (low, medium, and high). We performed the statistical analyses using Stata software V.16.0 (Stata, College Station, Texas, USA). The extracted data was cleaned and all the missing observations were dropped while subcategories of variables with small observations were merged. Percentages were used to present the prevalence of mother and newborn SSC. Later, we examined the distribution of mother and newborn SSC across the explanatory variables using a cross-tabulation. We adopted a binary logistic regression to select significant variables for the multivariable multilevel logistic regression. All the variables that had a p value <0.05 were considered statistically significant and included in the multilevel regression model. Four models of the multilevel regression analysis were built to examine the predictors of mother and newborn SSC. Model O (empty model) was built to examine the variation of the outcome variable (maternal and newborn SSC) attributed to the clustering of the primary sample units. Models I and II were fitted to include variables at the individual and community levels, respectively. The last model (model III) was fitted to include all the statistically significant explanatory variables from the binary logistic regression. The result of the multilevel binary logistic regression analysis was presented using the adjusted ORs (aORs), with their corresponding 95% confidence intervals (CIs). We used the Akaike Information Criterion (AIC) to assess the fitness of each model and for comparing the fitness across the models. All the analyses were weighted. The Stata command ‘svyset’ was employed in all the analyses to adjust for over-and-under sampling, non-response and to improve the generalisability of the findings. In writing the manuscript, we followed the guidelines from the Strengthening the Reporting of Observational Studies in Epidemiology statement (online supplemental table S1).23 bmjopen-2022-062422supp001.pdf Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.