rand()^ 2の分布がrand()* rand()と異なるのはなぜですか?


14

Libre Office Calcでは、rand()均一分布から0〜1のランダムな値を選択する関数が使用可能です。私は自分の確率で少しさびているので、次のような振る舞いを見たとき、私は困惑しました:

A = 200x1列 rand()^2

B = 200x1列 rand()*rand()

mean(A) = 1/3

mean(B) = 1/4

なぜmean(A)!=なの1/4か?


3
確率変数の平方の期待値は、その期待値の平方に等しくないためです。
マイケルM

3
rand()他の同様の演算子のように動作する場合、Aは同じ乱数の2乗であり、Bは乗算された2つの乱数です。
ピーターフロム-モニカの復職

わかります。数学が正確に記述されているか、これを行うリソースにリンクされているのを見ることができれば、非常に役立ちます。
ジェフトピア

4
状況を単純化すると、ポイントを理解するのに役立ちます。でRand()置き換えられたと仮定しInt(2*Rand())ます。これは、等しい確率で値とを取ります。その平方には2つの可能性があり、2つの(独立した)値の積には4つの可能性があります。期待値を計算するとどうなりますか?101
whuber

回答:


24

長方形を考えると役立つ場合があります。あなたが無料で土地を手に入れるチャンスがあると想像してください。土地のサイズは、(a)確率変数の1つの実現、または(b)同じ確率変数の2つの実現によって決定されます。最初のケース(a)では、面積は正方形であり、辺の長さはサンプリングされた値に等しくなります。2番目のケース(b)では、2つのサンプリング値は長方形の幅と長さを表します。どの選択肢を選びますか?

してみましょう正の確率変数の実現可能。U

a)1つの実現の期待値は、U 2に等しい正方形の面積を決定します。平均して、エリアのサイズは E [ U 2 ]になりますUU2

E[U2]

B)は、2つの独立した実現ある場合及びU 2は、エリアとなりU 1U 2。平均して、大きさが等しい E [ U 1U 2 ] = E 2 [ U ] 両方の実現が同じ分布と独立であるからです。U1U2U1U2

E[U1U2]=E2[U]

領域a)とb)のサイズの差を計算すると、次のようになります。

E[U2]E2[U]

上記の用語は本質的に0以上のと同じです。Var[うん]0

これは一般的な場合に当てはまります。

あなたの例では、一様分布からサンプリングされた。したがって、うん01

E2[U]=1

E[U]=12
Var[U]=1
E2[U]=14
Var[U]=112

我々が得 E [ U 2 ] = 1E[うん2]=Var[うん]+E2[うん]

E[うん2]=112+14=13

これらの値は分析的に得られたものですが、乱数ジェネレーターで取得した値と一致します。


任意の場合はB、私が手2 + B + bは2aba2+ab+b23
ジョン

それは分散の賢い使い方です。そして、ここで私は数学を直接打ち破ろうとしていました。
アフィン

これは理にかなっています。それはすべて、分散が非負であることにかかっています。また、ジョンがどのように答えを得たかについても興味があります。
ジェフトピア

基本的には、Svenが行ったことをそのまま踏襲しましたが、より一般的な均一分布の式に置き換えました。
ジョン

E [ U 2 ] E 2 [ U ]を読むべきではありませんか?E[U2]E[U2]E[U2]E2[U]
BoppreH

9

Svenの優れた答えに欠けているものがあることを示唆するものではありませんが、私はこの質問について比較的初歩的な見解を提示したかったのです。

ジョイント分布が非常に異なることを確認するために、各製品の2つのコンポーネントをプロットすることを検討してください。

u1対u2およびu1対u1のプロット

製品は両方のコンポーネントが大きい場合にのみ大きくなる傾向があることに注意してください(1に近い)。これは、2つのコンポーネントが独立ではなく完全に相関している場合により簡単に起こります。

1ϵϵϵ/2U2U1×U2ϵ2/2

かなり違います!

上記のようなグラフに等積等高線を描くと便利です。つまり、0.5、0.6、0.7、0.8、0.9のような値に対してxy = constantの曲線です。値が大きくなるにつれて、独立したケースの場合、輪郭の上部と右側のポイントの割合がはるかに速く低下します。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.