モーメント推定量

このQ&Aのポイント
  • モーメント推定量を求めるためには、X, Y, Zの分布からの無作為標本を用いて、a, b, c, dの値を推定する必要があります。
  • また、UとXの同時分布が正規分布である場合、UとXは独立であることが示されます。
回答を見る
  • ベストアンサー

モーメント推定量

統計の問題がわかりません。 問題 標準正規分布に従う確率変数X,Y,Z,U,Vについて Y=aX+bU X=cZ+dV が成り立ち、E[XU]=E[ZV]=E[ZU]=0 とする。a,b,c,dは正の定数とする。 1)(X、Y、Z)の定める分布からのn個の無作為標本(x1、y1、z1)、…、(xn,yn,zn)が得られたとき、a,b,c,dのモーメント推定量を求めよ。 2)(U,X)の同時分布も正規分布であるとき、UとXが独立であることを示せ。 1)はr次モーメントとr次標本モーメントをもとめて解くのかと思いつつ立式できません。。。 2)も何を示せばいいのかわかりません。。。 考え方を教えてください。

質問者が選んだベストアンサー

  • ベストアンサー
noname#227064
noname#227064
回答No.1

> 1)(X、Y、Z)の定める分布からのn個の無作為標本(x1、y1、z1)、…、(xn,yn,zn)が得られたとき、a,b,c,dのモーメント推定量を求めよ。 モーメント推定量が良くわかっていないので間違っているかもしれませんが… E[XU] = E[X(Y - aX)/b] = E[XY]/b - aE[X^2]/b = 0 から a = E[XY]/E[X^2] となるので、(Σxiyi)/(Σxi^2)がaの推定量となる。 また、Y=aX+bUでX,Y,Uが標準正規分布に従うことから、 E[Y^2] = E[(aX + bU)^2] = a^2 E[X^2] + 2ab E[XU] + b^2 E[U^2] = a^2 + b^2 = 1 となり、b > 0であることから√{1 - (Σxiyi)^2/(Σxi^2)^2}がbの推定量となる。 c, dも同様。 > 2)(U,X)の同時分布も正規分布であるとき、UとXが独立であることを示せ。 「(U,X)の同時分布も正規分布であるとき」というのは、二変量正規分布であるということでしょうか? そうであるなら、XとUの共分散が Cov(X, U) = E[XU] - E[X]E[U] = 0 となり、XとUの同時密度関数をXとUの周辺確率関数の積で表せることから独立であることを示せます。

Xepher13
質問者

お礼

わかりました! ありがとうございます!

