El arte de programar en R Un leguaje para la estadística - page 113

CAPÍTULO 5. ESCRITURA DE FUNCIONES
110
Figura 5.8: Ajuste del histograma a una curva continua
σ
. Estos valores, para el conjunto de datos que se ha venido manejando, son:
( mu
<-
mean
(pp
$
Precip) )
# media
## [1] 115.5
( sigma
<-
sd
(pp
$
Precip) )
# desviación estandar
## [1] 56.01
En R, la función de densidad normal de probabilidades es
dnorm()
, y si se
quiere saber su valor para varias precipitaciones, se puede hacer como sigue:
valores
<-
c
(
20
,
50
,
115
,
150
,
200
)
# precipitaciones
dnorm
(valores, mu, sigma)
## [1] 0.001666 0.003597 0.007122 0.005890 0.002281
Como
µ
y
σ
no varían para un conjunto de observaciones dado, conviene
hacer una función que solamente dependa de un solo argumento,
x
.
dnormX
<-
function
(
x
)
dnorm
(x, mu, sigma)
dnormX
(valores)
## [1] 0.001666 0.003597 0.007122 0.005890 0.002281
Graficar una función como ésta, en R es sencillo y se hace mediante la fun-
ción
curve()
. Para obtener una buena gráfica de la función de densidad normal
de probabilidades, conviene saber que, más o menos, el rango de valores de las
1...,103,104,105,106,107,108,109,110,111,112 114,115,116,117,118,119,120,121,122,123,...198
Powered by FlippingBook