3

当我在 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 ...
4

2 回答 2

3

cbsnivcat3是一个有序的因素，所以relevel不会工作。您可以做的是将其更改为普通因子：

cbsnivcat3 <- factor(cbsnivcat3, levels=c("Lower","Middle","Higher"))
于 2013-07-21T13:52:15.317 回答
1

试试这个

df <- data.frame(y=sample(c(0,1), size=100,replace=T),x1=rep(c('a','b','c'),100), x2=1:100)
summary(glm(y~x1+x2,data=df))

df$x1 <- relevel(df$x1, ref='b')
summary(glm(y~x1+x2,data=df))

df$x1 <- relevel(df$x1, ref='c')
summary(glm(y~x1+x2,data=df))
于 2014-12-13T00:03:45.840 回答