fit_binary.m

function results = fit_binary(data)
    
    % Fit two-alternative DDM model to data from a value-based choice task.
    %
    % USAGE: results = fit_binary(data)
    %
    % INPUTS:
    %   data - [S x 1] data structure, where S is the number of subjects; see likfun_binary for more details
    %
    % OUTPUTS:
    %   results - see mfit_optimize for more details
    %
    % Sam Gershman, Dec 2016
    
    % create parameter structure
    
    % drift rate differential action value weight
    param(1).name = 'b';
    param(1).logpdf = @(x) 0;  % uniorm prior
    param(1).lb = -20; % lower bound
    param(1).ub = 20;   % upper bound
    
    % decision threshold
    param(2).name = 'a';
    param(2).logpdf = @(x) 0;
    param(2).lb = 1e-3;
    param(2).ub = 40;
    
    % non-decision time
    param(3).name = 'T';
    param(3).logpdf = @(x) 0;
    param(3).lb = 0;
    param(3).ub = 1;
    
    % OPTIONAL: drift bias
    param(4).name = 'b0';
    param(4).logpdf = @(x) 0;
    param(4).lb = -20;
    param(4).ub = 20;
    
    % fit model
    f = @(x,data) likfun_binary(x,data);    % log-likelihood function
    results = mfit_optimize(f,param,data);