CAPÍTULO 6. GRAFICACIÓN CON R
155
F
(
y
n
) =
n
∑
i
=
1
f
(
y
i
)
(6.1)
En R, entonces, las dos funciones se pueden definir por medio del código
que se presenta a continuación.
# Función de densidad de probabilidades
ff
<-
rep
(
1
/
6
,
6
)
ff
## [1] 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
# En particular,
ff[
3
]
## [1] 0.1667
# Función de distribución de probabilidades
FF
<-
Reduce
(
"+"
, ff,
accumulate
=T)
FF
## [1] 0.1667 0.3333 0.5000 0.6667 0.8333 1.0000
# En particular,
FF[
3
]
## [1] 0.5
Nótese que, en este caso, en vez de utilizar la sintaxis para definición y uso
de funciones, se utiliza directamente la sintaxis de vectores, debido a que se
está tratando con variables discretas con valores de un conjunto finito. Ade-
más, para definir la distribución de probabilidades, se ha empleado la función
Reduce()
, con primer argumento
“+”
, que toma el vector
ff
y va aplicando
consecutivamente la función dada como primer argumento, la suma, en este
caso, entre el resultado anterior de la aplicación y cada elemento del vector.
Además, como se ha visto en la descripción de esta función en la sección 4.3.1
en la página 72, el argumento
accumulate=T
, obliga a entregar cada uno de los
resultados intermedios en el resultado final; de otro modo, el resultado sería
sólo la suma final y el llamado equivaldría al usar la función
sum(ff)
. Los grá-
ficos correspondientes a estas funciones se pueden hacer con el código que se
presenta a continuación.
# Para la función de densidad de probabilidades,
# se utiliza el tipo de de gráfico de agujas
plot
(ff,
type
=
"h"
,
col
=
"red"
,
lwd
=
2
,
xlab
=
"y"
,
ylab
=
"f(y)"
)
