function RD_err_stats(locs,nei)
%

 [a,b,step,at,err]...
     = RD_linear_nei(60,0.00000001,1500,4,4,0.25,0.0625,0.05,1,nei,0.15);
figure(2)
clf
RD_display_cells(a,locs,nei,0.15);
nei_bin = nei<0.15&nei>0;
a_15_m = mean(sum(nei_bin));

[a,b,step,at,err]...
    = RD_linear_nei(60,0.00000001,1500,4,4,0.25,0.0625,0.05,1,nei,0.21);
figure(2)
clf
RD_display_cells(a,locs,nei,0.21);
nei_bin = nei<0.21&nei>0;
a_21_m = mean(sum(nei_bin));

[a,b,step,at,err]...
    = RD_linear_nei(60,0.00000001,1500,4,4,0.25,0.0625,0.05,1,nei,0.25);
figure(2)
clf
RD_display_cells(a,locs,nei,0.25);
nei_bin = nei<0.25&nei>0;
a_25_m = mean(sum(nei_bin));

[a,b,step,at,err]...
    = RD_linear_nei(60,0.00000001,1500,4,4,0.25,0.0625,0.05,1,nei,0.31);
figure(2)
clf
RD_display_cells(a,locs,nei,0.31);
nei_bin = nei<0.31&nei>0;
a_31_m = mean(sum(nei_bin));

[a,b,step,at,err]...
    = RD_linear_nei(60,0.00000001,1500,4,4,0.25,0.0625,0.05,1,nei,0.35);
figure(2)
clf
RD_display_cells(a,locs,nei,0.35);
nei_bin = nei<0.35&nei>0;
a_35_m = mean(sum(nei_bin));

tch = 0;
