percentiles.fnc

Copia, Pega y Adapta

percentiles.fnc(datos, variable='tr', que.factor='zona')

percentiles.fnc(datos, variable=2:12, que.factor='zona')

percentiles.fnc(datos, variable='tr', que.factor='zona',

percentiles=c(0.25,0.50,0.75,0.95,0.99))

OBJETIVO

Estima los percentiles, definidos por el usuario, tanto para toda la base de datos como por nivel del factor definido por el usuario.

Percentiles

Partiremos de una base de datos ya utilizada en la página de uniones de bases de datos (fundir.objetos.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.

percentiles.fnc(MathAchieve, variable='MathAch')

*** Percentiles de la variable: MathAch ***

P.0% P.25% P.50% P.75% P.100%

-2.832 7.275 13.131 18.317 24.993

percentiles.fnc(MathAchieve, variable='MathAch',

que.factor='Minority')

*** Percentiles de la variable: MathAch en cada nivel de: 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

Podemos pedir los percentiles de múltiples variables de una sola vez:

percentiles.fnc(MathAchieve, variable=4:6, que.factor='Minority')

$SES

P.0. P.25. P.50. P.75. P.100.

Minority.No -3.758 -0.378 0.142 0.692 1.652

Minority.Yes -2.838 -0.918 -0.388 0.222 2.692

$MathAch

P.0. P.25. P.50. P.75. P.100.

Minority.No -2.832 8.795 14.498 19.3255 24.993

Minority.Yes -2.832 4.519 9.404 14.8042 24.993

$MEANSES

P.0. P.25. P.50. P.75. P.100.

Minority.No -1.188 -0.171 0.111 0.360 0.831

Minority.Yes -1.188 -0.588 -0.296 0.159 0.831

percentiles.fnc(MathAchieve, variable='MathAch',

percentiles=c(0.90, 0.95),

que.factor='Minority')


*** Percentiles de la variable: MathAch en cada nivel de: Minority ***

P.90. P.95.

Minority.No 22.455 23.465

Minority.Yes 18.938 21.253

En este punto puedes tener interés en saber como discretizar una variable en cortes percentílicos.