Print

---

 
Hi List,

I have a large data set with variables year, area, item, value (response
variable). I would like to identify where most of the variance in the
data comes from. This is a good chance to get familiar with create
doodles, but I am not succeeding. For ease of demonstration I am using
only year and value in the following example:

model
{
 for (i in 1:N) {
		mu[i]<-a[year[i]] +  error
    value[i]~dnorm(mu[i], tau.within)
 }

                          
 # temporal variation
 for(i in 1:Nyear){
 	a[i]~dnorm(0.0,tau.year)
 }
 a[Nyear+1]<-sum(a[1:Nyear])

#precision
 tau.within <- 1/v.within
 tau.year <- 1/v.year

 v.year <- sd.year*sd.year
 v.within <- sd.within*sd.within

#priors
 error~dnorm(0,0.00001)
 
 # put random uniform distribution on the SDs
 sd.year~dunif(0,1)
 sd.within~dunif(0,1)
}

I am not having much luck recreating above model with a doodle. What I
am doing wrong? I have tried to attach the doodle, alas the listserver
rejects attachments. So if you are able to help, please email and I will
send you the doodle trial.
 
Any help appreciated
Herry

Dr Alexander Herr - Herry
CSIRO, Sustainable Ecosystems
Gungahlin Homestead
Bellenden Street
GPO Box 284
Crace, ACT 2601
 
Phone/www
(02) 6242 1542; 6242 1705(fax)
0408679811 (mob)

home: www.csiro.au/people/Alexander.Herr
Webadmin ABS: http://ausbats.org.au
Sustainable Ecosystems: www.cse.csiro.au
--------------------------------------------

-------------------------------------------------------------------
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