error in using WinBUGS from R -


the code following:

require(brugs) require(r2winbugs) model<-function(){   for(i in 1:n){   y[i] ~ dnorm(x[i], sigma.y)   }   x[1] ~ dnorm(theta[1], sigma.y)   theta[1] <- 0   for(j in 2:n){     x[j] ~ dnorm(theta[j], sigma.x)     theta[j] <- b*x[j-1] # row wrong,  # right when set theta[j] <- 1*x[j-1]   }   ~ dunif(0, 1)   b ~ dunif(-1, 1)   sigma.y ~ dgamma(0.1, 0.1)   sigma.x ~ dgamma(0.1, 0.1) } data <- list( n <- 100, y <- rnorm(100)) data=list(n=100,y=rnorm(100))  inits=function(){   list(sigma.x = rgamma(1,0.1,0.1), sigma.y = rgamma(1, 0.1, 0.1), = dnorm(1, 0, 1), b = dnorm(1, -1, 1)) } parameters=c("a", "b", "x")  write.model(model, con = "model.bug") modelcheck("model.bug") # model syntactically correct ret.sim <- bugs(data, inits, parameters, "model.bug",                 n.chains = 1, n.iter = 1000,                 program= "winbugs",                 working.directory = null,                 debug = t)

i don't know why, program correct when replace theta[j] <- b*x[j-1] theta[j] <- 1*x[j-1], have defined b ~ dunif(-1, 1). indeed, need set theta[j] <- - b*x[j-1] in final model, , turns out wrong when try add a , b it. find problem ?

the problems in priors b (and a). don't know data perhaps range of current priors not include true values of , b. think if use continuous distribution(s):

a ~ dnorm(0,1) b ~ dnorm(0,1)

your problem might solved?

n.b. if trying create ar(1) model winbugs might want check out tsbugs package.


Comments

Post a Comment

Popular posts from this blog

Why does Ruby on Rails generate add a blank line to the end of a file? -

running rails generate controller foo home results in: class foocontroller < application controller def home end end (there's nothing on line in actual file, blank line) is there purpose blank line? it's common use newline indicate end of line, last line. many editors (like vi, use) add in newline silently. consensus text files (especially in unix world) should end in newline, , there have been problems historically if wasn't present. why should text files end newline? the tool used count lines in file "wc" counts newlines in file, without trailing newline show 3 instead of 4. it improves readability of templates used in generator. consider: https://github.com/rails/rails/blob/master/railties/lib/rails/generators/rails/controller/templates/controller.rb to remove trailing newline, template have last line of: end<% end -%> instead of: end <% end -%> that seems less readable me. the discussion below refers ...
Read more

keyboard - Smiles and long press feature in Android -

i've created android keyboard. how add features inserting smilies, , using long keypresses display symbols? to create long-press behavior allows select several options pop-up menu, follow guide (specifically pop-up menu section): http://developer.android.com/guide/topics/ui/menus.html smileys (also known emoticons) combination of keys, understood programs , parsed same way. have do, when clicks smiley face on keyboard, simulate pressing of key combination matching emoticon. example, when clicks smiley, should simulate key presses ':' , ')' :). here's list of emoticons: http://en.wikipedia.org/wiki/list_of_emoticons
Read more

node.js - Bad Request - node js ajax post -

i'm using express 3.0 , node v0.11.0. have button submits form , i'm using ajax post data json object node server. client code is: $('#contact').submit(function(e) { console.log('submit called'); var formdata = json.stringify($('form').serializeobject()); console.log(formdata); $.ajax({ url: "http://localhost:3000/savedata/", type: "post", datatype: 'json', data: {objectdata: formdata}, contenttype: "application/json", complete: function() { console.log('process complete'); }, success: function(data) { console.log('process sucess'); }, error: function() { console.log('process error'); }, }); return false; }); then on server: var express = require('express'); var app = express(); app.set('views', __dirname + '/views'); app.set('view engine', 'jade'); app.use(express.bodyparser());...
Read more