• 検索結果がありません。

補足資料 数理統計 2016 S1・S2 Kengo Kato

N/A
N/A
Protected

Academic year: 2018

シェア "補足資料 数理統計 2016 S1・S2 Kengo Kato"

Copied!
2
0
0

読み込み中.... (全文を見る)

全文

(1)

2016.6.27. 加藤賢悟 6/27の講義中に補足した証明ですが,いろいろタイポがあったので,訂正版を載せてお きます.

ΘはRの開区間とし,X = (X1, . . . , Xn)の同時確率(密度)関数は pn(x; θ) = H(x) exp{θT (x) − C(θ)}

の形に表せるとしよう.

Theorem 1. 次の形の検定はH0 : θ = θ0 vs. H1: θ ̸= θ0に対する水準αUMPU検定に なる:

δ(x) =





1 if T (x) < c1 or T (x) > c2

γi if T (x) = ci, i= 1, 2 0 otherwise

.

ただし,ci, γi, i= 1, 2は

Eθ0(X)] = α, Eθ0[T (X)δ(X)] = Eθ0[T (X)]α (*) をみたすように選ぶ.

Proof. 略証のみ与える.この証明では,pn(x; θ)は密度関数とし,微分と積分の順序交換を

自由に行う(ちゃんと正当化できる).さらに,(*)をみたすci, γi, i = 1, 2の存在は認めて, δがUMPU検定であることを示す.δを水準αの任意の不偏検定とすると,

βδ0) ≤ α, βδ(θ) ≥ α, ∀θ ̸= θ0

となる.ここで,βδ(θ) = Eθ[δ(X)]はθについて微分可能であることが示せる.βδ(θ)の連 続性より,βδ0) = αであり,βδ(θ)はθ= θ0で最小になる.よって,

βδ0) = α, 0 = βδ0) =

δ(x)

∂θpn(x; θ0)dx = Eθ0[T (X)δ(X)] − C0

となる.ここで,eC(θ) = eθT(x)H(x)dxの両辺をθについて微分して,C(θ) = Eθ[T (X)] を得る.従って,

δ(x)

∂θpn(x; θ0)dx = 0 ⇔ Eθ0[T (X)δ(X)] = Eθ0[T (X)]α である.

次に,θ1 ̸= θ0を任意に固定して,a1, a2∈ Rに対して,

ra1,a2(x) = pn(x; θ1) − a1pn(x; θ0) − a2

∂θpn(x; θ0)

= pn(x; θ0)eC(θ0)−C(θ1){e1−θ0)T (x)− ea1− ea2T(x)}

1

(2)

とおく.ここで,ea1 = eC(θ1)−C(θ0){a1 − C0)}, ea2 = eC(θ1)−C(θ0)a2である.c1, c2に対し て,a1, a2をc1, c2, θ0, θ1に依存させて適当に選べば,H(x) > 0なるxに対して,

δ(x) = 1 ⇔ ra1,a2(x) > 0, δ(x) = 0 ⇔ ra1,a2(x) < 0 が成り立つ.よって,

(x) − δ(x)}ra1,a2(x) ≥ 0 となるから,両辺を積分して,

βδ1) − βδ1) =

(x) − δ(x)}pn(x; θ1)dx

≥ a1

(x) − δ(x)}pn(x; θ0)dx + a2

(x) − δ(x)}

∂θpn(x; θ0)dx = 0

を得る.さらに,δ(x) ≡ αなる検定と比較して,βδ1) ≥ αを得る.θ1 ̸= θ0は任意だった から,δがUMPU検定であることが示された.

2

参照

関連したドキュメント

出典:第40回 広域系統整備委員会 資料1 出典:第50回 広域系統整備委員会 資料1.

マニピュレータで、プール 内のがれきの撤去や燃料取 り出しをサポートする テンシルトラスには,2本 のマニピュレータが設置さ

マニピュレータで、プール 内のがれきの撤去や燃料取 り出しをサポートする テンシルトラスには,2本 のマニピュレータが設置さ

[r]

章番号 ページ番号 変更後 変更前 変更理由.. 1 補足説明資

KK7 補足-028-08 「浸水 防護施設の耐震性に関す る説明書の補足説明資料 1.2 海水貯留堰における 津波波力の設定方針につ

2月 1月 12月 11月 10月 9月. 8月

3R・適正処理の促進と「持続可能な資源利用」の推進 自然豊かで多様な生きものと 共生できる都市環境の継承 快適な大気環境、良質な土壌と 水循環の確保 環 境 施 策 の 横 断 的 ・ 総