我有一个我正在研究的重复测量数据集。数据如下所示：
ID=c('X1', 'X1', 'X1', 'X1', 'X2', 'X2', 'X2', 'X3', 'X3', 'X3', 'X3', 'X4', 'X4', 'X4', 'X4', 'X5', 'X5', 'X5', 'X6', 'X6', 'X6', 'X6')
Diag=c('Con', 'Con', 'Con', 'Con', 'Con', 'Con', 'Con', 'AD', 'AD', 'AD', 'AD', 'AD', 'AD', 'AD', 'AD', 'FD', 'FD', 'FD', 'FD', 'FD','FD', 'FD')
Score=c(10, 9, 8, 8, 10, 9, 9, 5, 4, 3, 3, 5, 4, 3, 2, 5, 4, 3, 6, 5, 4, 3)
Time=c(1, 2, 3, 4, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 1, 2, 3, 4)
dat <- data.frame(ID, Time, Diag, Score)
其中ID =参与者ID，诊断=诊断，时间=重复评估，得分=得分
我想计算AD和FD组分数的z分数，与＆＃39; Con＆＃39;相关。小组得分（对照组）和时间＆＃39;。例如，时间2的AD和FD组的z分数应该与时间2的对照组相关（忽略缺失的变量）
我能想到的唯一方法是将数据框分成不同的“时间”，计算z分数并将数据框合并在一起，但这非常痛苦且耗时（并且正在给予我有时候尝试合并数据框时的NAs。
有更好的方法吗？