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
|