Funkcje łączące: probit a logit

Uogólnione modele liniowe w R

Richard Erickson

Instructor

Po co funkcje łączące?

  • Rozumienie i symulowanie GLM
  • Probit vs logit jako przykład
Uogólnione modele liniowe w R

Po co probit?

  • Ilustruje funkcję łączącą
  • Stosowany w niektórych dziedzinach (np. toksykologia)
  • Preferowany przez część badaczy
Uogólnione modele liniowe w R

Czym jest probit?

  • Jednostka probabilistyczna
  • Toksykologia — Chester Bliss, 1934
  • Obliczeniowo prostszy niż logit
  • Znany jako analiza probitowa, regresja probitowa lub model probitowy
Uogólnione modele liniowe w R

Równanie probitowe

Model danych binarnych

$Y = \text{Binomial}(p)$

Powiązany z równaniem liniowym

$\Phi^{-1} (p) = \beta_0 + \beta_1 x + \epsilon $

Uogólnione modele liniowe w R

Funkcja probitowa

Oparty na dystrybuancie rozkładu normalnego

$\Phi(z) = \frac{1}{\sqrt{2 \pi}} \int_{-\infty}^z e^{- \frac{1}{2} z^2} \, dz$

Uogólnione modele liniowe w R

Skumulowane prawdopodobieństwo probitu i logitu

Uogólnione modele liniowe w R

Dopasowywanie probitu w R

  • Opcja family w glm()
    • Znak: glm(..., family = "binomial")
    • Funkcja: glm(..., family = binomial())
  • Domyślnie: binomial(link = "logit")
  • Probit: binomial(link = "probit")
  • Stosować zgodnie z instrukcjami DataCamp
Uogólnione modele liniowe w R

Symulacja z probitem

Przeliczenie ze skali probitu na skalę prawdopodobieństwa:

p = pnorm(-0.2)

Użycie prawdopodobieństwa z rozkładem dwumianowym

rbinom(n = 10, size = 1, prob = p)
Uogólnione modele liniowe w R

Symulacja z logitem

Przeliczenie ze skali logitu na skalę prawdopodobieństwa:

p = plogis(-.2)

Użycie prawdopodobieństwa z rozkładem dwumianowym

rbinom(n = 10, size = 1, prob = p)
Uogólnione modele liniowe w R

Kiedy używać probitu, a kiedy logitu?

  • W dużej mierze zależne od dziedziny
  • Grubsze ogony logitu
  • Oba są dopuszczalne
Uogólnione modele liniowe w R

Czas na ćwiczenia!

Uogólnione modele liniowe w R

Preparing Video For Download...