VP_RPE_2020/MatlabScripts/2_Modeling/threeOutcomes/fit_model/RWmodels_RPE/ott_RW_RPE_mean.m

function [LH, probSpike, V, mean_predictedSpikes, RPE] = ott_RW_RPE_mean(startValues, spikeCounts, rewards, timeLocked)


probSpike = poisspdf(spikeCounts, mean(spikeCounts)); % mask rateParam to exclude trials where the animal didn't lick fast enough

mean_predictedSpikes = mean(spikeCounts);

if any(isinf(log(probSpike)))
    LH = 1e9;
else
    LH = -1 * sum(log(probSpike));
end

V = NaN;
RPE = NaN;