### klaerschlamm.S ### Read in the klaerschlamm data -- and write it to ### the datasets.Data collection ### ~~~~~~ ### attach("/u/sfs/S/datasets.Data", pos = 1) klaerschlamm <-read.table("/u/sfs/ueb/datasets/klaerschlamm.dat",header = T) attach(klaerschlamm) motif() boxplot(Pr1,Pr2,Pr3,Pr4,Pr5,Pr6,Pr7,Pr8,Pr9,names=names(klaerschlamm)) ## arithm. Mittel fuer jede Probe apply(klaerschlamm,2,mean) ## Median fuer jede Probe apply(klaerschlamm,2,median) ## streuung fuer jede Probe sqrt(apply(klaerschlamm,2,var)) ## Quartilsdifferenzen fuer jede Probe quant_apply(klaerschlamm,2,quantile) quant[4,]-quant[2,] ##system. Fehler in den einzelnen Labors Fehler_klaerschlamm-matrix(rep(med,21),byrow=T,ncol=9) par(mfrow=(c(2,2))) boxplot(Fehler[1,],Fehler[2,],Fehler[3,],Fehler[4,],Fehler[5,],Fehler[6,],Fehler[7,],Fehler[8,],Fehler[9,],Fehler[10,],names=row.names(klaerschlamm)[1:10],ylim=range(c(Fehler))) boxplot(Fehler[11,],Fehler[12,],Fehler[13,],Fehler[14,],Fehler[15,],Fehler[16,],Fehler[17,],Fehler[18,],Fehler[19,],Fehler[20,],Fehler[21,],names=row.names(klaerschlamm)[11:21],ylim=range(c(Fehler))) mtext("Fehler der 21 Labors",outer=T) boxplot(Fehler[1,],Fehler[2,],Fehler[3,],Fehler[4,],Fehler[5,],Fehler[6,],Fehler[7,],Fehler[8,],Fehler[9,],Fehler[10,],names=row.names(klaerschlamm)[1:10]) boxplot(Fehler[11,],Fehler[12,],Fehler[13,],Fehler[14,],Fehler[15,],Fehler[16,],Fehler[17,],Fehler[18,],Fehler[19,],Fehler[20,],Fehler[21,],names=row.names(klaerschlamm)[11:21])