宇宙論における線形モデル

100年も経たない昔、宇宙は静止した1つの銀河（おそらく数百万個の恒星を含む）だけで構成されていると考えられていました。現在では、数千億もの銀河が観測されており、それぞれが数千億の恒星を持ち、すべてが運動しています。

現代宇宙論のはじまりは、線形モデルの利用を含んだ1929年のEdwin Hubbleによる発表にあります。

この演習では、銀河の速度を地球からの距離の線形関数として表す定数であるハッブル定数（Hubble's Constant）を、傾きとして与えるモデルを構築します。

Hubbleのプロット

事前に読み込まれた names、distances、velocities の列を持つ DataFrame を使います。
formula="velocities ~ distances"、data=df を指定して、ols().fit() でモデルを構築して当てはめます。
model_fit.params を使って切片と傾きの推定値を取り出し、それぞれ a0 と a1 に代入します。
同様に、対応する不確かさについても model_fit.bse を使って値を取り出し、同じ手順で取得します。