
Hi:

I am trying to run your SCHOOLS example as following and I get an error message stated that : Error: could not find function "bugs"
Would you please help me since I would like to run some WINBUG programs from your R by reading external data and I am new to R?
Do I missing some function?  Should I install some extra files?


Thanks a lot in advance for your help!!


>  schools <- read.table ("schools.dat", header=TRUE)
>    J <- nrow(schools)
>    y <- schools$estimate
>    sigma.y <- schools$sd
>    data <- list ("J", "y", "sigma.y")
>    inits <- function() {list (theta=rnorm(J,0,100), mu.theta=rnorm(1,0,100), sigma.theta=runif(1,0,100))}
>    parameters <- c("theta", "mu.theta", "sigma.theta")
>    schools.sim <- bugs (data, inits, parameters, "schools.bug", n.chains=3, n.iter=1000)
Error: could not find function "bugs"


********************************************************************************************

Running WinBugs and OpenBugs from R
Example
*

8 schools analysis from Section 5.5 of "Bayesian Data Analysis"
*       Put the following<schools.bug> in the file "schools.bug" in your R working directory (it is convenient to edit this and other .bug and .R files in XEmacs<http://www.xemacs.org>):
model {
for (j in 1:J){
y[j] ~ dnorm (theta[j], tau.y[j])
theta[j] ~ dnorm (mu.theta, tau.theta)
tau.y[j] <- pow(sigma.y[j], -2)
}
mu.theta ~ dnorm (0.0, 1.0E-6)
tau.theta <- pow(sigma.theta, -2)
sigma.theta ~ dunif (0, 1000)
}
*       Save this file<schools.dat> as "schools.dat" in your R working directory.
*       Do the following<schools.R> commands in R:
schools <- read.table ("schools.dat", header=TRUE)
J <- nrow(schools)
y <- schools$estimate
sigma.y <- schools$sd
data <- list ("J", "y", "sigma.y")
inits <- function() {list (theta=rnorm(J,0,100), mu.theta=rnorm(1,0,100), sigma.theta=runif(1,0,100))}
parameters <- c("theta", "mu.theta", "sigma.theta")
schools.sim <- bugs (data, inits, parameters, "schools.bug", n.chains=3, n.iter=1000)
******************************************************************************


Aijun Gao
Senior Biostatistician



"Secure Server BSK" made the following
 annotations on 09/17/2008 01:39:59 PM
------------------------------"This e-mail, including attachments, may include confidential and/or proprietary information, and may be used only by the person or entity to which it is addressed. If the reader of this e-mail is not the intended recipient or his or her authorized agent, the reader is hereby notified that any dissemination, distribution or copying of this e-mail is prohibited. If you have received this e-mail in error, please notify the sender by replying to this message and delete this e-mail immediately."
==============================
	[[alternative HTML version deleted]]

