function BR_DT_KB
%

num_samples = 1000;
samples = zeros(num_samples,2);
actual = zeros(num_samples,1);
for s = 1:num_samples
    if rand<0.8
        s_size = 500 + floor(1000*rand);
        s_hood = 1;
    else
        s_size = 1500 + floor(2500*rand);
        if rand<0.7
            s_hood = 0;
        else
            s_hood = 1;
        end
    end
    samples(s,1) = s_size;
    samples(s,2) = s_hood;
    if s_size<800
        actual(s) = 1;
    elseif s_size<1500
        actual(s) = 2;
    elseif s_hood==0
        actual(s) = 3;
    else
        actual(s) = 4;
    end
end

indexes = randperm(num_samples);
indexes = indexes(1:200);
targets = actual;
for s  = 1:length(indexes)
    index = indexes(s);
    p = randperm(4);
    if p(1)~=actual(index)
        targets(index) = p(1);
    else
        targets(index) = p(2);
    end
end

y = zeros(num_samples,1);
for s = 1:num_samples
    if samples(s,1)<800
        y(s) = 1;
    elseif samples(s,1)<1500
        y(s) = 2;
    elseif samples(s,2)==0
        y(s) = 3;
    else
        y(s) = 4;
    end
end

pc = sum(y==targets)/num_samples;
tch = 0;
