帮助设计/分析重复测量方案你好，
我正在准备为我的研究方案设计实验，并试图确定一个模型来分析我的数据。以下是我的设计设置：
响应变量：心率
随机因素：主题
固定因素：治疗（对照、治疗A、治疗B）
目标：确定偏离基线心率的统计显着因素/交互作用
协议:
- 治疗前 60 分钟记录基线心率
- 在 20 分钟治疗期间记录 20 分钟心率
- 记录治疗后 60 分钟的心率
我很困惑，因为文献中使用了多种方法，而且网上关于重复测量数据分析的建议似乎相互矛盾。我已经了解了使用广义线性混合模型 (GLMM) 相对于重复测量方差分析的主要好处。我只是对使用我的协议收集的数据的 GLMM 的实际应用有点困惑。
我正在考虑将记录的心率平均到 5 分钟的区间中，这将为每个受试者提供 12 个基线、4 个治疗期间数据点和 12 个治疗后数据点。没有真正的韵律或理由为什么我选择了 5 分钟。大多数文献在分析之前都会按前 5 分钟的区间 (%) 对受试者心率进行标准化，但我厌倦了这样做，因为它不允许对受试者之间的心率水平进行比较。
现在，我有点陷入困境，不知道如何继续或理解我拥有什么样的模型。
希望熟悉此类协议的人可以提供帮助。
任何帮助深表感谢。
谢谢，
J.P
高级统计建模消费者和市场研究回复：帮助设计/分析重复测量协议谢谢你的帮助！
回复：帮助设计/分析重复测量协议不，不是的。RE和重复结构代表两种不同的变化来源。但是，对于某些重复结构（如UN），RE不再可估计 - 如果您尝试这样做，则会收到错误消息。AR（1）具有随机主题效果。
时间顺序对联合国无关紧要。
回复：帮助设计/分析重复测量协议有两种治疗水平（A =低剂量，B =高剂量）。
基于Cholesterol Stacked.jmp示例，我认为我需要设置数据集和分析：
- 平台 ：适合模特
- 个性 ：混合模型
- Y，回应 ：心率（bpm） - >连续
- 固定效果 ：治疗。，时区，Reg。研究所。- >所有分类
在重复结构选项卡下，如果我想使用'非结构化'：
- 学科 ： 学科
- 重复 ：Time Reg./Inst。的毗连。- >分类
如果我想使用'AR（1）'，请在重复结构选项卡下：
| 子。 | 对待。 | 时区 | 注册。研究所。 | Time Reg./Inst。的毗连。 | 时间（分钟） | 心率（bpm） |
| 1 | 一个 | 前处理 | 1 | 预Treatment_1 | 0 | 62 |
| 1 | 一个 | 前处理 | 2 | 预Treatment_2 | 五 | 60 |
| 1 | 一个 | 治疗期间 | 1 | 在治疗期间_1 | 10 | 64 |
| 1 | 一个 | 治疗期间 | 2 | During_Treatment_2 | 15 | 67 |
| 1 | 一个 | 后处理 | 1 | 后Treatment_1 | 20 | 65 |
| 1 | 一个 | 后处理 | 2 | 后Treatment_2 | 25 | 64 |
| 2 | 乙 | 前处理 | 1 | 预Treatment_1 | 0 | 71 |
| 2 | 乙 | 前处理 | 2 | 预Treatment_2 | 五 | 72 |
| 2 | 乙 | 治疗期间 | 1 | 在治疗期间_1 | 10 | 78 |
| 2 | 乙 | 治疗期间 | 2 | 治疗期间_2 | 15 | 76 |
| 2 | 乙 | 后处理 | 1 | 后Treatment_1 | 20 | 76 |
| 2 | 乙 | 后处理 | 2 | 后Treatment_2 | 25 | 75 |
| ... | ... | ... | ... | ... | ... | ... |
“重复结构”主题字段是否构成我的“子”。列随机效果还是需要在“随机效果”选项卡中指定？此外，使用“非结构化”选项，如果“时间（分钟）”列未用于该情况，模型如何知道每个时间点的顺序？
谢谢！
回复：帮助设计/分析重复测量协议多少治疗水平？对于时间组级别，它们应该像pre-1，pre_2，during_1，_2，...等。查看示例数据，胆固醇Stacked.jmp，用于JMP文档的重复测量示例。该数据的格式和布局与您的数据类似。
回复：帮助设计/分析重复测量协议谢谢，我一定会看看文字。
总而言之，我的堆叠数据的简化示例版本如下所示：
| 子。 | 对待。 | 时间组 | 从第一个数据点开始的时间（分钟） | 心率（bpm） |
| 1 | 一个 | 前处理 | 0 | 62 |
| 1 | 一个 | 前处理 | 五 | 60 |
| 1 | 一个 | 治疗期间 | 10 | 64 |
| 1 | 一个 | 治疗期间 | 15 | 67 |
| 1 | 一个 | 后处理 | 20 | 65 |
| 1 | 一个 | 后处理 | 25 | 64 |
| 2 | 乙 | 前处理 | 0 | 71 |
| 2 | 乙 | 前处理 | 五 | 72 |
| 2 | 乙 | 治疗期间 | 10 | 78 |
| 2 | 乙 | 治疗期间 | 15 | 76 |
| 2 | 乙 | 后处理 | 20 | 76 |
| 2 | 乙 | 后处理 | 25 | 75 |
| ... | ... | ... | ... | ... |
回复：帮助设计/分析重复测量协议您需要创建一个时间列（例如，0,10,20,30等）。 与主题ID列一起，两者唯一地标识每个心率测量值。指定重复结构时，需要这两个ID列。（注意：对于AR（1），时间列应定义为连续。对于其他人来说，这是绝对的。请参阅“重复协方差结构要求”部分 http://www.jmp.com/support/help/13/Fit_Model_Launch_Window_2.shtml#999895 。
您还需要创建一个单独的三级治疗柱，以用作固定效果。
最后，除了JMP文档，我强烈推荐这本书 https://www.sas.com/storefront/aux/en/spmixedmodel/59882_toc.pdf 。第5章是关于重复措施的分析。
回复：帮助设计/分析重复测量协议由于30个响应测量中的每一个实际上是每个5分钟窗口的摘要（平均值），我可以将其浓缩为总结为10分钟窗口以减少数据点的数量，我想我将研究Exchange和AR1结构。
作为一个更尖锐的问题，我是否需要为“时间组”创建一个列，其中有三个可能的级别：“预处理”，“治疗期间”和“治疗后”，并将其用作我的混合模特？
我问的原因是因为这些类型的研究通常在施用治疗之前具有零时间响应数据点，然后所有数据点被认为是治疗后的时间。在我的情况下，我的治疗在停止前连续给药30分钟，因此还需要定义“治疗期间”。
谢谢！
回复：帮助设计/分析重复测量协议1。 您可以将所有测量结果转换为“时间为零”的差异（或％变化）。是否转换应部分由您的学习目标指导。
2。重复测量的信息用于估计您将在模型中包含的治疗效果以及重复结构。需要注意的是，由于您已经记录了30个响应点，因此可能会有太多的方差 - 协方差参数来估算UN等结构。 我建议你从Exchangeable或AR（1）开始。
回复：帮助设计/分析重复测量协议你好，健，
谢谢您的回复。在回顾这个讨论之前，我刚刚完成了对重复测量教程的JMP Pro混合模型的观察。他们非常有帮助。
作为后续行动，我有两个问题：
1）最好将响应变量（心率）保持其原始形式以执行此混合模型分析，或者可以将心率（受试者内部随时间变化）标准化为从零时间或差异的百分比变化从运行分析之前的零时间开始？我的直觉告诉我，正常化可能不是一个好主意，因为它不允许进行主体间比较，但我很好奇你的意见。
2）在我的情况下，每个受试者在三个有序时间组内及时测量30个响应数据点，即：输注治疗前（12个数据点），输注治疗期间（6个数据点）和输注后治疗（12个数据点）。如何将此信息添加到混合模型平台？
谢谢！
回复：帮助设计/分析重复测量协议鉴于你的反应变量，心率，是连续的，我会用线性混合模型开始分析，这在JMP Pro中可用（ 分析=>拟合模型并选择混合模型如个性 ）。数据必须采用高格式，即每个观察都由“主题”和“时间”列标识。
在JMP混合模型中，可以使用“随机效果”选项卡指定随机主题效果，而重复测量之间的相关性由“重复结构”选项卡处理。当然，使用“固定效果”选项卡指定固定效果，例如治疗等。
使用JMP Pro 13，有许多不同的方差 - 协方差结构可供选择。在决定“最佳”结构之前，您可以尝试使用几个并评估模型。有关详细信息，请参阅文档 http://www.jmp.com/support/help/13/Fit_Model_Launch_Window_2.shtml#999895
