как генерировать случайные данные в winbugs - PullRequest
0 голосов
/ 09 июня 2019

Я дополнил свою модель приорами и функцией правдоподобия, но у меня нет реальных данных.Я хочу моделировать и генерировать случайные данные для оценки параметров.

model
{
#likelihood 
   for(i in 1:n){
      for(j in 1:m-1){
         x[i,j]~ dnorm( mu[i,j],  tau0[i])
         mu[i,j]<- mu1[i]*exp(-mu10/s)*(t[j+1]-t[j])*a+mu1[i]*exp(-mu10/s)*(tau[i]-t[j])*b+mu2[i]*exp(-mu20/s)*(t[j+1]-tau[i])*b+mu2[i]*exp(-mu20/s)*(t[j+1]-t[j])*c
         sigmatau0<-1/sqrt(tau0[i]*(t[j+1]-t[j]))      
         s~ dunif(1,100)                       
         a<-step( tau[i]-t[ j+1] )
         b<-(step(tau[i]- t[j])-(step(tau[i]- t[j+1])))
         c<-(1-step( tau[i]-t[j])

Я не знаю, как симулировать в winbugs и с нетерпением жду вашей помощи. Большое спасибо.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...