% E L I M O - U S A U T O
% ***********************
 dataauto=xlsread('datausagasm');
 close all
 format 
 USAauto=dataauto;
 N=max(size(USAauto(:,:)));
%
% PROMENNE EKONOMETRICKEHO MODELU
%
%     promenne bez zpozdeni [t]
 K   =USAauto(2:N,2);  % Vzdálenost ujetá osobním autem za rok    [mile]
 C   =USAauto(2:N,3);  % Počet automobilů v USA                  [mil.ks]
 M   =USAauto(2:N,4);  % Vzdálenost ujetá autem na galon benzínu  [mile]
 Pg  =USAauto(2:N,5);  % Cena benzínu                             [index 1953]
 Pt  =USAauto(2:N,6);  % Cena veřejné dopravy                    [index 1967]
 Pop =USAauto(2:N,7);  % Obyvatelsto USA                          [mil.osob]
 L   =USAauto(2:N,8);  % Zaměstnanost USA                        [mil.osob]
 Y   =USAauto(2:N,9);  % Příjem na obyvatele USA                 [index 1958]
 t   =USAauto(2:N,1);
 %
%     promenne se zpozdenim 1 [t-1]
 C_1 =USAauto(1:N-1,3);  % C(t-1)
%     podílová proměnná [LL]
 LL  =L./Pop;          % podíl ekonomicky aktivního obyvatelstva
%     aditivni konstanta
 const=ones(size(K)); 
%
 T=t;
%
% SPECIFIKACE 1.REGRESNI ROVNICE
%
  y1=K;
  X1=[const C M Pg Y];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  NOBD=length(y1);
  KPRO=kp;
  disp('  K=f(const      C        M        Pg        Y)')
  disp('  *********************************************')
  EST=1;
  outputOLST
  pause
% 
% SPECIFIKACE 2.REGRESNI ROVNICE
%
  y1=C;
  X1=[const C_1 Pg Pt Pop LL Y];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  NOBD=length(y1);
  KPRO=kp;
  disp(' C=f(const    C-1       Pg         Pt        Pop         LL       Y)')
  disp(' *******************************************************************')
  EST=1;
  outputOLST
  pause
%