Estou tentando calcular esta distribuição posterior:
O problema é que o numerador, que é o produto de um monte de probabilidades é muito pequeno. (Meu é grande, cerca de 1500).n
Portanto, os valores posteriores para todos são calculados como 0 (estou fazendo cálculos em R).
Para esclarecer, cada tem seu próprio , juntos esses formam um vetor de elementos para . Cada tem seu próprio vetor de elemento de .
EDIT: Adicionando um exemplo de reprodução (para o numerador)
p <- sample(seq(0,1,by=0.01), 1500, replace=T)
y <- sample(c(0,1), 1500, replace=T)
dbern(y, p) # 1500-element vector, each element is < 1
prod(dbern(y, p)) # produce 0
exp(sum(log(dbern(y, p)))) # produce 0 since the sum is very negative