function [features1 features2] = simMean(mean1, mean2, p1, p2)
% Generate feature matrix based on 
% Inputs:
% 1. mean1: mean vector of class1
% 2. mean2: mean vector of class2
% 3. p1: probability vector ( P(F_i != mean_1i) = p1_i)
% 4. p2: probability vector ( P(F_i != mean_2i) = p2_i)

% Output:
% 1. features1 - feature matrix using class1
% 2. features2 - feature matrix using class2

% ASSUMPTION: length(mean1) = length(mean2)
dim = length(mean1);
for i= 1:dim
    % if rand number <= p1, feature(i) != mean1(i)
    tmp = rand();
    if(tmp <= p1(i))
        features1(i) = abs(1-mean1(i));
    else
        features1(i) = mean1(i);
    end
    
    tmp = rand();
    if(tmp <= p2(i))
        features2(i)  = abs(1-mean2(i));
    else
        features2(i) = mean2(i);
    end
end