data( "Kmenta" )
eqDemand <- consump ~ price + income
eqSupply <- consump ~ price + farmPrice + trend
system <- list( demand = eqDemand, supply = eqSupply )
## perform OLS of the system
fitols <- systemfit( system, data = Kmenta )
## design matrix of the entire system
model.matrix( fitols )
## design matrix of the first equation
model.matrix( fitols$eq[[ 1 ]] )
Run the code above in your browser using DataLab