Hi,
I suggest to add the right loop the of the node "for(t in 2:(n-1)){"
at the end of the line alpha[t]<-exp(eta[t]).
try to
Quoting Bob O'Hara <[log in to unmask]>:
> On 01/08/16 11:41, Lecomte Jean-Baptiste wrote:
>> Hi,
>>
>> I think you have to give an initial value for every stochastic
>> nodes. For example, you have to give initial values for epsilon,
>> delta, gamma ....
>>
> No you don't: any stochastic node that doesn't have an initial value
> given will have one generated from its prior.
>
> I can't even compile the code (multiple definitions of eta, because
> it's indexed by t but in loops for both i and t), but I think the
> problem is mu[1]. It isn't defined in the model, but still has to be
> created by BUGS (because the indexing of mu starts at 2). In the
> inits it is given a value, but BUGS doesn't know what the node is,
> so can't give it a value. The solution is to either re-write the
> model or replace the first .1 by NA (so BUGS doesn't try to give
> mu[1] a value).
>
> Bob
>
>> Jean-Baptiste
>>
>> 2016-07-28 10:57 GMT+02:00 Abongo <[log in to unmask]
>> <mailto:[log in to unmask]>>:
>>
>> Hello, List!
>> This message came up when I tried to load initial values. This
>> seems to be an uncommon error since my search of Internet did not
>> get any hits. Any comments and suggests would be very much
>> appreciated. Thank you!
>>
>> model{
>> c<-10000
>> for(t in 2:(n-1)){
>> zero[t]<-0
>> zero[t]~dpois(mean.zero[t])
>> mean.zero[t]<- -l[t]+c
>>
>> l[t]<-0.5*(log(lambda.error)-log(2*22/7)-3*log(epsilon[t]))-0.5*lambda.error*pow((epsilon[t]-mu[t])/mu[t],2)/epsilon[t]
>> }
>> for(t in 2:(n-1)){
>> y[t]~dpois(lambda[t])
>> lambda[t]<-alpha[t]+lambda_1*y[t-1]+epsilon[t]
>> alpha[t]<-exp(eta[t])
>>
>> for(i in 1:s){
>>
>> eta[t]<-alpha_1+beta*t+(gamma[i]*sin(omega[i]*t)+delta[i]*cos(omega[i]*t))
>> omega[i]<-i*3.14/26
>> }
>> }
>> for(t in 2:(n-1)){mu[t]~dnorm(0.0,0.001) epsilon[t]~dgamma(0.01,0.01)}
>> for(i in 1:s){gamma[i]~dnorm(0.0,0.0001) delta[i]~dnorm(0.0,0.0001)}
>> lambda.error~dgamma(0.01,0.01)
>> dispersion.p<-1/lambda.error
>> lambda_1~dgamma(0.01,0.001)
>> alpha_1~dnorm(0.0,0.001)
>> beta~dnorm(0.0,0.0001)
>> }
>> list(y=c(1,0,5,2,1,2,0,4,0,0,0,3,1,1,0,0,2,2,0,2,6,3,2,1,2),s=10,n=25)
>>
>> list(mu=c(.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1),lambda.error=1,lambda_1=1,alpha_1=1,beta=1)
>>
>> -------------------------------------------------------------------
>> This list is for discussion of modelling issues and the BUGS software.
>> For help with crashes and error messages, first mail
>> [log in to unmask] <mailto:[log in to unmask]>
>> To mail the BUGS list, mail to [log in to unmask]
>> <mailto:[log in to unmask]>
>> Before mailing, please check the archive at
>> www.jiscmail.ac.uk/lists/bugs.html
>> <http://www.jiscmail.ac.uk/lists/bugs.html>
>> Please do not mail attachments to the list.
>> To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
>> <mailto:[log in to unmask]>
>> If this fails, mail [log in to unmask]
>> <mailto:[log in to unmask]>, NOT the whole list
>>
>>
>>
>>
>> --
>>
>> Dr. Jean-Baptiste Lecomte
>>
>> Postdoctoral Fellow, Groundfish Stock Assessment
>>
>> Pacific Biological Station
>>
>> Fisheries and Oceans Canada
>>
>> 3190 Hammond Bay Road
>>
>> Nanaimo, BC V9T 6N7
>>
>> PH: 250-756-7134
>>
>> -------------------------------------------------------------------
>> This list is for discussion of modelling issues and the BUGS
>> software. For help with crashes and error messages, first mail
>> [log in to unmask] <mailto:[log in to unmask]> To mail the
>> BUGS list, mail to [log in to unmask] <mailto:[log in to unmask]>
>> Before mailing, please check the archive at
>> www.jiscmail.ac.uk/lists/bugs.html
>> <http://www.jiscmail.ac.uk/lists/bugs.html> Please do not mail
>> attachments to the list. To leave the BUGS list, send LEAVE BUGS to
>> [log in to unmask] <mailto:[log in to unmask]> If this
>> fails, mail [log in to unmask]
>> <mailto:[log in to unmask]>, NOT the whole list
>
>
> --
>
> Bob O'Hara
>
> Biodiversity and Climate Research Centre
> Senckenberganlage 25
> D-60325 Frankfurt am Main,
> Germany
>
> Tel: +49 69 7542 1863
> Mobile: +49 1515 888 5440
> WWW: http://www.bik-f.de/root/index.php?page_id=219
> Blog: http://blogs.nature.com/boboh
> Journal of Negative Results - EEB: www.jnr-eeb.org
>
>
> -------------------------------------------------------------------
> This list is for discussion of modelling issues and the BUGS software.
> For help with crashes and error messages, first mail [log in to unmask]
> To mail the BUGS list, mail to [log in to unmask]
> Before mailing, please check the archive at
> www.jiscmail.ac.uk/lists/bugs.html
> Please do not mail attachments to the list.
> To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
> If this fails, mail [log in to unmask], NOT the whole list
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|