Метрополис Гастингс моделирование в R

Я пытаюсь внедрить Metropolis Hastings в R с помощью пакета MHadaptive. Но у меня есть некоторые проблемы. Почему эта ошибка появляется? объект 'z' не найден. И в результате должно появиться значение оптимизации для z.

q<-as.matrix(dataset1)
{
li_reg<-function(pars,data)
 z<-pars[1]
 pred<-(1/(1+exp((-z*(data[,2] +1)*(data[,2]-1)))/0.3208))
log_likelihood<-sum(dnorm(data[,2],pred,log = TRUE))
prior<-prior_reg(pars)
return(log_likelihood+prior)
}
prior<-function(pars)
{
 z<-pars[1]
return(prior_z)
}
mcmc_r<-Metro_Hastings(li_func = li_reg,pars = c(30),prop_sigma = NULL,par_names = c('z'),data=q,iterations = 2000,burn_in = 1000,adapt_par = c(100,20,0.5,0.75),quiet = FALSE)
mcmc_r<-mcmc_thin(mcmc_r)

0 ответов

Другие вопросы по тегам