% E L I M O - U S F O O D
% *******************************
  datafood=xlsread('datausfood');
  close all
  format 
  usaFOOD=datafood;
  N=max(size(usaFOOD(:,:)));
%
% PROMENNE EKONOMETRICKEHO MODELU
%
 QD   =usaFOOD(1:N,2);  % Spotreba potravin na hlavu                           [index]
 PD   =usaFOOD(1:N,3);  % Maloobchodni ceny potravin/index zivotnich nákladu   [index]
 Y    =usaFOOD(1:N,4);  % disponibilní duchod/index zivotních nákladu          [index]
 QS   =usaFOOD(1:N,5);  % Produkce potravin na hlavu                           [index]
 PS   =usaFOOD(1:N,6);  % Farmárské ceny potravin/index zivotních nákladu      [index]
 t    =usaFOOD(1:N,7);  % Promenná vyjadrující cas
%     aditivni konstanta
  const=ones(size(QD)); 
%     casova promenna
  T=t;
  for I = 1:N,
  T2T(I)=T(I)*T(I);
  T3T(I)=T(I)*T(I)*T(I);
  end
  T2=T2T'; 
  T3=T3T'; 
%     logaritmované veličiny
  LQD=log(QD);
  LT=log(T);
  LPD=log(PD);
  LY=log(Y);
  LQS=log(QS);
  LPS=log(PS);
%  
%  SPECIFIKACE 1.REGRESNI ROVNICE (Lineární poptávková funkce s cenou a příjmem)
%
  y1=QD;
  X1=[const PD Y];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  NOBD=length(y1);
  NOB1=NOBD-1;
  KPRO=kp;
  disp(' QD=f(const    PD        Y     )')
  disp(' *******************************')
  outputOLST
  pause
%
%  SPECIFIKACE 2.REGRESNI ROVNICE (Lineární poptávková funkce s cenou, příjmem a nabídkou )
%
  y1=QD;
  X1=[const PD Y QS];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' QD=f(const    PD        Y        QS   )')
  disp(' ***************************************')
  outputOLST
  pause
%
% SPECIFIKACE 3.REGRESNI ROVNICE ( Cobb-Douglasova poptávková funkce s cenou, příjmem a nabídkou )
%
  y1=LQD;
  X1=[const   LPD   LY    LQS] ;
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' lnQD=f(const   lnPD      lnY       lnQS   )')
  disp(' ********************************************')
  outputOLST
  pause
  %
% SPECIFIKACE 4.REGRESNI ROVNICE (Lineární poptávková funkce s kubickou parabolou)
%
  y1=QD;
  X1=[const PD Y QS T T2 T3];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' QD=f(const      PD        Y        QS         T        T2        T3 )')
  disp(' ********************************************************************')
  outputOLST
  pause
%
% SPECIFIKACE 5.REGRESNI ROVNICE (Lineární nabídková funkce s cenou a příjmem )
%
  y1=QS;
  X1=[const PS Y];
  [n kp]=size(X1);
  KPRO=kp;
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' QS=f(const       PS       Y     )')
  disp(' *********************************')
  outputOLST
  pause
%
% SPECIFIKACE 6.REGRESNI ROVNICE ( Lineární nabídková funkce s cenou, příjmem a nabídkou )
%
  y1=QS;
  X1=[const PS Y QD ];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' QS=f(const       PS       Y        QD       )')
  disp(' *********************************************')
  outputOLST
  pause
%
% SPECIFIKACE 7.REGRESNI ROVNICE ( Cobb-Douglasova nabídková funkce s cenou, příjmem a poptávkou )
%
  y1=LQS;
  X1=[  const   LPS    LY    LQD] ;
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' lnQS=f(const  lnPS      lnY       lnQD      )')
  disp(' *********************************************')
  outputOLST
  pause
%  
% SPECIFIKACE 8.REGRESNI ROVNICE (Lineární nabídková funkce s kubickou parabolou)
%
  y1=QS;
  X1=[const PS Y QD T T2 T3];
  [n kp]=size(X1);
  YPOZ = y1;
  XPOZ = X1;
  KPRO=kp;
  disp(' QS=f(const       PS         Y        QD         T         T2        T3 )')
  disp(' ************************************************************************')
  outputOLST
  pause