% E L I M O - U S A U T O
% ***********************
  close all
  clear all 
  datamodel=xlsread('dataUSAutoGR');
  format 
  USAauto=datamodel;
  N=max(size(USAauto(:,:)));
%
% PROMENNE EKONOMETRICKEHO MODELU
%
%     promenne bez zpozdeni [t]
  G   =USAauto(1:N,2);  % Celková spotřeba benzínu v USA
  Pg  =USAauto(1:N,3);  % Cenový index benzínu
  Y   =USAauto(1:N,4);  % Disponibilní příjem na hlavu
  PNC =USAauto(1:N,5);  % Cenový index nových aut                
  PUC =USAauto(1:N,6);  % Cenový index ojetých aut
  PPT =USAauto(1:N,7);  % Cenový index veřejné dopravy
  PD  =USAauto(1:N,8);  % agregátní index trvanlivých komodit
  PN  =USAauto(1:N,9);  % agregátní index netrvanlivých komodit
  PS  =USAauto(1:N,10); % agregátní index služeb
  Pop =USAauto(1:N,11); % počet obyvatel v USA
  t   =USAauto(1:N,12); % pořadová časová proměnná
%
%     podílové a transformované proměnné [GP,YP,PDNS]
%
  GP  =G./Pop ;             % spotřeba benzínu na obyvatele
  YP  =GP.*Pg ;             % výdaje na benzín na obyvatele
  PDNS=(PD+PN+PS)./3 ;      % součtovyY agregátní cenový index souhrnny
% 
%    aditivni konstanta
  const=ones(size(G)); 
%
  T=t;
%
% SPECIFIKACE 1.REGRESNI ROVNICE
%
  y1=G;
  X1=[const Y  Pg PNC PUC  PPT  Pop];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  nobs=length(y1);
  disp('  G=f(const   Y     Pg     PNC     PUC      PPT      Pop )')
  disp('  ********************************************************')
  EST=1;
  outputOLST
  pause
  EST=2;
  outputWLS
  pause
  EST=3;
  outputGLS
  pause
% 
% SPECIFIKACE 2.REGRESNI ROVNICE
%
  y1=GP;
  X1=[const Y  Pg  PNC PD PN PS ];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  nobs=length(y1);
  disp(' GP=f(const    Y     Pg     PNC      PD       PN       PS)')
  disp(' *********************************************************')
  EST=1;
  outputOLST
  EST=2;
  outputWLS
  EST=3;
  outputGLS
  pause
% 
% SPECIFIKACE 3.REGRESNI ROVNICE
%
  y1=YP;
  X1=[const PUC PPT PD PN PS];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  nobs=length(y1);
  disp(' YP=f(const    Y     PUC      PPT       PD       PN       PS )')
  disp(' *************************************************************')
  EST=1;
  outputOLST
  EST=2;
  outputWLS
  EST=3;
  outputGLS
  pause
%
% SPECIFIKACE 4.REGRESNI ROVNICE
%
  y1=YP;
  X1=[const  Y  PNC  PUC  PPT  PDNS  t];
  [nobs nvar]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  nobs=length(y1);
  disp(' YP=f(const    Y    PNC     PUC      PPT       PDNS     t)')
  disp(' *********************************************************')
  EST=1;
  outputOLST
  EST=2;
  outputWLS
  EST=3;
  outputGLS
  pause

