/*
Lesson 16.2a: SUR Panel Data Analysis
Cost of Production for Airline Services III
*/
use gpe2;
output file = gpe\output16.3 reset;
load data[91,6] = gpe\airline.txt;
panel = data[2:91,1:2];   @ panel definition @
n=6;
t=15;

/* stacked data series, by sections */
cs = ln(data[2:91,3]);    @ log cost (stacked) @
qs = ln(data[2:91,4]);    @ log output (stacked) @
pfs = ln(data[2:91,5]);   @ log fuel price (stacked) @
lfs = data[2:91,6];       @ load factor (stacked) @
xs = qs~pfs~lfs;

call reset;
_names = {"c","q","pf","lf"};
_const=0;
one=eye(n).*.ones(t,1);

_iter=200;
call system1(cs,one*~(xs~ones(n*t,1)),n,t);

end;

#include gpe\system1.gpe;