## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Interplay of approximate planning strategies

Edited by Larry F. Abbott, Columbia University, New York, NY, and approved December 31, 2014 (received for review July 25, 2014)

## Significance

Many problems, particularly sequential planning problems, are computationally very demanding. How humans combine strategies to approximate and simplify these problems is not understood. Using modelling to unpick performance in a planning task, we find that humans are able to exploit the structure of the task to subdivide it and reduce processing requirements nearly optimally. Subtasks are combined in a simple, greedy manner, however, and within subtasks there is evidence of inhibitory reflexes in response to losses.

## Abstract

Humans routinely formulate plans in domains so complex that even the most powerful computers are taxed. To do so, they seem to avail themselves of many strategies and heuristics that efficiently simplify, approximate, and hierarchically decompose hard tasks into simpler subtasks. Theoretical and cognitive research has revealed several such strategies; however, little is known about their establishment, interaction, and efficiency. Here, we use model-based behavioral analysis to provide a detailed examination of the performance of human subjects in a moderately deep planning task. We find that subjects exploit the structure of the domain to establish subgoals in a way that achieves a nearly maximal reduction in the cost of computing values of choices, but then combine partial searches with greedy local steps to solve subtasks, and maladaptively prune the decision trees of subtasks in a reflexive manner upon encountering salient losses. Subjects come idiosyncratically to favor particular sequences of actions to achieve subgoals, creating novel complex actions or “options.”

Humans and other animals often face complex tasks and environments in which they have to plan and execute long sequences of appropriate actions to achieve distant goals. One can represent the space of future actions and outcomes as a tree; such trees grow inordinately (often exponentially) large as a function of the length of the sequence (i.e., the depth of the tree). Rather little is definitively known about how this computational complexity is addressed. Work in the fields of reinforcement learning and artificial intelligence has suggested a number of heuristics that we describe below, namely, hacking, hierarchies, hoarding, and habitization (1⇓⇓–4). Various tasks have been designed to highlight individual heuristics; though how subjects generate and combine them without clear instruction has not been well characterized (however, see refs. 5 and 6).

We previously designed a moderately deep planning problem to elicit a specific heuristic, in this case hacking or pruning of the decision tree (4). However, the task contains many of the elements that make choosing appropriately tricky in general. Thus, we closely examined the nature of, and individual differences between, the performance of subjects, shedding light on the interaction of heuristics in the self-generation of adaptive control when faced with a complex planning problem.

Subjects had to plan a path through a maze so as to maximize their cumulative earnings. On each trial, they were placed in a random state and were asked to plan to a depth of 3, 4, or 5 moves (Fig. 1 *A* and *B*). Because each depth involved a binary choice, planning to depths 3, 4, and 5 corresponded to choosing among a set of 8, 16, or 32 possible sequences. We previously found that the large immediate losses at particular branch points in the tree (the red transitions) encouraged subjects to eliminate possibly lucrative subbranches beneath those points (4). This corresponds to suboptimal pruning or “hacking” of the decision tree (Fig. 1*C*). The analyses presented below show that this was by no means the only strategy subjects used.

Hierarchical task decomposition licenses strategies for reducing computational burdens based on divide and conquer (2, 7) or “chunking” (5, 8⇓⇓⇓–12). The resulting divided problems are easier to conquer because chunks are smaller and ignore aspects of the environment that do not impinge on their domains. The solutions to the subproblems can then be treated as larger-scale actions, often called macroactions or options (1). These simplify solving complex tasks by providing a way of building large decision trees out of smaller numbers of intermediate-sized parts (the macroactions) rather than larger numbers of smaller parts (each individual action). The downside is potential suboptimality. We use the precise form of suboptimality that our subjects exhibited to argue that they hierarchically fragmented the planning problems: Deep problems were solved by concatenating solutions to sequences of shallower problems (e.g., greedily adopting as a depth-5 solution the best depth-3 solution followed by the best remaining depth-2 solution; Fig. 1*D*). We then asked a critical question that has eluded previous approaches to hierarchical control, namely the degree to which the fragmentation is actually computationally advantageous—is the benefit of divide and conquer appropriately realized?

