clear all;
numIterations = 5;
maxFeatures = 4;

oneNN = zeros(1,maxFeatures-1);
threeNN = zeros(1,maxFeatures-1);
fiveNN = zeros(1,maxFeatures-1);
cent = zeros(1,maxFeatures-1);
avgcent = zeros(1,maxFeatures-1);
centparam = zeros(1,maxFeatures-1);
for w=1:maxFeatures-1
    q = w+1;
    q
    numFeatures = q;
    for numRuns = 1:numIterations
        numRuns
        run1overNsim;
        % save data in file name (whichiteration-howmanyfeatures)
        fname = [int2str(numRuns) '-' int2str(q)]
        save(fname);
        oneNN(1,w) = oneNN(1,w) + oneNNError;
        threeNN(1,w) = threeNN(1,w) + threeNNError;
        fiveNN(1,w) = fiveNN(1,w) + fiveNNError;
        cent(1,w) = cent(1,w) + centroidClassError;
        avgcent(1,w) = avgcent(1,w) + AvgcentroidClassError;
        centparam(1,w) = centparam(1,w) + centroidParamClassError;
    end
    oneNN(1,w) = oneNN(1,w)/numIterations;
    threeNN(1,w) = threeNN(1,w)/numIterations;
    fiveNN(1,w) = fiveNN(1,w)/numIterations;
    cent(1,w) = cent(1,w)/numIterations;
    avgcent(1,w) = avgcent(1,w)/numIterations;
    centparam(1,w) = centparam(1,w)/numIterations;
    save avgErrors;
end
save plotData;
w = 2:1:maxFeatures;
plot(w,oneNN,'-',w,cent,'+',w,fiveNN,'--',w,threeNN,'*',w,avgcent,'^',w,centparam, '-.');

% Details
% number of iterations = 1