関連するQ&A

  • 慣性モーメントの問題

    密度が一定(ρ_0とする)の次の指定された軸の周りの慣性モーメントを計算せよ。 (1)一辺の長さaの立方体。立方体の一辺を軸にした場合の慣性モーメントおよび、立方体の中心を通り、相対する2面に垂直な軸に関する慣性モーメント。 (2)xyz空間でx^2+y^2≦z^2,0≦z≦hで定義される物体。z軸に関する慣性モーメント。 この問題で、(1)の前者=∫∫∫_D ρ_0(x^2+y^2)dxdydz,D=(0≦x≦a,0≦y≦a,0≦z≦a)として解いてあり、(1)の後者=8∫∫∫_D ρ_0(x^2+y^2)dxdydz,D=(0≦x≦a/2,0≦y≦a/2,0≦z≦a/2)として解いてあり、(2)=∫∫∫_V ρ_0(x^2+y^2)dxdydz,v=(x^2+y^2≦z^2,0≦z≦h)として解いてあったのですが。 ここで質問です。 (1)確かに慣性モーメントの公式はI=∫∫∫ρ(x、y、z)r^2(x,y,z)dxdydzなのですが、なぜ、(1)(2)ともにr^2のところがx^2+y^2なのでしょうか?zが入ってもおかしくないと思うのですが。r^2は軸からの距離ですし。これは、自分の都合よく平面で考えて良いということなのですか? (2)(1)の後者で8倍してありますが、何をもって8倍なのでしょうか? (3)そもそも慣性モーメントとはなんなのか? この3つの私の質問の答えを教えてください。よろしくお願いします。

  • 二次形式のモーメント母関数について

    二次形式のモーメント母関数について 閲覧ありがとうございます。 確率統計について分からない問題があるので教えてください。 問 X[1],X[2],...,X[n]は独立で共通な標準正規分布N(0,1)に従う確率変数である。 X = (X[1],X[2],...,X[n])'とおく。 (’はベクトル・行列の転置を表す。) A=(a[i][j])をn次の実対称行列とし、固有値全体を{λ[1],λ[2],...,λ[n]} (λ[1]≧λ[2]≧...≧λ[n]) とする。 Aを係数行列にもつXに関する二次形式 Z = X' A X =Σ[n,i,j=1] a[i,j]X[i]X[j] このとき以下について答えよ (1)Mz(t)のモーメント関数を求めよ。 (2)Zの確率分布が自由度k (1≦k≦n)のカイ二乗分布ならば λ[1]=...=λ[n]=1,  λ[k+1]=...=λ[n]=0であることを証明せよ。 (1)については全然わかりません。。。すいません。 (2)について n次の直交行列をU,固有値を対角成分とする行列をLとしてA=U' L Uと分解します。 Y=UXとして、 Z=X' A X =Y' L Yとかける。 条件よりYは確率分布N(0,1)に従う確率分布である。 ここまで考えましたが後がわかりません。 よろしくお願いします。

  • 合成関数の偏微分法を用いた解き方

    いつもお世話になっています。 以下の問題を解いてみたのですが、あっているのか自信がもてません。 (特に、(4)(5)のsinθ,cosθが含まれるケース) 間違いなど、あればご指導のほど、よろしくお願いいたします。 【問題】 「合成関数の偏微分法」を用いて、継ぐの合成関数についてZu,Zv(またはZθ,Zr)を求めよ。 (2) z=x^2-y, x=u+v, y=uv Zu = Zx・Xu + Zy・Yu = 2x・1+(-1)・v=2x-v Zv = Zx・Xv + Zy・Yv = 2x・1+(-1)・u=2x-u (3) z=e^x・sin(y), x=u-v, y=uv Zu = Zx・Xu + Zy・Yu = e^x・sin(y)・1+e^x・cos(y)・v = e^x・sin(y)+v・e^x・cos(y) Zv = Zx・Xv + Zy・Yv = x^x・sin(y)・(-1)+e^x・cos(y)・u = -e^x・sin(y)+u・e^x・cos(y) (4) z=x+y, x=r・cosθ, y=r・sinθ Zθ= Zx・Xθ + Zy・Yθ = (1)・(-r・sinθ)+(1)・(r・cosθ) = (-r・sinθ)+(r・cosθ) = -r(sinθ-cosθ) Zr = Zx・Xr + Zy・Yr = (1)・(cosθ)+(1)・(sinθ) = sinθ+cosθ (5) z=x^2+2xy, x=r・cosθ, y=r・sinθ Zθ= Zx・Xθ + Zy・Yθ = (2x+2y)・(-r・sinθ)+(2x)・(r・cosθ) = 2{(x+y)(-r・sinθ)+x(r・cosθ)} = -2r{(x+y)(sinθ)-x(cosθ)} = -2r(x・sinθ+y・sinθ-x・cosθ) = -2r(r・cosθ・sinθ+r・sinθ・sinθ-r・cosθ・cosθ) = -2r^2(cosθ・sinθ+sinθ・sinθ-cosθ・cosθ) = -2r^2(sin^2θ+cosθsinθ-cos^2θ) Zr = Zx・Xr + Zy・Yr = (2x+2y)・(cosθ)+(2x)・(sinθ) = 2{(x+y)・(cosθ)+(x)・(sinθ)} = 2{x・cosθ+y・cosθ+x・sinθ} = 2{r・cosθ・cosθ+r・sinθ・cosθ+r・cosθ・sinθ} = 2r{cosθ・cosθ+sinθ・cosθ+cosθ・sinθ} = 2r{cos^2θ+2・sinθ・cosθ} = 2r・cosθ{cosθ+2・sinθ} 以上、よろしくお願いします。

  • 3次元混合正規分布のパラメータ推定

    3次元混合正規分布のパラメータ推定 標本(X,Y,Z,値)が大量にあり、これを3次元の混合正規分布で近似しようとしています。 現在、平均や分散などのパラメータの算出をEMアルゴリズムを用いて行おうと考えています。 いろいろ調べて勉強してみましたが、標本をどのように使えばいいか理解することができません。 この標本はどのような形で式につかえばよいでしょうか? そもそもEMアルゴリズムを用いてこの問題を解くことができるのでしょうか? また、この標本をもとにパラメータを決定する代替方法があれば、教えていただけると助かります。 説明不足な点があるかもしれませんが、よろしくお願いします。

  • 確率論の問題について

    (1)「確率変数Xが()(0,a)上の一様分布U(0,a)に従うとき、また()正規分布N(m,v)に従うとき、その標本化Zの分布密度関数を求めよ」 (2)「Xを標準正規分布N(0,1)に従う確率変数であるとする。Y=|X|の密度関数を求めよ。Yの平均と分散を求めよ」 というものなのですが(1)(2)ともにまったく手をつけることができません(泣)アドバイスなどお願いします(泣)

  • 円の断面二次モーメントの算出方法について

    円の断面二次モーメントの算出方法について 私なりに計算した方法を以下に示します。 x^2+y^2=r^2の円に原点よりy離れた場所に単位面積ΔA((2×√r^2-y^2)×dy)に応力σが作用している時 のモーメントM(x)=2×y×(√r^2-y^2)×σとなり、y=r×sinθ、σ=y・E/Rとすると(E:ヤング率、R:半径曲線)、 M(x)=2×E/R×r^4×sin^2θ×cos^2θとなります。 円全体のモーメントM=∫(-1/2π~1/2π)(2×E/R×r^4×sin^2θ×cos^2θ)dθ (sin^2θ×cos^2θ=1/4sin2θ) これは不定積分の公式より(sin^2θ=(x/2)-(1/4)sin 2x) 与式=π×E/R×r^4=π/32d×r^4×E/R (r=1/2d) となってしまいました。 公式ではπ/64d×r^4となります。 大変恐縮ですが、計算、考え方等間違いがありましたが、ご教授下さい。 どうかよろしくお願いします。

  • 大学の統計学です 確率母関数、ベルヌーイ分布、モーメント母関数

    明日試験なのですが、勉強不足で全然わかりません・・・・ ・2項分布B(n,p)の確率母関数を計算せよ ・幾何分布Ge(p)の確率母関数を計算せよ ・X1,X2....Xnを互いに独立でベルヌーイ分布に従うn個の確率変数とするとき、Sn=X1+X2+...+Xnの分布が2項分布となることを示せ またSn/nの平均値と分散を求めよ ・指数分布Exp(θ)のモーメント母関数、平均値(期待値)、分散を計算せよ ・2回のサイコロ投げにおいて、Xを最初の目、Yを2回目の目とするとき、Z=X+Y,W=X-Yとおく (1)ZとWの平均値を求めよ(2)ZとWの分散をもとめよ(c)ZとWの共分散を 求めよ ・X1,X2....Xnを互いに独立で同一の分布に従う確率変数とする。 E(Xi)=μ、V(Xi)=σ^2、i=1,....,nとしX1,X2....Xnの標本平均をZ=1/n(X1,X2....Xn)とおく。 E(Z)とV(Z)を計算せよ わかる方教えていただけたら嬉しいです!!!! よろしくお願いします。

  • 【統計学基礎】自由度とは?

    【統計学基礎】自由度とは? 以下の問題で自由度について答えさせられるのですが自由度とはなんでしょうか?自分はただ自由度=標本の大きさ-パラメータと考えていたのですが色々パラメータがどれを指すかによってかわってきて混乱してきています。 平均u=5、分散σ^2=4 の正規母集団より大きさn=16の無作為標本を抽出する。このときy=(n-1)s^2/σ^2は以下省略…また‐xを標準化したZは正規標準分布にしたがい、‐Xとs^2は統計的に有意であることから、z/√y/15は自由度???のt分布に従う。 という問題です。まずそもそもz/√y/15の式の意味もよくわからないのですがご回答お願いします。

  • 慣性モーメントがわかりません

    長半径a(x軸上)短半径b(y軸上)の楕円体の慣性モーメントを求めたいです(原点が楕円の中心でx軸y軸z軸それぞれの慣性モーメントを求めたい) ∫0からa ρ4x二乗b√1-(x二乗/a二乗)dx でy軸に関する慣性モーメントが求まるときいたんですがなぜなのかわかりません 慣性モーメントは∫r^2dmでもとまるので、rが上のxに対応してると思うのですがなぜそのままx^2とできるのかわからず、またb√1-(x二乗/a二乗)つまり楕円の式をy=にした形にρをかけたものがなぜdmになるのが想像つきません 急ぎではないので暇なときに答えてくれたらうれしいです

  • 偏微分方程式

    偏微分方程式の問題についてです。 (1)u_xx+4u_xx+4u_yy=0を正規形に変換して解け。 v=x,z=2x-yと置いて解くのですが、この後どうやって解くのでしょうか? (2)xu_xy=yu_yy+u_yを解け。 を(1)と同じでv=x,z=xyと置いて解くのですが、この後どうやってとくのでしょうか?