Karpdiagramm (Boxplot)

Näide 1. Horisontaalne

# Algandmete genereerimine:

x=rexp(15)

# Karpdiagrammi joonistamine:

boxplot(x, horizontal=TRUE,
main="Karpdiagramm", col="orange", xlab="Tunnuse X väärtused")

Näide 2. Mitme grupi jaoks korraga

# Algandmete genereerimine:

x=rexp(100, rate=rep(c(1,2,5),c(15,25,60)))
grupp=factor(rep(c(1,2,5),c(15,25,60)), labels=c("Mõmmik","Mesikäpp","Mörrr"))

# Karpdiagrammi joonistamine:

boxplot(x~grupp, main="Kolme karu võrdlus", col="green")

Näide 3. Palju vaatlusi - paksuke, vähe vaatluseid - kõhnake

# Algandmete genereerimine:

x=rexp(203, rate=rep(c(1,2,5,1),c(8,25,80,90)))
grupp=factor(rep(c(1,2,5,11),c(8,25,80,90)),
labels=c("Mõmmik","Mesikäpp","Mörrr","Kaisukaru"))

# Karpdiagrammi joonistamine:

boxplot(x~grupp, main="Kolme karu võrdlus", col="green", varwidth=TRUE)

Näide 4. Tähendusi

# Algandmed:

x=c(1:10,5,5,5,5,11,20)

# Karpdiagrammi enda joonistamine:

boxplot(x, notch=TRUE, col="orange", at=0.8, main="karpdiagramm - sisuseletus")

# Joonise tähenduse lahtimõistatamine:

text(1.2, 1, "miinimum")
text(1.2, 20, "maksimum \n (erind - kahtlaselt suur)")
text(1.3, 12, "suurim 'usutav' väärtus")
abi=""
text(1.3, quantile(x,0.25)-0.3,bquote(italic(lower-hinge) %~~% .(abi)))
text(1.4, quantile(x,0.25)-1,"alumine kvantiil")

text(1.3, quantile(x,0.75)+0.4, bquote(italic(upper-hinge) %~~% .(abi)))
text(1.4, quantile(x,0.75)-0.3,"ülemine kvantiil")

text(1.15, 6.2, "mediaan")

text(0.5,9,"ligikaudne usaldusintervall mediaanile", srt=90)
stat=boxplot.stats(x)
arrows(0.52,8, 0.59, stat$conf)

arrows( c(1.1,1.14,1.05,1.14,1.09, 1.1), c(stat$stats, stat$out)-c(0,0.1,-1,-0.2,-1,0), c(0.91, 1.01, 0.91, 1.01, 0.91, 0.82), c(stat$stats, stat$out) )