Ausgabe des Logistikmodells in R.

Ich versuche, die folgende Art von Logistikmodell zu interpretieren:

mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial)

Ist die Ausgabe predict(mdl)der erwarteten Erfolgsaussichten für jeden Datenpunkt? Gibt es eine einfache Möglichkeit, die Quoten für jede Faktorstufe des Modells und nicht für alle Datenpunkte zu tabellieren?

Die Hilfeseiten für

predict.glm

state: "Somit sind für ein Standard-Binomialmodell die Standardvorhersagen logarithmische Quoten (Wahrscheinlichkeiten auf der Logit-Skala) und 'type =" response "' gibt die vorhergesagten Wahrscheinlichkeiten an". So predict(mdl)gibt das Protokoll (Odds) und mit „type =‚Antwort‘liefert die vorhergesagten Wahrscheinlichkeiten Sie könnten dieses Spielzeug Beispiel lehrreich finden.:

> y <- c(0,0,0,1,1,1,1,1,1,1)
> prop.table(table(y))
y
  0   1 
0.3 0.7 
> glm.y <- glm(y~1, family = "binomial")
> ## predicted log(odds)
> predict(glm.y)
        1         2         3         4         5         6         7         8 
0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 
        9        10 
0.8472979 0.8472979 
> ## predicted probabilities (p = odds/(1+odds))
> exp(predict(glm.y))/(1+exp(predict(glm.y)))
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 
> predict(glm.y, type = "response")
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 

In Bezug auf Ihre zweite Frage sollten Sie sich das Effektpaket http://socserv.socsci.mcmaster.ca/jfox/Misc/effects/index.html von John Fox ansehen. siehe auch seinen JSS-Artikel "Effektanzeigen in R für verallgemeinerte lineare Modelle" (S. 8-10).