data( "Kmenta" )
eqDemand <- consump ~ price + income
eqSupply <- consump ~ price + farmPrice + trend
system <- list( demand = eqDemand, supply = eqSupply )
## perform a SUR estimation
fitsur <- systemfit( system, "SUR", data = Kmenta )
## residuals of all equations
logLik( fitsur )
Run the code above in your browser using DataLab