If I wanted to fit a logit model and account for clustering of observations, I would do something like:
library(Design) f <- lrm(Y1 ~ X1 + X2, x=TRUE, y=TRUE, data=d) g <- robcov(f, d$st.year) What would I do if I wanted to do the same thing with a probit model? ?robcov says the input model must come from the Design package, but the Design package appears not to do probit? Thanks very much! ____________________________________________________________ FREE 3D MARINE AQUARIUM SCREENSAVER - Watch dolphins, sharks & orcas on your desktop! Check it out at http://www.inbox.com/marineaquarium ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.