# NOT RUN {
data(heart)
heart.nb <- irls(death ~ anterior + hcabg + factor(killip),
a = 0.0001,
offset = log(heart$cases),
family = "negBinomial", link = "log",
data = heart)
# }
Run the code above in your browser using DataLab