MLE du paramètre de localisation dans une distribution de Cauchy

Après centrage, les deux mesures x et −x peuvent être supposées être des observations indépendantes d'une distribution de Cauchy avec fonction de densité de probabilité:

$f(x :\theta) =$ $1\over\pi (1+(x-\theta)^2)$ $, -∞ < x < ∞$

Montrer que si le MLE de est 0, mais si $x^2≤ 1$ $\theta$ $x^2>1$ il y a deux MLE de , égaux à ± $\theta$ $\sqrt {x^2-1}$

Je pense que pour trouver le MLE, je dois différencier la probabilité de journal:

$dl\over d\theta$ $=\sum$ $2(x_i-\theta)\over 1+(x_i-\theta)^2$ $=$ + $2(-x-\theta)\over 1+(-x-\theta)^2$ $2(x-\theta)\over 1+(x-\theta)^2$ $=0$

Donc,

$2(x-\theta)\over 1+(x-\theta)^2$ $=$ $2(x+\theta)\over 1+(x-\theta)^2$

que j'ai ensuite simplifié

$5x^2 = 3\theta^2+2\theta x+3$

Maintenant, j'ai heurté un mur. Je me suis probablement trompé à un moment donné, mais de toute façon, je ne sais pas comment répondre à la question. Quelqu'un peut-il aider?

— user123965
source

Veuillez expliquer pourquoi vous avez divisé x en -x et + x? C'est mes devoirs et je suis coincé à cette étape. Je suppose que vous y avez appliqué la méthode Raphson de Newton. Mais je ne sais pas comment l'appliquer. Pourriez-vous me le dire?

— user89929

Il y a une faute de frappe dans vos calculs. La condition de premier ordre pour un maximum est:

\begin{aligned} \frac{\partial L}{\partial θ} = 0 & \Rightarrow \frac{2 (x + θ)}{1 + (x + θ)^{2}} - \frac{2 (x - θ)}{1 + (x - θ)^{2}} & = 0 \\ \Rightarrow (x + θ) + (x + θ) (x - θ)^{2} - (x - θ) - (x - θ) (x + θ)^{2} & = 0 \\ \Rightarrow 2 θ + (x + θ) (x - θ) [x - θ - (x + θ] & = 0 \\ \Rightarrow 2 θ - 2 θ (x + θ) (x - θ) = 0 \Rightarrow 2 θ - 2 θ (x^{2} - θ^{2}) & = 0 \\ \Rightarrow 2 θ (1 - x^{2} + θ^{2}) = 0 \Rightarrow 2 θ (θ^{2} + (1 - x^{2})) & = 0 \end{aligned}

$\begin{align} \frac {\partial L}{\partial \theta}= 0 &\Rightarrow \frac {2(x+\theta)}{ 1+(x+\theta)^2} - \frac{2(x-\theta)}{ 1+(x-\theta)^2}&=0 \\[5pt] &\Rightarrow (x+\theta)+(x+\theta)(x-\theta)^2 - (x-\theta)-(x-\theta)(x+\theta)^2&=0 \\[3pt] &\Rightarrow 2\theta +(x+\theta)(x-\theta)\left[x-\theta-(x+\theta\right]&=0 \\[3pt] &\Rightarrow2\theta -2\theta(x+\theta)(x-\theta) =0\Rightarrow 2\theta -2\theta(x^2-\theta^2)&=0 \\[3pt] &\Rightarrow2\theta(1-x^2+\theta^2)=0 \Rightarrow 2\theta\big(\theta^2+(1-x^2)\big)&=0 \end{align}$

If $x^2\leq 1$ then the term in the parenthesis cannot be zero (for real solutions of course), so you are left only with the solution $\hat \theta =0$ .

If $x^2 >1$ you have $2\theta\big[\theta^2-(x^2-1)\big]=0$ so, apart from the candidate point $\theta =0$ you also get

\frac{\partial L}{\partial θ} = 0, for \hat{θ} = \pm \sqrt{x^{2} - 1}

$\frac {\partial L}{\partial \theta}= 0,\;\; \text{for}\;\;\hat \theta = \pm\sqrt {x^2-1}$

You also have to justify why in this case $\hat \theta =0$ is no longer an MLE.

ADDENDUM

For $x =\pm 0.5$ the graph of the log-likelihood is enter image description here

while for $x =\pm 1.5$ the graph of the log-likelihood is, enter image description here

Now all you have to do is to prove it algebraically and then wonder "fine -now which of the two should I choose?"

— Alecos Papadopoulos
source

Thanks! I can't see why

θ = 0

$\theta=0$ would no longer be an MLE though

— user123965

Work the 2nd order condition for a maximum, or evaluate the likelihood at the candidate solutions

— Alecos Papadopoulos

+1 great answer. Also, this might be interesting: wolframalpha.com/share/… wolframalpha.com/share/…

— random_user

@random_user Thanks! - I took the liberty to incorporate the plot in the answer.

— Alecos Papadopoulos

2nd derivative positive so indeed a local minimum

— Alecos Papadopoulos