data(birthwt)
# ADDITIVE MODEL FOR LOW BIRTHWEIGHT
# ESTIMATE RISK DIFFERENCE ASSOCIATED WITH MOTHER'S SMOKING STATUS AND AGE
fit <- blm(low~smoke+I(scale(age)),data=birthwt)
# ESTIMATE RISK DIFFERENCE ASSOCIATED WITH SMOKING
ci(fit,c(0,1,0),method="hessian")
Run the code above in your browser using DataLab