[R] HELP: BRUGS/WinBUGS/RBUGS Response is a combination of random variables
Haoda Fu
fuhds at yahoo.com.cn
Mon Jul 20 05:05:29 CEST 2009
Hi,
Is there anyone know if BUGS language allows the combination of variables as response
such as
Y[i] <- a*X1[i]+b*X2[i]
Y[i] ~ dnorm(c,d)
It seems doesn't work in my model. The problem is between two ######.
The error message is
> modelCheck("BayesBioMarker.BUGS")
model is syntactically correct
> modelData(paste("BUGS_data.txt",sep=""))
data loaded
> modelCompile(numChains=1)
multiple definitions of node bm[1]
===================================
model
{
for(iter in 1:numSubj){
bmd[iter] ~ dnorm(u[trt[iter]+1],sdbmd.precision);
###############
bm[iter] <- inprod(biomarker[iter,1:4], bta[1:4])
bm[iter] ~ dnorm(bmu[iter], sdbm.precision);
###############
bmu[iter] <- (1-delta[trt[iter]+1])*u[trt[iter]+1]+ delta[trt[iter]+1]*(u[1]+u[2])/2;
}
for(iter in 1:4){
bta[iter] ~ dnorm(0, 0.000001);
}
for (iter in 1:2){
delta[iter] ~ dbeta(0.0001,0.0001);
u[iter] ~ dnorm(0,000001);
}
sdbmd.precision ~ dgamma(0.00001,0.00001);
sdbm.precision ~ dgamma(0.00001,0.00001);
}
===================
Data
list(biomarker= structure(.Data= c(-1.6332505600E-01, -9.9820335000E-02, 4.0692602900E-01, -8.2650000000E-03, -6.6739468400E-01, -1.9444291700E-01, -2.5035248500E-01, -1.5926000000E-02, 2.5951119500E-01, 3.1943077100E-01, 2.1915456300E-01, -7.0600000000E-04, -1.1407237740E+00, -4.7155522000E-01, -1.4495155440E+00, 1.7036000000E-02, -1.9105523700E-01, -5.6980210000E-03, -8.4883829700E-01, 4.9390000000E-03), .Dim=c(5, 4)), bmd=c(-4.9440000000E-03, 3.5207000000E-02, 2.5975000000E-02, 2.7394000000E-02, 7.5010000000E-03), trt=c(1.0000000000E+00, 0.0000000000E+00, 0.0000000000E+00, 1.0000000000E+00, 1.0000000000E+00), numSubj=5.0000000000E+00)
More information about the R-help
mailing list