On 29.01.2011 15:06, martanair wrote:
Why this?
I write:
radon.data<- list ("n", "J", "x", "y", "county")
radon.inits<- function (){
list (a=rnorm(J), b=rnorm(1), mu.a=rnorm(1),
sigma.y=runif(1), sigma.a=runif(1))
}
radon.parameters<- c ("a", "b", "mu.a", "sigma.y", "sigma.a")
# with 10 iterations
radon.1<- bugs (radon.data, radon.inits, radon.parameters, "radon.1.bug",
n.chains=3, n.iter=10)
1. Let me assume your are talking about R2WinBUGS, right? Please say so
and indicate in your subject line.
2. We do not have radon1.bug. Hence we cannot reproduce anything. Please
read the posting guide.
plot (radon.bugs.1) # to get a plot similar to Figure 16.1
3. What is Figure 16.1?
print (radon.bugs.1) # to display the results in the R console
# with 500 iterations
radon.bugs.1<- bugs (radon.data, radon.inits, radon.parameters,
"radon.1.bug",
n.chains=3, n.iter=500)
We still do not have radon1.bug.
plot (radon.bugs.1) # to get Figure 16.1
print (radon.bugs.1) # to display the results in the R console
But after this
# with 10 iterations
radon.1<- bugs (radon.data, radon.inits, radon.parameters, "radon.1.bug",
n.chains=3, n.iter=10)
R programe give me an error:
Errore in FUN(X[[1L]], ...) : oggetto "n" non trovato
This is rather unlikely unless you have changed some of your objects in
the meantime. In any case, we can only help if you send reproducible code.
Best,
Uwe Ligges
Anybody can help me?
Thanks
Marta
______________________________________________
R-help@r-project.org mailing list
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.
