Abe, You should probably look at the systemfit package for this, which can be used to estimate systems of simultaneous equations. http://cran.r-project.org/web/packages/systemfit/vignettes/systemfit.pdf Mit besten Grüßen / With kind regards Hans Radtke