library("foreign") library("lme4") library("RLRsim") # read data setwd("c:/monographregression/speed") speed <- read.dta("data/SpeedDatingData.dta") # fit models fit1<-lmer(dec ~ gender + attr + shar + gender_attr + gender_shar + (1 | iid) , data=speed , family=binomial(link="logit")) print(fit1) fit2<-lmer(dec ~ gender + attr + shar + gender_attr + gender_shar + (1 | iid) + (0+ gender | iid) + (0+ attr | iid) + (0+ shar | iid) , data=speed , family=binomial(link="logit")) print(fit2)