% *********************** METODA GLS **********************
%  disp ('         ESTIMATOR GLS : Generalized Least Squares Method with AR(1)')
%  disp ('         ***********************************************************')

% VYPOCET KOEFICIENTU AUTOKORELACE 1. RADU
   UREZ=UREZOLS;
   clear UREZ1
   for I = 1:nob1,
   UREZ1(I)=UREZ(I+1);
   UREZZ(I)=UREZ(I);
   end
    
   UINV=UREZ1*UREZ1';
   UVEK=UREZ1*UREZZ';
   A1=UVEK./UINV;
%  TRANSFORMACE PROMENNYCH
    YTR(1)=YPOZ(1)*sqrt(1-A1*A1);
   for J = 1:nvar,
    XTR(1,J)=XPOZ(1,J)*sqrt(1-A1*A1);        
   end    
%      
   for I = 2:nobs,
    YTR(I)=YPOZ(I)-A1*YPOZ(I-1);
   for J = 1:nvar;
    XTR(I,J)=XPOZ(I,J)-A1*XPOZ(I-1,J);
   end     
   end
   YDEP=YTR';
   XIND=XTR;
   EST=3;
%  outputGLS
   
%   BPARGLS=BETA;
%   YVYRGLS=YVYR;
%   UREZGLS=URTR;
       
% GRAFICKE SROVNANI REZIDUI OLS, WLS a GLS
%  figure
%  plot(T,UREZOLS,T,UREZWLS,T,UREZGLS)
%  legend('rezidua OLS','rezidua WLS','rezidua GLS')
%  title('Srovnani rezidui vypoctenych pomoci OLS, WLS a GLS')
%  pause 
  return