% E L I M O - U S A U T O
% ***********************
% Výpočet 4 rovnic modelu Auto/Maddala(pouze metodou)OLS
%
 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);
%
%     promenna 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
%     Transformovaná proměnná [GL]
 GPop=C.*K./(M.*Pop);    % spotřeba benzínu na obyvatele
% 
 GL  =C.*K./(M.*L);      % spotřeba benzínu na ekon.aktivního obyvatele 
%     Transformovaná proměnná [GPop]
%
%     aditivni konstanta
 const=ones(size(K)); 
%
 T=t;
%
% SPECIFIKACE 1.REGRESNI ROVNICE
%
  y1=K;
  X1=[const C M Pg Y];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  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];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  disp(' C=f(const    C-1       Pg         Pt        Pop         LL       Y)')
  disp(' *******************************************************************')
  EST=1;
  outputOLST
  pause
% 
% SPECIFIKACE 3.REGRESNI ROVNICE
%
  y1=GPop;
  X1=[const K Pg Pt Pop Y];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  disp(' GPop=f(const     K        Pg         Pt        Pop        Y)')
  disp(' ************************************************************')
  EST=1;
  outputOLST
  pause
%
% SPECIFIKACE 4.REGRESNI ROVNICE
%
  y1=GL;
  X1=[const C_1 Pg Pt LL Y];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  disp(' GL=f(const     C-1       Pg         Pt        LL       Y)')
  disp(' *********************************************************')
  EST=1;
  outputOLST
  pause

