Skip to content

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).

time_to_first_vs_pma_at_full_oral_feeding_unadjusted.png

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).

time_to_first_vs_pma_at_full_oral_feeding_adjusted.png

Relationship Analysis

Adjusted plot shows relationship.