Dear R-users, On the basis of ZeligChoice Manual I wrote several commands to create a plot of probabilities for a bivariate logit model. However, after sim() command I receive information about non-conformable arguments. Could you interpret the message? R is, as far as I know, the only statistical programme that enable to estimate not only bivariate probit, but also bivariate logit model... Regards, Thank you in advance library(ZeligChoice) dane<-read.csv("dlaer.csv", header=TRUE, sep=";") z.out2<-zelig(list(mu1=c1~age+male+company+educ_4+educ_3+educ_2,mu2=c2~age+male+eduk_4+eduk_3+eduk_2),model="blogit",data=dane) x.low2 <- setx(z.out2, male = 1) x.high2 <- setx(z.out2, male = 0) s.out2 <- sim(z.out2, x1 = x.low2, x = x.high2, num = 500) AFTER THIS COMMAND I RECEIVE THE ANSWER ABOUT AN EROR: Error in all.coef[[1]] %*% as.matrix(xm[[1]]) : non-conformable arguments plot(s.out2) summary(z.out2) coef(z.out2) vcov(z.out2)
-- View this message in context: http://r.789695.n4.nabble.com/ZeligChoice-error-tp4686285.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.