（嵌套？）重复测量 ANOVA 的 BUGS 模型

Question

我想知道是否有人有用于重复测量方差分析的 BUGS/JAGS 模型的代码？基本上，我有一个响应 (y)，我想针对 一天中的时间 、天进行建模，以及治疗。我还想包括两个交互项，治疗 x 一天中的时间 和 治疗 x 天。研究中大约有 20 个人，他们在大约 1 周内每天接受 4 次测量。我不完全确定从哪里开始，我担心时间协变量也应该嵌套在日协变量中？如果有人有 BUGS/JAGS 模型的似然部分的代码，将不胜感激。我可以照顾先验。只是似乎无法摆脱这个。

Answer 1

你的问题有些含糊不清。

您希望一天中的时间和天作为连续协变量还是作为离散因子输入？
您希望个体身份作为固定效应还是随机效应进入模型？
如果天或 一天中的时间 是一个因素，您是否要将其作为固定因素包括在内还是随机效应？
您询问 Time of Day 是否应该嵌套在 Day 中。如果不进一步了解您的数据和目标，就无法回答这个问题。

这是一个代码示例，它假定您要将个体视为随机效应。

还假设：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)
}

（嵌套？）重复测量 ANOVA 的 BUGS 模型

BUGS model for (nested?) repeated measures ANOVA

bayesian

anova

winbugs

jags