3

ワイブル比例ハザード モデルを R のデータに当てはめようとしているときに、データの順序を変更すると出力が変わるgehanことに気付きました。gehan

library(MASS)
library(survival)
attach(gehan)

survreg(Surv(time, cens) ~ treat, gehan, dist = "weibull")

以下の出力が生成されます。

Call:
survreg(formula = Surv(time, cens) ~ treat, data = gehan, dist = "weibull")

Coefficients:
 (Intercept) treatcontrol 
 3.515687    -1.267335 

Scale= 0.7321944 

Loglik(model)= -106.6   Loglik(intercept only)= -116.4
Chisq= 19.65 on 1 degrees of freedom, p= 9.3e-06 
n= 42

ただし、Rに手動でデータを入力すると、

control <- c(1,1,2,2,3,4,4,5,5,8,8,8,8,11,11,12,12,15,17,22,23)
control.status <- c(rep(1,21))
treatment <- c(6,6,6,6,7,9,10,10,11,13,16,17,19,20,22,23,25,32,32,34,35)
treatment.status <- c(1,1,1,0,1,0,1,0,0,1,1,0,0,0,1,1,0,0,0,0,0)

new.time <- c(control,treatment)
new.status <- c(control.status,treatment.status)
group <- as.factor(c(rep("control",21),rep("treatment",21)))
parametrik.surv <- Surv(new.time, new.status)
parametrik.fit <- survreg( parametrik.surv  ~ group, dist="w")

以下は、わずかに異なる出力です。

 Call:
survreg(formula = parametrik.surv ~ group, dist = "w")

Coefficients:
(Intercept) grouptreatment 
  2.248352       1.267335 

Scale= 0.7321944 

Loglik(model)= -106.6   Loglik(intercept only)= -116.4
Chisq= 19.65 on 1 degrees of freedom, p= 9.3e-06 
n= 42

さらに、ある本では beta_0、beta_1 および scale はそれぞれ 0.98 、1.27 および 1.37 として指定されています。これの原因は何ですか？

4

1 に答える 1

2

処理の名前を変更したので、一方のモデルでは参照グループとして「コントロール」を使用し、もう一方のモデルでは参照グループとして「6-MP」を使用しています。

于 2015-04-06T16:21:03.817 に答える