差分
このページの2つのバージョン間の差分を表示します。
両方とも前のリビジョン前のリビジョン次のリビジョン | 前のリビジョン次のリビジョン両方とも次のリビジョン | ||
informatics:2011 [2011/06/16 13:43] – watalu | informatics:2011 [2011/06/16 17:54] – [他の先生の課題] watalu | ||
---|---|---|---|
行 1: | 行 1: | ||
+ | |||
+ | |||
+ | ==== 経営情報学コース分の課題 ==== | ||
+ | |||
+ | |西先生|これから40年どんな手段で世の中を良くするのか?技術や考え方をビジョン化し、具体的に説明しなさい| | ||
+ | |山田先生|講義内での課題を、A4一枚ぐらいまで膨らませて| | ||
+ | |水戸先生|標示資料の一番最後のページ| | ||
+ | |由良先生|今日の配付資料のとおり| | ||
+ | |山本|画鋲を30回、投げてみよう ([[http:// | ||
+ | |||
+ | |提出期限|2011年6月30日(木) 0500pmまで| | ||
+ | |提出場所|西五号館3階 総合情報学専攻事務室| | ||
+ | |||
+ | 宜しくお願いします。 | ||
+ | |||
==== 画鋲を投げて針が上向く確率の推定 ==== | ==== 画鋲を投げて針が上向く確率の推定 ==== | ||
行 24: | 行 39: | ||
=== ベイズ法 === | === ベイズ法 === | ||
- | この計算は、WinBUGSというソフトウェアを用いるのが簡単である。 | + | ベイズ法に関する計算、特に事後分布の算出は、WinBUGSというソフトウェアを用いるのが簡単である。 |
== データの与え方 == | == データの与え方 == | ||
< | < | ||
- | DATA list {n=32, x=c(1, | + | DATA list(n=32, x=c(1, |
</ | </ | ||
行 35: | 行 50: | ||
モデルと事前分布と初期値の指定。 | モデルと事前分布と初期値の指定。 | ||
< | < | ||
- | MODEL model { | + | MODEL model{ |
for ( i in 1:n ) { | for ( i in 1:n ) { | ||
x[i] ~ dbin(p,1) | x[i] ~ dbin(p,1) | ||
行 47: | 行 62: | ||
ここでベータ分布はパラメータを変えると、次のように動く。 | ここでベータ分布はパラメータを変えると、次のように動く。 | ||
- | {{: | + | {{: |
- | + | ここから先の手順は、多少複雑で[[http:// | |
- | ここから先の手順は、多少複雑で | + | |
行 68: | 行 82: | ||
- 履歴(history)、直近の履歴(trace)、密度関数(density)、分位点の推移(quantiles)、などを眺めて、定常分布となっていることを確認 (〔Sample Monitor Tool〕の〔node〕欄に ``*'' | - 履歴(history)、直近の履歴(trace)、密度関数(density)、分位点の推移(quantiles)、などを眺めて、定常分布となっていることを確認 (〔Sample Monitor Tool〕の〔node〕欄に ``*'' | ||
- 定常と判断した範囲(beg〜end)までについてstatsを表示させる (〔Sample Monitor Tool〕でbegの値を1ではなく、定常と判断した歩数に変更) | - 定常と判断した範囲(beg〜end)までについてstatsを表示させる (〔Sample Monitor Tool〕でbegの値を1ではなく、定常と判断した歩数に変更) | ||
- | | ||
- | 詳細は[[http:// | ||
=== ブートストラップ法 === | === ブートストラップ法 === | ||
行 93: | 行 105: | ||
これは点推定の方法は正しくても、分布の仮定がそれほど正しくないかもしれない時に、一番最初の信頼区間とは異なった範囲を示す。 | これは点推定の方法は正しくても、分布の仮定がそれほど正しくないかもしれない時に、一番最初の信頼区間とは異なった範囲を示す。 | ||
- | {{: | + | {{: |
ただし、ブートスラップ反復回数 B の設定には要注意。 | ただし、ブートスラップ反復回数 B の設定には要注意。 | ||
+ | === ベータ分布の図 === | ||
+ | < | ||
+ | # Beta(1,1) | ||
+ | plot(c(1: | ||
+ | # Beta(3,2) | ||
+ | lines(c(1: | ||
+ | # Beta(2,3) | ||
+ | lines(c(1: | ||
+ | # Beta(10,5) | ||
+ | lines(c(1: | ||
+ | # Beta(5,10) | ||
+ | lines(c(1: | ||
+ | # Beta(15,5) | ||
+ | lines(c(1: | ||
+ | # Beta(5,15) | ||
+ | lines(c(1: | ||
+ | |||
+ | </ |