MLE của

Đặt là một mẫu ngẫu nhiên từ một bản phân phối với pdf $X_{1},X_{2},X_{3},...,X_{n}$

$$f(x;\alpha,\theta)=\frac{e^{-x/\theta}}{\theta^{\alpha}\Gamma(\alpha)}x^{\alpha-1}I_{(0,\infty)}(x ),\alpha,\theta>0$$

Tìm công cụ ước tính khả năng tối đa của và . Đặt$\alpha$$\theta$$\Psi(\alpha)=\frac{d\Gamma(\alpha)}{d\alpha}$

Nỗ lực của tôi,

$$\begin{eqnarray*} \mathcal{L}(\alpha,\theta)&=&\prod_{i=1}^{n}f(x_i)\\ &=&\prod_{i=1}^{n}\frac{e^{-x_i/\theta}}{\theta^{\alpha}\Gamma(\alpha)}x_i^{\alpha-1}\\ &=&\frac{1}{\Gamma^{n}(\alpha)\cdot \theta^{n \alpha}}(\prod_{i=1}^{n}x_i)^{\alpha-1}\exp(-\sum_{i=1}^{n}\frac{x_i}{\theta}) \end{eqnarray*}$$ $$\begin{eqnarray*} \ell(\alpha,\theta)&=&-n\log(\Gamma(\alpha))-n\alpha\log(\theta)+(\alpha-1)\sum_{i=1}^{n}\log(x_i)-\frac{1}{\theta}\sum_{i=1}^{n}x_i\\ \frac{\delta \ell(\alpha,\theta)}{\delta \theta}&=&-\frac{n\alpha}{\theta}+\frac{1}{\theta^2}\sum_{i=1}^{n}x_i=0\\ \frac{1}{\theta^2}\sum_{i=1}^{n}x_i&=&\frac{n\alpha}{\theta}\\ \hat{\theta}&=&\frac{\sum_{i=1}^{n}x_i}{n\alpha}\\ &=&\frac{1}{\alpha}\bar{x}\\ \end{eqnarray*}$$ $$\begin{eqnarray*} \frac{d \ell(\alpha,\hat{\theta})}{d\alpha}&=&\frac{-n \cdot \Gamma'(\alpha)}{\Gamma(\alpha)}-n\log(\frac{1}{\alpha}\bar{x})+\sum_{i=1}^{n}\log(x_i)=0\\ &=&\frac{-n \cdot \Gamma'(\alpha)}{\Gamma(\alpha)}+n\log(\alpha)-n\log(\bar{x})+\sum_{i=1}^{n}\log(x_i)=0\\ \log(\alpha)-\frac{\Gamma'(\alpha)}{\Gamma(\alpha)}&=&\log(\bar{x})-\frac{\sum_{i=1}^{n}\log(x_i)}{n} \end{eqnarray*}$$

Tôi không thể tiếp tục tìm kiếm . Thứ hai, tôi không biết cách sử dụng \ Psi (\ alpha) = \ frac {d \ Gamma (\ alpha)} {d \ alpha} như được nêu trong câu hỏi. Hy vọng ai đó có thể giải thích cho tôi.$\alpha$$\Psi(\alpha)=\frac{d\Gamma(\alpha)}{d\alpha}$

Cảm ơn trước.

Đặt vì vậy là hàm digamma (Tôi đang sử dụng chứ không phải của bạn ).$\psi(\alpha) = \frac{\Gamma'(\alpha)}{\Gamma(\alpha)}$$\psi$$\psi$$\Psi$

Theo bất đẳng thức AM-GM nên (trong đó và được xác định gần như chắc chắn). Hơn nữa, đẳng thức chỉ giữ cho là sự kiện có xác suất , vì vậy gần như chắc chắn.

$$\bar x \geq \left(\prod_i x_i\right)^{1/n}$$ $$\log \bar x - \overline{\log x} \geq 0$$ $\log \bar x$$\log x_i$$x_1=\dots=x_n$$0$$\log \bar x - \overline{\log x} > 0$

Để đơn giản, tôi sẽ lấy .$y = \log \bar x - \overline{\log x}$

Xem xét trên . Đây là liên tục và do đó theo định lý giá trị trung gian đạt mọi số thực trong . Cụ thể, điều này có nghĩa là tức là có ít nhất một điểm trong được ánh xạ tới , vì .$f(\alpha) = \log(\alpha) - \psi(\alpha)$$(0,\infty)$

$$\lim_{\alpha\to 0} f(\alpha) = \infty$$ $$\lim_{\alpha\to\infty} f(\alpha) = 0$$ $f$$(0,\infty)$ $$f^{-1}\left(\left\{y\right\}\right) \neq \emptyset$$ $(0,\infty)$$y$$y > 0$

Hơn nữa, hóa ra được tiêm vào là nên thực sự có một duy nhất với .$f$$(0,\infty)$$f' < 0$$\hat \alpha$$f(\hat\alpha) = y$

Trên thực tế, việc tìm kiếm này sẽ yêu cầu các phương thức số, như @StubbornAtom nói.$\hat \alpha$