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... 阅读全帖 |
|
E****2
发帖数: 247 | 2 请教大家一个问题:
At first, I used,
proc phreg data = pfsdays;
class trt;
model ept*eptcen(1) = trt / rl;
assess var = (trt) ph;
run;
I got : ERROR: Cannot assess the functional form of ATRSORT since it is not
a continuous explanatory variable.
Then I created a numeric var, trtpn, 0, 1 correspoding to two treatment
groups. No error. And I got
two plots for the assess statemetns, Cumulative Martingale function form,
and Standardized Score Process Plot
ods graphics on;
proc ... 阅读全帖 |
|
t****x
发帖数: 1650 | 3 这两句有问题
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
发帖数: 102 | 4 多谢多谢！！！
我也是感觉这两句不是很确定，实际上我觉得如果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; |
|
