%Calculation of Sampling Distributions

k = 100 ;

n = length(SS) ;

OUT = zeros(k,11) ;

X(:,1) = ones(n,1) ;

X(:,2:3) = SS ;

i = 1 ;

u = SIM(:,28) ;   
   
   x1 = SS(:,1) ;
   
   x2 = SS(:,2) ;
   
   y = ones(n,1) + .4*x1 +.8*x2 + u ;
   
   DAT = geo_regr(y,SS,SS) ;
   
   b = DAT{1} ;
   
   p = DAT{2} ;   
    
   e = y - XX*b ;
   
   var = (1/n)*e'*e ;
   
   %OLS calculations
   
   b0 = X\y ;   
      
   e0 = y - XX*b0 ;
   
   var0 = (1/n)*e0'*e0 ;
   
   %construct output
   
   OUT(i,1:3) = b' ;
   OUT(i,4:6) = p' ;
   OUT(i,7:9) = b0' ;
   OUT(i,10) = var0 ;
   OUT(i,11) = var ;