rats_model.txt

matlab中实现openbugs或winbugs的功能调用

model
{
	for( i in 1 : N ) {
		for( j in 1 : T ) {
			Y[i , j] ~ dnorm(mu[i , j],tau.c)
			mu[i , j] <- alpha[i] + beta[i] * (x[j] - xbar)
		}
		alpha[i] ~ dnorm(alpha.c,alpha.tau)
		beta[i] ~ dnorm(beta.c,beta.tau)
	}
	tau.c ~ dgamma(0.001,0.001)
	sigma <- 1 / sqrt(tau.c)
	alpha.c ~ dnorm(0.0,1.0E-6)	   
	alpha.tau ~ dgamma(0.001,0.001)
	beta.c ~ dnorm(0.0,1.0E-6)
	beta.tau ~ dgamma(0.001,0.001)
	alpha0 <- alpha.c - xbar * beta.c	
}