注文統計の変換


9

確率変数およびは独立しており、分散されていると仮定します。ことを示し有する\ text {Exp}(1)分布。X1,...,XnY1,...,YnU(0,a)Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))Exp(1)

私は設定して、この問題を始めました{X1,...,Xn,Y1,...Yn}={Z1,...,Zn}次にmax(Yn,Xn)=Z(2n)(\ frac {z} {a})^ {2n}として配布され(za)2nmin(Yn,Xn)=Z(1)1(1za)2n 密度は、f_ {Z_ {1}}(z)=(2n)(1- \ frac {z} {a})^ {2n-1} \ fracとして簡単に見つけることができます{1} {a}fZ1(z)=(2n)(1za)2n11aおよびfZ(2n)(z)=(2n)(za)2n11a

これは、これらが計算されたため、次にどこに行くべきかを知るのに苦労しているところです。変革で何かをしなければならないと思っていますが、よくわかりません...


確かに、XiY_i iidYiだけでなく、X_iY_jからXi独立していると想定する必要があります。それで\ log(Z_i)を直接操作することを考えましたか?Yjlog(Zi)
whuber

@whuberコメントからの私の考えは、n * log(Z)の密度を解決する変換を設定することでしょうか?i
スーザン

私は少し再フォーマットしました(特にとをと変換しました)が気に入らない場合は、以前のバージョンにロールバックできます( "edited <x> ago"リンクをクリックして)投稿の下部にある私の墓地の上にあります)、次に以前のバージョンの上にある[ロールバック]リンクをクリックします。logminlogmin
Glen_b-2015

3
スーザン、あなたは質問を誤解/誤解しているようです。質問は、の比率を求めます。分母は、:はの最大次数統計、は最大次数統計s。つまり、 は、すべてのとの最小値ではなく、min(maxX、maxY)をシークするため、Zトリックを使用できません。すべてのX値とY値を平坦化/結合します。.......
max(Y(n),X(n))min(Y(n),X(n))
min(Y(n),X(n))Y(n)YX(n)Xmin(Y(n),X(n))XY
wolfies

2
いずれにしても、別の問題として、の密度を計算することは(これまでに行ったように)意味がなく、の密度を個別に計算することはできません。一般的に独立しています。の比率を見つけるには、最初にの結合pdfを見つける必要がありますそれが問題だった場合手元にあります(そうではありません)。Z(1)Z(2n)Z(2n)/Z(1)(Z(1),Z(2n))
wolfies

回答:


2

この問題は定義だけで解決できます。唯一の高度な計算は単項式の積分です。


予備観察

変数と全体的に操作してみましょう。これは変更しませんが、 iidを均一分布にして、計算における気が散るすべての外観を排除します。したがって、一般性を失うことなくと仮定できます。Xi/aYi/aZn(X1,,Yn)(0,1)aa=1

とその一様分布の独立性は、である任意の数値に対して、Yiy0y1

Pr(yY(n))=Pr(yY1,,yYn)=Pr(yY1)Pr(yYn)=yn,

について同じ結果が保持されます。将来の参考のために、これにより、X(n)

E(2X(n)n)=012xnd(xn)=012nx2n1dx=1.

解決

してみましょう正の実数とします。の分布を見つけるには、その定義を置き換えて、結果の不等式を単純化します。tZn

Pr(Zn>t)=Pr(Zn/n>t/n)=Pr(exp(Zn/n)>et/n)=Pr(max(X(n),Y(n))min(X(n),Y(n))>et/n)=Pr(et/nmax(X(n),Y(n))>min(X(n),Y(n))).

このイベントは、またはのどちらが小さいかによって、確率が2つのケースに分かれます(確率がゼロの交差は無視できます)。したがって、これらのケースの1つ(たとえば、の方が小さいの可能性を計算し、それを2倍にするだけで済みます。以来、私たち(せる時にできるように、役割を果たしてof)予備セクションの計算を適用するには:X(n)Y(n)Y(n)t00et/nX(n)1et/nX(n)y

Pr(Zn>t)=2Pr(et/nX(n)>Y(n))=2E[(et/nX(n))n]=etE[2X(n)n]=et.

これが、がExp分布を持つことの意味です。Zn(1)


4

解決策をスケッチします。ここでは、コンピュータ代数システムを使用して、重要なことを行います...

解決

場合サイズのサンプルである親上の 、次いで試料最大のPDFがある: およびについても同様です。X1,...,XnnXUniform(0,a)

fn(x)=nanxn1
Y

1アプローチ:の関節PDFを探す(X(n),Y(n))

とは独立しているため、2つのサンプルの最大値の結合pdfは、2つのpdfの積です。たとえば、:XY(X(n),Y(n))f(n)(x,y)

ここに画像の説明を入力してください

与えられ 。次に、の累積分布関数はです。Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))ZnP(Zn<z)

ここに画像の説明を入力してください

自動化するためにMathematicaProbmathStaticaパッケージの関数を使用しています。cdf wrt微分すると、のpdfが標準の指数として得られます。zZn


アプローチ2:注文統計

次数統計を使用して、Max関数とMin関数を処理する必要があるメカニズムを「バイパス」できます。

もう一度:が親 サイズサンプルである場合である場合、サンプルの最大値のは、言う、: X1,...,XnnXUniform(0,a)W=X(n)fn(w)

ここに画像の説明を入力してください

サンプルの最大値およびは、この分布からの2つの独立した図面にすぎません。つまり、のおよびオーダーの統計(サイズ2のサンプル)は、まさに私たちが探しているものです。X(n)Y(n)W1st2ndW

  • W(1)=min(Y(n),X(n))

  • W(2)=max(Y(n),X(n))

の結合pdfは、サイズ2のサンプル、たとえば、次のようになります。(W(1),W(2))g(.,.)

ここに画像の説明を入力してください

与えられ 。次に、の累積分布関数はです。Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))ZnP(Zn<z)

ここに画像の説明を入力してください

このアプローチの利点は、確率計算にmax / min関数が含まれなくなるため、(特に手動で)導出を多少簡単に表現できるようになることです。

その他の

上記の私のコメントによると、あなたは質問を誤解しているようです...

私たちは見つけることを求められます:

Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))

ここで、分母はmin(xMax、yMax)、...すべてのとの最小値ではありません。XY


あなたのスケッチに続いて、私が質問を誤って解釈したことを理解しました。2つのサンプルの最大値の結合pdfを計算する方法は理解していますが、max / minの比率をどのように解釈するかはまだわかりません。
スーザン

最大/最小を「回避」する注文統計を使用して、別の派生を追加しました。
wolfies

データのログであるスーザンから始めた場合は、比率ではなく注文統計の違いを確認することになります。
whuber

比率がExp(1)確率変数である理由を説明するには、コンピューター形式の計算を使用するのが最善の方法だとは思いません。
西安

1
良い点... OPが理由を尋ねないことを除いて...それがExp [1]であることを示すため。これが宿題(または課題)であるかどうかもわかりません...これは、実際にはコンピュータを使用する優れた利点の1つです。手順を実行し、結果を確認して、正しいアプローチをとることができます。 、しかしメカニックはまだOPに任されています。誰かが最初にログを取るという@whuberの提案を探索するのは良いことです。
wolfies
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.