gmmclf.m

% Author: Baihan Lin (doerlbh@gmail.com)
% Date: Jan 2020

function result = gmmclf(data,opts)

if opts.oracle == 1
    nArms = opts.nOptions;
    labls = [];
    labln = [];
else
    nArms = 1; % 1 for the new arm
    labls = ['new'];
    labln = [1];
end

y = data.y;
y_true = data.full_y;
x = data.rec;

rew = 0;
r = [];
a = [];

% kms = zeros(nArms,data.dim);
kms = x(:,1:nArms)';
kmv = ones(nArms,data.dim);
kmn = zeros(nArms,1);
idx = 1:nArms;

hbk = cell(nArms,1);
for i = 1:nArms
    hbk{i} = [];
end

for t = 1:data.t
    disp(strcat('gmm - ',num2str(t)))
    
    feat = x(:,t)';
    labl = y(t);
    labl_true = y_true(t);
    
    stillWrong = 0;
    stillCorrect = 0;
    
    nlls = [];
    for i = 1:nArms
        nll = 0;
        for j = 1:data.dim
            nll = nll + normlike([kms(i,j),kmv(i,j)],feat(j));
        end
        nlls = [nlls; nll];
    end
    [snlls,scores] = sort(nll);
    pred = scores(1);
    
    if opts.oracle == 1
        
        if ~any(labln == pred)
            if labl ~= "-1" && (isempty(labls) || ~any(labls == labl))
                labls = [labls;labl];
                labln = [labln;length(labls)];
                pred = length(labls);
            else
                stillWrong = 1;
            end
        else
            if labl ~= "-1" && (isempty(labls) || ~any(labls == labl))
                labls = [labls;labl];
                labln = [labln;length(labls)];
            end
        end
        
    else
        
        if ~any(labls == labl) && labl ~= "-1"
            nArms = nArms + 1;
            labls = [labls;labl];
            labln = [labln;nArms];
            kms = [kms;feat];
            kmv = [kmv;ones(1,data.dim)];
            kmn = [kmn; 0];
            idx = 1:nArms;
            hbk{nArms} = [];
            if pred == 1
                stillCorrect = 1;
            end
        end
    end
    
    if ~stillWrong && (stillCorrect || labls(labln == pred) == labl_true)
        rew = rew + 1;
    end
    
    r = [r;rew];
    acc = rew / t;
    a = [a;acc];
    
    if labl ~= "-1"
        assignment = labln(labls == labl);
        kmn(assignment) = kmn(assignment) + 1;
        hbk{assignment} = [hbk{assignment};feat];
        kms(assignment,:) = mean(hbk{assignment},1);
        kmv(assignment,:) = std(hbk{assignment},1);
    end
    
end
acc = rew / data.t;

result.acc = acc;
result.rew = rew;
result.a = a;
result.r = r;

end