免費(fèi)的域名解析楓樹seo網(wǎng)
文章目錄
- 點(diǎn)估計(jì)和估計(jì)量的求法
- 點(diǎn)估計(jì)概念
- 矩估計(jì)法
- 極大似然估計(jì)法
- 參考文獻(xiàn)
參數(shù)估計(jì)是數(shù)理統(tǒng)計(jì)中重要的基本問題之一。通常,稱參數(shù)的可容許值的全體為參數(shù)空間,并記為 Θ \Theta Θ。所謂參數(shù)估計(jì)就是由樣本對總體分布所含的未知參數(shù)做出估計(jì)。另外,在有些實(shí)際問題中,由于事先并不知道總體 X X X 的分布類型,而要對其某些數(shù)字特征,如均值、方差等做出估計(jì),習(xí)慣上也把這些數(shù)字特征稱為參數(shù),對它們進(jìn)行估計(jì)也屬于參數(shù)估計(jì)范疇。
點(diǎn)估計(jì)和估計(jì)量的求法
點(diǎn)估計(jì)概念
設(shè)總體 X X X 的分布函數(shù)是 F ( x ; θ 1 , . . . , θ l ) F(x;\theta_1,...,\theta_l) F(x;θ1?,...,θl?),其中 θ 1 , . . . , θ l \theta_1,...,\theta_l θ1?,...,θl? 是未知參數(shù), X 1 , . . . , X n X_1,...,X_n X1?,...,Xn? 是來自總體 X X X 的樣本, x 1 , . . . , x n x_1,...,x_n x1?,...,xn? 是相應(yīng)的樣本值,參數(shù)點(diǎn)估計(jì)就是研究如何構(gòu)造適當(dāng)?shù)慕y(tǒng)計(jì)量 θ ^ i ( X 1 , . . . , X n ) \hat{\theta}_i(X_1,...,X_n) θ^i?(X1?,...,Xn?),并分別用觀察值 θ ^ i ( x 1 , . . . , x n ) \hat{\theta}_i(x_1,...,x_n) θ^i?(x1?,...,xn?) 作為未知參數(shù) θ i \theta_i θi? 的估計(jì)。
通常,稱用作估計(jì)的統(tǒng)計(jì)量 θ ^ i ( X 1 , . . . , X n ) \hat{\theta}_i(X_1,...,X_n) θ^i?(X1?,...,Xn?) 為估計(jì)量,稱其觀察值 θ ^ i ( x 1 , . . . , x n ) \hat{\theta}_i(x_1,...,x_n) θ^i?(x1?,...,xn?) 為估計(jì)值。
由于對不同的樣本值,得到的參數(shù)估計(jì)值往往不同,因此,點(diǎn)估計(jì)問題的關(guān)鍵在于構(gòu)造估計(jì)量的方法。下面介紹求估計(jì)量的一些常用方法。
矩估計(jì)法
設(shè)總體 X X X 的分布中含有 l l l 個(gè)未知參數(shù) θ 1 , . . . , θ l \theta_1,...,\theta_l θ1?,...,θl?,又設(shè)總體 X X X 的前 l l l 階原點(diǎn)矩 α k = E ( X k ) ( k = 1 , . . . , l ) \alpha_k=E(X^k)(k=1,...,l) αk?=E(Xk)(k=1,...,l) 存在,且是 θ 1 , . . . , θ l \theta_1,...,\theta_l θ1?,...,θl? 的函數(shù),即 α k = α k ( θ 1 , . . . , θ l ) \alpha_k=\alpha_k(\theta_1,...,\theta_l) αk?=αk?(θ1?,...,θl?),令
α k ( θ ^ 1 , . . . , θ ^ l ) = A k , k = 1 , . . . , l \alpha_k(\hat{\theta}_1,...,\hat{\theta}_l)=A_k,\quad k=1,...,l αk?(θ^1?,...,θ^l?)=Ak?,k=1,...,l
解此方程組可得 θ ^ 1 , . . . , θ ^ l \hat{\theta}_1,...,\hat{\theta}_l θ^1?,...,θ^l?,并將它們分別作為 θ 1 , . . . , θ l \theta_1,...,\theta_l θ1?,...,θl? 的估計(jì)量。這種求估計(jì)量的方法稱為矩估計(jì)法,用矩估計(jì)法求得的估計(jì)量稱為矩估計(jì)量。
例:設(shè)總體 X X X 的二階矩存在, X 1 , . . . , X n X_1,...,X_n X1?,...,Xn? 為總體 X X X 的樣本,求總體均值 μ \mu μ 與總體方差 σ 2 \sigma^2 σ2 的矩估計(jì)。
解:因 α 1 = μ , α 2 = σ 2 + μ 2 \alpha_1=\mu, \alpha_2=\sigma^2+\mu^2 α1?=μ,α2?=σ2+μ2,令 { μ ^ = A 1 = X ˉ σ ^ 2 + μ ^ 2 = A 2 = 1 n ∑ i = 1 n X i 2 \begin{cases} \hat{\mu}=A_1=\bar{X} \\ \hat{\sigma}^2+\hat{\mu}^2=A_2=\frac{1}{n}\sum_{i=1}^n X_i^2 \end{cases} {μ^?=A1?=Xˉσ^2+μ^?2=A2?=n1?∑i=1n?Xi2??
解得 μ \mu μ 與 σ 2 \sigma^2 σ2 的矩估計(jì)分別為
μ ^ = X ˉ \hat{\mu}=\bar{X} μ^?=Xˉ σ ^ 2 = A 2 ? X ˉ 2 = S 2 \hat{\sigma}^2=A_2-\bar{X}^2=S^2 σ^2=A2??Xˉ2=S2
極大似然估計(jì)法
以下用 X = ( X 1 , . . . , X n ) T \boldsymbol{X}=(X_1,...,X_n)^T X=(X1?,...,Xn?)T 表示樣本, x = ( x 1 , . . . , x n ) T \boldsymbol{x}=(x_1,...,x_n)^T x=(x1?,...,xn?)T 表示樣本點(diǎn), f ( x ; θ ) f(\boldsymbol{x};\theta) f(x;θ) 表示樣本分布。
極大似然法的提出是基于如下的想法:
當(dāng)給定 θ \theta θ 時(shí), f ( x ; θ ) f(\boldsymbol{x};\theta) f(x;θ) 度量樣本 X \boldsymbol{X} X 在 x \boldsymbol{x} x 點(diǎn)發(fā)生的可能性。對于樣本空間中的兩個(gè)不同樣本點(diǎn) x 1 , x 2 ∈ X \boldsymbol{x}_1, \boldsymbol{x}_2 \in \mathcal{X} x1?,x2?∈X,如果有 f ( x 1 ; θ ) > f ( x 2 ; θ ) f(\boldsymbol{x}_1;\theta) > f(\boldsymbol{x}_2;\theta) f(x1?;θ)>f(x2?;θ),自然會認(rèn)為樣本 X \boldsymbol{X} X 更可能在 x 1 \boldsymbol{x}_1 x1? 點(diǎn)發(fā)生。
現(xiàn)在換個(gè)角度來看待 f ( x ; θ ) f(\boldsymbol{x};\theta) f(x;θ)。當(dāng)給定樣本點(diǎn) x \boldsymbol{x} x 時(shí),對參數(shù)空間中的兩個(gè)不同參數(shù) θ 1 , θ 2 ∈ Θ \theta_1,\theta_2 \in \Theta θ1?,θ2?∈Θ,如果有 f ( x ; θ 1 ) > f ( x ; θ 2 ) f(\boldsymbol{x};\theta_1) > f(\boldsymbol{x};\theta_2) f(x;θ1?)>f(x;θ2?),那么會認(rèn)為樣本點(diǎn) x \boldsymbol{x} x 更像是來自總體 f ( X ; θ 1 ) f(\boldsymbol{X};\theta_1) f(X;θ1?),所以,數(shù) f ( x ; θ ) f(\boldsymbol{x};\theta) f(x;θ) 的大小可作為參數(shù) θ \theta θ 對產(chǎn)生樣本觀察值 x \boldsymbol{x} x 有多大似然性的一種度量。
當(dāng)給定樣本點(diǎn) x \boldsymbol{x} x 時(shí),稱 f ( x ; θ ) f(\boldsymbol{x};\theta) f(x;θ) 為 θ \theta θ 的似然函數(shù),記為 L ( θ ; x ) L(\theta;\boldsymbol{x}) L(θ;x),即
L ( θ ; x ) = f ( x ; θ ) = { ∏ i = 1 n p ( x i ; θ ) , 總體 X 為離散型隨機(jī)變量 ∏ i = 1 n f ( x i ; θ ) , 總體 X 為連續(xù)型隨機(jī)變量 L(\theta;\boldsymbol{x})=f(\boldsymbol{x};\theta)=\begin{cases} \prod_{i=1}^np(x_i;\theta), & 總體 X 為離散型隨機(jī)變量 \\ \prod_{i=1}^nf(x_i;\theta), & 總體 X 為連續(xù)型隨機(jī)變量 \end{cases} L(θ;x)=f(x;θ)={∏i=1n?p(xi?;θ),∏i=1n?f(xi?;θ),?總體X為離散型隨機(jī)變量總體X為連續(xù)型隨機(jī)變量?
而稱 ln ? f ( x ; θ ) \ln f(\boldsymbol{x};\theta) lnf(x;θ) 為對數(shù)似然函數(shù),記為 ln ? L ( θ ; x ) \ln L(\theta;\boldsymbol{x}) lnL(θ;x)。
若有統(tǒng)計(jì)量 θ ^ ? θ ^ ( X ) \hat{\theta}\bumpeq \hat{\theta}(\boldsymbol{X}) θ^?θ^(X),使得
L ( θ ^ ( x ) ; x ) = sup ? θ ∈ Θ { L ( θ ; x ) } L(\hat{\theta}(\boldsymbol{x});\boldsymbol{x})=\sup_{\theta \in \Theta}\{L(\theta;\boldsymbol{x})\} L(θ^(x);x)=θ∈Θsup?{L(θ;x)}
或等價(jià)的,使得
ln ? L ( θ ^ ( x ) ; x ) = sup ? θ ∈ Θ { ln ? L ( θ ; x ) } \ln L(\hat{\theta}(\boldsymbol{x});\boldsymbol{x})=\sup_{\theta \in \Theta}\{\ln L(\theta;\boldsymbol{x})\} lnL(θ^(x);x)=θ∈Θsup?{lnL(θ;x)}
則稱 θ ^ ( X ) \hat{\theta}(\boldsymbol{X}) θ^(X) 為參數(shù) θ \theta θ 的極大似然估計(jì)量(Maximum Likelihood Estimators, MLE)。
例:設(shè)總體 X ~ P ( λ ) , λ > 0 X \sim P(\lambda),\lambda>0 X~P(λ),λ>0,試求參數(shù) λ \lambda λ 的極大似然估計(jì)量。
解: X X X 的概率函數(shù)為
P { X = x } = λ x x ! e ? λ , x = 0 , 1 , 2 , . . . P\{X=x\}=\frac{\lambda^x}{x!}e^{-\lambda},\quad x=0,1,2,... P{X=x}=x!λx?e?λ,x=0,1,2,...
故 λ \lambda λ 的似然函數(shù)為
L ( λ ) = ∏ i = 1 n ( λ x i x i ! e ? λ ) = e ? n λ λ ∑ i = 1 n x i ∏ i = 1 n ( x i ! ) L(\lambda)=\prod_{i=1}^n (\frac{\lambda^{x_i}}{x_i!}e^{-\lambda})=e^{-n\lambda}\frac{\lambda^{\sum_{i=1}^nx_i}}{\prod_{i=1}^n(x_i!)} L(λ)=i=1∏n?(xi?!λxi??e?λ)=e?nλ∏i=1n?(xi?!)λ∑i=1n?xi??
對數(shù)似然函數(shù)為
ln ? L ( λ ) = ? n λ + ln ? λ ∑ i = 1 n x i ? ∑ i = 1 n ln ? ( x i ! ) \ln L(\lambda)=-n\lambda+\ln \lambda \sum_{i=1}^nx_i-\sum_{i=1}^n \ln(x_i!) lnL(λ)=?nλ+lnλi=1∑n?xi??i=1∑n?ln(xi?!)
令
? ln ? L ( λ ) ? λ = ? n + 1 λ ∑ i = 1 n x i = 0 \frac{\partial \ln L(\lambda)}{\partial \lambda}=-n+\frac{1}{\lambda}\sum_{i=1}^nx_i=0 ?λ?lnL(λ)?=?n+λ1?i=1∑n?xi?=0
該似然方程有唯一解 λ ^ = 1 n ∑ i = 1 n x i = x ˉ \hat{\lambda}=\frac{1}{n}\sum_{i=1}^nx_i=\bar{x} λ^=n1?∑i=1n?xi?=xˉ,又因
? 2 ln ? L ( λ ) ? λ 2 ∣ λ = x ˉ < 0 \frac{\partial^2 \ln L(\lambda)}{\partial \lambda^2}|_{\lambda=\bar{x}}<0 ?λ2?2lnL(λ)?∣λ=xˉ?<0
故 λ \lambda λ 的極大似然估計(jì)量為 λ ^ = X ˉ \hat{\lambda}=\bar{X} λ^=Xˉ。
參考文獻(xiàn)
[1] 《應(yīng)用數(shù)理統(tǒng)計(jì)》,施雨,西安交通大學(xué)出版社。