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