t1=2.4390e-02; t2=2.4994e-02; t3=1.2748e-01; t4=9.3777e-02;

% ********* lc mine ************* 5: filter 1 s=1, 6: filter 1 s=2
%	7: filter 1 s=4
x=[[	2.8997e-02   3.0366e-02   1.1713e-01   9.1178e-02   6.3987e-04 ...
	7.6710e-04   2.4986e-03   2.2285e-04]
[	2.8865e-02   3.0277e-02   1.1694e-01   9.1224e-02   6.4969e-04 ...
	8.1283e-04   2.7745e-03   2.1529e-04]
[	2.8877e-02   3.0395e-02   1.1673e-01   9.1205e-02   6.9958e-04 ...
	8.8605e-04   3.0168e-03   2.1496e-04]
[	2.7545e-02   2.8878e-02   1.1785e-01   9.2329e-02   7.7841e-04 ...
	9.9055e-04   3.2918e-03   2.3216e-04]
[	3.0030e-02   3.1634e-02   1.1645e-01   9.0015e-02   6.1790e-04 ...
	7.2806e-04   2.4160e-03   2.1693e-04]
[	3.2480e-02   3.4635e-02   1.1475e-01   8.7492e-02   5.8475e-04 ...
	6.5753e-04   2.2201e-03   2.0611e-04]
[	3.8946e-02   4.2221e-02   1.0983e-01   8.1270e-02   5.1498e-04 ...
	5.3217e-04   1.7254e-03   1.8445e-04]];
x(:,1)=x(:,1)/t1*100; x(:,5)=x(:,5)/t1*100;
x(:,2)=x(:,2)/t2*100; x(:,6)=x(:,6)/t2*100;

sel=[1 5 6];
figure(3),p1=plot(x(sel,1)-100,x(sel,5),'s'); hold on
        plot(x(sel,1)-100,x(sel,5)); hold off; plerr(x(sel,1)-100,x(sel,5))
figure(4),p2=plot(x(sel,2)-100,x(sel,6),'s'); hold on
        plot(x(sel,2)-100,x(sel,6)); hold off; plerr(x(sel,2)-100,x(sel,6))
%subplot(223),p3=plot(abs(x(sel,3)-t3),x(sel,7)-2.9e-4,'s'); hold on
%        plot(abs(x(sel,3)-t3),x(sel,7)-2.6e-4); hold off
%	plerr(abs(x(sel,3)-t3),x(sel,7)-2.6e-4)
%subplot(224),p4=plot(abs(x(sel,4)-t4),x(sel,8),'s'); hold on
%        plot(abs(x(sel,4)-t4),x(sel,8)); hold off
%	plerr(abs(x(sel,4)-t4),x(sel,8))

% ********* huber delta=0.004 *************
x=[[	6.0159e-02   6.1436e-02   1.5352e-01   1.1846e-01   2.8677e-03 ...
	3.1587e-03   7.6374e-03   8.6796e-04]
[	2.6922e-02   2.7578e-02   1.2302e-01   9.4014e-02   9.2443e-04 ...
	1.0940e-03   3.3598e-03   3.2921e-04]
[	2.6827e-02   2.7649e-02   1.1857e-01   9.3393e-02   7.8877e-04 ...
	9.4385e-04   2.7944e-03   2.9861e-04]
[	2.8834e-02   2.9832e-02   1.1455e-01   9.2522e-02   7.5390e-04 ...
	8.5511e-04   1.9315e-03   2.7846e-04]
[	3.1864e-02   3.3290e-02   1.1086e-01   9.0958e-02   7.0709e-04 ...
	8.0397e-04   1.4905e-03   2.6513e-04]
[	3.6380e-02   3.8614e-02   1.0789e-01   8.8395e-02   6.2722e-04 ...
	7.0940e-04   1.1342e-03   2.3775e-04]
[	4.6761e-02   5.0766e-02   1.0527e-01   8.1663e-02   3.6911e-04 ...
	4.2781e-04   4.4311e-04   1.4414e-04]
[	5.2659e-02   5.6699e-02   1.0185e-01   7.4692e-02   1.1541e-04 ...
	1.3456e-04   9.9207e-05   5.9333e-05]];
x(:,1)=x(:,1)/t1*100; x(:,5)=x(:,5)/t1*100;
x(:,2)=x(:,2)/t2*100; x(:,6)=x(:,6)/t2*100;

sel=[2:7];
figure(3),hold on; q1=plot(x(sel,1)-100,x(sel,5),'o');
        plot(x(sel,1)-100,x(sel,5),'--'); hold off; plerr(x(sel,1)-100,x(sel,5))
figure(4),hold on; q2=plot(x(sel,2)-100,x(sel,6),'o');
        plot(x(sel,2)-100,x(sel,6),'--'); hold off; plerr(x(sel,2)-100,x(sel,6))	
%subplot(223),hold on; q3=plot(abs(x(sel,3)-t3),x(sel,7),'o');
%        plot(abs(x(sel,3)-t3),x(sel,7),'--'); hold off
%	plerr(abs(x(sel,3)-t3),x(sel,7))
%subplot(224),hold on; q4=plot(abs(x(sel,4)-t4),x(sel,8),'o');
%        plot(abs(x(sel,4)-t4),x(sel,8),'--'); hold off
%	plerr(abs(x(sel,4)-t4),x(sel,8))

% ********* huber delta=0.002 *************
x=[[	3.2663e-02   3.3415e-02   1.2677e-01   9.6873e-02   1.3443e-03 ...
	1.5066e-03   3.9583e-03   4.0885e-04]
[	2.6807e-02   2.7471e-02   1.2301e-01   9.4045e-02   9.3481e-04 ...
	1.0819e-03   3.3953e-03   3.2987e-04]
[	2.8249e-02   2.9143e-02   1.1200e-01   9.2856e-02   7.7336e-04 ...
	8.8753e-04   2.3134e-03   2.8368e-04]
[	3.1215e-02   3.2507e-02   1.0665e-01   9.1361e-02   7.2662e-04 ...
	8.4508e-04   1.2894e-03   2.6894e-04]
[	4.3792e-02   4.7506e-02   9.7810e-02   8.4315e-02   6.4300e-04 ...
	6.8716e-04   7.3977e-04   2.2103e-04]];
x(:,1)=x(:,1)/t1*100; x(:,5)=x(:,5)/t1*100;
x(:,2)=x(:,2)/t2*100; x(:,6)=x(:,6)/t2*100;

sel=[1:5];
figure(3),hold on; r1=plot(x(sel,1)-100,x(sel,5),'+');
        plot(x(sel,1)-100,x(sel,5),'-.'); hold off; %plerr(x(sel,1)-100,x(sel,5))
figure(4),hold on; r2=plot(x(sel,2)-100,x(sel,6),'+');
        plot(x(sel,2)-100,x(sel,6),'-.'); hold off; %plerr(x(sel,2)-100,x(sel,6))
%subplot(223),hold on; r3=plot(abs(x(sel,3)-t3),x(sel,7),'+');
%        plot(abs(x(sel,3)-t3),x(sel,7),'-.'); hold off
%	plerr(abs(x(sel,3)-t3),x(sel,7))
%subplot(224),hold on; r4=plot(abs(x(sel,4)-t4),x(sel,8),'+');
%        plot(abs(x(sel,4)-t4),x(sel,8),'-.'); hold off
%	plerr(abs(x(sel,4)-t4),x(sel,8))

% ********* huber delta=0.008 *************
x=[[	2.7053e-02   2.7865e-02   1.1985e-01   9.3241e-02   7.8497e-04 ...
	9.3605e-04   2.7914e-03   3.0057e-04]
[	2.9373e-02   3.0384e-02   1.1696e-01   9.2177e-02   7.5715e-04 ...
	8.3748e-04   2.2303e-03   2.8270e-04]
[	3.7100e-02   3.9204e-02   1.1085e-01   8.7800e-02   5.9820e-04 ...
	6.3890e-04   1.4920e-03   2.3451e-04]
[	4.2285e-02   4.5286e-02   1.0672e-01   8.4129e-02   5.2043e-04 ...
	5.5363e-04   1.1676e-03   2.0750e-04]
[	5.4811e-02   5.8575e-02   9.5393e-02   7.3729e-02   4.1500e-04 ...
	4.5143e-04   6.4406e-04   1.4590e-04]
[	5.9453e-02   6.1496e-02   8.1321e-02   6.1429e-02   1.5634e-04 ...
	1.6787e-04   1.6228e-04   6.0026e-05]];
x(:,1)=x(:,1)/t1*100; x(:,5)=x(:,5)/t1*100;
x(:,2)=x(:,2)/t2*100; x(:,6)=x(:,6)/t2*100;

sel=[1:5];
figure(3),hold on; s1=plot(x(sel,1)-100,x(sel,5),'x');
        plot(x(sel,1)-100,x(sel,5),':'); hold off; %plerr(x(sel,1)-100,x(sel,5))
figure(4),hold on; s2=plot(x(sel,2)-100,x(sel,6),'x');
        plot(x(sel,2)-100,x(sel,6),':'); hold off; %plerr(x(sel,2)-100,x(sel,6))
%subplot(223),hold on; s3=plot(abs(x(sel,3)-t3),x(sel,7),'x');
%        plot(abs(x(sel,3)-t3),x(sel,7),':'); hold off
%        plerr(abs(x(sel,3)-t3),x(sel,7))
%subplot(224),hold on; s4=plot(abs(x(sel,4)-t4),x(sel,8),'x');
%        plot(abs(x(sel,4)-t4),x(sel,8),':'); hold off
%        plerr(abs(x(sel,4)-t4),x(sel,8))

figure(3),legend([q1 r1 s1 p1],'Huber:  \delta=0.004',...
        'Huber:  \delta=0.002','Huber:  \delta=0.008','Proposed method');
axis([0 80 1 6])
%xlabel('normalized bias (%)'); ylabel('normalized standard deviation (%)');
figure(4),legend([q2 r2 s2 p2],'Huber:  \delta=0.004',...
        'Huber:  \delta=0.002','Huber:  \delta=0.008','Proposed method');
axis([0 80 1 6])
