ggplot2 a binomický GLM

Zobecněné lineární modely v R

Richard Erickson

Instructor

Co vidím ve svých datech?

Ovlivňuje vzdálenost dojíždění pravděpodobnost jízdy autobusem?

ggplot(bus, aes(x = MilesOneWay, y = Bus)) + geom_point()

Graf pouze s body

geom_jitter()

gg_jitter <- ggplot(bus, aes(x = MilesOneWay, y = Bus)) +
             geom_jitter(width = 0, height = 0.05)
print(gg_jitter)

Graf bodů s jitterem

geom_smooth()

gg_jitter + geom_smooth()

Předchozí jitter graf s vyhlazeným geomem.

Faktor na číslo

str(bus)
bus$Bus2 <- as.numeric(bus$Bus) - 1

geom_smooth()

gg_jitter + geom_smooth()

Změna dat z charakteru na číslo

Lineární modely

gg_jitter + geom_smooth(method = 'glm')

Přidání lm do ggJitter

Logistická regrese

ggJitter + 
geom_smooth(method = 'glm',
            method.args = list(family = "binomial"))

ggJitter s logistickou křivkou

gg_jitter + 
    geom_smooth(method = 'glm',
                method.args = list(family = binomial(link = 'logit')),
                se = FALSE, color = 'red') +
    geom_smooth(method = 'glm',
                method.args = list(family = binomial(link = 'probit')),
                se = FALSE, color = 'blue')

ggJitter s křivkami logit a probit

Shrnutí postupu

Vykreslete jako jitter, aby se body nepřekrývaly
Přidejte vyhlazený geom
Zadejte správnou metodu a rodinu
Upravte graf (tento kurz se tím nezabývá)

Pojďme procvičovat!

Zobecněné lineární modely v R