[R] Binary data on JAGS
Bert Gunter
bgunter@4567 @end|ng |rom gm@||@com
Tue Jun 1 19:39:51 CEST 2021
Unless you have got reason not to, always reply to the list (included in
this response). I cannot help, but someone else may be able to.
Bert Gunter
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Tue, Jun 1, 2021 at 10:34 AM Massimiliano Tripoli <
massimiliano.tripoli using gmail.com> wrote:
> I tried to write
> p[i]<y[i]/n[i] but I received a runtime error
> (Attempt to redefine node p[1]).
> Thanks
>
> Il mar 1 giu 2021, 17:22 Bert Gunter <bgunter.4567 using gmail.com> ha scritto:
>
>> Where is p defined before it is used? (Is this part of what jags provides
>> somehow?)
>>
>> Bert Gunter
>>
>> "The trouble with having an open mind is that people keep coming along
>> and sticking things into it."
>> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>>
>>
>> On Tue, Jun 1, 2021 at 7:57 AM Massimiliano Tripoli <
>> massimiliano.tripoli using gmail.com> wrote:
>>
>>> Dear R users,
>>> I'm trying to reproduce the example 6.5.1 (Dobson (1983)) in BUGS book in
>>> linux using JAGS.
>>> Below the code as
>>>
>>>
>>>
>>> https://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-the-bugs-book/bugs-book-examples/the-bugs-book-examples-chapter-6-6-5-1/
>>>
>>> # By R code:
>>> library('rjags')
>>>
>>> jags <- jags.model(
>>> file = "651.bug",
>>> list (x = c(1.6907, 1.7242, 1.7552, 1.7842, 1.8113, 1.8369, 1.8610,
>>> 1.8839),
>>> n = c(59, 60, 62, 56, 63, 59, 62, 60),
>>> y = c(6, 13, 18, 28, 52, 53, 61, 60)) ,
>>> inits = list(beta = 0, alpha = 50),
>>> n.chains = 3
>>> )
>>>
>>> # By JAGS code: (651.bug)
>>>
>>> model {
>>> for (i in 1:8) {
>>> y[i] ~ dbin(p[i], n[i])
>>> logit(p[i]) <- alpha + beta*(x[i] - mean(x[]))
>>> phat[i] <- y[i]/n[i]
>>> yhat[i] <- n[i]*p[i]
>>> }
>>> alpha ~ dnorm(0, 0.0001)
>>> beta ~ dnorm(0, 0.0001)
>>> }
>>> # I received this error in R:
>>>
>>> # Resolving undeclared variables
>>>
>>> # Allocating nodes
>>>
>>> #Graph information:
>>>
>>> # Observed stochastic nodes: 8
>>>
>>> # Unobserved stochastic nodes: 2
>>>
>>> # Total graph size: 78
>>>
>>>
>>>
>>> #Initializing model
>>>
>>> #Deleting model
>>>
>>>
>>>
>>> #Error in jags.model(file = "651.bug", list(x = c(1.6907, 1.7242, 1.7552,
>>> :
>>> # Error in node y[1]
>>>
>>> #Node inconsistent with parents
>>>
>>> Thanks in advance
>>> M. Tripoli
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________________________
>>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
[[alternative HTML version deleted]]
More information about the R-help
mailing list