SunkCostModelsAndData2022/Test_DecisionTimeWZ.m at main · adredish/SunkCostModelsAndData2022 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
function SCout = Test_DecisionTimeWZ(decisionTimeToTest, qTdecayThroughDT, varargin)

% KPTest Delayed Start
% if delay start by Xs, what happens?

SCout.decisionTimeWZ = 0:9;
nBoot = 3;
process_varargin(varargin);

nD = length(SCout.decisionTimeWZ);

SCout.sc = nan(nBoot, nD, 30);
SCout.scSum = nan(nBoot, nD);

assert(islogical(qTdecayThroughDT));
switch (decisionTimeToTest)
    case 'decisionTimeWZ_sigmaN'
        SCout.T = '\sigma_N=0 for a few s';
        SCout.Lstr = '\sigma_N=0 0 for %s s';
    case 'decisionTimeWZ_hardset'
        SCout.T = 'cannot quit for a few s';
        SCout.Lstr = 'cannot quit for %s s';
end
if qTdecayThroughDT
    SCout.T = [SCout.T, '; T_{WZ} decays at WZ entry (through DZ)'];
else
    SCout.T = [SCout.T, '; T_{WZ} decays only after DZ'];
end

hWait = waitbar(0);
for iB = 1:nBoot
    waitbar(0, hWait, sprintf('%d/%d', iB, nBoot));
    for iD = 1:nD
        waitbar(iD/nD, hWait);

        R = GenerateSimulation(decisionTimeToTest, SCout.decisionTimeWZ(iD), ...
            'offerEnterFunction', @(x)ones(size(x)), ...
            'qTdecayThroughDT', qTdecayThroughDT);
        pEarn = Calculate_pEarn(R);
        [S,C] = Calculate_SunkCostMeans(pEarn,R);
        SCout.sc(iB, iD, :) = S(:,1) - C(:,1);
        SCout.scSum(iB, iD) = nansum(S(:,1) - C(:,1));

        if iB==1
            Show_SunkCostBubble(pEarn, R);
            title(sprintf('%d s', iD-1));
            axis square; xticks([0 30]); yticks([0 1]);

            FigureLayout('scaling', 5);
        end
    end
end