平均は分散と等しいですか？

Poisson モデルでは、平均が分散と等しいという仮定がありました。講義で学んだとおり、この仮定が成り立たない場合は過分散（overdispersion）が生じます。過分散を調整しないと、モデルの標準誤差を誤って解釈してしまいます。

この演習では、まずメスのカブトガニにおけるサテライト（衛星オス）の数の平均と分散を計算します。

crab データセットはワークスペースに読み込まれています。

この演習はコースの一部です

Pythonで学ぶ一般化線形モデル

コースを見る

演習の手順

numpy の mean() を使って、サテライト数の変数 sat の標本平均を sat_mean として計算し、表示します。
numpy の var() を使って、サテライト数の変数 sat の標本分散を sat_var として計算し、表示します。
sat_var を sat_mean で割った比を計算し、表示します。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Compute and print sample mean of the number of satellites: sat_mean
____ = np.____(____.____)

print('Sample mean:', round(____, 3))

# Compute and print sample variance of the number of satellites: sat_var
____ = np.____(____.____)
print('Sample variance:', round(____, 3))

# Compute ratio of variance to mean
print('Ratio:', round(____/____, 3))

コードを編集して実行