CAPÍTULO 6. GRAFICACIÓN CON R
162
# Igual, aquí se incluye la expresión 'qq'
legend
(
"bottomright"
,
legend
=qq,
lty
=
1
,
pch
=
21
,
pt.bg
=cc)
El resultado de este código se puede ver el la Fig. 6.21. Nótese que, para
representar las funciones de densidad de probabilidades no se ha hecho con un
gráfico tipo agujas, ya que al sobreponer las tres funciones, se hubiera prestado
a confusión. En vez de ello, los puntos contiguos se han unido con una línea,
sólo como una guía, ya que entre dos puntos consecutivos la función no tiene
valor alguno.
Finalmente, queda por responder la pregunta de cuál es la probabilidad
de que en cada una de las líneas de cajas se presenten 14 o más clientes por
minuto. Para responderla, basta con usar la función
ppois()
, como se muestra
a continuación.
# En orden son para lambda=6, 10, 13
lapply
(lambdas,
function
(
ll
)
1
-
ppois
(
14
, ll))
## [[1]]
## [1] 0.0014
##
## [[2]]
## [1] 0.08346
##
## [[3]]
## [1] 0.3249
# En este caso se ha restado la probabilidad resultante
# de la unidad (1) porque el resultado directo es que el
# valor sea menor o igual a 14.
# El mismo efecto se logra cambiando el argumento
# lower.tail, como sigue
lapply
(lambdas,
function
(
ll
)
ppois
(
14
, ll,
lower.tail
=F))
## [[1]]
## [1] 0.0014
##
## [[2]]
## [1] 0.08346
##
## [[3]]
## [1] 0.3249
Así, las probabilidades de que se presenten 14 o más clientes por minuto
son 0.14%, 8.35% y 32.49% para las lineas de cajas de minusválidos, normales
y rápidas, respectivamente.