A third, “hoarding” heuristic is known as memoization. If subjects are repeatedly faced with the task of finding a good policy at a state, then rather than building and searching the decision tree each time it is sensible to recall a previous solution and use that. If the previous solution cannot be guaranteed to be correct, then storing and deciding among several past solutions could be wise. When such storage and recall are probabilistic, the heuristic is stochastic memoization. It has been most extensively investigated in computational linguistics (13–15), and recently imported into decision making (16). Hoarding and hierarchies interact closely: Stored solutions can exactly be considered as macroactions or chunks, and so stochastic memoization can be seen as an answer to another poorly explored question, namely, how hierarchical decompositions arise. In particular, we will see that different subjects fragmented the task in idiosyncratic ways, putatively based on the way that they memoized.

We used both flexible and constrained statistical analyses to examine the use of these heuristics. For the constrained analyses, we stipulated a particular mathematical form for each cognitive process and implemented it in a model that, after its parameters had been fit, reported the likelihood of the subjects’ choices. More complex models should provide better fits and were penalized by computing integrated Bayesian information criterion (iBIC) scores, which approximate Bayes factors (4, 17). We tested models including and excluding particular cognitive processes. Those processes present in the model with the best iBIC score were taken as putatively present in subjects’ decision making. Importantly, this approach always tests the ability of the hypothesized set of cognitive processes to account for the entire dataset, rather than only hand-selected aspects of the data.

## Results

Subjects were trained extensively on both the transition matrix and the rewards associated with each of the transitions. As previously reported (4), subjects “pruned” (Fig. 1*C*) extensively, with outcomes distal to large losses being discounted at a faster rate than outcomes distal to non-large loss outcomes (*Supporting Information*, *Pruning* and Fig. S1). We previously found this to be insensitive to the actual size of the large loss and to reduce earnings overall, and therefore interpreted it as a reflexive “Pavlovian” influence on goal-directed decisions. All findings below incorporate and extend these findings, that is, we correct for them throughout, and ultimately reassess pruning in the light of the more refined analyses of other heuristics. The time pressures imposed in this version of task, and the fact that planning had to be completed before any move was registered, led to few apparent differences in choices relative to the original version (4).

### Fragmenting Decisions.

We first looked for a very simple type of hierarchy, whereby a smaller initial subproblem of a difficult task is addressed greedily, leaving whatever remains to be solved in turn. It seemed that subjects indeed adopted this approach. For instance, starting from state 1, the optimal path for a depth-3 problem is 1–2–3–4, whereas for a depth-4 problem it is 1–2–5–1–2. However, subjects had a strong predilection for the path 1–2–3–4–2 (Fig. S2, first column, middle row). A similar pattern was qualitatively observed throughout (other panels of Fig. S2). To quantify this statistically, we computed each subject’s distribution over depth-3 sequences (cf. Fig. S2) in the depth-3 problems and at the beginning of depth-4 problems and then performed Spearman rank correlation between these distributions over choice sequences. The correlation was on average 0.44. Comparison with correlations obtained from permutations of the distributions revealed that it was individually significant (*P* < .05) for 35/37 subjects (94%). Repeating the analysis by comparing the initial depth-4 sequences of depth-5 choices to the depth-4 choices, the correlation was on average 0.47 and individually significant in all subjects.

Subjects thus seemed to solve harder problems hierarchically, by exploiting the solutions to fragments, which are themselves smaller problems. However, just as there are many different ways of subdividing a complex task into simpler subtasks, there are many ways simpler problem solutions could be substituted into harder problems. To measure directly the hierarchical decomposition and avoid the potential biases in the above simple analysis, we fitted an exploratory model with sufficient flexibility to capture all possible fragmentations (i.e., for each trial, we asked which fragmented decomposition best fitted subjects’ actual choices).

