make_stopos_file.m

function [] = make_stopos_file()

% Code to fit the history-dependent drift diffusion models as described in
% Urai AE, de Gee JW, Tsetsos K, Donner TH (2019) Choice history biases subsequent evidence accumulation. eLife, in press.
%
% MIT License
% Copyright (c) Anne Urai, 2019
% anne.urai@gmail.com

% ============================================ #
% parameter file for HDDM
% ============================================ #s

nsmp       = [5000];
datasets   = [0]; % dataset number, 0-5 for all main ones
models     = [0:3]; % the nr of the models
nrTraces   = 30; % nr of chains, 15 cores/node (so make sure this is a multiple of 15)

alldat = [];
for n = nsmp,
    for b = models,
        for a = datasets,
            for c = 0:nrTraces-1, % put all chains of same model together on a node
                alldat = [alldat; a b c n];
            end
        end
    end
end

% write to a file
dlmwrite('hddmparams', alldat, 'delimiter', ' ');
%size(alldat)

fprintf('to submit all these jobs, need an array of %d jobs \n', length(datasets)*length(models)*2);

% PPC
alldat = [];
for a = datasets,
    for v = models,
        alldat = [alldat; a v];
    end
end

% write to a file
dlmwrite('hddmparams_PPC', alldat, 'delimiter', ' ');
% size(alldat)
%% for PPC & chiSquare