中心極限定理のシミュレーション問題: 解法とプログラムの作り方
- 乱数表から数を抽出し、それに応じた変換を行う問題です。
- プログラミング知識がなくてもFORTRAN77を使ったプログラムで解けます。
- 50回の繰り返しを行い、変換した数25個の平均値をグラフで表示します。
- ベストアンサー
この問題、どのように解けばよいですか?
中心極限定理のシュミレーションとして出された問題なのですが、乱数表(00から99まで)から無作為に数を抽出し、その数が00から49なら1、50から79なら2.80から99なら3に変換して、数の抽出を25回行い、変換した数25個の平均値を出すという作業を50回繰り返したものをグラフにするという問題なのですが、どのように解けばよいでしょうか?プログラムを使って作るようにいわれたのですが、プログラムに関してはほとんど知識はないので、(FORTRAN77ならかろうじて学校でやっているので少しはわかるのですが、もしできればFORTRAN77を使ったプログラムでお願いします)できるだけ丁寧に教えてください
- furenntirappu
- お礼率55% (5/9)
- 数学・算数
- 回答数2
- ありがとう数1
- みんなの回答 (2)
- 専門家の回答
質問者が選んだベストアンサー
プログラムの概要は ------------------------------------------------------------ for i=1,50 #25回の抽出作業 #平均値の計算 #平均値を出力 next i ------------------------------------------------------------ のようになります.#部分は数行ずつ必要になります. まず,「#25回の抽出作業」の部分は, ------------------------------------------------------------ sum=0 for j=1,25 #数値の00から99の数値を無作為に抽出→xとする. #xを変換しyとする sum=sum+y next j ------------------------------------------------------------ となります.「#平均値の計算」と「#平均値の出力」は, ------------------------------------------------------------ ave=sum/25 write "F16.8",ave ------------------------------------------------------------ のようになります.write文の書式は記憶違いしているかもしれないので, furenntirappuさんの知識をもとにして,必要なら修正してください. さて, #数値の00から99の数値を無作為に抽出→xとする. の部分ですが, 0から1の乱数を取り出す関数ってありませんか?忘れてしまいましたが. その値を百倍し,整数部分だけとりだせば,00から99のうちの一つになります. あとは, if文で0<=x かつx<=49ならy=1 if文で50<=x かつx<=79ならy=2 if文で80<=x かつx<=99ならy=3 をそのまま書いてyに変換すれはよいです. 以上まとめると, ------------------------------------------------------------ for i=1,50 sum=0 for j=1,25 0から1の乱数を取り出し百倍して整数化する. if文で0<=x かつx<=49ならy=1 if文で50<=x かつx<=79ならy=2 if文で80<=x かつx<=99ならy=3 sum=sum+y next j ave=sum/25 write "F16.8",ave next i ------------------------------------------------------------ fortranの文法をよく覚えていないので,適当に修正してください. グラフ化もfortranでやるのなら,「write "F16.8",ave」 の部分に少し工夫が必要です.
その他の回答 (1)
- wolv
- ベストアンサー率37% (376/1001)
グラフの表示は,例えば 0 1 2 3 1: * 2: * 3: * : 25: * というように,折れ線グラフの線のないものを縦にしたようなグラフを書くのが 一番簡単です. ------------------------------------------------------------ write "F16.8",ave ------------------------------------------------------------ のかわりに, ------------------------------------------------------------ 文字列のクリア(string=""など) for k=1,int(ave*10)-1 文字列=文字列+” ” next k 文字列=文字列+”*” write 文字列 ------------------------------------------------------------ のようなことをやればグラフになります. for文を使って,ave(の10倍)の長さに相当する文字列を作り, それを表示してます. 特定の長さの文字列を生成する関数があったかもしれませんが, 自信がないので,あえて自分で作ってみました. ------------------------------------------------------------ この他に, 配列などを利用して度数分布を計算し,50回の作業が終わったところで 一気にグラフとして出力することも考えられます.
関連するQ&A
- PERT と中心極限定理 2
中心極限定理は、次のような定理だと思います。 平均μ、分散σ2の母集団から無作為にn個の標本を抽出してその平均値mを求めることを繰り返すと、母集団がどのような分布を示す集団であるかに拘わらず、nが充分大きいとき、mの分布は平均μ、分散σ2/nの正規分布で近似される。 次に、PERTにおいて、n個の作業から成るプロジェクトの全体工程Tを求める方法は、一般に次のように説明されています。 作業iの所用時間がベータ分布に従うと仮定すると、その期待値ei、楽観値oi、最可能値mi、悲観値pi、分散σi2の間には次の関係がある。 ei=(oi+4mi+pi)÷6 (式-1) σi2=(pi-oi)2÷36 (式-2) 一般に平均と分散については加法定理が成り立つので、クリティカルパス上のn個の作業の総所要時間(n個の作業の所要時間の合計)Tの期待値eと分散σ2は次のように表される。 e=Σei (式-3) σ2=Σσi2 (式-4) 中心極限定理により、Tは期待値e、分散σ2の正規分布で近似されるので、今、e=20、σ2=25であるとすると、95%の確率でTが完了する工期は、標準正規分布表の95%点=1.960から、20-1.960×5≦T≦20+1.960×5となる。 それで、次の(1)、(2)が分かりません。 (1)中心極限定理は、「平均μ、分散σ2の母集団から無作為にn個の標本を抽出してその平均値mを求めること」から始まる定理なのに、上記Tを求めた過程には、「平均μ、分散σ2の母集団」も「n個の標本の抽出」も「その平均値m」も、一切何もありません。母集団、抽出、平均値にあたるものは、上記Tを求めた過程のどの値または計算なのでしょうか。 (2)「中心極限定理により、Tは期待値e、分散σ2の正規分布」とありますが、いったいどう考えればTを正規分布であるとみなせるのでしょうか。
- ベストアンサー
- 数学・算数
- PERT と中心極限定理
中心極限定理は、次のような定理だと思います。 平均μ、分散σ2の母集団から無作為にn個の標本を抽出してその平均値mを求めることを繰り返すと、母集団がどのような分布を示す集団であるかに拘わらず、nが充分大きいとき、mの分布は平均μ、分散σ2/nの正規分布で近似される。 次に、PERTにおいて、n個の作業から成るプロジェクトの全体工程Tを求める方法は、一般に次のように説明されています。 作業iの所用時間の期待値ei、楽観値oi、最可能値mi、悲観値pi、分散σi2の間には次の関係がある。 ei=(oi+4mi+pi)÷6 (式-1) σi2=(pi-oi)2÷36 (式-2) 一般に平均と分散については加法定理が成り立つので、クリティカルパス上のn個の作業の総所要時間(n個の作業の所要時間の合計)Tの期待値eと分散σ2は次のように表される。 e=Σei (式-3) σ2=Σσi2 (式-4) 中心極限定理により、Tは期待値e、分散σ2の正規分布で近似されるので、今、e=20、σ2=25であるとすると、・・・(と来て、Tの確率を求めるのですが、長くなるので以下省略します)。 ここで分からないのは、「中心極限定理により、Tは期待値e、分散σ2の正規分布で近似される」というところです。なぜ、いきなりこんなことが言えるのでしょうか。具体的に分からない点は次の(1)です。 (1)中心極限定理は、「平均μ、分散σ2の母集団から無作為にn個の標本を抽出してその平均値mを求めること」から始まる定理なのに、上記Tを求めた過程には、「平均μ、分散σ2の母集団」も「n個の標本の抽出」も「その平均値m」も、一切何もありません。一体これらは、上記Tにおいては、どこへ行ってしまったのでしょうか。
- ベストアンサー
- 数学・算数
- 中心極限定理について
中心極限定理についていろいろ調べたんですが、よくわからないことがあります。 (ほとんどの)任意の母集団(平均μ、分散σ^2)からn個の確率変数x1, x2, x3 .... xn を無作為抽出すると、平均値X を求めると、その平均値の分布は、nが大きくなると正規分布(平均μ、分散σ^2 / nの平方根)に近づく と書いています。 ある母集団分布をおいて、n = 6 として、 サンプル1: x1, x2, x3, x4, x5, x6 を抽出し、平均値 X1 を求める サンプル2: また、x1, x2, x3, x5, x6 を抽出し、平均値 X2 を求める サンプル3: また、x1, x2, x3, x5, x6 を抽出し、平均値 X3 を求める と同様に、やっていくのですよね? で、この「1回に抽出するデータがn」(上記では n = 6)であり、この nが大きくなると正規分布に近づくということなんですが、 サンプル数(平均値Xの数)はいくつを想定しているのでしょうか? サンプルが無限だったら、n = 2 だろうが、n = 100 だろうが正規分布のような気がするのですが。 いろいろなサイトで、 n = 1 や、n = 2 のケースでやったときの分布図が掲載されているのですが、これはサンプル自体の数が多いのでしょうか? n の意味と、サンプルの数(平均値Xの数) が混乱しているようです。 教えてください。
- ベストアンサー
- 数学・算数
- Rにおいて、中心極限定理を確認するために、乱数の標本平均のヒストグラム
Rにおいて、中心極限定理を確認するために、乱数の標本平均のヒストグラムと正規分布のグラフを重ね合わせたいのですが、軸ラベルがおかしくなってしまいます。例えば、1000個の標本平均のヒストグラム hist(x)と平均0,分散1の正規分布のグラフを上手く重ね合わせるにはどうしたらいいですか?
- ベストアンサー
- 数学・算数
- 有限母集団の中心極限定理
訳あって確率統計の勉強をしておりますが、中心極限定理について質問があります。 「母集団の平均をm、分散をvとすると、そこから抽出したn個の標本の平均の分布は、平均m、分散v/nという分布になり、標本数をn→無限大とすると、分布は母集団の分布によらず正規分布に近づく。」 とありますが、母集団が有限個(N個)の集合ならどうなるでしょうか。 その場合、標本数をnがNに等しくなった時点で平均はm、分散0、つまり標本から母集団の平均が完全に推定(決定)することになります。 ●有限母集団の場合も中心極限定理は成り立つのか? ●成り立つならn→Nで分散が0になるという点はどう表現されるのか? このあたりを教えてください。
- 締切済み
- 数学・算数
- 中心極限定理 一様分布 平均値の意味
統計学からの疑問です。 以前も中心極限定理について質問したことがあったのですが、また疑問が発生しました。 中心極限定理は、母集団が何分布であれ、そこからサンプリングされた標本の平均値は正規分布に従って分布するということで、この定理があるからこそ色んな統計手法を使うことができますよね。 例えば、プレス機の設定荷重を5.0トンに設定し、実荷重を複数回測定するとします。 この場合、実荷重は、4.9トンとか5.1トンとか、誤差的にばらつきますよね。 なので複数回測定した実荷重の平均値=プレス機の実力であると言えると思います。 では母集団が一様分布、例えばサイコロの場合、 サイコロを2個以上振って出た目の平均値は3.5ですが、この値は一体何を意味してるのでしょうか?サイコロは3.5の目がもっともよく出るというわけはありませんし・・・ 中心極限定理を解説するために、サイコロを使った説明があり、 母集団が一様分布であっても中心極限定理がちゃんと現れることは確認できたのですが、 母集団が一様分布の場合の統計量(平均値)の意味がわかりません。 変な質問かもしれませんが、疑問を解決したいのでどなたか教えて下さい。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- 分からなくて困ってます
投げると3/5の確立で表が出るコインがあるとします。 (1)このコインを600回投げるとき、表が380回以上出る確率を中心極限定理を用いて 近似的に求めよ。 (2)このコインを600回投げるとき、表が出る回数350回以下となる確率を、中心極限定理を用いて近似的に求めよ。 組み合わせを使ってやるやりかたしか分かりません。 中心極限定理を用いるのが分からないんで教えてください。 お願いします
- ベストアンサー
- 数学・算数
- 確率統計の問題
確率密度関数f(X)= 1(0<=X<=1) 0(その他) で与えられている。 Z=Σ(i=1,12)Xi-6 として新しい確率変数Zを定義するとき、Z~N(0,1)と近似できることを中心極限定理によって説明せよ。 という問題です。 私のやり方としては、まずf(X)の平均と分散を求める。 μ=1/2 σ^2=1/12 と求まります。 ここで、Zを変形すると、 Z=(1/√12)*(Σ(i=1,12)(Xi-1/2)/(√1/12)) となり、ここの12をnと置き換えればn→∞のときに Z~N(0,1)と近似できるとわかりますが、 中心極限定理ではn が大きいときに、 Z=(1/√n)*(Σ(i=1,n)(Xi-μ)/(σ)) の分布は標準正規分布でよく近似できるとありますが、 問題文ではn=12であって、これが何を持って大きいといえるのです。 確率統計の初心者ですので、馬鹿な質問になり、申し訳ありません。 わかる方がいらっしゃいましたら、ご教授お願いします!
- ベストアンサー
- 数学・算数
お礼
質問してから早い時間にありがとうございました。おかげでできました、さすがにグラフまではできませんでしたが。ちなみに少し修正させていただくとwrite(*,*)でforはdoになります。