In this model, we first examined fragment endpoints (red asterisk in Fig. 1*D*). To do so, we extracted all fragmented choices, that is, we extracted the maximum a posteriori decomposition of each subject’s choices, determined the start and endpoints of each fragment, and constructed histograms of endpoints as a function of the fragment origin. Fig. 2*A*, *Left* shows that fragments tended to terminate in state 2 irrespective of the start state (seen as a dark horizontal band), or in the next state around the outer ring of the maze (seen as a dark band below the diagonal). This pattern accounted for 90% of target end states on average (range 0.83–0.99) and was present in all subjects above chance (*P* < 10^{−10} for each; *Supporting Information*, *Optimal Fragmentation* and Fig. S3). Fragments of length 1 tended to end in the next state along the circle (Fig. 2*A*, *Middle*), whereas fragments of greater length terminated in state 2 irrespective of where they started (Fig. 2*A*, *Right*) and typically reaped the big reward when transitioning from state 1. Importantly, this structure was stable through the hierarchical decomposition: It was present when considering only the first fragment at the root of the decision tree and also when considering only the second fragment inside the decision tree (Fig. S4*A*). Finally, we examined the delays between button presses. Button press delays preceding depth-1 fragments were on average 9.7ms longer than the decision time before longer fragments [*t*(36) = 2.17, *P* = 0.037]. Because the computational requirement for shorter fragments is smaller than that for long fragments, this suggests that emitting a short fragment is the result of failed deeper searches. Thus, subjects seemed to build decision trees repeatedly until they found a path through the salient reward (transitioning from state 1 to 2); failing this, they would move one step along the circle (the off-diagonal band) and then reattempt to build a decision tree that would lead them to state 2. This suggests that subjects might have treated state 2 as a subgoal, which makes intuitive sense because this usually corresponds to earning the large reward of +140 when reached from state 1.

The driving force behind fragmentation should be the reduction in computational cost, a simple measure of which is the total number of computations—here largely additions—required to sum the rewards along all paths in a tree. A full evaluation of a depth-3 tree would require a sum of three outcomes for each of the eight arms, that is, *d* tree. Using this measure of computational cost, the fragmentation in Fig. 1*D* would require evaluation of a depth-2 and then a depth-1 tree, resulting in *D*). (The results remain the same when quantifying the cost of using dynamic programming, which is in general more efficient; see *Supporting Information*, *Optimal Fragmentation* for further details.) Given these costs, one can identify, for each state–depth problem, the fragmentation with the least computational cost that still selects the optimal choice. The distribution over fragment endpoints for the optimal decomposition is shown in Fig. 2*B* and matches the empirically inferred distributions well (Fig. 2*A*). Fig. 3 *A* and *B* show that the frequency of fragment depths and the overall distribution over end states also match the optimal distribution well. This suggests that participants chose a hierarchical task decomposition that nearly optimally reduced computational cost. For this task, the near-optimal fragmentation can be interpreted in terms of two organizing principles: one aligned with the perceptually salient transition structure (the ring of states) and the other with the reinforcement structure (the large rewards).

The distribution of fragment endpoints showed substantial structure that was present in every single subject. However, precisely because of the exploratory intention to allow for the discovery of any target pattern, this model (called “baseline + unrestricted fragmentation”) enjoyed very many degrees of freedom. It was able to reproduce the detailed choice patterns with high fidelity and accounted for 64% of the variability (Fig. S2, blue line). However, for the same reason, it overfitted the data. When integrating over all fragmentations it had a worse iBIC score than the baseline model without fragmentation (ΔiBIC = 682).

