CAPÍTULO 6. GRAFICACIÓN CON R
159
## [1] 0.0000454 0.0004994 0.0027694 0.0103361 0.0292527
## [6] 0.0670860 0.1301414 0.2202206 0.3328197 0.4579297
## [11] 0.5830398 0.6967761 0.7915565 0.8644644 0.9165415
## [16] 0.9512596 0.9729584 0.9857224 0.9928135 0.9965457
## [21] 0.9984117
# Con la función nativa del lenguaje, se usa 'x'
Ps10
<-
ppois
(x,
lambda
=
10
)
Ps10
## [1] 0.0000454 0.0004994 0.0027694 0.0103361 0.0292527
## [6] 0.0670860 0.1301414 0.2202206 0.3328197 0.4579297
## [11] 0.5830398 0.6967761 0.7915565 0.8644644 0.9165415
## [16] 0.9512596 0.9729584 0.9857224 0.9928135 0.9965457
## [21] 0.9984117
# Los resultados son "semejantes"
Al igual que para el caso de la Fig. 6.19, se pueden producir los gráficos
correspondientes a las dos funciones por medio de los datos obtenidos,
ds10
y
Ps10
, con el siguiente código.
# La función de densidad:
plot
(
x
=x,
y
=ds10,
type
=
"h"
,
lwd
=
2
,
xlab
=
"y"
,
ylab
=
"ds10(y)"
)
# Por claridad se agregan puntos al gráfico:
points
(x, ds10,
pch
=
21
,
bg
=
"red"
)
# La función de distribución:
plot
(x, Ps10,
type
=
"s"
,
col
=
"red"
,
lwd
=
2
,
xlab
=
"y"
,
ylab
=
"Ps10(y)"
)
points
(x, Ps10,
pch
=
16
)
El resultado de este código se muestra en la Fig. 6.20. Nótese que en todos
los llamados a la funciones
plot()
y
points()
, se ha tenido que especificar
tanto las X como las Y. Esto se debe a que ahora el rango de las X comienza con
0 y no con 1, como anteriormente.
El problema original, enunciado al principio de esta sección, en la pági-
na 156, se mencionan tres líneas de cajas, con razones promedio de 13, 10 y
6 clientes por minuto. Si se quiere ver los tres gráficos correspondientes a las
funciones de densidad de probabilidades y los tres, correspondientes a las fun-
ciones de distribución de probabilidades, se pueden generar con el código que
se presenta a continuación.
x
<-
0
:
20
# vector con los datos
# Colores para cada función:
cc
<-
c
(
"yellow"
,
"springgreen"
,
"navyblue"
)
# Las razones promedio para cada función
