Вычислить апостериорную вероятность с учетом вероятности, распределенной по Бернулли

Предложите прочитать это:

https://css-tricks.com/just-one-of-those-weird-things-about-css-background-on-body/

По существу, в отсутствие элемента background на элементе html, body background будет охватывать страницу. Если в элементе html есть background, body background ведет себя точно так же, как и любой другой элемент.

1
задан merv 16 January 2019 в 21:08
поделиться

1 ответ

Этот ответ был предоставлен через раздел комментариев @JosephClarkMcIntyre.

Вот краткое изложение:

  • В испытании Бернулли, когда N - общее количество испытаний - и z - общее количество успеха велико, а основной параметр тета равен маленький, лучше работать только в пространстве журнала и никогда не брать экспоненциальный.
  • Кроме того, поскольку функция логарифма увеличивается, сравнение логарифмических постеров двух распределений аналогично сравнению апостериорного.
  • Вышеуказанная реализация была неправильной, потому что формула для вычисления доказательств не верна. p (доказательство) = сумма (вероятность * предшествующая), p (log_evidence) = сумма (log_likelihood + log_prior)

Это окончательный код, в котором предшествующее значение, вероятность и свидетельство находятся в пространстве журнала. :

  a = 1  # a and b are the beta distribution's paramteres
  b= 1
  num_steps = 1e5
  z= 17220 #Number of heads
  N= 143293 #Total number of flips

  Theta = seq(from=0.07,to=0.12, length.out= num_steps)
  lprior = dbeta(Theta, a,b,log=TRUE) #Compute the log prior at each value 

  llikelihood = log(Theta)*z + log1p(-Theta)*(N-z) #log likelihood

  lpData = sum(llikelihood + lprior) #compute log of the evidence

  lposterior = llikelihood+lprior - lpData
  plot(Theta,log(dbeta(Theta,a+z,N-z+b)))
  plot(Theta, lposterior, type="l")

Тем не менее, аналитические и вычисленные логарифмические аппостеры не совпадают с показанными на графике.

Не стесняйтесь комментировать, если считаете, что в этом ответе есть изъян или объясните, почему аналитический и вычисленный логарифмический апостериор не совпадают. ^^

0
ответ дан Salma Bouzid 16 January 2019 в 21:08
поделиться
Другие вопросы по тегам:

Похожие вопросы: