仮説の複雑性を算出する関数を作ってみた

情報仮説に基づく仮説の複雑性を算出する関数を作ってみました。ここではとりあえずコードだけですが公開しておきます。

詳しくは11月28日に開催される第3回 R勉強会@比治山大(広島) #HijiyamaRの発表でふれていきます(https://atnd.org/events/68737)。情報仮説のことについては,あとで記事にしようと思います。

基本的には,v_numには仮説に含まれる変数の数を,ie_numには仮説に含まれる不等号の数を,e_numには仮説に含まれる等号の数を入力すれば,仮説の複雑性を評価します。以下の使用例の場合,3つの変数があって,不等号は1つ,等号は1つというようになります(e.g., H1: mu1 < (mu2 = mu3))。

このコードを走らせるとこんな結果が返ってきます。 複雑性結果

差がないと考える変数について,積極的に同じであるかを検討したい場合はnullの結果を,特に同じであるという仮定をしない場合はunconstraintの結果を仮説の複雑性として使うとよいでしょう。Hoijtink, Klugkist, Boelen (2010)のハンドブックの中では,両方やってみて比較したらよいようなことが書いてあったと思います。

Leave a Comment

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です