quantile.fnc
Copy, Paste and Adapts
quantile.fnc(dat, variable='rt', which.factor='zone')
quantile.fnc(dat, variable='rt', which.factor='zone', quantiles=c(0.25,0.50,0.75,0.95,0.99))
Objetive
Estima los percentiles, definidos por el usuario, tanto para toda la base de datos como por nivel del factor definido por el usuario.
Quantiles
Partiremos de una base de datos ya utilizada en la página de uniones de bases de datos (merge.data.fnc). Hablamos de MathAchieve de la librería nlme que contiene la medida del rendimiento en matemáticas de 7185 alumnos. Activamos la base de datos MathAchieve mediante esta llamada:
data(package='nlme','MathAchieve')
head(MathAchieve)
School Minority Sex SES MathAch MEANSES
1 1224 No Female -1.528 5.876 -0.428
2 1224 No Female -0.588 19.708 -0.428
3 1224 No Male -0.528 20.349 -0.428
4 1224 No Male -0.668 8.781 -0.428
5 1224 No Male -0.158 17.898 -0.428
6 1224 No Male 0.022 4.583 -0.428
Vamos a solicitar los percentiles (cuartiles) de la variable rendimiento en matemáticas (MathAch) primero con caracter general para toda la base de datos y posteriormente en cada nivel del factor Minority.
quantile.fnc(MathAchieve, variable='MathAch')
*** Quantiles of variable: MathAch ***
P.0% P.25% P.50% P.75% P.100%
-2.832 7.275 13.131 18.317 24.993
quantile.fnc(MathAchieve, variable='MathAch',
which.factor='Minority')
*** Quantiles of variable: MathAch in each level of: Minority
P.0. P.25. P.50. P.75. P.100.
Minority.No -2.832 8.795 14.498 19.326 24.993
Minority.Yes -2.832 4.519 9.404 14.804 24.993
quantile.fnc(MathAchieve, variable='MathAch',
quantiles=c(0.90, 0.95),
which.factor='Minority')
*** Quantiles of variable: MathAch in each level of: Minority
P.90. P.95.
Minority.No 22.455 23.465
Minority.Yes 18.938 21.253
En este punto puedes tener interes en saber como cuts una variable en cortes percentílicos.