跳到主要内容

法政大学 理工学研究科 システム理工学専攻 経営システム系 2024年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

[3] 確率密度関数 fX(x)f_X(x)

fX(x)=βxβ1(0<x<1)f_X(x) = \beta x^{\beta-1} \quad (0 < x < 1)

で与えられている確率分布からの無作為標本を (X1,,Xn)(X_1, \dots, X_n) とする。ただし、 0<β<0 < \beta < \infty である。 このとき、以下の問いに答えよ。

(1) E[ln(X)]E[\ln(X)]β\beta を用いて表せ。

(2) θ=1β\theta = \frac{1}{\beta} としたとき、 θ\theta の最尤推定量 θ^\hat{\theta} を示せ。

(3) θ^\hat{\theta}θ\theta の不偏推定量となるか否か示せ。

题目描述

【3】设 (X1,,Xn)(X_1,\dots,X_n) 是来自下列概率分布的随机样本:

fX(x)=βxβ1(0<x<1),f_X(x)=\beta x^{\beta-1}\qquad(0<x<1),

其中 0<β<0<\beta<\infty。回答下列问题。

(1)用 β\beta 表示 E[ln(X)]E[\ln(X)]

(2)令

θ=1β.\theta=\frac1\beta.

θ\theta 的最大似然估计量 θ^\hat{\theta}

(3)判断并说明 θ^\hat{\theta} 是否为 θ\theta 的无偏估计量。

Kai

(1)

E[ln(X)]=01ln(x)fX(x)dx=01ln(x)βxβ1dxE[\ln(X)] = \int_0^1 \ln(x) f_X(x) dx = \int_0^1 \ln(x) \beta x^{\beta - 1} dx

u=ln(x)u = \ln(x) , dv=βxβ1dxdv = \beta x^{\beta - 1} dx , 则 du=1xdxdu = \frac{1}{x} dx , v=xβv = x^{\beta} . 所以,

E[ln(X)]=[ln(x)xβ]0101xβ1xdx=[ln(x)xβ]0101xβ1dxE[\ln(X)] = [\ln(x) x^{\beta}]_0^1 - \int_0^1 x^{\beta} \frac{1}{x} dx = [\ln(x) x^{\beta}]_0^1 - \int_0^1 x^{\beta - 1} dx
=(ln(1)1βlimx0ln(x)xβ)[xββ]01= (\ln(1) \cdot 1^{\beta} - \lim_{x \to 0} \ln(x) x^{\beta}) - [\frac{x^{\beta}}{\beta}]_0^1
=00(1β0)=1β= 0 - 0 - (\frac{1}{\beta} - 0) = -\frac{1}{\beta}

因此, E[ln(X)]=1βE[\ln(X)] = -\frac{1}{\beta} .

(2) 似然函数为:

L(β)=i=1nfX(xi)=i=1nβxiβ1=βn(i=1nxi)β1L(\beta) = \prod_{i=1}^n f_X(x_i) = \prod_{i=1}^n \beta x_i^{\beta - 1} = \beta^n (\prod_{i=1}^n x_i)^{\beta - 1}

取对数:

lnL(β)=nln(β)+(β1)i=1nln(xi)\ln L(\beta) = n \ln(\beta) + (\beta - 1) \sum_{i=1}^n \ln(x_i)

β\beta 求导数:

lnL(β)β=nβ+i=1nln(xi)=0\frac{\partial \ln L(\beta)}{\partial \beta} = \frac{n}{\beta} + \sum_{i=1}^n \ln(x_i) = 0

因此,

nβ=i=1nln(xi)    β^=ni=1nln(xi)\frac{n}{\beta} = -\sum_{i=1}^n \ln(x_i) \implies \hat{\beta} = -\frac{n}{\sum_{i=1}^n \ln(x_i)}

因为 θ=1β\theta = \frac{1}{\beta} , 所以

θ^=1β^=i=1nln(xi)n=1ni=1nln(xi)\hat{\theta} = \frac{1}{\hat{\beta}} = -\frac{\sum_{i=1}^n \ln(x_i)}{n} = -\frac{1}{n} \sum_{i=1}^n \ln(x_i)

(3)

E[θ^]=E[1ni=1nln(Xi)]=1ni=1nE[ln(Xi)]=1ni=1n(1β)=1nnβ=1β=θE[\hat{\theta}] = E[-\frac{1}{n} \sum_{i=1}^n \ln(X_i)] = -\frac{1}{n} \sum_{i=1}^n E[\ln(X_i)] = -\frac{1}{n} \sum_{i=1}^n (-\frac{1}{\beta}) = \frac{1}{n} \frac{n}{\beta} = \frac{1}{\beta} = \theta

所以 θ^\hat{\theta}θ\theta 的不偏推定量.