Este código sirve para observar si dada una colección de datos, se puede ajustar una distribución paramétrica. Como por ejemplo una distribución exponencial o una weibull La librería “survival” contiene varias funciones para trabajar con datos de tupo supervivencia, es decir, tiempos de falla que pueden ser censurados por la derecha, por la izquierda o por intervalos y truncados.

require("survival")
## Loading required package: survival

Se genera una función de ditribución con ciertos tiempos de fallo y se grafica, des esta forma si esta se parece a alguna distribución paramétrica entonces se utiliza esa para el estudio.

Se hace la prueba con una distribución weibull y una exponencial

#Grafica Q-Q
p<-c(.01,.1,.2,.3,.4,.5,.6,.7,.8,.9)
Qw<-qweibull(p, shape=1, scale = 1/0.06, lower.tail = TRUE, log.p = FALSE)
Qe<-quantile(tiempo,probs = p)
plot(Qw,Qe)
abline(a=0,b=1)

#Grafica de linearizacion
#Exponencial
plot(f3$time,log(f3$surv),type = "l",main = "Diagnóstico modelo exponencial",xlab = "tiempo", ylab="Log de supervivencia estimada")

#Weibull
plot(log(f3$time),log(-log(f3$surv)),type = "l",main="Diagnóstico modelo Weibull")