clear all

namef1='lin_inv_prb_data1.mat';
load(namef1)

betatrue=beta_0;

betay=inv(X'*X)*X'*y;
betav=inv(X'*X)*X'*v;
betaw=inv(X'*X)*X'*w;

figure;
set(0,'DefaultAxesFontSize',18)
title('True \beta_0 versus three estimates')

subplot(311)
    plot(betay,betatrue,'o')
    xlabel('Estimate \beta_y')
    ylabel('True \beta_0')

subplot(312)
    plot(betav,betatrue,'o')
    xlabel('Estimate \beta_v')
    ylabel('True \beta_0')

subplot(313)
    plot(betaw,betatrue,'o')
    xlabel('Estimate \beta_w')
    ylabel('True \beta_0')
  
    
    
clear all    
namef2='lin_inv_prb_data2.mat';

load(namef2)

betatrue=beta_0;

betay=inv(X'*X)*X'*y;
betav=inv(X'*X)*X'*v;
betaw=inv(X'*X)*X'*w;

figure;
set(0,'DefaultAxesFontSize',18)
title('True \beta_0 versus three estimates')


subplot(311)
    plot(betay,betatrue,'o')
    xlabel('Estimate \beta_y')
    ylabel('True \beta_0')

subplot(312)
    plot(betav,betatrue,'o')
    xlabel('Estimate \beta_v')
    ylabel('True \beta_0')

subplot(313)
    plot(betaw,betatrue,'o')
    xlabel('Estimate \beta_w')
    ylabel('True \beta_0')