Simulación de Variables Aleatorias

El siguiente código es para simular una muestra de variables aleatoras exponencial con media lambda. Se muestra un ejemplo con las funciones de densidad, distribución y supervivencia. También se hace un ejemplo con variables normal estándar.

library(ggplot2)
library(survival)

Distribución Exponencial

Se simulan 1000 exponenciales con tasa 2

lambda <- 2
exp <- rexp(10000, rate = lambda)

Función de densidad con la muestra de exponenciales

f1 <- dexp(exp, rate = lambda)

Función de distribución de la muestra

Fx1 <- pexp(exp, rate = lambda)

Función de sobrevivencia de la muestra

S1 <- (1-Fx1)
###S1 <- qexp(exp, rate = lambda)

Tasa de riesgo/Fuerza de mortalidad/ Tasa Hazard

h1 <- f1/S1

Gráfica función de densidad

plot(exp,f1)

Gráfica función de distribución

plot(exp,Fx1)

Gráfica función de supervivencia

plot(exp,S1)

Distribución normal

mu <- 0
sigma <-1
nor <- rnorm(10000,mean = mu, sd = sigma)

f2 <- dnorm(nor, mean = mu, sd = sigma )
Fx2 <- pnorm(nor, mean = mu, sd = sigma)
S2 <- (1- Fx2)
h2 <- f2/S2

Gráfica función de densidad

plot(nor,f2)

Gráfica función de distribución

plot(nor,Fx2)

Gráfica función de supervivencia

plot(nor,S2)