# Invoca la libreria de supervivencia
library(survival)

# Pon en un objeto los datos
larynx <- read.table("C:/Mis documentos/larynx.dat",
        col.names=c("estado", "tiempo", "edad", "fecha", "estatus"))
# Crea el factor ETAPA en los datos larynx:
larynx$ETAPA <- factor(larynx$estado, labels = c("I", "II", "III", "IV"))

# Ajusta el modelo Weibull:
larynx.regwei <- survreg(Surv(tiempo, estatus)~ ETAPA + edad,
                 data=larynx, dist='weibull')

#Los coefficientes y el estimador de log(sigma) son:
coeficientes <- c(larynx.regwei$coefficients, log(larynx.regwei$scale))

# los errores estandares son:
errores <- sqrt(diag(larynx.regwei$var))

# el estadistico Wald de Ji cuadrada es:
ji.wald <- (coeficientes/errores)^2

# obten el valor p para cada parametro:
valor.p <- 1 - pchisq(ji.wald, 1)