両方とも前のリビジョン前のリビジョン次のリビジョン | 前のリビジョン次のリビジョン両方とも次のリビジョン |
prob:2011 [2011/06/23 11:32] – [#11] watalu | prob:2011 [2011/08/04 12:56] – [#16 2011.07.28] watalu |
---|
==== 授業計画 ==== | ==== 授業計画 ==== |
| |
^回^テーマ^トピック^実施日^レポート課題^ | ^回^テーマ^トピック^予定日^実施日^レポート課題^ |
|#01|ガイダンス、確率の基礎概念(1)|事象,確率|2011.04.07| | | |#01|ガイダンス、確率の基礎概念(1)|事象,確率|2011.04.07| | | |
|#02|確率の基礎概念(2)|事象,確率|2011.04.14| | | |#02|確率の基礎概念(2)|事象,確率|2011.04.14| | | |
|#03|確率の基礎概念(3)|条件つき確率と独立性,ベイズの定理|2011.04.21| | | |#03|確率の基礎概念(3)|条件つき確率と独立性,ベイズの定理|2011.04.21| | | |
|#04|確率変数と分布関数(1)|確率変数,確率分布,分布関数|2011.04.28|2011.05.11〆切で出題(第1章分)| | |#04|確率変数と分布関数(1)|確率変数,確率分布,分布関数|2011.04.28| |2011.05.11〆切で出題(第1章分)| |
|#05|確率変数と分布関数(2)|確率変数のモーメント,分散|2011.05.12| | | |#05|確率変数と分布関数(2)|確率変数のモーメント,分散|2011.05.12| | | |
|#06|モーメント母関数とその応用(1),離散型確率モデル(1)|モーメント母艦数,ベルヌーイ分布,二項分布 |2011.05.19| | | |#06|モーメント母関数とその応用(1),離散型確率モデル(1)|モーメント母艦数,ベルヌーイ分布,二項分布 |2011.05.19| | | |
|#07|離散型確率モデル(2),連続型確率モデル(1)|幾何分布,負の二項分布,ポアソン分布,指数分布|2011.05.26|2011.06.01〆切で出題予定(第2章分)| | |#07|離散型確率モデル(2),連続型確率モデル(1)|幾何分布,負の二項分布,ポアソン分布,指数分布|2011.05.26| |2011.06.01〆切で出題予定(第2章分)| |
|#08|連続型確率モデル(1)|ポアソン分布と指数分布の関係,正規分布|2011.06.02| | | |#08|連続型確率モデル(1)|ポアソン分布と指数分布の関係,正規分布|2011.06.02|2011.06.23| | |
|<del>#09</del>|<del>中間試験および解説</del>| | | | | |<del>#09</del>|<del>中間試験および解説</del>| | | | | |
|#09|確率ベクトルと分布関数(1)|確率ベクトル,同時分布,周辺分布|2011.06.09| | | |#09|確率ベクトルと分布関数(1)|確率ベクトル,同時分布,周辺分布|2011.06.09|2011.06.23, 2011.06.30| | |
|#10|確率ベクトルと分布関数(2)|確率変数の独立性,同時モーメント,共分散,相関係数|2011.06.16|2011.06.22〆切で出題予定(第5,6章分)| | |#10|確率ベクトルと分布関数(2)|確率変数の独立性,同時モーメント,共分散,相関係数|2011.06.16|2011.06.30, 2011.07.07|2011.06.22〆切で出題予定(第5,6章分)| |
|#11|モーメント母関数とその応用(2)| |2011.06.23| | | |#11|モーメント母関数とその応用(2)| |2011.06.23| | | |
|#12|連続型確率モデル(3)|2変量正規分布|2011.06.30| | | |#12|連続型確率モデル(3)|2変量正規分布|2011.06.30|2011.07.07| | |
|#13|大数の法則| |2011.07.07|2011.07.14〆切で出題予定(第3章分)| | |#13|大数の法則| |2011.07.07|2011.07.14|2011.07.14〆切で出題予定(第3章分)| |
|#14|中心極限定理| |2011.07.14| | | |#14|中心極限定理| |2011.07.14| | | |
|#15|標本分布論| |2011.07.21| | | |#15|標本分布論| |2011.07.21|パス| | |
|#16|期末試験| |試験期間中| | | |#16|期末試験| |試験期間中| | | |
| |
==== #1 2011.04.07 ==== | ==== #1 2011.04.07 ==== |
==== #4 2011.04.28 ==== | ==== #4 2011.04.28 ==== |
| |
[[http://stat.inf.uec.ac.jp/library/prob.2011/prob-4-note-and-quiz-20110428.pdf|課題 #1]] | [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-4-note-and-quiz-20110428.pdf|課題 #1]]とその[[http://stat.inf.uec.ac.jp/library/prob.2011/prob-4-solutions-20110627.pdf|略解]]。 |
| |
* 確率変数 <jsm>X</jsm> | * 確率変数 <jsm>X</jsm> |
</jsmath> | </jsmath> |
です。この訂正、土曜日以降、掲示もします。ご免なさい。 | です。この訂正、土曜日以降、掲示もします。ご免なさい。 |
| |
| ==== #8 ==== |
| |
| ==== #9 ==== |
| |
==== #10 ==== | ==== #10 ==== |
==== #11 ==== | ==== #11 ==== |
| |
[[http://stat.inf.uec.ac.jp/library/prob.2011/prob-a-solutions-20110623-draft.pdf|課題#3の解答・手書き版]] (これからタイプに回します) | [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-a-solutions-20110623-draft.pdf|課題#3の解答・手書き版]] (昨日、別の仕事のトラブル対応に追われて、タイプが間に合いませんでした。これからタイプに回しますが、とりあえず手書きのノートを暫定公開します。) |
| |
| [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-a-solutions-20110629.pdf|課題#3の解答・暫定公開版]] (タイプは終わりましたが、未推敲のため、タイプミスがありそうな気がします。) |
| |
| * 2次元連続確率変数の |
| * 同時累積分布関数と同時密度関数 |
| * 周辺累積分布関数と周辺密度関数 |
| * 条件付き累積分布関数と条件付き密度関数 |
| |
| [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-b-note-20110623.pdf|グラフ]] |
| |
| <code> |
| B1 := plot3d(exp(-(1/2)*x^2)/sqrt(2*Pi)*(exp(-(1/2)*y^2)/sqrt(2*Pi)), x = -3 .. 3, y = -3 .. 3, axes = boxed); |
| B2 := plots[pointplot3d]({seq([3, (1/100)*y, exp(-(1/2)*((1/100)*y)^2)/sqrt(2*Pi)], y = -300 .. 300)}); |
| B3 := plots[pointplot3d]({seq([(1/100)*x, 3, exp(-(1/2)*((1/100)*x)^2)/sqrt(2*Pi)], x = -300 .. 300)}); |
| A1 := plot(eval(exp(-(1/2)*x^2)/sqrt(2*Pi)*(exp(-(1/2)*y^2)/sqrt(2*Pi)), x = -3), y = -3 .. 3, axes = boxed); |
| A2 := plot(eval(exp(-(1/2)*x^2)/sqrt(2*Pi)*(exp(-(1/2)*y^2)/sqrt(2*Pi)), x = -2), y = -3 .. 3, axes = boxed); |
| A3 := plot(eval(exp(-(1/2)*x^2)/sqrt(2*Pi)*(exp(-(1/2)*y^2)/sqrt(2*Pi)), x = -1), y = -3 .. 3, axes = boxed); |
| A4 := plot(eval(exp(-(1/2)*x^2)/sqrt(2*Pi)*(exp(-(1/2)*y^2)/sqrt(2*Pi)), x = 0), y = -3 .. 3, axes = boxed); |
| plots[display]({B1, B2, B3}); |
| plots[display]({A1, A2, A3, A4}) |
| </code> |
| |
| |
| <code> |
| library(mvtnorm) |
| |
| postscript("prob-b-correlated-bivariate-normal-distribution-scatterplots-positive.eps", width=6, height=6) |
| par(mfrow=c(2,2)) |
| par(cex=0.5) |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,0,0,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: 0", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,0.5,0.5,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: 0.5", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,0.8,0.8,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: 0.8", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,0.95,0.95,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: 0.95", xlab="X.1", ylab="X.2") |
| |
| postscript("prob-b-correlated-bivariate-normal-distribution-scatterplots-negative.eps", width=6, height=6) |
| par(mfrow=c(2,2)) |
| par(cex=0.5) |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,0,0,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: 0", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,-0.5,-0.5,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: -0.5", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,-0.8,-0.8,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: -0.8", xlab="X.1", ylab="X.2") |
| plot(rmvnorm(n=3000, mean=c(0,0), sigma=matrix(c(1,-0.95,-0.95,1),ncol=2)),pch=20, |
| main ="Correlation Coefficient: -0.95", xlab="X.1", ylab="X.2") |
| graphics.off() |
| </code> |
| |
| ==== #12 ==== |
| |
| * 同時分布と条件付き分布と周辺分布の関係 |
| * 期待値ベクトル |
| * 分散共分散行列 |
| |
| ==== #13 ==== |
| |
| * 二変量正規分布 |
| * 期待値ベクトル |
| * 分散共分散行列 |
| * 相関係数 |
| * 条件付き分布 |
| * 周辺分布 |
| |
| ==== #14 2011.07.14 ==== |
| |
| [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-e-note-and-quiz-20110714.pdf|課題4]] ([[http://stat.inf.uec.ac.jp/library/prob.2011/prob-e-solutions-20110802.pdf|解答例]] (手書きの汚いノートで済みません。レポートを8/2と8/3の二日間は、西5号館3階の事務室で返却しています。8/4には試験会場に持って行くので、事務室からは引き上げます。), [[http://stat.inf.uec.ac.jp/library/prob.2011/prob-e-solutions-20110802-last-page.pdf|解答例の最後のページ]] (スキャンミスしました。重ね重ね済みません。)) |
| |
| * マルコフの不等式 |
| * チェビシェフの不等式 |
| * 大数の法則 |
| |
| 参考書 |
| |
| * 宮川雅巳(1998)「統計技法」, 共立出版. |
| |
| === 課題4への補足 === |
| == 離散分布の共分散と相関係数 == |
| |
| <jsm>\left(X_1, X_2\right)</jsm>が有限個の値の組み合わせしかとらない場合、それぞれの値の組み合わせをとる確率を <jsm>p\left(x_1, x_2\right)</jsm> と置くと、 |
| <jsmath> |
| E\left[X_1\right] = \sum_{x_1} \sum_{x_2} x_1 p\left(x_1, x_2\right) = \sum_{x_1} x_1 \sum_{x_2} p\left(x_1, x_2\right) = \sum_{x_1} x_1 p_1\left(x_1\right) |
| </jsmath> |
| と |
| <jsmath> |
| E\left[X_2\right] = \sum_{x_1} \sum_{x_2} x_2 p\left(x_1, x_2\right) = \sum_{x_2} x_2 \sum_{x_1} p\left(x_1, x_2\right) = \sum_{x_2} x_2 p_2\left(x_2\right) |
| </jsmath> |
| と |
| <jsmath> |
| E\left[X_1X_2\right] = \sum_{x_1} \sum_{x_2} x_1 x_2 p\left(x_1, x_2\right) |
| </jsmath> |
| から、今回の参考書からの出題で証明させられる |
| <jsmath> |
| Cov\left[X_1, X_2\right] = E\left[X_1X_2\right]-E\left[X_1\right] E\left[X_2\right] |
| </jsmath> |
| から共分散を求めることができる。 |
| |
| 相関係数を計算するにも、 |
| <jsmath> |
| E\left[X_1^2\right] = \sum_{x_1} \sum_{x_2} x_1^2 p\left(x_1, x_2\right) = \sum_{x_1} x_1^2 \sum_{x_2} p\left(x_1, x_2\right) = \sum_{x_1} x_1^2 p_1\left(x_1\right) |
| </jsmath> |
| と |
| <jsmath> |
| E\left[X_2^2\right] = \sum_{x_1} \sum_{x_2} x_2^2 p\left(x_1, x_2\right) = \sum_{x_2} x_2^2 \sum_{x_1} p\left(x_1, x_2\right) = \sum_{x_2} x_2^2 p_2\left(x_2\right) |
| </jsmath> |
| を用いて、周辺分散を |
| <jsmath> |
| V\left[X_1\right] = E\left[X_1^2\right]-\left\{E\left[X_1\right]\right\}^2 |
| </jsmath> |
| と |
| <jsmath> |
| V\left[X_2\right] = E\left[X_2^2\right]-\left\{E\left[X_2\right]\right\}^2 |
| </jsmath> |
| のように求め、これらと先に求めた共分散とを合わせて、 |
| <jsmath> |
| \rho\left[X_1, X_2\right] = \frac{Cov\left[X_1, X_2\right]}{\sqrt{V\left[X_1\right]V\left[X_2\right]}} |
| </jsmath> |
| を得る。この手順が一番、計算間違いしにくいんじゃないかと思う。 |
| |
| == 条件付き期待値 == |
| |
| 定数としての条件付き期待値と確率変数としての条件付き期待値の区別。<jsm>\left(X_1, X_2\right)\sim F\left(x_1, x_2\right)</jsm> とする。 |
| <jsm>X_2</jsm>の定義域を<jsm>\Omega_2</jsm>と置くと、 |
| |
| <jsmath> |
| E\left[\phi\left(X_2\right)|X_1=x_1\right]=E_{X_2|X_1}\left[\phi\left(X_2\right)|X_1=x_1\right]=\int_{x_2\in\Omega_2} \phi\left(v\right)dF_{X_2|X_1}\left(v|x_1\right) |
| </jsmath> |
| |
| 最後の積分は、連続分布の場合には、 |
| |
| <jsmath> |
| \int_{x_2\in\Omega_2} \phi\left(v\right)f_{X_2|X_1}\left(v|x_1\right)dv = \int_{x_2\in\Omega_2} \phi\left(v\right) \frac{f_{X_1,X_2}\left(x_1, v\right)}{f_{X_1}\left(x_1\right)}dv |
| </jsmath> |
| |
| と書ける。離散分布の場合にも同様に、総和記号と条件付き確率の公式(あるいはベイズの定理)を用いて、表せる。いずれにせよ、右辺に大文字は残らないので、これは定数。これを |
| |
| <jsmath> |
| \mu_\phi\left(x_1\right)=E\left[\phi\left(X_2\right)|X_1=x_1\right] |
| </jsmath> |
| |
| と置く。 |
| |
| 次に、 |
| |
| <jsmath> |
| E\left[\phi\left(X_2\right)|X_1\right] |
| </jsmath> |
| |
| の方だが、手続きとしてはまず、上の<jsm>X_1=x_1</jsm>を与えた条件付き期待値を計算してから、改めて、<jsm>x_1</jsm>を確率変数<jsm>X_1</jsm>で置き換えることになる。 |
| これはすなわち、<jsm>\mu_\phi\left(x_1\right)</jsm>の<jsm>x_1</jsm>を確率変数とみなせ、という意味で、<jsm>\mu_\phi\left(X_1\right)</jsm> を考えよ、ということだから、これは確率変数 <jsm>X_1</jsm> の関数なので、確率変数。 |
| |
| ==== #15 2011.07.21 ==== |
| |
| * 中心極限定理 |
| * レポート一斉返却 |
| |
| |
| 参考書: |
| * 清水良一(1976)「中心極限定理」, 教育出版. |
| * 竹内啓(1975)「確率分布の近似」, 教育出版. |
| * 竹内啓(1974)「統計的推定の漸近理論」, 教育出版. |
| * D. Williams(1991, 赤堀・原・山田・訳, 2004)「マルチンゲールによる確率論」, 培風館. |
| |
| |
| ==== #16 2011.07.28 ==== |
| |
| 16回目なので休講。 |
| 期末試験が16週目になる。 |
| |
| ==== #Exam 2011.08.04 ==== |
| |
| 期末試験 |
| |
| |日時|2011.08.04 0240pm-0410pm| |
| |場所|C-301| |
| |
| ルール |
| |
| * 通信機能を持たない電卓の持ち込みは可とする |
| * 出席をとるので学生証を持参のこと |
| * 退室の願い出は、試験開始の30分後から許可する |
| |
| お願いごと |
| |
| * 回答用紙は、可能な限り1ページ単位で使用してほしい |
| |