Displaying 5 results from an estimated 5 matches for "bugsdata".
Did you mean:
bugdata
2008 Apr 08
2
Error: expected the collection operator c error pos 98 (error on line 1)
Greetings,
I implemented BRugs to run the EWOC model with a cohort size N=1. I output the simulation data using bugsdata(data), where data is the following list.
> data
$Dose
[1] 140
$Y
[1] 0
bugsdata(data) puts out the file data.txt as follows.
list(Dose=1.40000E+02, Y=0.00000E+00)
after I type the following line
modelData("data.txt") # read data file
I get th...
2008 Dec 09
1
expected variable name error pos 98349 WInBUGS in R
...t;
> The following object(s) are masked from anal.data ( position 4 ) :
> NAUSEA SITEID
> > data.bugs <-
> list("NAUSEA"=NAUSEA,"SITEID"=SITEID,"n.samples"=n.samples,"n.sites"=n.sites,"n.params"=n.params)
> > bugsData(data.bugs, fileName = "nauseadata.txt")
> > inits.bugs <- list("alpha"=rep(0,n.sites), "tau"=1)
> > bugsInits(list(inits.bugs), fileName = "nauseainit.txt")
> > modelCheck("nausea_random.txt") # check model file
> model...
2008 Dec 07
1
Reading txt file in R
...h(unique(SITEID))
data.bugs <- list("NAUSEA"=NAUSEA,"SITEID"=SITEID,"BMI"=BMI, "USUBJID"=USUBJID,"n.samples"=n.samples,"n.sites"=n.sites,"n.params"=n.params)
bugsData(data.bugs, fileName = "nauseadata.txt")
inits.bugs <- list("b"=rep(0,n.params),"alpha"=rep(0,n.sites), "tau"=1)
bugsInits(list(inits.bugs), fileName = "nauseainit.txt")
> modelCheck("nausea_random.txt") # check model file
model...
2008 Dec 08
1
Reading txt file in R to run Random Intercept Model
...) :
NAUSEA SITEID
The following object(s) are masked from anal.data ( position 4 ) :
NAUSEA SITEID
> data.bugs <-
list("NAUSEA"=NAUSEA,"SITEID"=SITEID,"n.samples"=n.samples,"n.sites"=n.sites,"n.params"=n.params)
> bugsData(data.bugs, fileName = "nauseadata.txt")
> inits.bugs <- list("alpha"=rep(0,n.sites), "tau"=1)
> bugsInits(list(inits.bugs), fileName = "nauseainit.txt")
> modelCheck("nausea_random.txt") # check model file
model is syntactically co...
2013 Jan 18
0
problem that arises after using the new version of "BRugs"
...p("currently only working for chains longer than 100")
if ((100 * round(nBurnin/100) != nBurnin) | (100 * round(nIter/100) !=
nIter))
warning("chain lengths not multiples of 100; truncation may occur.")
modelCheck(modelFile)
modelData(bugsData(data))
modelCompile(numChains = 1)
modelInits(bugsInits(inits))
burninIncSize <- round(nBurnin/(nThin * 100))
for (i in 1:100) {
modelUpdate(burninIncSize, thin = nThin)
}
samplesSet(parametersToSave)
iterIncSize <- round(nIter/(nThi...