Question

我试图将我的半参数比例风险模型和加速失效时间模型与指数基线进行比较。系数应该“相同”，但标志不同。

我在R中使用coxph和survreg函数.Coxph给我结果没有拦截，surverg给我带截距的结果。据我所知，我无法打开对coxph基线危害的拦截。有什么方法可以关闭幸存的拦截？我在下面附上了我的结果。

另外，当我的目的是了解治疗效果时，对协变量进行拦截是否很重要？

> coxph(Surv(Time,Reason=="I")~AgeG, data = PTARdata)
Call:
coxph(formula = Surv(Time, Reason == "I") ~ AgeG, data = PTARdata)


           coef exp(coef) se(coef)      z       p
AgeG0   -3.4241    0.0326    1.005 -3.406 6.6e-04
AgeG1    0.1232    1.1311    0.154  0.803 4.2e-01
AgeG2   -0.2774    0.7578    0.137 -2.023 4.3e-02
AgeG3   -0.0503    0.9509    0.146 -0.345 7.3e-01
AgeG4    0.2063    1.2292    0.203  1.014 3.1e-01
AgeG5   -0.1531    0.8581    0.147 -1.039 3.0e-01
AgeG6   -0.7706    0.4627    0.182 -4.238 2.3e-05

Likelihood ratio test=70.5  on 7 df, p=1.15e-12  

> PTAR.ph.age2<-survreg(Surv(Time,Reason=="I")~AgeG,data = PTARdata, dist='weibull',scale=1)
> summary(PTAR.ph.age2)

Call:
survreg(formula = Surv(Time, Reason == "I") ~ AgeG, data = PTARdata, 
    dist = "weibull", scale = 1)
              Value Std. Error       z       p
(Intercept)  8.9611     0.0891 100.588 0.00000
AgeG0        3.1638     1.0040   3.151 0.00163
AgeG1       -0.1332     0.1535  -0.868 0.38537
AgeG2        0.2740     0.1371   1.998 0.04571
AgeG3        0.0423     0.1458   0.290 0.77176
AgeG4       -0.2398     0.2031  -1.181 0.23778
AgeG5        0.1042     0.1465   0.712 0.47673
AgeG6        0.5469     0.1751   3.123 0.00179

Scale fixed at 1 

Weibull distribution
Loglik(model)= -5087.9   Loglik(intercept only)= -5115.2
        Chisq= 54.51 on 7 degrees of freedom, p= 1.9e-09 
Number of Newton-Raphson Iterations: 11

在幸存和coxph中拦截

0 个答案: