Omezení lineárních modelů

Zobecněné lineární modely v R

Richard Erickson

Instructor

Přehled kurzu

Kapitola 1: Přehled a limity lineárního modelu a Poissonovy regrese
Kapitola 2: Logistická (binomická) regrese
Kapitola 3: Interpretace a vizualizace GLM
Kapitola 4: Vícenásobná regrese s GLM

Tahoun datové vědy

Lineární modely: Tahoun datové vědy

¹ US Department of Agriculture https://www.nal.usda.gov/exhibits/ipd/localfoods/exhibits/show/farmto–table/the–roads–of–rural–america

Lineární modely

Jak lineární koeficienty vysvětlují data?
Intercept pro základní efekt
Sklon pro lineární prediktor
$y = \beta_0 + \beta_1 x + \epsilon $

Lineární modely v R

lm(y ~ x, data = dat)

Předpoklad linearity

Příklad lineárních a nelineárních dat

Předpoklad normality

Normální vs. nenormální rezidua

Předpoklad spojitých proměnných

Spojitá vs. nespojitá data .

Krabicový diagram hmotností kuřat na konci studie

Vliv diety na hmotnost kuřat

Data ChickWeight z balíčku datasets
ChickWeightsEnd – poslední pozorování ze studie
Jak se diety 2, 3 a 4 liší od diety 1?

lm(formula = weight ~ Diet, data = ChickWeightEnd)

    Call:
    lm(formula = weight ~ Diet, data = ChickWeightEnd)
    Coefficients:
    (Intercept)        Diet2        Diet3        Diet4  
         177.75        36.95        92.55        60.81

Co přežití nebo počty?

Co přežití nebo počty kuřat?
Ani jedno není spojité!
Potřebujeme nový nástroj
Zobecněný lineární model

Zobecněný lineární model

Podobné lineárním modelům
Nenormální rozdělení chyb
Spojovací funkce: $y = \psi(b_0 + b_1 x + \epsilon)$

GLM v R

glm( y ~ x, data = data, family = "gaussian")

lm() je totéž jako glm( ..., family = "gaussian")

Lass uns üben!

Zobecněné lineární modely v R