Proyecto 3

Explicación de la función

Esta función la utilizamos para automatizar el análisis de los datos de supervivencia, queremos ver cuando se resulta de la primera falla dependiendo del caso.

library(survival)  
library(KMsurv) #se cargan las librerias Survival y KMsurv para poder analizar los datos de genfan y aids

library(dplyr)
library(ggplot2)

data("genfan")
genfan %>% View()

data("aids")
aids %>% View()

data = read.csv("C:\Users\acer1\Documents\septimo semestre\estadistica\HR_comma_sep.csv")

f <- function(columna, columna_f, muestra){ #La funcion llevara tres parametros la columna de tiempo, la de tiempo de falla y el numero de observaciones que queremos.
  m <- data.frame(columna,columna_f) %>% #hacemos un data frame para poder usar las funciones de ggplot2
    sample_n(muestra)
  
  ggplot(m) + 
    geom_point(aes(y = m$columna,
                   x = 1:muestra,
                   shape = m$columna_f %>% as.factor(),
                   colour = m$columna_f %>% as.factor()),
               size = 3) +
    geom_linerange(aes(x = 1:muestra,
                       ymin = 0,
                       ymax = m$columna), #la linea s昼㸳lo ira hasta el punto de donde sea la falla
                   linetype = "dashed") +
    labs(x = "Tiempo",
         y = "Observaciones",
         shape = "Ocurrencia del evento",
         colour = "Ocurrencia del evento") +
    coord_flip() #este ultimo es para voltear el gr攼㸱fico y as攼㹤 observar de manera mas practica el tiempo de falla
}

f(aids$infect,aids$adult,60)
f(genfan$hours,genfan$status,60)
f(HR_comma_sep$time_spend_company,HR_comma_sep$left,60)

Aids

Genfan

HR