Posterior for Pareto distribution with Jeffreys prior

bayesianprobabilityprobability distributions

For some known parameter m, the data is IID pareto distribution: $X_1,..,X_n \sim \text{Pareto}(\theta, m)$

$f(x | \theta) = \theta m ^\theta x^{-(\theta + 1)} \textbf{1}{\{m < x \}}$

I need to find a posterior for the Jeffreys prior:
$$\pi(\theta) \propto I(\theta)^{1/2}$$

I calculated the prior as follows:

Fisher information:
$$I(\theta) = -E\bigg[\frac{\partial^2\log f(X|\theta)}{\partial\theta^2} \bigg]$$

$$\log f(X|\theta) = \log \theta + \theta \log m – (\theta + 1)\log x$$

$$ I(\theta) = -\frac{1}{\theta^2}$$

Now, in order to obtain the posterior, I'm trying to use the fact that the conjugate prior for Pareto is Gamma, and the Jeffreys prior is the limiting case of the conjugate one, but I can't arrive at an appropriate limit. Does this strategy make any sense? Simply multiplying the Likelihood with the obtained Jeffreys prior doesn't seem to work. Any hints highly aprreciated!

Best Answer

it is not necessary to calucalte the prior with the normalizing constant, thus you can leave

$$\pi(\theta)\propto \frac{1}{\theta}$$

(the minus sign is an obvious error given that in your Fisher's information you have $-E\{\dots\}$)

Noe let's focus on the likelihood (any term not depending on $\theta$ can be canceled)

$$p(\mathbf{x}|\theta)\propto \theta^n\cdot\left( \frac{\Pi_ix_i}{m^n}\right)^{-\theta}=\theta^n\cdot\exp\{-\theta[\Sigma_i\log x_i-n\log m]\}$$

thus the posterior is

$$\pi(\theta|\mathbf{x})\propto \theta^{n-1}\cdot\exp\{-\theta[\Sigma_i\log x_i-n\log m]\}$$

that is a Gamma

$$\pi(\theta|\mathbf{x})\sim\text{Gamma}[n;\Sigma_i\log x_i-n\log m]$$

Related Question