plot-käsk

Plot-käsk

Näide 1. Joone joonistamine

# Algandmed

y=c(450969, 572836, 579418, 763479, 871488)/1000
x=1998:2002

# Joonistamine

plot(x, y, type="l",
    xlab="aasta", ylab="miljonites kroonides",
    main="Kulud teadus- ja arendustegevusele \n Eesti Vabariigis (eraettevõtted+riik)",
    lwd=3, ylim=c(0,900))

# Märkused:
# type="l" - ühendab x ja y-koordinaatige antud punktid joonega,
# proovi ka type="b" (joonistatakse jooned ja punktid)
# lwd=3 - tavalisest paksem joon

Näide 2. Lihtsad punktid

# Algandmete tekitamine

x=seq(0,10,0.1)
y=log((x)/12)+rnorm(101, sd=0.5)

# Punktide joonistamine

plot(x,y)

Näide 3. Trepi-tüüpi joonis

plot(ecdf(rnorm(15, sd=0.5)),
   main="Jõevee kõrgus. Jaotusfunktsiooni hinnang",
   ylab="Tõenäosus", xlab="Veetase (m)")

# Kommentaar:
#     funktsioon ecdf (empirical cumulative distribution function) hindab vaatluste põhjal jaotusfunktsiooni.
# NB! plot - käsu käitumine sõltub ka sellest, mida kästakse joonistada.
# Vaata ka kahte järgmist näidet!

Näide 4. Erinevaid võimalusi treppide joonistamiseks

# Algandmete tekitamine:

x=c(1,2,4,7,10)

# Joonistamine:
par(mfrow=c(2,2)) # Neli joonist korraga ekraanile

plot(ecdf(x), main="plot(ecdf(x))")
plot(ecdf(x), verticals=TRUE, main="plot(ecdf(x), verticals=TRUE)")
plot(x,1:5, type="s", main="plot(x,1:5, type='s')")
plot(x,1:5, type="S", main="plot(x,1:5, type='S')")

par(mfrow=c(1,1)) # Taastame olukorra - üks joonis korraga

Näide 5. Plot käsu erinev käitumine sõltuvalt argumendi tüübist

x=rnorm(15, sd=0.5)
y=10+2*x+rnorm(15)

par(mfrow=c(2,3)) # Mitme graafiku mahutamiseks

plot(ecdf(x))
plot(lm(y~x), main="lm(y~x)")
plot(density(x))

par(mfrow=c(1,1)) # Taastame algseisu - 1 graafik korraga

# Pane tähele! Käsk plot(lm(y~x)) joonistab koguni 4 graafikut korraga (kõik seotud regressioonmudeli sobivuse kontrolliga)!