added Hmisc to avoid error in snippet 30