Argument family dla rozkładu dwumianowego

Główna różnica między regresją liniową wykonywaną za pomocą lm() a regresją logistyczną wykonywaną za pomocą glm() polega na tym, że w glm() trzeba ustawić argument family na binomial. binomial() to funkcja zwracająca listę innych funkcji, które wskazują glm(), jak przeprowadzać obliczenia podczas regresji. Dwie najciekawsze z nich to linkinv i linkfun – służą do przekształcania zmiennych z całej osi liczbowej (od minus nieskończoności do nieskończoności) na prawdopodobieństwa (od zera do jedynki) i z powrotem.

Dostępny jest wektor wartości x oraz wektor prawdopodobieństw p.

To ćwiczenie jest częścią kursu

Regresja średnio zaawansowana w R

Zobacz kurs

Instrukcje do ćwiczenia

Zbadaj strukturę funkcji binomial(). Zwróć uwagę, że zawiera ona dwa elementy będące funkcjami: binomial()$linkinv* oraz* binomial()$linkfun.
Wywołaj binomial()$linkinv() na x, przypisując wynik do linkinv_x.
Sprawdź za pomocą all.equal(), czy linkinv_x i plogis() zastosowane na x dają takie same wyniki.
Wywołaj binomial()$linkfun() na p, przypisując wynik do linkfun_p.
Sprawdź, czy linkfun_p i qlogis() zastosowane na p dają takie same wyniki.

Interaktywne ćwiczenie praktyczne

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

# Look at the structure of binomial() function
___

# Call the link inverse on x
linkinv_x <- ___

# Check linkinv_x and plogis() of x give same results 
___

# Call the link fun on p
linkfun_p <- ___

# Check linkfun_p and qlogis() of p give same results  
___

Edytuj i uruchom kod