model
{
# likelihood
for(subject in 1:N){ # loop in subject
intercept[subject]~dnorm(0,tau.subject);
for(time in 1:T){ # loop in time
mu[subject,time]<-intercept[subject] +
alpha.c[group[subject]] + (beta[group[subject]]*time);
ability[subject,time]~dnorm(mu[subject,time],tau.resid);
} # end of time loop
} # end of subject loop
# priors
var.subject~dunif(0,1.0E4);
var.resid~dunif(0,1.0E4);
beta[1]~dnorm(0,1.0E-4); # Linear effect of time (group=A)
beta[2]~dnorm(0,1.0E-4); # Linear effect of time (group=B)
beta[3]~dnorm(0,1.0E-4); # Linear effect of time (group=C)
alpha.c[1]~dnorm(0,1.0E-4); # Centered intercept (group=A)
alpha.c[2]~dnorm(0,1.0E-4); # Centered intercept (group=B)
alpha.c[3]~dnorm(0,1.0E-4); # Centered intercept (group=C)
# scalars
tau.subject<-1/var.subject;
tau.resid<-1/var.resid;
rho<-var.subject/(var.resid+var.subject); # Within-sub corr.
b.diff[1]<-beta[2]-beta[1];
b.diff[2]<-beta[3]-beta[1];
alpha[1]<-alpha.c[1]+50; # re-scale intercepts
alpha[2]<-alpha.c[2]+50;
alpha[3]<-alpha.c[3]+50;
a.diff[1]<-alpha[2]-alpha[1];
a.diff[2]<-alpha[3]-alpha[1];
}
这是model代码!不知道怎么载入数据



雷达卡



京公网安备 11010802022788号







