fit1 <- glm(Survived ~ Age + Class + Sex, data = as.data.frame(Titanic),
weights = Freq, family = binomial("probit"))
getdata.DN(fit1)
library(survival)
fit2 <- coxph(Surv(time, status) ~ age + strata(sex) + ph.ecog, data = lung)
getdata.DN(fit2)
Run the code above in your browser using DataLab