Continuo lendo em revistas de economia sobre um resultado específico usado em modelos de utilidade aleatórios. Uma versão do resultado é: se Gumbel ( , então:
onde é a constante de Euler-Mascheroni. Eu verifiquei que isso faz sentido usando R, e faz. O CDF para a distribuição Gumbel é:
Estou tentando encontrar uma prova disso e não tive sucesso. Eu tentei provar isso sozinho, mas não consigo passar por uma etapa específica.
Alguém pode me indicar uma prova disso? Caso contrário, talvez eu possa postar minha prova de tentativa até onde estou preso.
expected-value
gumbel
Jason
fonte
fonte
Respostas:
Agradeço o trabalho exibido na sua resposta: obrigado por essa contribuição. O objetivo deste post é fornecer uma demonstração mais simples. O valor da simplicidade é a revelação: podemos obter facilmente toda a distribuição do máximo, não apenas sua expectativa.
Ignore absorvendo-o no e assumindo que o tenha uma distribuição Gumbel . (Ou seja, substitua cada por e altere para .) Isso não altera a variável aleatóriaμ δi ϵi (0,1) ϵi ϵi−μ δi δi+μ
A independência do implica para todo real que é o produto das chances individuais . Pegar logs e aplicar propriedades básicas de exponenciais produzϵi x Pr(X≤x) Pr(δi+ϵi≤x)
Este é o logaritmo do CDF de uma distribuição Gumbel com o parâmetro de localização Isso é,λ=log∑ieδi.
Isso é muito mais informações do que o solicitado. A média dessa distribuição é implicaγ+λ,
QED.
fonte
Acontece que um artigo da Econometrica de Kenneth Small e Harvey Rosen mostrou isso em 1981, mas em um contexto muito especializado, de modo que o resultado requer muita escavação, sem mencionar algum treinamento em economia. Decidi provar isso da maneira que achar mais acessível.
Prova : SejaJ o número de alternativas. Dependendo dos valores do vetorϵ={ϵ1,...,ϵJ} , a função assume valores diferentes. Primeiro, foque nos valores de modo que . Ou seja, integraremos no conjunto :maxi(δi+ϵi) ϵ maxi(δi+ϵi)=δ1+ϵ1 δ1+ϵ1 M1≡{ϵ:δ1+ϵ1>δj+ϵj,j≠1}
O termo acima é o primeiro de desses termos em . Especificamente,J E[maxi(δi+ϵi)]
Agora aplicamos a forma funcional da distribuição Gumbel. Isto dá
onde a segunda etapa vem da coleta de um dos termos exponenciados no produto, junto com o fato de queδj−δi=0 se .i=j
Agora definimos e fazemos a substituição , de modo que e . Observe que, conforme aproxima do infinito, aproxima de 0, e como aproxima do infinito negativo, aproxima do infinito.Di≡∑jeδj−δi x=Dieμ−ϵi dx=−Dieμ−ϵidϵi⇒−dxDi=eμ−ϵidϵi ϵi=μ−log(xDi) ϵi x ϵi x
A Função Gamma é definida como . Para valores de que são números inteiros positivos, isso é equivalente a, então . Além disso, sabe-se que a constante de Euler – Mascheroni satisfazΓ(t)=∫∞0xt−1e−xdx t Γ(t)=(t−1)! Γ(1)=0!=1 γ≈0.57722
A aplicação desses fatos dá
Então somamos para obteri
Lembre-se de que . Observe que as probabilidades familiares de escolha do logit são inversas dos 's, ou seja, . Observe também que . Então nós temosDi=∑jeδj−δi=∑jeδjeδi Pi=eδi∑jδj Di Pi=1/Di ∑iPi=1
fonte