s*****e
发帖数: 157 | 1 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
time2测量了一次.所以一个id可能有多个repeated measure.
我想建个mixed model比较tr对response的影响
proc mixed;
class id tr time;
model response=tr|time;
repeated /subject=id;
run;
程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
呢？请大家支招，多谢！
dataset:
id week tr time response
01 1 A 1 xx
01 1 A 2 xx
01 2 A 1 xx
01 2 A 2 xx
01 3 B 1 xx
01 3 B 2 xx
02 1 A 1 xx
02 1 A 2 xx
02 3 B 1 xx
02 3 B 2 xx
03 3 B 1 xx
03 3 B 2 xx
04 3 B 1 xx
04 3 B 2 xx
04 4 A 1 xx
04 4 A 2 xx |
j*******y
发帖数: 58 | 2 这好像应该用random id instead of repeated / subjiect=id
【在 s*****e 的大作中提到】
: 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
: 与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
: time2测量了一次.所以一个id可能有多个repeated measure.
: 我想建个mixed model比较tr对response的影响
: proc mixed;
: class id tr time;
: model response=tr|time;
: repeated /subject=id;
: run;
: 程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
|
k****i
发帖数: 347 | 3 程序无法运行的error message是什么？
你更care time而不是week的effect?
repeated statement没有给structure |
j*******y
发帖数: 58 | 4 ft这很明显每个id每个time都有multiple points，不出错才怪。除非specify
repeated time*week/subj=id,不过也有可能有些id在某个time*week有multiple
points。
【在 k****i 的大作中提到】
: 程序无法运行的error message是什么？
: 你更care time而不是week的effect?
: repeated statement没有给structure
|
s*****e
发帖数: 157 | 5 time1是pre, time2是post,所以要比较。week是随机对应两种tr的，我不知道还要不要
考虑？
structrure随便给了一个。error具体不记得了，但似乎是由于多个repeated measure
导致的
【在 k****i 的大作中提到】
: 程序无法运行的error message是什么？
: 你更care time而不是week的effect?
: repeated statement没有给structure
|
s*****e
发帖数: 157 | 6 谢谢指点，一个id一个week只测两次，这样应该不会有multiple points?
这种情况可以用genmod吗？
【在 j*******y 的大作中提到】
: ft这很明显每个id每个time都有multiple points，不出错才怪。除非specify
: repeated time*week/subj=id,不过也有可能有些id在某个time*week有multiple
: points。
|
j*******y
发帖数: 58 | 7 2=multiple
【在 s*****e 的大作中提到】
: 谢谢指点，一个id一个week只测两次，这样应该不会有multiple points?
: 这种情况可以用genmod吗？
|
C*******1
发帖数: 2134 | 8 id 和 time 不要放在class里面。
【在 s*****e 的大作中提到】
: 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
: 与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
: time2测量了一次.所以一个id可能有多个repeated measure.
: 我想建个mixed model比较tr对response的影响
: proc mixed;
: class id tr time;
: model response=tr|time;
: repeated /subject=id;
: run;
: 程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
|
s*****e
发帖数: 157 | 9 time1 和 time2 各一次，这个不是导致错误的重复阿
【在 j*******y 的大作中提到】
: 2=multiple
|
C*******1
发帖数: 2134 | 10 你这个时间太乱了，建议生成一个新的列，里面的时间是distinct的。
week*time
【在 s*****e 的大作中提到】
: 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
: 与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
: time2测量了一次.所以一个id可能有多个repeated measure.
: 我想建个mixed model比较tr对response的影响
: proc mixed;
: class id tr time;
: model response=tr|time;
: repeated /subject=id;
: run;
: 程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
|
|
|
j*******y
发帖数: 58 | 11 you have the hell 2 weeks at time1 and 2 weeks at time2. look at the hell
subject 4.
【在 s*****e 的大作中提到】
: time1 和 time2 各一次，这个不是导致错误的重复阿
|
s*****e
发帖数: 157 | 12 但是我要比较time1(pre)和time2(post)response的不同，如果time和week混起来，怎
么比较呢？
【在 C*******1 的大作中提到】
: 你这个时间太乱了，建议生成一个新的列，里面的时间是distinct的。
: week*time
|
C*******1
发帖数: 2134 | 13 我试了一下你的程序，没有问题，可以运行。
结果是否正确就不知道了。
你是不是data 都没有读进去啊？
【在 s*****e 的大作中提到】
: 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
: 与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
: time2测量了一次.所以一个id可能有多个repeated measure.
: 我想建个mixed model比较tr对response的影响
: proc mixed;
: class id tr time;
: model response=tr|time;
: repeated /subject=id;
: run;
: 程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
|
j*******y
发帖数: 58 | 14 直接看anova output不就行了吗，ft
【在 s*****e 的大作中提到】
: 但是我要比较time1(pre)和time2(post)response的不同，如果time和week混起来，怎
: 么比较呢？
|
k****i
发帖数: 347 | 15 握手。我早就说了，把error message贴出来就知道了
【在 C*******1 的大作中提到】
: 我试了一下你的程序，没有问题，可以运行。
: 结果是否正确就不知道了。
: 你是不是data 都没有读进去啊？
|
s*****e
发帖数: 157 | 16 Note:
An infinite likelihood is assumed in iteration 0 because of a nonpositive
definite
estimated R matrix for ID ***
Thanks!
【在 k****i 的大作中提到】
: 握手。我早就说了，把error message贴出来就知道了
|
k****i
发帖数: 347 | 17 你的model里没有week，会造成同一个id同样的time下有多个点
这样并没有关系
然而如果同一个id同样的time里面的多个点都是同一个数值，就是identical rows，有
可能因此造成infinite likelihood
这个error里已经给出一个有问题的id了 |
j*******y
发帖数: 58 | 18 are you sure? repeated 是specify一个covariance matrix，跟response是不是一样
没有关系吧？我觉得可能是因为这个data是unbalanced，每个id观察值个数不一样，所
以sas连这个covariance matrix的维数都确定不了。
【在 k****i 的大作中提到】
: 你的model里没有week，会造成同一个id同样的time下有多个点
: 这样并没有关系
: 然而如果同一个id同样的time里面的多个点都是同一个数值，就是identical rows，有
: 可能因此造成infinite likelihood
: 这个error里已经给出一个有问题的id了
|
k****i
发帖数: 347 | 19 不太确定，猜的。呵呵
【在 j*******y 的大作中提到】
: are you sure? repeated 是specify一个covariance matrix，跟response是不是一样
: 没有关系吧？我觉得可能是因为这个data是unbalanced，每个id观察值个数不一样，所
: 以sas连这个covariance matrix的维数都确定不了。
|
s*****e
发帖数: 157 | 20 The responses per id*time are different.
【在 k****i 的大作中提到】
: 你的model里没有week，会造成同一个id同样的time下有多个点
: 这样并没有关系
: 然而如果同一个id同样的time里面的多个点都是同一个数值，就是identical rows，有
: 可能因此造成infinite likelihood
: 这个error里已经给出一个有问题的id了
|
|
|
P****D
发帖数: 11146 | 21 一个id一个时间下有多个点，我以前干过这事，虽然是无意的，是数据错误。运行原则
上无问题，但老是convergence failure。
建议把week加入model中，如果一周内有多个treatment，那就用treatment开始的那天
calendar date。 |
s*****e
发帖数: 157 | 22 One week corresponds to a treatment. If I useed time*week in repeated
statement, there was WARNING: Unable to make hessian positive definite.
【在 P****D 的大作中提到】
: 一个id一个时间下有多个点，我以前干过这事，虽然是无意的，是数据错误。运行原则
: 上无问题，但老是convergence failure。
: 建议把week加入model中，如果一周内有多个treatment，那就用treatment开始的那天
: calendar date。
|
P****D
发帖数: 11146 | 23 just model response=week tr time tr*time; I cannot see any scientific reason
to put time*week in the model.
【在 s*****e 的大作中提到】
: One week corresponds to a treatment. If I useed time*week in repeated
: statement, there was WARNING: Unable to make hessian positive definite.
|
p*******r
发帖数: 1951 | 24 首先要明确你要做什么比较？看你说的意思你的null hypothesis是不是
H0: treatment effect = time2 - time1 = 0 ?
如果是，那你的dependent variable 完全可以用 effect = time2 - time1。而不是单
独的time1, time2。
其次，你的数据显示不同的subject有完全不同的时间点， 看你的SAS程序里面也根本
没有用到时间变量week。所以你想做的不是longitudinal data analysis？你只是想比
较不同treatment, pre-, post- response的差别。这样的话用repeated measure
ANOVA就成了。
很久没有用过mixed model了。说的不对的地方大家指教。
【在 s*****e 的大作中提到】
: 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
: 与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
: time2测量了一次.所以一个id可能有多个repeated measure.
: 我想建个mixed model比较tr对response的影响
: proc mixed;
: class id tr time;
: model response=tr|time;
: repeated /subject=id;
: run;
: 程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
|
j*******y
发帖数: 58 | 25 what is the difference between
repeated / subj=id
and
random id?
I imagined they are the same and tested in SAS and they are different, but I
don't know why they are different.
【在 p*******r 的大作中提到】
: 首先要明确你要做什么比较？看你说的意思你的null hypothesis是不是
: H0: treatment effect = time2 - time1 = 0 ?
: 如果是，那你的dependent variable 完全可以用 effect = time2 - time1。而不是单
: 独的time1, time2。
: 其次，你的数据显示不同的subject有完全不同的时间点， 看你的SAS程序里面也根本
: 没有用到时间变量week。所以你想做的不是longitudinal data analysis？你只是想比
: 较不同treatment, pre-, post- response的差别。这样的话用repeated measure
: ANOVA就成了。
: 很久没有用过mixed model了。说的不对的地方大家指教。
|
