data(birthwt)
# ESTIMATE RISK DIFFERENCE ASSOCIATED WITH MOTHER'S SMOKING STATUS, AGE, WEIGHT
fit <- lexpit(low~smoke,low~I(scale(age))+I(scale(lwt)),data=birthwt)
# RISK OF LOW BIRTH WEIGHT OF CHILD BORN TO MOTHER OF AVERAGE AGE, WEIGHT
ci(fit,C=list(linear=1,expit=c(1,0,0)),coef=FALSE,
average=list(linear=FALSE,expit=rep(FALSE,3)),
method="logit")
Run the code above in your browser using DataLab