r - 使用带有二项式族的 glm 更改引用组

Question

当我在 R 中运行二项式回归时，其独立因子变量由三个级别“Higher”、“Middle”和“Lower”组成，我想使用 relevel 更改引用类别我得到这个错误:

“Error in relevel.ordered(cbsnivcat3, "Lower") : 
  'relevel' only for factors”

我检查过cbsnivcat3是否是一个因素

> is.factor(data$cbsnivcat3)
[1] TRUE

> levels(data$cbsnivcat3)
[1] "Higher" "Middle" "Lower" 

> t1m4=glm(tertiary ~ relevel(cbsnivcat3, "Lower") , family = binomial, data = data)
    Error in relevel.ordered(cbsnivcat3, "Lower") : 
      'relevel' only for factors

但与 factor glm(tertiary ~ factor(cbsnivcat3) , family = binomial, data = data) 相同。

MQ1:是什么原因导致错误，我做错了什么？

提前感谢您提供的提示、命令和解决方案

海峡(数据)

$ cbsnivcat3     : Ord.factor w/ 3 levels "Higher"<"Middle"<..: 1 1 1 1 1 1 1 1 1 1 ...
$ tertiary       : num  0 0 1 1 1 1 0 1 0 0 ...

Answer 1

cbsnivcat3 是一个有序因子，因此 relevel 将不起作用。您可以做的是将其更改为普通因子:

cbsnivcat3 <- factor(cbsnivcat3, levels=c("Lower","Middle","Higher"))