One way to capture the fragmentation structure would be to build a process model (i.e., an account of how subjects actually solve the metacontrol problem of decomposing a problem). This is beyond the scope of the current work. However, in an effort to constrain possible theoretical frameworks for this, we built a reduced model in which the full tree was replaced by the most frequently used fragmented version of the tree for each particular subject. For instance, if the tree starting in state 3 was most frequently decomposed into a tree of depth 2 followed by a depth-1 tree, as displayed in Fig. 1*D*, then this unique fragmentation was assumed to be fixed for the entire experiment for that state and depth. Fixing fragmentations for individuals did not substantially alter the fragment endpoint distributions (Fig. S4*B*) and the optimal pattern was still significantly present in every subject (*P* < 10^{−10}, binomial test). This model (“baseline + restricted fragmentation”) improved the fit over the baseline model (Fig. S5 *A* and *B*) and led to a substantial improvement in formal model comparison (Fig. 3*C*). Importantly, this model was as identifiable as the previous simpler models without fragmentation: We fitted each model to surrogate data generated from each of the models and were always able to recover the correct model (*Supporting Information*, *Robustness of Inference* and Table S1).

However, repeating this procedure, but now forcing all subjects to decompose problems in the same manner, produced a worse model (ΔiBIC = 116 compared with baseline model). Thus, fragmentation strategies were stable within individuals, but varied across different subjects. These results suggest that subgoals (the salient reward) and the decomposition of plans are intimately related. They constrain the processes that generate the fragmentation in the first place by identifying salient reinforcements as one central influence.

### Stochastic Memoization.

The fact that fragmentation is consistent for a given subject but varies between different subjects suggests the possibility that each subject generates one or a few possible decompositions and then sticks with this limited collection. Different subjects could generate and stick with different decompositions. For instance, the subject might initially generate action sequences through a laborious and computationally expensive tree search, but later on simply reuse a past solution. This memoization process (Fig. S6) would make later choices duplicate early fragmentations. One sign of this could be in the temporal evolution of the use of fragments, with those ultimately used most frequently coming to dominate the distribution of fragments slowly. Thus, the distribution over fragments should become more strongly concentrated on a few fragments as the task progresses. We computed the probability of fragments over time in the model “baseline + unrestricted fragmentation.” Fig. 3*D* shows that the frequency of the most commonly used fragment increased gradually over time. In marked contrast, the frequency with which all other choice fragments were chosen decayed over time. Fig. 3*E* shows that this results in the entropy over the fragment distribution decaying steeply.

Precisely because the task is too hard to solve perfectly, however, subjects cannot be sure that their previously computed choice sequence really represents the best option. Stochastic memoization refers to probabilistic, as opposed to deterministic, reuse. It is more appropriate when the result of the computation might change if it were recomputed, for instance due to incomplete or error-prone computation. One formalization of such a process is inspired by a method invented in computational linguistics (14, 16) that employs a distribution known as the “Chinese restaurant process” (CRP) (18). A CRP defines a probability distribution consisting of two terms. The first term is proportional to the frequency of past samples, whereas the second is the “base distribution” from which samples are drawn in the first place. Applied to the current problem, this model assumes that the probability of emitting a particular fragment is a weighted sum of two probabilities: the frequency of that particular fragment in the collection of previous choices and the probability that the fragment would be chosen anew if the solution to the problem were recomputed (i.e., the probability under the model baseline + restricted fragmentation). One critical feature of the CRP statistical model is a gradual change, with the choice probability being initially mainly driven by the base model (implying recomputation), and later by subjects’ past choices (implying stronger reuse later; see Eq. **2**). This transfer from flexible but costly computation to inflexible reliance on past experience is reminiscent of arguments about the transfer from goal-directed to habitual controllers (3, 19). However, by relying only on which choice was emitted rather than on how good it was, it also differs from certain formalizations of habits (3).

Fig. 3*C* and Fig. S5 show that the CRP addition in the augmented model “baseline + restricted fragmentation + stochastic memoization” outperformed the other models. Fitted reward sensitivity parameters correlated closely with the true reward sizes (mean of 0.994), and the target structure of the fragments was again not substantially altered by the inclusion of stochastic memoization (Fig. S4*C*). Adding stochastic memoization to the model with unrestricted fragmentation also improved all measures of model fit drastically (log likelihood improved for every subject by 43 ± 22; 14% more variance explained; ΔiBIC = 4,162). The same was true when controlling for an increase in the scaling of the reward sensitivities [i.e., an increase in exploitation (20)] over time (ΔiBIC = 86), and the model was also clearly identifiable on surrogate datasets (*Supporting Information*, *Robustness of Inference* and Table S1). Finally, corresponding parameters were highly correlated between all models tested (0.81 ± 0.06), suggesting that parameters captured similar variability in different models.

### Pruning.

Finally, we considered whether the pruning that we had previously seen using a similar task (4) might have been an artifact of the incomplete analysis of fragmentation and memoization. The baseline model that we fit included two pruning parameters: one that discounted outcomes distal to large losses *Supporting Information*, *Pruning*). However, the relationship between these two parameters was not constrained by the models. We examined the pruning parameters in the (overfitting) model “baseline + unrestricted fragmentation + stochastic memoization” because this model captured 78% of the variance, and hence controlled most strongly for all other processes. Fig. 3*F* shows that the continuation rate after outcomes other than large losses was indistinguishable from 1 (and hence *B*; however, note that owing to the length of the fragments the actual number of choices being part of longer fragments is higher), this also underscores our previous contention that pruning is a Pavlovian and reflexive response to aversive outcomes (4).

### Intelligence Quotient.

It has been suggested that subjects’ ability to decompose problems into larger chunks is a key ingredient of intelligence (21). The correlation between mean fragment length and verbal intelligence quotient (IQ) measured by a reading test was not significant (*ρ* = 0.08, *P* = 0.64).

## Discussion

Our results suggest that humans naturally decompose problems in a way that efficiently trades computational cost for performance; that the fragmentation of the task and the search strategy are shaped both by salient rewards and by salient perceptual features of the task; that subjects initially generate action sequences using a tree-search process, but later rely on a stored representation thereof; and that whenever subjects perform a mental search of a decision tree they have a tendency to prune the tree upon encountering salient losses. These features operate cumulatively: Subjects concurrently use multiple approximations when solving planning tasks. We were able to find sufficient evidence in favor of relatively complex models, potentially because of the high accuracy exhibited by the subjects—Fig. S2, bottom row, shows that even at depth 5 subjects rarely made very poor choices, even given the relatively tight time constraints.

### Fragmentation.

The fact that the decomposition achieved in this task is so close to the optimum (Fig. 3 *A* and *B*) is striking, particularly because finding an optimal fragmentation is typically more difficult than finding an optimal path. We were not able to address directly how this was achieved (i.e., to build a process model) but rather fitted all possible fragmentations to identify the relevant features.

One possibility is that the fragmentation arises from the exact sequence in which parts of the decision tree are searched. A hint comes from the suggestion that subjects search for the large reward, accept whatever path they find that leads them there, and then start recomputing from there onward. If they fail to find a path to the large reward, they move one step along the circle and try again. One immediate prediction from this is borne out: Fragments that only move one step along the circle often represent a failure of the search and should take longer to produce than longer fragments, even though they involve less computational cost. Thus, if it is this strategy that drives subjects’ internal search through the decision tree, then the optimality of the decomposition hinges directly on the relationship between the subgoal subjects aim for and the likelihood that this subgoal is on the optimal path. In the present task, the large reward subgoal was often on the optimal path. The fact that subjects rely so strongly on the salient reward in defining the subgoal bears some resemblance to the impact the large losses have. We have previously argued that the losses induce pruning in a reflexive, approximate, Pavlovian manner (4), and it may be that the selection of subgoals follows similar rules, as opposed to being derived in an adaptive manner from a clever insight into the task structure and optimality of various decompositions (see also ref. 21). It would certainly be highly instructive to alter the reward matrix such that this is not the case any more—one would expect the optimality of the decomposition to then break (see Fig. S7 for a predicted fragment endpoint distribution for a simple alteration to the reward matrix).

One might compare the hierarchical decomposition that we observed with those that have been studied in frameworks that explicitly set out to study hierarchical control (rather than to study pruning, which was our original target). A central construct in those tasks is that of a functional bottleneck—a state that makes a worthy subgoal because many paths have to flow through it by virtue of its position in state space (22⇓⇓⇓⇓–27). In our task, no state has this status—all states are equally connected—it is the reward structure that licenses the particular fragmentation.

Against our expectations (21) we did not find a correlation between verbal IQ and the average length of fragments. This is possibly because the Wechsler Test of Adult Reading (WTAR) used here is more a measure of verbal than fluid IQ. Alternatively, optimality in this task implies efficient selection of paths at the least computational cost. It is conceivable that the prediction should have been the opposite: that participants with higher IQ should have smaller average fragment length but achieve similar outcomes. A mixture between these two effects may explain the current null finding.

### Memoization and Option Generation.

One important insight that has engendered extensive research in decision making over the past decade is the distinction between goal-directed decisions and cached habits. Whereas the former suffer the sort of computational complexities that justify fragmentation, the latter suffer from a requirement for sampling from the world, or a model thereof: Instead of thinking through the future, the consequences of choices are experienced, and these experiences are cached to determine future choices (3, 28). Memoization is explicitly a form of caching, exhibiting the signature characteristic that if the environment suddenly changes (for instance via outcome devaluation or contingency degradation) the cached values will remain the same, and so control based on them will look maladaptive.

One common formalization of cached habit is in terms of state-action or **2**, because the probability of performing a whole new tree search (i.e., sampling from the base measure) decreases with the number of relevant trials so far. However, note that our form of stochastic memoization was independent of reward (i.e., memoization did not depend on the actual quality of the solution produced; it is therefore more like a refined form of choice kernel); indeed, it is known that nonhuman primate choices, for instance, depend substantially on their own past choices, above and beyond the rewards associated with the decisions (30, 31). Similar arguments have been made for human choices in a variety of tasks and settings (32, 33) and have been argued to be under dopaminergic (34) and serotonergic (35) control.

## Materials and Methods

### Participants.

We recruited 41 healthy volunteers (21 female; 23.3 ± 3.7 y) via the University College London psychology subject pool. They were screened for past and present psychiatric disorders (including drug and alcohol abuse) with the Mini International Neuropsychiatric Inventory (36). Subjects with past or current axis I diagnosis were excluded (one participant was excluded owing to previous substance dependence). Subjects completed the WTAR (37) (mean = 111, SD = 4.2) to assess IQ. The study was approved by the University College London Graduate School Ethics Committee. Subjects provided written, informed consent and were remunerated based on performance, up to a maximum of £40.

### Task.

The task is described in Fig. 1 and was adapted for functional MRI (fMRI) from one described in detail elsewhere (4) and programmed in Cogent 2000 (www.vislab.ucl.ac.uk/Cogent), a stimulus presentation toolbox for MATLAB (version 7.1). The fMRI results will be reported elsewhere, and we here report analyses only of the behavior of the same 37 subjects included there. Subjects were first extensively trained on the transition and reward matrix and all passed a test. Each of 90 trials of the main experiment began in a random starting state, but the combination of starting state and depth were biased such that in 60 trials it was optimal to transition through large losses, whereas in 30 trials the optimal path did not involve transition through a large loss. As part of the training, subjects had performed 32 trials that matched those of the main experiment, but in 18 of which there was no time restriction. The analyses presented here include these training trials. The experiment contained additional “restricted plan” trials, where subjects chose between two predefined paths, as a control condition for the fMRI analysis. These trials were not analyzed here.

### Fragmentation.

This model subdivided action sequences. The probability of a sequence **a** was represented as the product of the probability of *K* fragment action sequences:

where *s* and *d* denote start state and overall depth. The probabilities depend on the value *Supporting Information*, *Pruning*). The *d* can be subdivided into fragments of lengths 1 to *d* in

### Stochastic Memoization.

This model allowed for the reuse of fragments. Each entire action sequence was again subdivided as above into *k* segments. The probability of generating the particular segment *Supporting Information*, *Pruning*). The second component was proportional to how frequently that particular fragment action *s* so far, and *α*:

As *n* grows with time, the second term vanishes whereas the first term remains *n* updated accordingly.

### Model Fitting and Model Comparison.

We applied a nested model comparison strategy. We start from the simplest model and always evaluate whether additional model complexity is warranted by computing approximate Bayes factors (the integrated group-level BIC scores; see ref. 4) for each model. All models were fitted using MATLAB version 8.0 (MathWorks). We used the parallel processing toolbox and the function fminunc. All parameters were transformed to lie on the real line for inference. Models had the following number of parameters: lookahead, 1; discount, 2; pruning, 3; pruning + loss, 6; pruning + loss + restricted fragmentation, 24; and pruning + loss + restricted fragmentation + reuse, 25. See *Supporting Information*, *Robustness* for an assessment of the robustness of this approach.

## Acknowledgments

This study has appeared in abstract form and was funded by a British Academy grant to J.P.R. Q.J.M.H. received funds from the German Research Foundation and N.L. from a Wellcome Trust–National Institutes of Health (NIH) studentship. N.E. is supported by a Sackler Fellowship in Psychobiology and NIH Grants T32GM007753 and F30MH100729. E.S. is supported by the Swiss National Science Foundation, the University of Zürich, and the Neuroscience Centre Zürich and P.D. by the Gatsby Charitable Foundation.

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: qhuys{at}cantab.net. ↵

^{2}P.D. and J.P.R. contributed equally to this work.

Author contributions: Q.J.M.H., N.L., N.E., P.D., and J.P.R. designed research; Q.J.M.H., N.L., and P.F. performed research; Q.J.M.H. and S.J.G. contributed new analytic tools; Q.J.M.H., P.D., and J.P.R. analyzed data; and Q.J.M.H., N.L., P.F., N.E., E.S., S.J.G., P.D., and J.P.R. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

See Commentary on page 2929.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1414219112/-/DCSupplemental.

## References

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Dietterich TG

- ↵.
- Jog MS,
- Kubota Y,
- Connolly CI,
- Hillegaart V,
- Graybiel AM

- ↵.
- Koechlin E,
- Ody C,
- Kouneiher F

- ↵
- ↵.
- Ostlund SB,
- Winterbauer NE,
- Balleine BW

- ↵.
- Logan GD,
- Crump MJC

- ↵
- ↵.
- O’Donnell TJ,
- Goodman ND,
- Tenenbaum JB

- ↵.
- O'Donnell TJ

*Productivity and Reuse in Language: A Theory of Linguistic Computation and Storage*(MIT Press, Cambridge, MA) - ↵.
- Wingate D,
- Diuk C,
- O’Donnell T,
- Tenenbaum J,
- Gershman S

- ↵
- ↵
- ↵
- ↵.
- Sutton RS,
- Barto AG

- ↵
- ↵
- ↵.
- McGovern A,
- Barto AG

- ↵
- ↵.
- Şimşek Ö,
- Wolfe AP,
- Barto AG

- ↵
- ↵
- ↵.
- Huys QJM,
- Guitart-Masip M,
- Dolan RJ,
- Dayan P

- ↵.
- Watkins CJCH

- ↵
- ↵.
- Seo H,
- Barraclough DJ,
- Lee D

- ↵
- ↵
- ↵.
- Rutledge RB, et al.

- ↵.
- Seymour B,
- Daw ND,
- Roiser JP,
- Dayan P,
- Dolan R

- ↵
- ↵.
- Wechsler D

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Biological Sciences
- Neuroscience

## Jump to section

## You May Also be Interested in

_{2}into liquid methanol using artificial marine floating islands.

## See related content: