% Procedura OUTPUT3SLS
% ********************
%
% YDEP ... vektor zavisle promenne
% XIND ... matice nezavisle promennych
% XPRE ... matice predeterminovaných proměnných
% BETA ... vektor regresnich koeficientu z procedury 3SLS
%  
  NOB1=NOBD-1;
  YDEP=YC(:,kk);

% SOUČET ČTVERCU 3SLS-REZIDUI
% URTR=YDEP-YVYR;
  UR3S2=0;
  i=1;
  for i = 1 : NOBD,
  UR3S2=UR3S2 + URE3(i)*URE3(i);
  end
  
% ODHAD REZIDUALNIHO ROZPTYLU 3SLS-REZIDUÍ  [SIGMA3]
  SIG3=UR3S2/NOB1;
  SIGT=sqrt(SIG3);
  
% SMERODATNE ODCHYLKY PARAMETRU [SBET] A T-STATISTIKY [TBET]
  DIAX=diag(INAM);
  PB3S=zeros(1,KPRO);
  SB3S=zeros(1,KPRO);
  TB3S=zeros(1,KPRO);

  j=1;
  for j = 1 : KPRO;
    PB3S(j)=BE3S(j);  
    SB3S(j)=SO3S(j); 
%   Stávající výpočet operuje se směr odchylkami násobenými SIGT ! ! ! :
%   SB3S(j)=SIGT*SO3S(j);
    TB3S(j)=PB3S(j)/SB3S(j);
  end
   
% ZOBRAZENI HLAVNICH VYSLEDKU KVANTITATIVNI ANALYZY 

  disp([' regresni parametry ' int2str(kk) '.rovnice : '])
  disp([PB3S])
  disp([' smérodatné odchylky parametru ' int2str(kk) '.rovnice : '])
  disp([SB3S])
  disp([' t-statistiky parametru ' int2str(kk) '.rovnice : '])
  disp([TB3S]) 

% KOEFICIENT DETERMINACE [R2]
  YPR=mean(YDEP);
  YDEV=YDEP-YPR;
  YDEV2=YDEV'*YDEV;
  R2=1.-UR3S2/YDEV2;

% DURBIN-WATSONUV KOEFICIENT [DW]
  I=1;
  for I = 1 : NOB1,
   URTD(I)=URE3(I+1)-URE3(I);
  end
  URTD2=URTD*URTD';
  DW=URTD2/UR3S2;

% von NEUMANNUV POMER [VN]
  VN=DW*NOBD/(NOBD-KPRO);

% ZOBRAZENI MEZIVYSLEDKU [vystiznost modelu]
  disp('   rezidualni chyba  koef.determinace')
  disp([   SIGT                 R2])

% ZOBRAZENI CHARAKTERISTIK PRIPADNE AUTOKORELACE REZIDUI
  disp('   Durbin-Watson  von Neumannuv pomer')
  disp([   DW                   VN])
  pause
  T=t;
  plot(T,YDEP,T,YVYR3,T,URE3)
  legend ( 'pozorovane','vyrovnane','rezidua' )
  title( '3SLS : pozorovane a vyrovnane hodnoty , rezidua rovnice')
  return
