Punkti kuju muutmine
Joonistatava punkti kuju saab muuta käsuga pch=, punkti suurust
saab muuta käsuga cex= ja tema värvi käsuga col= (värvide nimesid vaata
siit)
Näide 1.
# Algandmete genereerimine
x=runif(100,min=0, max=10)
sugu=rbinom(100,1,0.5)
y=3+0.8*x-0.03*x*x+rnorm(100,sd=(x/10+0.5))+1*sugu+0.04*x*x*sugu
# Joonistame (vaatlus)punktid
plot(x,y, ylim=c(0,15), pch=sugu+1, col=c("blue","red")[sugu+1])
# Leiame ja lisame graafikule jooned
mudel1=lm(y~poly(x,2)+as.factor(sugu)+as.factor(sugu)*I(x^2))
xx=seq(0,10,0.1)
y1=predict(mudel1, data.frame(x=xx, sugu=1))
y2=predict(mudel1, data.frame(x=xx, sugu=0))
y3=predict(mudel2, data.frame(x=xx))
lines(xx,y1, lty=2, col="red")
lines(xx,y2, lty=3, col="blue")
# Lisame selgitused (legendi)
legend(1,13,c("isased","emased"), pch=c(1,2), col=c("blue","red"), lty=c(3,2))
Näide 2.
Samad andmed, mis eelmiseski näites. Veidi teistsugune programm punktide joonistamiseks:
# Algandmete genereerimine
x=runif(100,min=0, max=10)
sugu=rbinom(100,1,0.5)
y=3+0.8*x-0.03*x*x+rnorm(100,sd=(x/10+0.5))+1*sugu+0.04*x*x*sugu
# Joonistame naisi iseloomustavad punktid
plot(x[sugu==0],y[sugu==0], ylim=c(0,15),
pch=21, col="red", bg="yellow", cex=2,
xlab="x-telg",ylab="y-telg", main="Lihtsalt näide")
# Lisame mehi iselomustavad punktid:
points(x[sugu==1],y[sugu==1],
pch=21, col="slateblue2", bg="skyblue1", cex=2, lwd=2)
# Lisame legendi:
legend(1,14, c("mehed","naised"),
col=c("slateblue2","red"), pch=c(21,21), bty="n",
pt.bg=c("skyblue1","yellow"), pt.cex=2, pt.lwd=c(2,1))
