Gráficas de Tiempos de Falla y Censura

A continuación se presenta una función que recibe como parámetros el vector de tiempo antes de la falla dentro del estudio y el vector de censura de estos tiempos. También se recibe como parámetro el tamaño de la muestra que se graficará. Esta función realiza un gráfico que permite observar los tiempos de falla o censura en el estudio para una muestra de tamaño n. La función recibe los parámetros, genera una muestra y transforma a factor el vector de censuras. Finalmente, realiza la gráfica para las fallas y censuras.

Parámetros

#Funcion que realiza grafico de tiempos de falla en un estudio
graf<-function(datos,cen,n){
  df<-data.frame(datos, cen)
  m<-df%>%select(datos, cen)%>%sample_n(n)
  m$cen<-as.factor(m$cen); colnames(m)[colnames(m)=="cen"]<-"Censura"
  ggplot(m)+geom_point(aes(x=1:n,y=m$datos, shape=Censura, colour=Censura),size=5)+
    geom_linerange(aes(x=1:n, ymin=0, ymax=m$datos), linetype='dashed')+
    coord_cartesian(ylim = c(0,max(m$datos)+5))+
    coord_flip()+
    ggtitle("Tiempos de falla en el estudio")+ylab("Tiempo")+xlab("Observacion")
  #+scale_color_discrete(name='Censura')
}

Hagamos ahora varios ejemplos de la función

Ejemplo 1

La base de datos trabaja con pacientes que presentan cáncer en sus ovarios. Se analiza el tiempo en morir y se clasifican los datos en censurados y no censurados.

data(ovarian)
set.seed(123)

graf(ovarian$futime,ovarian$fustat,11)

Ejemplo 2

La base de datos trabaja con paciente que presentan cáncer de médula. Se analiza el tiempo en el que fallece después de entrar al estudio o su salida del mismo. De manera análoga, se clasifican los datos en censurados (0) y no censurados (1).

data(myeloid)
graf(myeloid$futime, myeloid$death, 15)

Ejemplo 3

La base de datos trabaja con el tiempo de contrato de empleados en una cierta empresa y si han dejado o no la empresa.

hr <- read_csv("~/R/Estadistica 3/HR_comma_sep.csv")
graf(hr$time_spend_company, hr$left, 20)