tmpar <- par(mar=c(3.1,3.1,0,0),pty="s")

mplot(leprosy$x,leprosy$y,group=leprosy$drug,xlab="",ylab="")
mtext("before",1,2)
mtext("after",2,2)
#abline(0,1,lty=2)
mlines(leprosy$x, predict(leprosy.fit), group = leprosy$drug,lty=1:3)

tmp <- mean(leprosy$x)
abline(v=tmp,lty=4)
points(rep(tmp,3),predict(leprosy.fit,data.frame(drug=levels(leprosy$drug),
	x=rep(tmp,3))),pch=18)
points(tapply(leprosy$x,leprosy$drug,mean),tapply(leprosy$y,leprosy$drug,mean),
	pch=5)

legend(3,23,c("A","B","C"),lty=1:3)
lsd.bar(leprosy.fit,leprosy,xpos=18,ypos=5)

par(tmpar)