best unbiased estimator of θ² for a Poisson model

A mostly traditional question on X validated about the “best” [minimum variance] unbiased estimator of θ² from a Poisson P(θ) sample leads to the Rao-Blackwell solution

\mathbb{E}[X_1X_2|\underbrace{\sum_{i=1}^n X_i}_S=s] = -\frac{s}{n^2}+\frac{s^2}{n^2}=\frac{s(s-1)}{n^2}

and a similar estimator could be constructed for θ³, θ⁴, … With the interesting limitation that this procedure stops at the power equal to the number of observations (minus one?). But,  since the expectation of a power of the sufficient statistics S [with distribution P(nθ)] is a polynomial in θ, there is de facto no limitation. More interestingly, there is no unbiased estimator of negative powers of θ in this context, while this neat comparison on Wikipedia (borrowed from the great book of counter-examples by Romano and Siegel, 1986, selling for a mere $180 on amazon!) shows why looking for an unbiased estimator of exp(-2θ) is particularly foolish: the only solution is (-1) to the power S [for a single observation]. (There is however a first way to circumvent the difficulty if having access to an arbitrary number of generations from the Poisson, since the Forsythe – von Neuman algorithm allows for an unbiased estimation of exp(-F(x)). And, as a second way, as remarked by Juho Kokkala below, a sample of at least two Poisson observations leads to a more coherent best unbiased estimator.)

2 Responses to “best unbiased estimator of θ² for a Poisson model”

  1. Unless I’m very confused about something, -1 to the power S is not an unbiased estimator of exp(-theta). exp(-theta) is P(X_1 = 0; theta) so I_{X_1 =0} would be an unbiased estimator, as well as its Rao-Blackwellization (1 – 1/N)^S (

    As far as I understand, the linked “neat comparison” in Wikipedia is actually about (-1)^S being the only unbiased estimator of exp(-2*theta) in the case N=1. For N>=2, there is I_{X_1=0, X_2=0} and its Rao-Blackwellization (1 – 2/N)^S.

    • Thanks for correcting my confused explanation. Indeed, (-1) to the power S is an unbiased estimator of exp(-2Nθ) when S is the sum and indeed there are much more efficient solutions in that case, since the Rao-Blackwellised solution is the UMVUE.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.

%d bloggers like this: