YA_full_optiontaskanalysis

Young adults - options task

Participants

The young adult analysis recruited 88 participants.

1 experimental dataset excluded as did experiment twice.

1 participant had incomplete questionnaire data for Dimensional Apathy Scale - not included in these analyses.

3 rows with RT <0.3s removed

1 participant had low accuracy but retained as seemingly normal RT

Analysis notes

New analysis uses a combined variable called ‘Difficulty’ which is the mean of the inverse value difference between 2 best options and the effort level of the best option. This is because…

  1. Value difference and bestoption cannot be considered independent of eachother and therefore violates assumptions for linear models and introduces problems of collinearity: The value of ‘bestoption’ is dependent on the value difference (i.e. one cannot have bestoption of 5 and value difference of 4) and the two values are negatively correlated. This introduces problems with interpretation of parameters - especially in the context of two- and three-way interactions between (inverse) value difference and bestoption, e.g.:

Here, the estimated marginal means from an accuracy model finding a significant two-way interaction between value difference and bestoption suggest a steeper decline in accuracy in ‘hard’ (high inverse value difference) choices with high bestoption values - however this is a non-existent condition in the data.

  1. However both value difference and bestoption are important attributes of choice difficulty. Both correlate with RT, accuracy and subjective choice difficulty. Although in the design choice difficulty was manipulated primarily by value difference, bestoption is an even stronger predictor (particularly of RT) and is necessary to fully capture the essence of choice difficulty.

  2. I have previously explored various ways to incorporate both bestoption and value difference, including inclusion of bestoption as a control variable to hold at its mean value, or as an interacting variable. However bestoption captures much of the variance and, in drift diffusion models, leads to difficulty in intepreting models with multiple terms and three-way interactions, produces divergences and inconsistent model findings (and weird effects, such as bimodal peaks for non-decision time).

Therefore I have found the most suitable way to address this is to combine both bestoption and (inverse) value difference into a new variable called difficulty - which is the mean of both (i.e. a linear combination, which appears a valid way to approach this and gets round all problems of model violations and interpretability). This new variable also correlates extremely well with RT, accuracy and subjective difficulty ratings:

Analyses

For linear models of RT, I could either use a GLMM with gamma log link to account for skewed data, or log transform RTs in LMM. The advantage of the latter is you can use buildmer to refine random effects, as I haven’t for the life of me been able to find a way it can apply the gamma log link function.

The format of mixed models is therefore:

Mixed models

Effort task only (n=86)

LMM: logRT ~ Difficulty * Options + (1 +Difficulty * Options | ID)

GLMM (binomial): accuracy ~ Difficulty * Options + (1 +Difficulty * Options | ID)

GLMM (gamma link): Rating ~ Difficulty * Options + (1 +Difficulty * Options | ID)

GLMM (negative binomial): Timeouts ~ Difficulty * Options + (1 + Difficulty * Options | ID)

LMM : Force (AUC) ~ Option worked * Apathy

Models of above also including main, two- and three way interactions with DAS total, DAS executive, DAS auto-activation, DAS emotional-affective.

Effort and reward task (n=62)

LMM: logRT ~ Difficulty * Options * Task + (1 +Difficulty * Options * Task | ID)

GLMM (binomial): accuracy ~ Difficulty * Options * Task + (1 +Difficulty * Options * Task | ID)

GLMM (gamma link): Rating ~ Difficulty * Options + (1 +Difficulty * Options * Task | ID)

GLMM (negative binomial): Timeouts ~ Difficulty * Options + (1 + Difficulty * Options * Task | ID)

Drift diffusion modelling

(a=threshold, v=drift, z=bias, t=non-decision-time)

Start with most complex model:

a ~ Difficulty + Options + Difficulty:Options + (1 +Difficulty + Options + Difficulty:Options | ID)

v ~ Difficulty + Options + Difficulty:Options + (1 +Difficulty + Options + Difficulty:Options | ID)

z (group-level)

t (group-level)

Refine sequentially on basis of variance contribution and test with loo-cv.

Also try models with collapsing bounds.

Results

Reaction times

Base models

The base model of reaction times in 87 participants doing effort task:

  RT
Predictors Estimates CI p
(Intercept) 1.00 0.95 – 1.04 0.871
Difficulty sc 1.16 1.14 – 1.18 <0.001
noopts s2z1 0.92 0.91 – 0.93 <0.001
Difficulty sc × noopts
s2z1
0.99 0.98 – 0.99 <0.001
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.Difficulty_sc 0.00
τ11 shortid.noopts_s2z1 0.00
τ11 shortid.Difficulty_sc:noopts_s2z1 0.00
ρ01 -0.38
0.17
0.30
ICC 0.10
N shortid 87
Observations 13852
Marginal R2 / Conditional R2 0.244 / 0.320

The base model of reaction times incorporating task as dummy main and interacting variable, in 61 participants who did both effort and reward task: -Faster in reward task -Faster particularly in easy choices in reward vs effort task -Fsster particularly in two option choices in reward vs effort task

  RT
Predictors Estimates CI p
(Intercept) 1.02 0.97 – 1.08 0.447
Difficulty sc 1.15 1.13 – 1.17 <0.001
noopts s2z1 0.92 0.92 – 0.93 <0.001
Task [Reward_task] 0.86 0.82 – 0.90 <0.001
Difficulty sc × noopts
s2z1
0.99 0.98 – 1.00 0.001
Difficulty sc × Task
[Reward_task]
1.03 1.01 – 1.05 0.009
noopts s2z1 × Task
[Reward_task]
0.98 0.97 – 0.99 <0.001
(Difficulty sc × noopts
s2z1) × Task
[Reward_task]
1.00 0.99 – 1.00 0.330
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.Difficulty_sc 0.00
τ11 shortid.noopts_s2z1 0.00
τ11 shortid.TaskReward_task 0.01
τ11 shortid.Difficulty_sc:noopts_s2z1 0.00
τ11 shortid.Difficulty_sc:TaskReward_task 0.00
τ11 shortid.noopts_s2z1:TaskReward_task 0.00
ρ01 -0.36
0.16
-0.66
0.24
0.30
-0.13
ICC 0.09
N shortid 61
Observations 19482
Marginal R2 / Conditional R2 0.302 / 0.362

Models incorporating apathy

Below shows the model incorporating DAS (Dimensional Apathy Score) total score:

  RT
Predictors Estimates CI p
(Intercept) 1.00 0.96 – 1.05 0.956
Difficulty sc 1.16 1.14 – 1.18 <0.001
noopts s2z1 0.92 0.91 – 0.93 <0.001
DAS total sc 1.04 1.00 – 1.09 0.060
Difficulty sc × noopts
s2z1
0.99 0.98 – 0.99 <0.001
Difficulty sc × DAS total
sc
0.99 0.97 – 1.00 0.148
noopts s2z1 × DAS total
sc
1.01 1.00 – 1.01 0.146
(Difficulty sc × noopts
s2z1) × DAS total sc
1.00 0.99 – 1.00 0.644
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.Difficulty_sc 0.00
τ11 shortid.noopts_s2z1 0.00
τ11 shortid.Difficulty_sc:noopts_s2z1 0.00
ρ01 -0.38
0.15
0.31
ICC 0.09
N shortid 86
Observations 13692
Marginal R2 / Conditional R2 0.257 / 0.327

And below, incorporating DAS executive score: -Slower responses with higher executive apathy -Slower responses especially in easier tasks in those with higher apathy

  RT
Predictors Estimates CI p
(Intercept) 1.00 0.96 – 1.05 0.947
Difficulty sc 1.16 1.14 – 1.18 <0.001
noopts s2z1 0.92 0.91 – 0.93 <0.001
DAS Exec sc 1.05 1.01 – 1.10 0.027
Difficulty sc × noopts
s2z1
0.99 0.98 – 0.99 <0.001
Difficulty sc × DAS Exec
sc
0.98 0.96 – 0.99 0.002
noopts s2z1 × DAS Exec sc 1.01 1.00 – 1.02 0.061
(Difficulty sc × noopts
s2z1) × DAS Exec sc
1.00 1.00 – 1.01 0.757
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.Difficulty_sc 0.00
τ11 shortid.noopts_s2z1 0.00
τ11 shortid.Difficulty_sc:noopts_s2z1 0.00
ρ01 -0.36
0.14
0.28
ICC 0.09
N shortid 86
Observations 13692
Marginal R2 / Conditional R2 0.264 / 0.332

There are no signfiicant main or interacting effects of DAS auto-activation or DAS emotional-affective dimensions on RT (not shown).

When effort and reward subtasks are analysed separately (for 61 participants that did both), significant difficulty:DAS executive interaction is only seen in effort task (Effort task: p=0.014; Reward task p=0.77).

Accuracy

Base models

Below shows the base model for accuracy:

  accuracy
Predictors Odds Ratios CI p
(Intercept) 36.26 27.17 – 48.39 <0.001
Difficulty sc 0.28 0.25 – 0.32 <0.001
noopts s2z1 1.22 1.12 – 1.34 <0.001
Difficulty sc × noopts
s2z1
0.99 0.93 – 1.07 0.868
Random Effects
σ2 3.29
τ00 shortid 1.42
τ11 shortid.Difficulty_sc 0.15
ρ01 shortid -0.55
ICC 0.32
N shortid 87
Observations 13852
Marginal R2 / Conditional R2 0.251 / 0.493

Below shows the base accuracy model incorporating Task: -Higher accuracy in more difficult choices in reward versus effort task

  accuracy
Predictors Odds Ratios CI p
(Intercept) 34.11 24.43 – 47.62 <0.001
Difficulty sc 0.28 0.24 – 0.33 <0.001
noopts s2z1 1.17 1.05 – 1.30 0.003
Task [Reward_task] 1.10 0.78 – 1.56 0.579
Difficulty sc × noopts
s2z1
1.00 0.92 – 1.09 0.933
Difficulty sc × Task
[Reward_task]
1.40 1.14 – 1.73 0.002
noopts s2z1 × Task
[Reward_task]
0.93 0.80 – 1.08 0.343
(Difficulty sc × noopts
s2z1) × Task
[Reward_task]
1.00 0.88 – 1.13 0.942
Random Effects
σ2 3.29
τ00 shortid 1.32
τ11 shortid.Difficulty_sc 0.15
τ11 shortid.TaskReward_task 1.07
τ11 shortid.Difficulty_sc:TaskReward_task 0.26
ρ01 -0.55
-0.72
0.38
ICC 0.26
N shortid 61
Observations 19482
Marginal R2 / Conditional R2 0.220 / 0.420

Models incorporating apathy

No models show significant effects of DAS total or any of the three subscales on accuracy.

Timeouts

Timeouts have been analysed in two ways:

  1. A simple negative binomial regression model specified as count of timeouts ~ apathy score. In 86 people doing effort task, hese find significant effects of DAS total (p=0.049), DAS executive (p=0.007), but neither DAS auto-activation or DAS emotional. A similar linear model with predictors of Task*Task order in 61 ppts that did both effort and reward task finds significant effects of Task (less timeouts in reward task, p=0.011) and no effect of task order or interaction. Incorporating apathy scores into this model shows main effects of DAS total and DAS executive but no interactions with task.
  2. A mixed effects binomial regression with timeout (yes or no) as the dependent variable and Difficulty + Options + Difficuly:Options as predictors. The base model shows significant main effects of difficulty (p<0.001), options (p<0.001) but no interaction. Incorporating apathy into models shows main effects of DAS executive (p=0.015). No interactions with DAS executive, and no main or interacting effects of DAS total or any other subscales.

Ratings

Below shows the mixed effect linear model for subjective ratings (z-scored to individual participants) in 87 participants who did effort task. Here a gamma log link is applied for skewed data.

  Dependent variable
Predictors Estimates CI p
(Intercept) 1.39 1.14 – 1.70 0.001
Difficulty 1.59 1.51 – 1.68 <0.001
noopts s2z1 0.96 0.91 – 1.00 0.068
Difficulty × noopts s2z1 0.98 0.94 – 1.03 0.469
Random Effects
σ2 0.70
τ00 shortid 0.84
ICC 0.55
N shortid 87
Observations 2953
Marginal R2 / Conditional R2 0.124 / 0.603

Interestingly, incorporating apathy scores into these subjective rating models produces some interesting results. There is a positive effect of DAS executive score on ratings, and also a highly significant Difficulty:Executive apathy score interaction, in accordance with the data - more apathetic people find particularly easier decisions more subjectively difficulty

  Dependent variable
Predictors Estimates CI p
(Intercept) 1.40 1.15 – 1.69 0.001
Difficulty 1.60 1.52 – 1.68 <0.001
noopts s2z1 0.96 0.91 – 1.00 0.078
DAS Exec sc 1.26 1.04 – 1.53 0.018
Difficulty × noopts s2z1 0.99 0.94 – 1.04 0.675
Difficulty × DAS Exec sc 0.88 0.84 – 0.92 <0.001
noopts s2z1 × DAS Exec sc 1.01 0.97 – 1.06 0.545
(Difficulty × noopts
s2z1) × DAS Exec sc
1.00 0.95 – 1.05 0.999
Random Effects
σ2 0.70
τ00 shortid 0.80
ICC 0.53
N shortid 86
Observations 2919
Marginal R2 / Conditional R2 0.163 / 0.609

There is also a significant interaction between difficulty and the auto-activation dimension. The emmeans plots suggest this is in the opposite direction - that those with greater behavioural apathy find less difficulty decisions easier and more difficult decisions harder. However the mean plots below don’t support this interaction with the auto-activation dimension. Of note there are no main effects or interactions with DAS total or DAS emotional scores.

  Dependent variable
Predictors Estimates CI p
(Intercept) 1.39 1.14 – 1.70 0.001
Difficulty 1.59 1.51 – 1.67 <0.001
noopts s2z1 0.96 0.91 – 1.00 0.074
DAS Behav sc 0.87 0.71 – 1.07 0.188
Difficulty × noopts s2z1 0.98 0.94 – 1.03 0.462
Difficulty × DAS Behav sc 1.11 1.05 – 1.17 <0.001
noopts s2z1 × DAS Behav
sc
1.00 0.95 – 1.05 0.930
(Difficulty × noopts
s2z1) × DAS Behav sc
0.99 0.94 – 1.03 0.561
Random Effects
σ2 0.71
τ00 shortid 0.89
ICC 0.56
N shortid 86
Observations 2919
Marginal R2 / Conditional R2 0.134 / 0.616

Control analyses

One question is whether observed effects of apathy are due to worse learning of stimulus-associations causing slower responses. This seems unlikely as no correlations with accuracy and effects are greater in easier tasks which are well-learned. It is also possible that unsatisfacory learning for certain stimuli may skew results (i.e. consistently get one stimulus wrong).

Therefore can perform subanalysis on participants of good learners by excluding those who consistently scored below 50% on accuracy on at least one of the different stimuli. This excludes 24 participants however the executive apathy results remain significant:

  RT
Predictors Estimates CI p
(Intercept) 0.98 0.93 – 1.03 0.357
Difficulty sc 1.17 1.15 – 1.19 <0.001
noopts s2z1 0.92 0.91 – 0.93 <0.001
DAS Exec sc 1.05 1.00 – 1.10 0.055
Difficulty sc × noopts
s2z1
0.99 0.98 – 0.99 <0.001
Difficulty sc × DAS Exec
sc
0.98 0.96 – 1.00 0.029
noopts s2z1 × DAS Exec sc 1.00 1.00 – 1.01 0.337
(Difficulty sc × noopts
s2z1) × DAS Exec sc
1.00 0.99 – 1.01 0.787
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.Difficulty_sc 0.00
τ11 shortid.noopts_s2z1 0.00
τ11 shortid.Difficulty_sc:noopts_s2z1 0.00
ρ01 -0.32
0.13
0.26
ICC 0.09
N shortid 62
Observations 9885
Marginal R2 / Conditional R2 0.294 / 0.358

Another question is whether the effects of apathy relate to adaptive speeding of responses over time, particularly after stimuli have been presented multiple times (the Gratton effect) - as opposed to static slowed responses, relating either to executive impairments or effort aversion.

To look into this you can separate choices into stimulus of the option chosen and plot RT against trial number of choosing that stimulus. Here trials are shown up to the point at which 50% of participants have selected that option.

Specifying a linear model of RT ~ Trial * DAS executive * Option chosen + (1 | ID) (more complex models don’t converge) suggests that both effects are present:

-Main effect of DAS executive - apathy associated with slower responses

-Interaction between DAS executive and trial - less adaptive speeding with time.

-Same DAS executive and Option chosen interaction as demonstrated previously - lower apathy individuals have quicker responses in low effort (easy) tasks

  RT
Predictors Estimates CI p
(Intercept) 0.97 0.93 – 1.01 0.176
Option chosen 1.14 1.12 – 1.16 <0.001
DAS Exec sc 1.06 1.01 – 1.11 0.012
it 0.98 0.97 – 0.98 <0.001
Option chosen × DAS Exec
sc
0.98 0.96 – 1.00 0.013
Option chosen × it 0.99 0.99 – 1.00 0.069
DAS Exec sc × it 1.01 1.01 – 1.02 <0.001
(Option chosen × DAS Exec
sc) × it
1.00 0.99 – 1.00 0.445
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.scale(Option_chosen) 0.00
ρ01 shortid -0.25
ICC 0.10
N shortid 86
Observations 12020
Marginal R2 / Conditional R2 0.211 / 0.288

Force

I have analysed force (AUC as a proportion of max AUC over all trials in effort task). I have specified a model of Force ~ Option worked * Apathy score + (1 + Option worked | ID).

These show that the executive apathy domain is the main driver of less force, although there are weaker correlations with DAS total and DAS behavioural.

# A tibble: 7 × 2
  Option_worked `mean(AUC_norm)`
          <int>            <dbl>
