gg_qqplot: Graphique des quantiles (qqplot) et droite interquartile

Description

Dessine le graphique des quantiles ou qqplot et la droite interquartile (passant par le premier et le troisi?me quartile ? la mani?re de la fonction qqline) avec la biblioth?que graphique ggplot2.

Usage

gg_qqplot(df,var,qdist=qnorm,params=list(),qq.line=TRUE,color="red",alpha=.5)

Arguments

Un jeu de donn?es (dataframe)

var

Le nom d'une variable de df

qdist

La fonction quantile d'une (famille de) distribution. Par d?faut celle de la famille des lois normales.

params

Une liste de param?tres pour sp?cifier la loi ? utiliser. Par d?faut la loi normale centr?e et r?duite. Les param?tres peuvent ?tre estim?s avec la fonction fitdistr de la biblioth?que MASS.

qq.line

Une valeur logique. Affiche ou masque la droite interquartile.

color

Le nom d'une couleur. Sp?cifie la couleur ? utiliser pour la droite interquartile.

alpha

Indice de transparence. Sp?cifie la transparence ? utiliser pour repr?senter les valeurs de l'?chantillon.

Value

ggplot

Un graphique utilisant la biblioth?que ggplot2. Affiche les valeurs des quartiles th?oriques par lesquels passe la droite ainsi que son ordonn?e ? l'origine et sa pente si la trac? de celle-ci est demand?.

References

F. Bertrand, M. Maumy-Bertrand, Initiation ? la Statistique avec R, Dunod, 2e, 2014.

Examples

Run this code

# NOT RUN {
glycine.blanche<-subset(Mesures,subset=(Mesures$espece=="glycine blanche"))
gg_qqplot(glycine.blanche,"taille")

#bonus ajustement avec une autre loi (ici Student (car dist = qt) dont on estime les ddl)
lauriers.roses<-subset(Mesures,subset=(Mesures$espece=="laurier rose"))
shapiro.test(lauriers.roses$taille) 
#pas issu d'une loi normal au risque alpha=5%
gg_qqplot(lauriers.roses,"taille")
gg_qqplot(lauriers.roses,"taille",qq.line=FALSE)
#essayons un qqplot avec une loi de Student
# }
# NOT RUN {
require(MASS)
params <- as.list(fitdistr(lauriers.roses$taille, "t")$estimate)
#avec la droite
gg_qqplot(lauriers.roses,"taille",qt,params)
#essayons un qqplot avec une loi gamma
params <- as.list(fitdistr(lauriers.roses$taille,"gamma")$estimate)
#avec la droite
gg_qqplot(lauriers.roses,"taille",qgamma,params)
#essayons un qqplot avec une loi du chi-deux
params <- list(df=fitdistr(lauriers.roses$taille,"chi-squared",start=list(df=5),
method="Brent",lower=1,upper=40)$estimate)
#avec la droite
gg_qqplot(lauriers.roses,"taille",qchisq,params)
# }