Estimador de polarização do momento da distribuição lognormal

Estou fazendo um experimento numérico que consiste em amostrar uma distribuição lognormal $X\sim\mathcal{LN}(\mu, \sigma)$ e tentar estimar os momentos $\mathbb{E}[X^n]$ por dois métodos:

1. Olhando para a média amostral do $X^n$
2. Estimando $\mu$ e $\sigma^2$ usando as médias da amostra para $\log(X), \log^2(X)$ e depois usando o fato de que, para uma distribuição lognormal, temos $\mathbb{E}[X^n]=\exp(n \mu + (n \sigma)^2/2)$ .

A questão é :

Descobri experimentalmente que o segundo método tem um desempenho muito melhor que o primeiro, quando mantenho o número de amostras fixo e aumento $\mu, \sigma^2$ por algum fator T. Existe alguma explicação simples para esse fato?

Estou anexando uma figura na qual o eixo x é T, enquanto o eixo y são os valores de $\mathbb{E}[X^2]$ comparando os valores reais de $\mathbb{E}[X^2] = \exp(2 \mu + 2 \sigma^2)$ (laranja linha), para os valores estimados. método 1 - pontos azuis, método 2 - pontos verdes. o eixo y está na escala do log

EDITAR:

Abaixo está um código Mathematica mínimo para produzir os resultados para um T, com a saída:

   ClearAll[n,numIterations,sigma,mu,totalTime,data,rmomentFromMuSigma,rmomentSample,rmomentSample]
(* Define variables *)
n=2; numIterations = 10^4; sigma = 0.5; mu=0.1; totalTime = 200;
(* Create log normal data*)
data=RandomVariate[LogNormalDistribution[mu*totalTime,sigma*Sqrt[totalTime]],numIterations];

(* the moment by theory:*)
rmomentTheory = Exp[(n*mu+(n*sigma)^2/2)*totalTime];

(*Calculate directly: *)
rmomentSample = Mean[data^n];

(*Calculate through estimated mu and sigma *)
muNumerical = Mean[Log[data]]; (*numerical \[Mu] (gaussian mean) *)
sigmaSqrNumerical = Mean[Log[data]^2]-(muNumerical)^2; (* numerical gaussian variance *)
rmomentFromMuSigma = Exp[ muNumerical*n + (n ^2sigmaSqrNumerical)/2];

(*output*)
Log@{rmomentTheory, rmomentSample,rmomentFromMuSigma}

Saída:

(*Log of {analytic, sample mean of r^2, using mu and sigma} *)
{140., 91.8953, 137.519}

acima, o segundo resultado é a média amostral de , abaixo dos outros dois resultados$r^2$

Há algo intrigante nesses resultados, já que

1. o primeiro método fornece um estimador imparcial de , a saber $\mathbb{E}[X^2]$ temE[X2]como sua média. Portanto, os pontos azuis devem estar em torno do valor esperado (curva laranja); $$\frac{1}{N}\sum_{i=1}^N X_i^2$$ $\mathbb{E}[X^2]$
2. O segundo método proporciona um estimador tendencioso de , isto é, E [ exp ( n μ + n 2 σ 2 / 2 ) ] > exp ( n μ + ( n σ ) 2 / 2 ) quando$\mathbb{E}[X^2]$ $$\mathbb{E}[\exp(n \hat\mu + n^2 \hat{\sigma}^2/2)]>\exp(n \mu + (n \sigma)^2/2)$$ e σ ²são estimadores imparciais deμeσ$\hat\mu$$\hat\sigma²$$\mu$$\sigma²$ respectivamente, e é estranho que os pontos verdes estejam alinhados com a curva laranja.

$\mu_T$$\sigma_T$

Aqui está o código R correspondente:

moy1=moy2=rep(0,200)
mus=0.14*(1:200)
sigs=sqrt(0.13*(1:200))
tru=exp(2*mus+2*sigs^2)
for (t in 1:200){
x=rnorm(1e5)
moy1[t]=mean(exp(2*sigs[t]*x+2*mus[t]))
moy2[t]=exp(2*mean(sigs[t]*x+mus[t])+2*var(sigs[t]*x+mus[t]))}

plot(moy1/tru,col="blue",ylab="relative mean",xlab="T",cex=.4,pch=19)
abline(h=1,col="orange")
lines((moy2/tru),col="green",cex=.4,pch=19)

$\mu$$\sigma$$\mu$$\sigma$

$\mathbb{E}[X^2]$$X^2$$X^2$$e^{2\mu}$$X^2$$\exp\{2\mu+2\sigma\epsilon\}$$\epsilon\sim\mathcal{N}(0,1)$$\sigma$$\sigma\epsilon$$\sigma^2$$X$$\mathcal{LN}(\mu,\sigma)$

$$\begin{align*}\mathbb{P}(X^2>\mathbb{E}[X^2])&=\mathbb{P}(\log\{X^2\}>2\mu+2\sigma^2)\\&=\mathbb{P}(\mu+\sigma\epsilon>\mu+\sigma^2)\\&=\mathbb{P}(\epsilon>\sigma)\\ &=1-\Phi(\sigma)\end{align*}$$

Eu pensei em vomitar algumas figuras mostrando que as tramas de user29918 e Xi'an são consistentes. A Figura 1 mostra o que o usuário29918 fez e a Figura 2 (com base nos mesmos dados) faz o que Xi'an fez por seu gráfico. Mesmo resultado, apresentação diferente.

$\frac{1}{n} \sum_i x_i^2$

Comentários adicionais:

1. Um estimador imparcial não significa que o estimador deve estar próximo! Os pontos azuis não precisam estar próximos da expectativa. Por exemplo. uma única observação escolhida aleatoriamente fornece uma estimativa imparcial da média da população, mas não se espera que esse estimador esteja próximo.
2. A questão está surgindo à medida que a variação está se tornando absolutamente astronômica. À medida que a variação se torna uma bobagem, a estimativa para o primeiro método está sendo conduzida com apenas algumas observações. Você também começa a ter uma probabilidade minúscula de um número insana, insana e insanamente grande ...
3. $P(X^2 > E[X^2]) = 1 - \Phi(\sigma)$$\sigma$$X^2 > E[X^2]$