1             1            0.196
2             2            0.277
3             3            0.365
4             4            0.467
5             5            0.523
6             6            0.569
7             7            0.581
  AUC_norm
Predictors Estimates CI p
(Intercept) 0.31 0.29 – 0.34 <0.001
Option worked 0.11 0.10 – 0.12 <0.001
DAS total 0.00 -0.02 – 0.03 0.796
Option worked × DAS total -0.01 -0.02 – 0.00 0.052
Random Effects
σ2 0.03
τ00 shortid 0.01
τ11 shortid.scale(Option_worked) 0.00
ρ01 shortid 0.57
ICC 0.33
N shortid 87
Observations 13900
Marginal R2 / Conditional R2 0.227 / 0.486
  AUC_norm
Predictors Estimates CI p
(Intercept) 0.31 0.29 – 0.34 <0.001
Option worked 0.11 0.10 – 0.12 <0.001
DAS Behav 0.00 -0.02 – 0.03 0.779
Option worked × DAS Behav -0.01 -0.02 – 0.00 0.088
Random Effects
σ2 0.03
τ00 shortid 0.01
τ11 shortid.scale(Option_worked) 0.00
ρ01 shortid 0.57
ICC 0.33
N shortid 87
Observations 13900
Marginal R2 / Conditional R2 0.227 / 0.486
  AUC_norm
Predictors Estimates CI p
(Intercept) 0.31 0.29 – 0.34 <0.001
Option worked 0.11 0.10 – 0.12 <0.001
DAS Exec -0.01 -0.03 – 0.02 0.522
Option worked × DAS Exec -0.01 -0.02 – -0.00 0.007
Random Effects
σ2 0.03
τ00 shortid 0.01
τ11 shortid.scale(Option_worked) 0.00
ρ01 shortid 0.55
ICC 0.33
N shortid 87
Observations 13900
Marginal R2 / Conditional R2 0.231 / 0.486

This does not seem to be related to the number of failures in a simple linear model of per participant failures count of failures being predicted by apathy scores (ps all >0.3) or cumulative force production (sum of AUC for all trials).

Distractor effects

Finally I have looked at effects of distracting options. Here I have only analysed 3 option tasks. I have split task on the basis of the best option and second best option, and divided choices into low distractor value (DV) and high distractor value. This equates to:

Bestoption Second best Low DV High DV Excluded
1 2 3 4 6 7 5
2 3 4 5 6 7
3 4 5 7 6
4 5 6 7
5 6 - - 7

In these plots the best option goes across the x axis at the top and second best option along the y axis. Most difficult choices (lowest value difference) therefore lie across the diagnoal. First I have plotted participants’ mean RTs:

Although the effect sizes look small, there is a significant effect of distractor value:

  RT
Predictors Estimates CI p
(Intercept) 1.02 0.97 – 1.08 0.377
bestoption 1.18 1.15 – 1.22 <0.001
value difference 0.98 0.96 – 1.01 0.169
dval [low] 0.97 0.95 – 0.99 0.008
bestoption × value
difference
1.01 0.99 – 1.03 0.271
bestoption × dval [low] 0.98 0.95 – 1.00 0.085
value difference × dval
[low]
0.99 0.97 – 1.02 0.504
(bestoption × value
difference) × dval [low]
0.99 0.97 – 1.01 0.360
Random Effects
σ2 0.08
τ00 shortid 0.01
τ11 shortid.scale(bestoption) 0.00
τ11 shortid.scale(value_difference) 0.00
ρ01 -0.22
0.16
ICC 0.17
N shortid 87
Observations 4161
Marginal R2 / Conditional R2 0.231 / 0.362

I have also looked at accuracy but this doesn’t show anything significant.

  accuracy
Predictors Odds Ratios CI p
(Intercept) 67.17 33.86 – 133.23 <0.001
bestoption 0.54 0.34 – 0.86 0.009
value difference 7.28 4.15 – 12.76 <0.001
dval [low] 1.00 0.66 – 1.50 0.988
bestoption × value
difference
0.66 0.45 – 0.97 0.032
bestoption × dval [low] 0.97 0.57 – 1.65 0.896
value difference × dval
[low]
0.97 0.65 – 1.45 0.869
(bestoption × value
difference) × dval [low]
1.08 0.68 – 1.71 0.757
Random Effects
σ2 3.29
τ00 shortid 2.74
τ11 shortid.scale(bestoption) 0.17
τ11 shortid.scale(value_difference) 0.78
ρ01 0.34
0.97
ICC 0.52
N shortid 87
Observations 4161
Marginal R2 / Conditional R2 0.460 / 0.739