Vous êtes ici : Accueil > Salle du prof > Il est libre Math > Les fonctions pour les boîtes à moustaches
Par : mathazay
Publié : 14 juin 2006

Les fonctions pour les boîtes à moustaches

Les paramètres « conformes » d’une série

fbam(serie)

function(serie){xt<-sort(serie)
xbam<-boxplot(xt,outliers=T,plot=F)
q1<-as.numeric(quantile(xt,0.25,type=1))
q3<-as.numeric(quantile(xt,0.75,type=1))
d1<-as.numeric(quantile(xt,0.1,type=1))
d9<-as.numeric(quantile(xt,0.9,type=1))
extr<-unique(c(subset(xt,xt<d1),subset(xt,xt>d9)))
nbextr<-rep(1,length(extr))
xbam$stats[1,1]<-d1
xbam$stats[2,1]<-q1
xbam$stats[4,1]<-q3
xbam$stats[5,1]<-d9
xbam$out<-extr
xbam$group<-nbextr
return(xbam)
}

La boîte à moustaches horizontale d’une série (serie), notée sur maxnote (10 ou 20)

mabam(serie)

function(serie,maxnote) {bxp(fbam(serie),h=T,axes=F,boxfill="gray",xlim=c(0,maxnote))
box()
abline(v=0:maxnote,lw=1,lty=2,col="light gray")
axis(1,0:maxnote)
bxp(fbam(serie),h=T,axes=F,boxfill="gray",xlim=c(0,maxnote),add=T)}

La boîte à moustaches horizontale d’une série (serie), notée sur maxnote (10 ou 20)

mabamquant(mini,d1,q1,med,q3,d9,maxi)

function(mini,d1,q1,med,q3,d9,maxi){
xbam<-boxplot(c(mini,d1,q1,med,q3,d9,maxi),plot=F)
xbam$stats[1,1]<-d1
xbam$stats[2,1]<-q1
xbam$stats[3,1]<-med
xbam$stats[4,1]<-q3
xbam$stats[5,1]<-d9
xbam$out<-c(mini,maxi)
xbam$group<-c(1,1)
return(xbam)
}
function(mini,d1,q1,med,q3,d9,maxi) {bxp(fbamquant(mini,d1,q1,med,q3,d9,maxi),h=T,axes=F,boxfill="gray",xlim=c(0,maxnote))
box()
abline(v=0:maxnote,lw=1,lty=2,col="light gray")
axis(1,0:maxnote)
bxp(fbamquant(mini,d1,q1,med,q3,d9,maxi),h=T,axes=F,boxfill="gray",xlim=c(0,maxnote),add=T)}