1. 学ぶ
    2. /
  2. コース
    3. /
  3. 統計的思考 ケーススタディ
Connected

演習

指数分布に対するK-S検定

パークフィールドの地震シーケンスにおける地震間隔が指数分布に従うという帰無仮説を検定します。つまり、地震は前回の発生時刻に関係なくランダムに起きるという仮説です。注：この計算は処理量が多く（1億を超える乱数を生成します）、完了までに約10秒かかります。

指示

100 XP
  • np.random.exponential() を使って指数分布から10,000件のレプリケートを生成してください。地震の平均時間間隔は mean_time_gap に格納されており、これは前の演習で計算したものです。結果を x_f に保存してください。
  • このサンプル x_f と、time_gap に格納されている実際の時間間隔を使って、dcst.ks_stat() でコルモゴロフ・スミルノフ統計量を計算してください。
  • 前の演習で作成した関数（dcst.draw_ks_reps() として利用可能）を使い、指数分布からK-Sレプリケートを10,000件生成してください。対象の指数分布からのサンプリングには size=10000 キーワード引数を使用してください。レプリケートを reps として保存してください。
  • p値を計算して出力してください。この場合、「少なくとも同程度に極端」とは、帰無仮説のもとでの検定統計量が観測値以上であることを意味します。