results = pcr_model.fit()
n_points = pca_df.shape[0]
y_output = pca_df['childHeight'].values.reshape(n_points, 1)
print(results.summary())