1. Programowanie Jak dodać linie do wykresu w R.

Autorstwa Andrie de Vries, Joris Meys

W R dodajesz linie do wykresu w bardzo podobny sposób, jak dodajesz punkty, z tym wyjątkiem, że do tego celu używasz funkcji lines ().

Ale najpierw użyj odrobiny R magii, aby utworzyć linię trendu w danych, zwaną modelem regresji. Użyj funkcji lm (), aby oszacować model regresji liniowej:

fit <- lm (oczekiwanie ~ erupcje, dane = wierny)

Wynikiem jest obiekt klasy lm. Użyj funkcji dopasowanej (), aby wyodrębnić dopasowane wartości z modelu regresji. Jest to przydatne, ponieważ można następnie wykreślić dopasowane wartości na wykresie. Zrób to potem.

Aby dodać tę linię regresji do istniejącego wykresu, wystarczy użyć linii funkcyjnych (). Możesz także określić kolor linii za pomocą argumentu col:

> fabuła (wierna)
> linie (wierne erupcje, dopasowanie (dopasowanie), col = „niebieski”)

Kolejną przydatną funkcją jest abline (). Umożliwia to rysowanie linii poziomych, pionowych lub nachylonych. Aby narysować pionową linię przy erupcjach pozycji == 3 w kolorze fioletowym, użyj następujących poleceń:

> abline (v = 3, col = „fioletowy”)

Wynikowa grafika powinna mieć pionową fioletową linię przy erupcjach == 3 i niebieską linię regresji.

image0.jpg

Aby utworzyć linię poziomą, użyj również abline (), ale tym razem określisz argument h. Na przykład utwórz poziomą linię w średnim czasie oczekiwania:

> abline (h = średnia (wierne oczekiwanie))

Możesz także użyć funkcji abline (), aby utworzyć nachyloną linię przez wykres. W rzeczywistości, podając argumenty aib, możesz narysować linię, która pasuje do równania matematycznego y = a + b * x. Innymi słowy, jeśli jako argumenty aib określisz współczynniki swojego modelu regresji, otrzymasz linię przez dane, które są identyczne z twoją linią predykcji:

> abline (a = coef (dopasowanie) [1], b = coef (dopasowanie) [2])

Co więcej, możesz po prostu przekazać obiekt lm do abline (), aby narysować linię bezpośrednio. (Działa, ponieważ istnieje metoda abline.lm ().) Dzięki temu kod jest bardzo łatwy:

> abline (fit, col = „red”)