Ändern der Referenz-Gruppe mit binomial-glm mit Familie

Wenn ich eine binomial-regression in R mit einem Faktor unabhängig variable, bestehend aus drei Ebenen der "Höheren" , "Mittleren" und "Unteren", den ich ändern möchten die referenzkategorie mit relevel bekomme ich diesen Fehler:

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

Habe ich geprüft, ob cbsnivcat3 ist ein Faktor 

> 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

aber das gleiche mit factor glm(tertiary ~ factor(cbsnivcat3) , family = binomial, data = data) arbeitet oke.

MQ1: Was verursacht die Fehler, was mache ich falsch?

Vielen Dank im Voraus für die Tipps, - Befehle und-Lösung(s)

str(Daten)

$ 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 ...
  • vermutlich cbsnivcat3 ist ein geordneter Faktor, nicht ein gewöhnlicher. Können wir sehen str(data) bitte?
InformationsquelleAutor Adam | 2013-07-21
  1. 3

    cbsnivcat3 ist ein geordneter Faktor, so relevel wird nicht funktionieren. Was Sie tun können, ist zu ändern, um einen ordentlichen Faktor:

    cbsnivcat3 <- factor(cbsnivcat3, levels=c("Lower","Middle","Higher"))
    • Danke Euch beiden @Ooi Hong und Ben Bolker. Nach dem Tipp von Ben Bolker, dachte ich "ja das ist es".
    InformationsquelleAutor Hong Ooi
  2. 1

    Versuchen, diese

    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))
    InformationsquelleAutor Diogo
Schreibe einen Kommentar