use “lab3.dta”, clear glm death i.hlamat c.age_don c.age i.txtype c.year, family(binomial) link(log) eform logistic death i.hlamat c.age_don c.age i.txtype c.year margins txtype, post nlcom log(_b[1.txtype])-log(_b[0.txtype]) logistic death i.hlamat c.age_don c.age i.txtype c.year replace txtype = 1 predict p1, p summ p1 replace txtype = 0 predict p0, p summ p0