pcr_model = sm.ols(formula = 'childHeight ~ pc1 + pc2', data=pca_df)