a***e
发帖数: 102 | 1 一写contrast就犯晕，请大牛帮忙给看一下对不对， group是1到6，10mg 是1， 以此
类推， EE 是5， PBO是6， 用的是
logistic reg。多谢了！
roc logistic data=adc descending;
class TRTPN;
model RFL=TRTPN /expb firth;
oddsratio TRTPN ;
contrast '10mg vs EE' trtpn 1 0 0 0 -1 / estimate=exp;
contrast '20mg vs EE' trtpn 0 1 0 0 -1 / estimate=exp;
contrast '35mg vs EE' trtpn 0 0 1 0 -1 / estimate=exp;
contrast '55mg vs EE' trtpn 0 0 0 1 -1 / estimate=exp;
contrast 'ALLTRT vs EE' trtpn 1 1 1 1 -1 / estimate=exp;
contrast '10mg vs PBO' trtpn 2 1 1 1 1 / estimate = exp;
contrast '20mg vs PBO' trtpn 1 2 1 1 1 / estimate = exp;
contrast '35mg vs PBO' trtpn 1 1 2 1 1 / estimate = exp;
contrast "55mg vs PBO" trtpn 1 1 1 2 1 / estimate = exp;
contrast 'EE vs PBO' trtpn 1 1 1 1 2 / estimate = exp;
contrast "ALLTRT vs PBO" trtpn 2 2 2 2 1 / estimate = exp;
run; |
t****x
发帖数: 1650 | 2 这两句有问题
contrast 'ALLTRT vs EE' trtpn 1 1 1 1 -1 / estimate=exp;
contrast "ALLTRT vs PBO" trtpn 2 2 2 2 1 / estimate = exp;
需要estimate的话还要加个divider之类，你要是看不清就加个PARAM=GLM吧
【在 a***e 的大作中提到】
: 一写contrast就犯晕，请大牛帮忙给看一下对不对， group是1到6，10mg 是1， 以此
: 类推， EE 是5， PBO是6， 用的是
: logistic reg。多谢了！
: roc logistic data=adc descending;
: class TRTPN;
: model RFL=TRTPN /expb firth;
: oddsratio TRTPN ;
: contrast '10mg vs EE' trtpn 1 0 0 0 -1 / estimate=exp;
: contrast '20mg vs EE' trtpn 0 1 0 0 -1 / estimate=exp;
: contrast '35mg vs EE' trtpn 0 0 1 0 -1 / estimate=exp;
|
a***e
发帖数: 102 | 3 多谢多谢！！！
我也是感觉这两句不是很确定，实际上我觉得如果compare 所有trt的 avg to EE的话
，是不是需要把ee的coe 变成-4呢？ 但是logistic 好像不需要所有的coe sum as 0。
这点和anova又不同.有点糊涂了。。。
contrast 'ALLTRT vs EE' trtpn 1 1 1 1 -4 / estimate=exp; |
t****x
发帖数: 1650 | 4 我的理解是的， -4，但是你要estimate的话得加个 DIVISOR=4
还是用class 里面加param=glm 吧，每个trt分开，简单明了
【在 a***e 的大作中提到】
: 多谢多谢！！！
: 我也是感觉这两句不是很确定，实际上我觉得如果compare 所有trt的 avg to EE的话
: ，是不是需要把ee的coe 变成-4呢？ 但是logistic 好像不需要所有的coe sum as 0。
: 这点和anova又不同.有点糊涂了。。。
: contrast 'ALLTRT vs EE' trtpn 1 1 1 1 -4 / estimate=exp;
|
