diagnostyka: Diagnostyka regresji.

In tzoltak/EWDogolny: Utility functions developed while working on Polish Educational Value-Added indicators

Description

Funkcja rysuje wykresy diagnostyczne dla modeli regresji i w tym sensie powiela funkcjonalno<c5><9b><c4><87> print(lm(...)), jednak jest zoptymalizowana do tego, <c5><bc>eby dawa<c4><87> <c5><82>adne wykresy dla du<c5><bc>ych danych (i dawa<c4><87> je szybko). Poza tym obs<c5><82>uguje automatyczne zapisywanie wykres<c3><b3>w do plik<c3><b3>w PNG.

Usage

diagnostyka(model, zapiszPng = NULL, tytul = NULL, smoothScatter = TRUE,
  loess = TRUE, span = 0.5, control = loess.control(surface =
  "interpolate", statistics = "approximate", trace.hat = "approximate", cell =
  1), paleta = colorRampPalette(c("white", blues9)), pch = 16, cex = 1,
  lwd = 2, alpha = NULL)

Arguments

model	model regresji (typowo wynik dzia<c5><82>ania funkcji lm, glm lub lmer)
zapiszPng	tekst - <c5><9b>cie<c5><bc>ka do katalogu, w kt<c3><b3>rym maj<c4><85> by<c4><87> zapisane pliki PNG z wykresami lub NULL, je<c5><9b>li wykresy maj<c4><85> nie by<c4><87> zapisywane
tytul	tekst - tytu<c5><82> wykres<c3><b3>w (od tytu<c5><82>u b<c4><99>d<c4><85> si<c4><99> te<c5><bc> rozpoczyna<c4><87> nazwy zapisywanych plik<c3><b3>w PNG)
smoothScatter	warto<c5><9b><c4><87> logiczna - czy do rysowania wykres<c3><b3>w rozrzutu u<c5><bc>ywa<c4><87> funkcji smoothScatter?
loess	warto<c5><9b><c4><87> logiczna - czy do wykres<c3><b3>w rozrzutu dodawa<c4><87> lini<c4><99> regresji nieparametrycznej (estymowanej funkcj<c4><85> loess)?
span	parametr steruj<c4><85>cy wyg<c5><82>adzaniem regresji nieparametrycznej (p. loess)
control	parametry regresji nieparametrycznej - lista zwr<c3><b3>cona przez funkcj<c4><99> loess.control
paleta	funkcja definiuj<c4><85>ca palet<c4><99> kolor<c3><b3>w do wykorzystania przy rysowaniu wykres<c3><b3>w, typowo wynik wywo<c5><82>ania funkcji colorRampPalette
pch	je<c5><9b>li smoothScatter=FALSE, znak punktu u<c5><bc>ywanego do rysowania obserwacji (p. points)
cex	je<c5><9b>li smoothScatter=FALSE, wielko<c5><9b><c4><87> punkt<c3><b3>w u<c5><bc>ywanych do rysowania obserwacji (p. points)
lwd	je<c5><9b>li loess=TRUE, grubo<c5><9b><c4><87> rysowanej linii (p. lines)
alpha	je<c5><9b>li smoothScatter=FALSE i nie odpowiada Ci wynik dzia<c5><82>ania domy<c5><9b>lnego argumentu ustalajacego, jak bardzo prze<c5><ba>roczyste s<c4><85> punkty, mo<c5><bc>esz zada<c4><87> prze<c5><ba>roczysto<c5><9b><c4><87> sam

Details

Rysowane s<c4><85> nast<c4><99>puj<c4><85>ce wykresy:

• reszty w funkcji przewidywania do diagnostyki niezale<c5><bc>no<c5><9b>ci reszt od przewidywania,

• wykres kwantyl-kwantyl do diagnostyki normalno<c5><9b>ci rozk<c5><82>adu reszt,

• pierwiastek z modu<c5><82>u reszt w funkcji przewidywania do diagnostyki homoscedatyczno<c5><9b>ci.

Dodatkowo dla modeli mieszanych efekt<c3><b3>w, dla ka<c5><bc>dego efektu losowego:

• wykres kwantyl-kwantyl do diagnostyki normalno<c5><9b>ci rozk<c5><82>adu BLUPs<c3><b3>w,

• BLUPsy w funkcji <c5><9b>redniej przewidywa<c5><84> z cz<c4><99><c5><9b>ci sta<c5><82>ej do diagnostyki niezale<c5><bc>no<c5><9b>ci cz<c4><99><c5><9b>ci sta<c5><82>ej i losowej modelu.

Value

funkcja nic nie zwraca

See Also

lm, glm, lmer, smoothScatter, loess

Examples

## Not run: 
  x = rnorm(100000)
  y = x + rnorm(100000)
  diagnostyka(lm(y ~ x), tytul = "y~x")

  library(lme4)
  fm1 = lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
  diagnostyka(fm1)

## End(Not run)

tzoltak/EWDogolny documentation built on May 3, 2019, 2:05 p.m.