Próbuję zinterpretować następujący typ modelu logistycznego:
mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial)
Czy wynik predict(mdl)
oczekiwanych szans powodzenia dla każdego punktu danych? Czy istnieje prosty sposób na zestawienie szans dla każdego poziomu czynnika w modelu, a nie dla wszystkich punktów danych?
fac1
ifac2
.Design
paczkę od Francka Harrella; ma bardzo fajne funkcjelrm()
dla GLM i pokrewnych rzeczy.Odpowiedzi:
Strony pomocy dla
state: „Zatem dla domyślnego modelu dwumianowego domyślnymi przewidywaniami są iloraz szans (prawdopodobieństwa w skali logit), a„ type = „response” „daje przewidywane prawdopodobieństwa”. Tak,
predict(mdl)
zwraca log (OR) i za pomocą „type =«odpowiedź»zwraca przewidywanych prawdopodobieństw Można znaleźć ten przykład zabawki pouczające.:Jeśli chodzi o twoje drugie pytanie, możesz sprawdzić pakiet efektów http://socserv.socsci.mcmaster.ca/jfox/Misc/effects/index.html autorstwa Johna Foxa; patrz także jego artykuł JSS „Wyświetlacze efektów w R dla uogólnionych modeli liniowych” (str. 8-10).
źródło