Метрополис-Гастингс R интерпретация кода - PullRequest
Новая
       16

Метрополис-Гастингс R интерпретация кода

0 голосов
/ 23 марта 2020

У меня есть этот R-код для реализации алгоритма Метрополиса Гастингса: 

install.packages('coda')
library(coda)

rmh <- function(lamb,L,x0,s){
x <- numeric(L)
x[1] <- xold <- x0
for(i in 2:L){
  xnew <- xold + sample(c(-1,1),1)
  rho <- min(1,dpois(xnew,lamb)/dpois(xold,lamb)*(xnew>=0))
  u <- runif(1)
  x[i] <- xold <- (u<rho)*xnew + (u>rho)*xold
}
x}


x0 <- 0
L <- 10000
x <- rmh(10,L,x0,s)

plot(table(x)/L)
points(1:19,dpois(1:19,10),pch=2,col=2)

Меня спрашивают о нескольких вещах:

- опишите гистограмму, эргоди c график, автокорреляции, плотность, статистика c резюме.

Tbh Я немного растерялся в этом классе, кто-то может мне помочь, как поступить?

Спасибо заранее,

Алекс

...