A distribuição normal bivariada com média e matriz de covariância pode ser reescrita em coordenadas polares com raio e ângulo . Minha pergunta é: Qual é a distribuição amostral de , isto é, a distância de um ponto ao centro estimado dada a matriz de covariância da amostra ?
Antecedentes: A verdadeira distância de um ponto à média segue uma distribuição de Hoyt . Com os valores próprios de e , seu parâmetro de forma é e seu parâmetro de escala é . Sabe-se que a função de distribuição cumulativa é a diferença simétrica entre duas funções Q Marcum.
A simulação sugere que a inserção de estimativas e para e no cdf verdadeiro funcione para amostras grandes, mas não para amostras pequenas. O diagrama a seguir mostra os resultados de 200 vezes
- simulando 20 vetores normais 2D para cada combinação de ( eixo ), (linhas) e quantil (colunas)
- para cada amostra, calculando o quantil fornecido do raio observado para
- Para cada amostra, o cálculo do quantil do Hoyt teórica (normal 2D) ED, e a partir do CDF teórico Rayleigh depois de ligar as estimativas de amostra e .
À medida que aproxima de 1 (a distribuição se torna circular), os quantis estimados de Hoyt se aproximam dos quantis estimados de Rayleigh que não são afetados por . À medida que o cresce, a diferença entre os quantis empíricos e os estimados aumenta, notadamente na cauda da distribuição.q ω
Respostas:
Como você mencionou em seu post, sabemos a distribuição da estimativa de se assim sabemos a distribuição da estimativa do valor verdadeiro . μ ^ r 2 t r u e r2rtrueˆ μ r2trueˆ r2
Queremos encontrar a distribuição de onde são expressos como vetores de coluna.xi
Agora fazemos o truque padrão
Observe que é o traço da matriz de covariância da amostra e depende apenas apenas da média da amostra . Assim, escrevemos como a soma de dois variáveis aleatórias independentes. Conhecemos as distribuições de e e, portanto, terminamos o truque padrão usando esse funções características são multiplicativas.r2ˆ S (x¯¯¯−μ)T(x¯¯¯−μ) x¯¯¯
Editado para adicionar:
Isso significa que o pdf de é||xi−μ||2
Para facilitar a notação definir , e .a=1−q44q2ω b=−(1+q2)24q2ω c=121+q2qω
A função geradora de momento de é||xi−μ||2
Assim, a função geradora de momento de é e a função de geração de momentos de ér2trueˆ
Isso implica que a função geradora de momento de ér2ˆ
A aplicação da conversão inversa de Laplace fornece que possui pdfr2ˆ
fonte