Nooled, usalduspiirid (arrows-käsk)


Noolte ja usalduspiiride kandmiseks joonisele sobib arrows-käsk:
    
plot(x,y, main="Erind", ylim=c(0,20))
arrows(8, 4, 10.8, 6.9)
text(8,4, "Veider\n vaatlus", adj=c(1, 0.5))
    
plot(ecdf(tunnus), xlim=c(-3,4), ylim=c(0,1), 
	verticals=T, do.points=F, 
	main="Emp. ja teoreetiline jaotusfn")

x1=seq(-3,4, length=100)
y1=pnorm(x1)
lines(x1,y1, col="gray50", lwd=2)

arrows(-0.45, ecdf(tunnus)(-0.45), -0.45, pnorm(-0.45), 
	code=3, angle=15, col="darkred", length=0.1)
text(-0.45, (ecdf(tunnus)(-0.45)+pnorm(-0.45))/2, 
	" Maksimaalne erinevus", col="darkred", adj=c(0, 0.5))
    
a=barplot(prop.table(table(y,x),2)[1,], 
	ylim=c(0,1), main="95%-Usaldusintervall", 
	ylab="Tõenäosus")
lim1 = binom.test(table(y[x==0]))$conf.int
lim2 = binom.test(table(y[x==1]))$conf.int
lim3 = binom.test(table(y[x==2]))$conf.int
arrows(a, c(lim1[1], lim2[1], lim3[1]), 
       a, c(lim1[2], lim2[2], lim3[2]), 
	code=3, angle=90, col="gray40", lwd=2)
Süntaks:

arrows(<kust X>,<kust Y>,<kuhu X>,<kuhu Y>, code=3, angle=90, length=0.1, ...)

kus:
   <kust X>, <kust Y> - alguspunkti koordinaadid, kust alata noole joonistamist
   <kuhu X>, <kuhu Y> - lõpp-punkti koordinaadid, kuhu nool suundub
    code=... - Kuhu otsa nool teha. Vaikimisi code=2 (nool tehakse kuhu-otsa); code=3 - nool mõlemasse otsa (kasulik usaldusintervallide joonistamisel)
    angle=... - nurk e. kui terav nool teha. Usaldusintervallide lisamisel joonisele vahel kaval kasutada angle=90.
    length=... - Kui pikka nooleotsa soovitakse (tollides).
    ... - muud joone värvi (col=...), laiust (lwd=...) või tüüpi (lty=...) määravad parameetrid.

Referentsjooned (abijooned)


Abijoonte lisamiseks sobivad grid(), abline(v=...), abline(h=...)käsud.
    
plot(x,y)
grid()
    
plot(x,y)
grid(lty=3, col="orange")
points(x,y, pch=21, col="green4", 
	bg="green", lwd=2, cex=1.5)
    
barplot(table(x))
abline(h=seq(0,50,10), col="gray50", lty=2)
barplot(table(x), add=T)

Vaatluste lisamine joonise serva


Vahel soovime vaatluspunkte lisada joonise servadesse. Selleks sobib rug()-käsk:
    
hist(x)
rug(x)
    
plot(density(x), main="rug(x)", 
	col="gray25", lwd=2, xlab="x",ylab="f(x)")
rug(x)

Graafiku lisamine joonisele


Graafiku lisamiseks olemasolevale joonisele saab kasutada käsku par(plt=c(...)). Täpsemalt vaata graafikute kombineerimise kohta siit.
    

barplot(table(Kvaliteet, Linn),  
	col=c("white", "lightblue", "mistyrose"), 
	main="Katuste kvaliteet")
par(plt=c(0.4, 1, 0.55, 0.8), new=T)
pie(table(Kvaliteet), cex=0.8)

# Algseisundi taastamine
par(plt=c(0, 1, 0, 1))




Programm, mis genereeris pildid siia leheküljele.