the Matching() package by Jasjeet Sekhon does propensity score matching in a very user friendly way. (as you said you don't want to reinvent the wheel...) just feed it with the fitted values from a glm model (fitted$myglmmodel). afaik, you may additionally match on some covariates directly. HTH marc _________________________________________________________________________ In 5 Schritten zur eigenen Homepage. Jetzt Domain sichern und gestalten!
______________________________________________ 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.