Ekraani jagamine võrdse suurusega tükkideks (mfrow, mfcol)

Graafikaparameetrite mfrow ja mfcol korral tuleb öelda, mitmeks reaks ja mitmeks veeruks graafikaaken tuleb jagada. Järjestikused joonistatavad pildid lähevad siis järgmistesse "pesadesse". See, millises järjekorras graafikuid joonistatakse, on mfrow-l ja mfcol-il erinev: ühel juhul joonistatakse graafikuid ridade, teisel juhul veergude kaupa.
       

Näide 1

# Algandmed

x=rnorm(100)

# Joonistamine

# Jagame ekraani 6 osaks - 2 rida ja 3 veergu graafikuid
par(mfrow=c(2,3))

# Jätame 6 pildi ühise tiitli jaoks rohkem ruumi
par(oma=c(1,1,3,1))

# Joonistame kuus pilti
plot(x, bg=rgb(0,0,runif(100)), cex=2, pch=21, col="blue", main="plot")
hist(x, col=heat.colors(10),main="hist")
barplot(table(cut(x, breaks=5)), col=terrain.colors(9), main="barplot")
pie(table(cut(x, breaks=5)), main="pie", labels=NA)
qqnorm(x, main="qqnorm"); qqline(x)
plot(ecdf(x))

# Lisame ühise tiitli
title("Kuus pilti - par(mfrow(2,3))", outer=TRUE, cex.main=1.75)


   

Näide 2

Sageli peaksid erinevate graafikute teljed olema teineteisega võrreldavad. Sellisel juhul on soovitav graafikute joonistamisel kasutada xlim= ja ylim= parameetreid. Toodud näites paigutatakse valimi tihedus- ja jaotusfunktsioon samale skaalale just tänu xlim-parameetri kasutamisele:
# Algandmete genereerimine:
x=rchisq(500, df=3)

# Ekraani jagamine kaheks
par(mfrow=c(2,1))

# Jätame ruumi ühisele tiitlile
par(oma=c(1,1,4,1))

# Joonistame kaks graafikut
plot(density(x), xlim=c(0,5), main="tihedusfunktsioon")
plot(ecdf(x), xlim=c(0,5), main="jaotusfunktsioon")

# Lisame tiitli
title("Tihedus- ja jaotusfunktsiooni hinnangud", outer=TRUE)


   

Näide 3

Kui soovitakse jooniste maatriksis üle joonistada mõnda konkreetset pilti, saab seda teha graafikaparameetri mfg=c(rida, veerg) abil:
# Jagame ekraani 9 osaks - 3 rida ja 3 veergu graafikuid
par(mfrow=c(3,3))

# Vähendame ühe pildi ümber jäetavat vaba ruumi
par(mar=c(3,3,1,1))

# Joonistame üheksa pilti (9 valimi histogrammid)
for (i in 1:9){
  x=rnorm(40)
  hist(x, col=heat.colors(10), main="", xlim=c(-3,3))
}

# Täiendame 2-rea 3. veeru pilti kakukesega:
par(mfg=c(2,3))
pie(1:10, col=terrain.colors(10), labels=NA)






Programm, mis genereeris siin leheküljel olevad pildid