function [data1, data2, NewVar1, NewVar2] = extr_by_varAA(data, AllVar, column)
% divides the data in a vector into 2 separate 
% ones, each corresponding to one value of a binary variable. 
% Also divides the variable matrix into 2. 
data1 = data;
data2 = data;
NewVar1 = AllVar;
NewVar2 = AllVar;
for i = length(AllVar(:,column)):-1:1
    if NewVar1(i,column) == 0
        NewVar1(i,:) = [];
        data1(i,:) = [];
    end
end
for i = length(AllVar(:,column)):-1:1
    if NewVar2(i,column) == 1
        NewVar2(i,:) = [];
        data2(i,:) = [];
    end
end
end