我想知道是否有人有重复测量ANOVA的BUGS / JAGS模型代码？基本上，我有一个回复（ y ），我想根据时间，日和治疗。我还想包括两个互动术语，治疗x时间和治疗x天。该研究中约有20个人，在约1周内每天测量4次。我不完全确定从哪里开始，我担心时间协变量也应该嵌入日间协变量？如果有人拥有BUGS / JAGS模型的似然部分的代码，那将非常感激。我可以照顾先辈。似乎无法用这个开始。
答案 0 :(得分：0)
你的问题有些含糊不清。
您希望时间和日作为连续协变量或离散因子输入？
您是否希望个人身份作为固定或随机效果进入模型？
如果日或时间是一个因素，是否要将其包含为固定或随机效果？
您询问日期是否应嵌套在日内。如果不了解更多有关您的数据和目标的话，就无法回答这个问题。
以下是代码的示例，假设您希望将个人视为随机效果。
还假设：Treatment，Time.of.day和Day在所有个体中都有不变的斜率。将该模型扩展到固定或随机斜率模型是直截了当的，其中不同的个体获得单独的模型化斜率。例如，对于随机斜率模型，您只需修改下面的beta参数，以类似于alpha参数的方式处理它们。
根据OP的要求，这只是可能性部分，不包括先验部分。
for(i in 1:n.observations){
y[i] ~ dnorm(alpha[individual[[i]] + beta1*Day[i] + beta2*Time.of.day[i] + beta3*Treatment[i] + beta4*Treatment[i]*Day[i] + beta5*Treatment[i]*Time.of.day[i], tau.obs)
}
# individual[i] contains the numerical index representing the individual that corresponds to observation i.
for(j in 1:n.individuals){
alpha[j] ~ dnorm(mu, tau)
}