Existe um conjugado antes da distribuição de Laplace ? Caso contrário, existe uma expressão de forma fechada conhecida que se aproxima do posterior para os parâmetros da distribuição de Laplace?
Eu pesquisei bastante por aí, sem sucesso, então meu palpite atual é "não" nas perguntas acima ...
bayesian
conjugate-prior
laplace-distribution
Rasmus Bååth
fonte
fonte
Respostas:
Vamos examiná-los um de cada vez primeiro (considerando o outro como determinado).
No link (com a modificação de seguir a convenção de usar símbolos gregos para parâmetros):
- parâmetro de escala :
para certos valores de e S . Essa é a probabilidade é de forma gama inversa.k S
Portanto, o parâmetro de escala possui um conjugado anterior - por inspeção, o conjugado anterior é gama inversa.
- parâmetro de localização
Esta é, de fato, mais complicado, porque não se simplifica em algo conveniente em µ ; Não acho que exista nenhuma maneira de "coletar os termos" (de certa forma existe, mas não precisamos).∑Eu| xEu- μ | μ
Um anterior uniforme simplesmente truncará o posterior, o que não é tão ruim de se trabalhar, se isso parece plausível como anterior.
Uma possibilidade interessante que pode ocasionalmente ser útil é que é bastante fácil incluir um Laplace anterior (um com a mesma escala dos dados) usando uma pseudo-observação. Pode-se também aproximar-se de outro (mais rígido) antes através de várias pseudo-observações)
Também é flexível o suficiente para ser usado para aproximar outros anteriores.
(De um modo mais geral, ainda é possível trabalhar na escala de toras e usar um côncavo de toras linear contínuo e por partes, a anterior e a posterior também seriam dessa forma; isso incluiria Laplace assimétrico como um caso especial)
Exemplo
Apenas para mostrar que é muito fácil lidar com - abaixo é um anterior (cinza pontilhado), probabilidade (tracejado, preto) e posterior (sólido, vermelho) para o parâmetro de localização para um Laplace ponderado (... isso foi com escalas conhecidas )
A abordagem ponderada de Laplace funcionaria bem no MCMC, eu acho.
-
Gostaria de saber se o modo posterior resultante é uma mediana ponderada?
- na verdade (para responder minha própria pergunta), parece que a resposta é sim. Isso torna bastante agradável trabalhar.
-
Articulação prévia
Sem dúvida, algo mais geral para o conjunto anterior é bem possível, mas acho que não vou prosseguir com o caso conjunto além disso aqui.
-
Eu nunca tinha visto ou ouvido falar dessa abordagem anterior ponderada, mas era bastante simples de elaborar, por isso provavelmente já foi feito. (As referências são bem-vindas, se alguém souber de alguma.)
Se ninguém souber de nenhuma referência, talvez eu deva escrever alguma coisa, mas isso seria surpreendente.
fonte