在我的输出中未显示的 cox 回归分类级别

如何解决在我的输出中未显示的 cox 回归分类级别

在我的数据中，Hemoglobin_group 是一个具有 6 个级别的分类变量。当我运行以下代码时，我在输出中看不到 Hemoglobin_group 级别。我该如何解决这个问题？

fit <- coxph(Surv(Time,Status)~Hemoglobin_group,data)
fit

我的输出：

                      coef   exp(coef) se(coef)  z        p
    Hemoglobin_group -0.06585  0.93627  0.01874 -3.514 0.000441

解决方法

您的变量 Hemoglobin_group 可能被视为一个数值。试试：

Hemoglobin_groupF <- factor(Hemoglobin_group)
fit <- coxph(Surv(Time,Status) ~ Hemoglobin_groupF,data)
fit

参考组将是第一个因素。您可以使用函数 relevel

轻松更改参考系数 ,

str(data)
tibble [2,021 x 21] (S3: tbl_df/tbl/data.frame)
$ status             : num [1:2021] 1 1 1 1 1 1 0 0 0 0 ...
$ id                 : num [1:2021] 1 1 1 1 1 1 2 2 2 2 ...
$ Time               : num [1:2021] 20.3 20.3 20.3 20.3 20.3 ...
$ t1                 : num [1:2021] 0 1 2 3 4 5 0 1 2 3 ...
$ t2                 : num [1:2021] 1 2 3 4 5 ...
$ sex_string         : chr [1:2021] "MALE" "MALE" "MALE" "MALE"...
$ sex                : num [1:2021] 1 1 1 1 1 1 1 1 1 1 ...
$ age                : num [1:2021] 89 89 89 89 89 89 77 77 77 77 ...
$ hemoglobin         : num [1:2021] 9.71 10.22 11.3 11.8 11.2 ...
$ Diabet             : num [1:2021] 1 1 1 1 1 1 2 2 2 2 ...
$ Hemoglobin_group   : num [1:2021] 4 4 4 4 4 4 5 5 5 5 ...
$ Kreatinin          : num [1:2021] 7.19 7.19 7.19 7.19 7.19 ...
$ fosfor             : num [1:2021] 4.14 4.14 4.14 4.14 4.14 ...
$ Kalsiyum           : num [1:2021] 8.5 8.5 8.5 8.5 8.5 ...
$ CRP                : num [1:2021] 1.33 1.33 1.33 1.33 1.33 ...
$ Albumin            : num [1:2021] 4.19 4.19 4.19 4.19 4.19 ...
$ Ferritin           : num [1:2021] 428 428 428 428 428 ...
$ months             : num [1:2021] 1 2 3 4 5 6 1 2 3 4 ...

看起来像我写str（数据）的时候。我认为通过在我的数据中执行以下代码，我正在转变为一个因子。我想我无法转型。没看懂？

我写的转换为factor的代码如下

sex<-as.factor(sex)
is.factor(sex)

Diabet<-as.factor(Diabet)
is.factor(Diabet)


Status<-as.factor(Status)
is.factor(Status)

months<-as.factor(months)
is.factor(months)

Hemoglobin_group<-as.factor(Hemoglobin_group)
is.factor(Hemoglobin_group)
                         

当我运行这段代码时，R 控制台看起来像：

> sex<-as.factor(sex)
> is.factor(sex)
 [1] TRUE
> 
> Diabet<-as.factor(Diabet)
> is.factor(Diabet)
[1] TRUE
> 
> 
> Status<-as.factor(Status)
> is.factor(Status)
[1] TRUE
> 
> months<-as.factor(months)
> is.factor(months)
[1] TRUE
> 
> Hemoglobin_group<-as.factor(Hemoglobin_group)
> is.factor(Hemoglobin_group)
[1] TRUE
                    

在这种情况下，我的数据中的分类变量不会变成因子吗？