Well spotted, the b is a muck up (what happens when you are basing stuff on
someone else's code). The bit you though might be a dimension mismatch seems
to work ok, but the bit that I was worried about from the start doesn't; the
reason I have sqrt of var y/ var x is because my posterior is for a
The comma *before* acc=1 ?
I also wondered whether (further up) this should work:
s2y[i,]=s2y[i-1] # would think this to result in a dimension mismatch
This looks sketchy as well:
s2y[i,j] = s2y[b[i-1,j] + rnorm(1,mean=0, sd=s2yscale[j])
^ ^ ^ # unmatched sqr-brackets
It w
2 matches
Mail list logo
