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)
}
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)}
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)}
2004-2021 © Mathazay - Tous droits réservés
Ce site est géré sous
SPIP 3.0.17 [21515]
et utilise le squelette
EVA-Web 4.2
Dernière mise à jour : lundi 7 septembre 2015