Secondary Analysis: time_to_first vs pma_at_full_oral_feeding¶
Analysis of time_to_first vs pma_at_full_oral_feeding.
Unadjusted Model Results¶
Simple linear regression examining the direct relationship between time_to_first and pma_at_full_oral_feeding.
Model diagnostics warning: High multicollinearity detected (condition number > 30). Residuals may not be normally distributed (p < 0.05).
Outcome: pma_at_full_oral_feeding Predictor: time_to_first
OLS Regression Results
====================================================================================
Dep. Variable: pma_at_full_oral_feeding R-squared: 0.363
Model: OLS Adj. R-squared: 0.360
Method: Least Squares F-statistic: 128.9
Date: Tue, 05 Aug 2025 Prob (F-statistic): 6.26e-24
Time: 04:28:00 Log-Likelihood: -470.98
No. Observations: 228 AIC: 946.0
Df Residuals: 226 BIC: 952.8
Df Model: 1
Covariance Type: nonrobust
=================================================================================
coef std err t P>|t| [0.025 0.975]
---------------------------------------------------------------------------------
const 33.7417 0.284 118.911 0.000 33.183 34.301
time_to_first 0.0704 0.006 11.352 0.000 0.058 0.083
==============================================================================
Omnibus: 89.496 Durbin-Watson: 1.813
Prob(Omnibus): 0.000 Jarque-Bera (JB): 333.155
Skew: 1.602 Prob(JB): 4.53e-73
Kurtosis: 7.980 Cond. No. 102.
==============================================================================
Notes:
[1] Standard Errors assume that the covariance matrix of the errors is correctly specified.
Unadjusted Plot¶
Plot showing relationship between time_to_first and pma_at_full_oral_feeding (LOESS with CI and OLS prediction).

Model Variables (Adjusted Analysis)¶
Outcome: pma_at_full_oral_feeding Indep: time_to_first Gestational Age Var: gestational_age_weeks Confounders: ['gestational_age_weeks', 'birth_weight_grams', 'mechanical_ventilation', 'o2_device_at_36_weeks'] Covariates: ['medicalcomplexity']
Adjusted Model Results¶
Multiple regression examining the relationship between time_to_first and pma_at_full_oral_feeding, controlling for gestational age and other potential confounders.
Model diagnostics warning: High multicollinearity detected (condition number > 30). Residuals may not be normally distributed (p < 0.05).
Outcome: pma_at_full_oral_feeding Predictors: time_to_first, gestational_age_weeks, birth_weight_grams, mechanical_ventilation_no history of mechanical ventilation, o2_device_at_36_weeks_did NOT require respiratory support at 36 weeks, medicalcomplexity_Children with noncomplex chronic disease, medicalcomplexity_Children without chronic disease
OLS Regression Results
====================================================================================
Dep. Variable: pma_at_full_oral_feeding R-squared: 0.470
Model: OLS Adj. R-squared: 0.454
Method: Least Squares F-statistic: 27.92
Date: Tue, 05 Aug 2025 Prob (F-statistic): 2.65e-27
Time: 04:28:02 Log-Likelihood: -449.95
No. Observations: 228 AIC: 915.9
Df Residuals: 220 BIC: 943.3
Df Model: 7
Covariance Type: nonrobust
=========================================================================================================================================
coef std err t P>|t| [0.025 0.975]
-----------------------------------------------------------------------------------------------------------------------------------------
const 12.9174 3.840 3.364 0.001 5.349 20.486
time_to_first 0.1162 0.013 9.192 0.000 0.091 0.141
gestational_age_weeks 0.7190 0.133 5.394 0.000 0.456 0.982
birth_weight_grams -0.0005 0.001 -0.776 0.439 -0.002 0.001
mechanical_ventilation_no history of mechanical ventilation -0.7477 0.292 -2.559 0.011 -1.324 -0.172
o2_device_at_36_weeks_did NOT require respiratory support at 36 weeks -0.3311 0.309 -1.073 0.285 -0.939 0.277
medicalcomplexity_Children with noncomplex chronic disease -0.1569 0.711 -0.221 0.825 -1.558 1.244
medicalcomplexity_Children without chronic disease -0.3238 0.559 -0.579 0.563 -1.426 0.779
==============================================================================
Omnibus: 85.387 Durbin-Watson: 1.649
Prob(Omnibus): 0.000 Jarque-Bera (JB): 245.761
Skew: 1.645 Prob(JB): 4.30e-54
Kurtosis: 6.878 Cond. No. 4.24e+04
==============================================================================
Notes:
[1] Standard Errors assume that the covariance matrix of the errors is correctly specified.
[2] The condition number is large, 4.24e+04. This might indicate that there are
strong multicollinearity or other numerical problems.
Adjusted Plot¶
Plot showing adjusted relationship between time_to_first and pma_at_full_oral_feeding (LOESS with CI).

Relationship Analysis¶
Adjusted plot shows relationship.