Dear Subscribers, On the basis of syntax from the ZeligChoice Manual I wanted to receive s.out results after BLOGIT (e.g. qi$ev, qi$fd). However, after sim() command I get information about non-conformable arguments. Could you interpret the message? Is it a problem with matrix dimensions? Actually, I wanted to compare results of BIPROBIT computed in Stata with BLOGIT (bivariate logit is not available in Stata). Unfortunately, I do not use R often enough :( Regards, Thank you in advance library(ZeligChoice) data<-read.csv("data.csv", header=TRUE, sep=";") z.out2<-zelig(list(mu1=y1~age+gender+company+educ_4+educ_3+educ_2+qual_3+qual_2+qual_1,mu2=y2~age+gender+educ_4+educ_3+educ_2),model="blogit",data=data) x.hi <- setx(z.out2, gender = 1) x.lo <- setx(z.out2, gender = 0) s.out2 <- sim(z.out2, x = x.lo, x1 = x.hi) AFTER THIS COMMAND I GET AN EROR: Error in all.coef[[1]] %*% as.matrix(xm[[1]]) : non-conformable arguments