# NOT RUN {
charitable <- dplyr::mutate(charitable,
logdon = log(donation) - log(25))
char_form <- logdon ~ log(donparents) + log(income) +
education + religion + married + south
ml <- tobit1(char_form, data = charitable)
scls <- update(ml, method = "trimmed")
haustest(scls, ml, omit = "(Intercept)")
# }
Run the code above in your browser using DataLab