ประมาณค่าเส้นโค้งการรอดชีพด้วย survreg()

# Weibull model
wb <- ___(___(___, ___) ~ ___, data = ___)

# Retrieve survival curve from model probabilities 
surv <- seq(.99, .01, by = -.01)

# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(1))

# Create data frame with the information
surv_wb <- data.frame(time = ___, surv = ___)

# Look at first few lines of the result
head(___)