Statistical analysis
REGRESSION
Analyses were performed using R v4.1 [We tested for associations between PPQ and EPDS scores (log-transformed) and intervention group using a linear regression model. We first tested for main effects of the intervention. We adjusted for additional covariates using a hybrid approach, forcing in site, and then using a backwards stepwise regression to identify demographic and clinical covariates that contributed p < .1 to the final model.We then tested for moderating effects on the intervention with four factors: infant gestational age at birth, infant chronological age at study enrollment (log transformed), whether the infant was discharged on a respiratory or feeding device, and NICU-related stress level at study enrollment (log-transformed), controlling for additional covariates as described above. Finally, for factors that had an interaction with the intervention group (denoted intervention moderators) with p < .05, we conducted an exploratory per protocol analysis testing for an interaction between the moderator and each of the following mFICare intervention components over the course of the NICU hospitalization: whether the parent had a mentor; participated in at least 2 clinical rounds; attended at least 1 parent class; or logged into the We3health™ app at least 4 times. We report results meeting a nominal p < .05, as none reached Bonferroni-corrected significance.