R help - Sep 2008 - A Simple Question

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

Gao, Aijun <Aijun.Gao <at> i3statprobe.com> writes:
> 
> 
> 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"
  I think you need 

library(R2WinBUGS)

before you try to run the bugs command

this is all assuming that you have WinBUGS installed ...

  Ben Bolker