SOBREVIVENCIA

Importamos las siguientes librerias

library(survival)
library(survminer)
library(flexsurv)
library(dplyr)

Cargamos datos

times = rexp(500)
event = rbinom(500,1,.7)
df = data.frame(times=times, event=event)
Surv(times, event)
surv_obj=Surv(times, event)
s_fit = flexsurvreg(surv_obj~1, data = df, dist = "exp")

s_fit$n.risk
s_fit$n.event
s_fit$n.censor
s_fit$surv
s_fit$cumhanz

Graficamos funciones de sobrevivencia

plot(surv_obj)

ggsurvplot(s_fit,df)

ggsurvplot(s_fit,cumevents= TRUE)

ggsurvplot(s_fit,fun ="cumhaz")

ggsurvplot(s_fit,surv.median.line = "hv")

ggsurvplot(s_fit,risk.table =TRUE)

Riesgo de distribuciones

Distribucion exponencial con lamda = 3

rate=3
exp=rexp(100, rate=rate)
f_exp=dexp(exp, rate = rate)
F_exp=pexp(exp, rate = rate)
S_exp=1-F_exp
h_exp=f_exp/S_exp
plot(h_exp)

Gamma con parametros (100,2)

g =rgamma(100, 2)
f_g=dgamma(g, 2)
F_g=pgamma(g, 2)
S_g=1-F_g
h_t=f_g/S_g
plot(sort(h_t),type="l")

Weibull con parametros (100,2,3)

a =rweibull(100,shape =  2,scale = 3)
f_w=dweibull(a,shape =  2,scale = 3)
F_w=pweibull(a,shape =  2,scale = 3)
S_w=1-F_w
h_w=f_w/S_w
plot(h_w %>% sort)

plot de una funcion

plot(1:100,18*(1:100))

Funcion comparar riesgo

Exponencial

riesgo_exp<- function(n1,rate1,n2,rate2) {
  exp1=rexp(n1, rate=rate1)
  f_exp1=dexp(exp1, rate = rate1)
  F_exp1=pexp(exp1, rate = rate1)
  S_exp1=1-F_exp1
  h_exp1=f_exp1/S_exp1
  
  exp2=rexp(n2, rate=rate2)
  f_exp2=dexp(exp2, rate = rate2)
  F_exp2=pexp(exp2, rate = rate2)
  S_exp2=1-F_exp2
  h_exp2=f_exp2/S_exp2
  
  plot(h_exp2)
  par(new=TRUE)
  plot(h_exp1)
  title(" compracion riesgo de exponenciales")
}
riesgo_exp(100,3,100,5)

Gamma

riesgo_gamma<- function(n1,rate1,n2,rate2) {
  g1 =rgamma(n1,rate1)
  f_g1=dgamma(g1,rate1)
  F_g1=pgamma(g1,rate1)
  S_g1=1-F_g1
  h_t1=f_g1/S_g1
  
  g2 =rgamma(n2,rate2)
  f_g2=dgamma(g2,rate2)
  F_g2=pgamma(g2,rate2)
  S_g2=1-F_g2
  h_t2=f_g2/S_g2
  
  plot(sort(h_t2),type="l")
  par(new=TRUE)
  plot(sort(h_t1),type="l")
  title(" compracion riesgo gammas")
}
riesgo_gamma(100,3,100,5)

Weibull

riesgo_weibull<- function(n1,sh1,sc1,n2,sh2,sc2) {
  w1 =rweibull(n1,shape =sh1,scale =sc1)
  f_w1=dweibull(w1,shape =sh1,scale =sc1)
  F_w1=pweibull(w1,shape =sh1,scale =sc1)
  S_w1=1-F_w1
  h_w1=f_w1/S_w1
  w2 =rweibull(n2,shape =sh2,scale =sc2)
  f_w2=dweibull(w2,shape =sh2,scale =sc2)
  F_w2=pweibull(w2,shape =sh2,scale =sc2)
  S_w2=1-F_w2
  h_w2=f_w2/S_w2
  
  plot(h_w2 %>% sort)
  par(new=TRUE)
  plot(h_w1 %>% sort)
  title(" compracion weibull")
}
riesgo_weibull(100,2,3,100,3,5)