Many statistical analyses are based on models containing systems of structurally related equations. In cases where cross-equation disturbances are correlated, full information methods are required (Zellner, 1962). If exogenous variables are stochastically dependent on the disturbances in the system, then instrumental variable estimation methods should be used (Zellner and Theil, 1962) The package systemfit provides the capability to estimate systems of linear equations within the R programming environment