Abstract
Encoding precision in visual working memory decreases with the number of encoded items. Here, we propose a normative theory for such set size effects: the brain minimizes a weighted sum of an errorbased behavioral cost and a neural encoding cost. We construct a model from this theory and find that it predicts set size effects. Notably, these effects are mediated by probing probability, which aligns with previous empirical findings. The model accounts well for effects of both set size and probing probability on encoding precision in nine delayedestimation experiments. Moreover, we find support for the prediction that the total amount of invested resource can vary nonmonotonically with set size. Finally, we show that it is sometimes optimal to encode only a subset or even none of the relevant items in a task. Our findings raise the possibility that cognitive "limitations" arise from rational cost minimization rather than from constraints.
https://doi.org/10.7554/eLife.34963.001eLife digest
You can read this sentence from beginning to end without losing track of its meaning thanks to your working memory. This system temporarily stores information relevant to whatever task you are currently performing. However, the more items you try to hold in working memory at once, the poorer the quality of each of the resulting memories.
It has long been argued that this phenomenon – known as the set size effect – occurs because the brain devotes a fixed amount of neural resources to working memory. But this theory struggles to account for certain experimental results. It also fails to explain why the brain would not simply recruit more resources whenever it has more items to remember. After all, your heart does something similar by beating faster whenever you increase your physical activity.
Van den Berg and Ma break with the idea that working memory resources are fixed. They propose that resource allocation is flexible and driven by two conflicting goals: maximize memory performance, but use as few neural resources as necessary. Indeed, a computer simulation that follows this strategy mimics the set size effects seen in healthy volunteers. In the model, the items most relevant for a task are stored more accurately than less important ones, a phenomenon also observed in participants. Lastly, the simulation predicts that the total amount of resources devoted to working memory will vary with the number of items to be remembered. This too is consistent with the results of previous experiments.
Working memory thus appears to be more flexible than previously thought. The amount of resources that the brain allocates to working memory is not fixed but could be the result of balancing resource cost against cognitive performance. If this is confirmed, it may be possible to improve working memory by offering rewards, or by increasing the perceived importance of a task.
https://doi.org/10.7554/eLife.34963.002Introduction
A wellestablished property of visual working memory (VWM) is that the precision with which items are encoded decreases with the number of encoded items (Ma et al., 2014; Luck and Vogel, 2013). A common way to explain this set size effect has been to assume that there is a fixed amount of resource available for encoding: the more items, the less resource per item and, therefore, the lower the precision per item. Different forms have been proposed for this encoding resource, such as samples (Palmer, 1994; Sewell et al., 2014), Fisher information (van den Berg et al., 2012; Keshvari et al., 2013), and neural firing rate (Bays, 2014). Models with a fixed amount of resource generally predict that the encoding precision per item (defined as inverse variance of the encoding error) is inversely proportional to set size. This prediction is often inconsistent with empirical data, which is the reason that more recent studies instead use a power law to describe set size effects (Bays et al., 2009; Bays and Husain, 2008; van den Berg et al., 2012; van den Berg et al., 2014; Devkar et al., 2015; Elmore et al., 2011; Mazyar et al., 2012; Wilken and Ma, 2004; Donkin et al., 2016; Keshvari et al., 2013). In these powerlaw models, the total amount of resource across all items is no longer fixed, but instead decreases or increases monotonically with set size. These models tend to provide excellent fits to experimental data, but they have been criticized for lacking a principled motivation (Oberauer et al., 2016; Oberauer and Lin, 2017): they accurately describe how memory precision depends on set size, but not why these effects are best described by a power law – or why they exist at all. In the present study, we seek a normative answer to these fundamental questions.
While previous studies have used normative theories to account for certain aspects of VWM, none of them has accounted for set size effects in a principled way. Examples include our own previous work on change detection (Keshvari et al., 2012; Keshvari et al., 2013), change localization (van den Berg et al., 2012), and visual search (Mazyar et al., 2012). In those studies, we modelled the decision stage using optimalobserver theory, but assumed an ad hoc power law to model the relation between encoding precision and set size. Another example is the work by Sims and colleagues, who developed a normative framework in which working memory is conceptualized as an optimally performing information channel (Sims, 2016; Sims et al., 2012). Their informationtheoretic framework offers parsimonious explanations for the relation between stimulus variability and encoding precision (Sims et al., 2012) and the nonGaussian shape of encoding noise (Sims, 2015). However, it does not offer a normative explanation of set size effects. In their early work (Sims et al., 2012), they accounted for these effects by assuming that total information capacity is fixed, which is similar to other fixedresource models and predicts an inverse proportionality between encoding precision and set size. In their later work (Orhan et al., 2014; Sims, 2016), they add to this the assumption that there is an inefficiency in distributing capacity across items and fit capacity as a free parameter at each set size. Neither of these assumptions has a normative motivation. Finally, Nassar and colleagues have proposed a normative model in which a strategic tradeoff is made between the number of encoded items and their precision: when two items are very similar, they are encoded as a single item, such that there is more resource available per encoded item (Nassar et al., 2018). They showed that this kind of "chunking" is rational from an informationtheoretical perspective, because it minimizes the observer’s expected estimation error. However, just as in much of the work discussed above, this theory assumes a fixed resource budget for item encoding, which is not necessarily optimal when resource usage is costly.
The approach that we take here aligns with the recent proposal that cognitive systems are "resourcerational," that is, trade off the cost of using resources against expected task performance (Griffiths et al., 2015). The starting point of our theory is the principle that neural coding is costly (Attwell and Laughlin, 2001; Lennie, 2003; Sterling and Laughlin, 2015), which may have pressured the brain to trade off the behavioral benefits of high precision against the cost of the resource invested in stimulus encoding (Pestilli and Carrasco, 2005; Lennie, 2003; Ma and Huang, 2009; Christie and Schrater, 2015). We hypothesize that set size effects – and limitations in VWM in general – may be the result of making this tradeoff nearoptimally. We next formalize this hypothesis in a general model that can be applied to a broad range of tasks, analyze the theoretical predictions of this model, and fit it to data from nine previous delayedestimation experiments.
Theory
General theoretical framework: tradeoff between behavioral and neural cost
We define a vector Q={Q_{1},…, Q_{N}} that specifies the amount of resource with which each of N taskrelevant items is encoded. We postulate that Q affects two types of cost: an expected behavioral cost ${\overline{C}}_{\text{neural}}\left(\mathbf{Q}\right)$ induced by task errors and an expected neural cost ${\overline{C}}_{\text{neural}}\left(\mathbf{Q}\right)$ induced by spending neural resources on encoding. The expected total cost is a weighted combination,
where the weight λ≥0 represents the importance of the neural cost relative to the behavioral cost. Generally, increasing the amount of resource spent on encoding will reduce the expected behavioral cost, but simultaneously increase the expected neural cost.
The key novelty of our theory is that instead of assuming that there is a fixed resource budget for stimulus encoding (a hard constraint), we postulate that the brain – possibly on a trialbytrial basis – chooses its resource vector Q in a manner that minimizes the expected total cost. We denote the vector that yields this minimum by Q_{optimal}:
Under this policy, the total amount of invested resource – the sum of the elements of Q_{optimal} – does not need to be fixed: when it is "worth it" (i.e. when investing more resource reduces the expected behavioral cost more than it increases the expected neural cost), more resource may be invested.
Equations (1) and (2) specify the theory at the most general level. To derive testable predictions, we next propose specific formalizations of resource and of the two expected cost functions.
Formalization of resource
As in our previous work (Keshvari et al., 2012; Keshvari et al., 2013; Mazyar et al., 2012; van den Berg et al., 2012; van den Berg et al., 2014), we quantify encoding precision as Fisher information, J. This measure provides a lower bound on the variance of any unbiased estimator (Cover and Thomas, 2005; Ly et al., 2017) and is a common tool in the study of theoretical limits on stimulus coding and discrimination (Abbott and Dayan, 1999). Moreover, we assume that there is itemtoitem and trialtotrial variation in precision (Fougnie et al., 2012; van den Berg et al., 2012; van den Berg et al., 2014; Keshvari et al., 2013; van den Berg et al., 2017). Following our previous work, we model this variability using a gamma distribution with a mean $\overline{J}$ and shape parameter τ ≥0 (larger τ means more variability); we denote this distribution by gamma $\left(J;\overline{J},\tau \right)$.
We specify resource vector Q as the vector with mean encoding precisions, $\overline{\mathbf{J}}$, such that the general theory specified by Equations (1) and (2) modifies to
and
In this formulation, it is assumed that the brain has control over resource vector $\overline{J}$, but not over the variability in how much resource is actually assigned to an item. It should be noted, however, that our choice to incorporate variability in J is empirically motivated and not central to the theory: parameter τ mainly affects the kurtosis of the predicted estimation error distributions, not their variance or the way that the variance depends on set size (which is the focus of this paper). We will show that the theory also predicts set size effects when there is no variability in J.
Formalization of expected neural cost
To formalize the neural cost function, we make two general assumptions. First, we assume that the expected neural cost induced by encoding a set of N items is the sum of the expected neural cost associated with each of the individual items. Second, we assume that each of these “local” neural costs has the same functional dependence on the amount of allocated resource: if two items are encoded with the same amount of resource, they induce equal amounts of neural cost. Combining these assumptions, the expected neural cost induced by encoding a set of N items with resource $\overline{J}=\left\{{{\displaystyle \overline{J}}}_{1},\dots ,{{\displaystyle \overline{J}}}_{N}\right\}$ takes the form
where we introduced the convention to denote local costs (associated with a single item) with small c, to distinguish them from the global costs (associated with the entire set of encoded items), which we denote with capital C.
We denote by c_{neural}(J) the neural cost induced by investing an amount of resource J. The expected neural cost induced by encoding an item with resource $\overline{J}$ is obtained by integrating over J,
The theory is agnostic about the exact nature of the cost function ${c}_{\mathrm{n}\mathrm{e}\mathrm{u}\mathrm{r}\mathrm{a}\mathrm{l}}(J)$: it could include spiking and nonspiking components (Lennie, 2003), be associated with activity in both sensory and nonsensory areas, and include other types of cost that are linked to “mental effort” in general (Shenhav et al., 2017).
To motivate a specific form of this function, we consider the case that the neural cost is incurred by spiking activity. For many choices of spike variability, including the common one of Poissonlike variability (Ma et al., 2006), Fisher information J of a stimulus encoded in a neural population is proportional to the trialaveraged neural spiking rate (Paradiso, 1988; Seung and Sompolinsky, 1993). If we further assume that each spike has a fixed cost, we find that the local neural cost induced by each item is proportional to J,
where α is the amount of neural cost incurred by a unit increase in resource. Combining Equations (5–7) yields
Hence, the global expected neural cost is proportional to the total amount of invested resource and independent of the amount of variability in J. Although we use this linear expected neural cost function throughout the paper, we show in Appendix 1 that the key model prediction – a decrease of the optimal resource per item with set size – generalizes to a broad range of choices.
Formalization of expected behavioral cost for local tasks
Before we specify the expected behavioral cost function, we introduce a distinction between two classes of tasks. First, we define a task as "local" if the observer’s response depends on only one of the encoded items. Examples of local tasks are singleprobe delayedestimation (Blake et al., 1997; Prinzmetal et al., 1998; Wilken and Ma, 2004), singleprobe change detection (Todd and Marois, 2004; Luck and Vogel, 1997), and singleprobe change discrimination (Klyszejko et al., 2014). By contrast, when the task response depends on all memorized items, we define the task as "global." Examples of global tasks are wholedisplay change detection (Luck and Vogel, 1997; Keshvari et al., 2013), change localization (van den Berg et al., 2012), and delayed visual search (Mazyar et al., 2012). The theory that we developed up to this point – Equations (1–8) – applies to both global and local tasks. However, from here on, we develop our theory in the context of local tasks only; we will come back to global tasks at the end of the Results.
As in local tasks only one item gets probed, the expected behavioral cost across all items is a weighted average,
where p_{i} is the experimentally determined probing probability of the i^{th} item and ${{\displaystyle \overline{c}}}_{\text{behavioral},i}\left({{\displaystyle \overline{J}}}_{i};\tau \right)$ is the local expected behavioral cost associated with reporting the i^{th} item. We will refer to the product ${p}_{i}{{\displaystyle \overline{c}}}_{\text{behavioral},i}\left({{\displaystyle \overline{J}}}_{i};\tau \right)$ as the 'expected behavioral cost per item'. The only remaining step is to specify ${{\displaystyle \overline{c}}}_{\text{behavioral},i}\left({{\displaystyle \overline{J}}}_{i};\tau \right)$. This function is taskspecific and we will specify it after we have described the task to which we apply the model.
A resourcerational model for local tasks
Combining Equations 3, 8, and 9 yields the following expected total cost function for local tasks:
As parameters α and λ have interchangeable effects on the model predictions, we will fix α = 1 and only treat λ as a free parameter.
We recognize that the righthand side of Equation 10 is a sum of independent terms. Therefore, each element of ${{\displaystyle \overline{J}}}_{\text{optimal}}$, Equation 4, can be computed independently of the other elements, by minimizing the expected total cost per item,
This completes the specification of the general form of our resourcerational model for local tasks. Its free parameters are λ and τ.
Set size effects result from cost minimization and are mediated by probing probability
To obtain an understanding of the model predictions, we analyze how ${{\displaystyle \overline{J}}}_{\text{optimal}}$ depends on probing probability and set size. We perform this analysis under two general assumptions about the local expected behavioral cost function: first, that it monotonically decreases with $\overline{J}$ (i.e. increasing resource reduces the expected behavioral cost) and, second, that it satisfies a law of diminishing returns (i.e. the reductions per unit increase of resource decrease with the total amount of already invested resource). It can be proven (see Appendix 1) that under these assumptions, the domain of probing probability p_{i} consists of three potential regimes, each with a different optimal encoding strategy (Figure 1A). First, there might exist a regime 0≤p_{i}<p_{0} in which it is optimal to not encode an item, ${{\displaystyle \overline{J}}}_{\text{optimal}}=0$. In this regime, the probing probability of an item is so low that investing any amount of resource can never reduce the local expected behavioral cost by more than it increases the expected neural cost. Second, there might exist a regime p_{0}≤p_{i}<p_{∞} in which it is optimal to encode an item with a finite amount of resource, ${{\displaystyle \overline{J}}}_{\text{optimal}}\in \left(0,\infty \right)$. In this regime, ${{\displaystyle \overline{J}}}_{\text{optimal}}$ increases as a function of p_{i}. Finally, there may be a regime p_{∞}≤p_{i}≤1 in which the optimal strategy is to encode the item with an infinite amount of resource, ${{\displaystyle \overline{J}}}_{\text{optimal}}=\infty $. This last regime will only exist in extreme cases, such as when there is no neural cost associated with encoding. The threshold p_{0} depends on the importance of the neural cost, λ, and on the derivative of the local expected behavioral cost evaluated at $\overline{J}=0$; specifically, $p}_{0}=\frac{\lambda}{{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(0\right)$ . The threshold p_{∞} depends on λ and on the derivative of the local expected behavioral cost evaluated at $\overline{J}\to \infty $; specifically, $p}_{\mathrm{\infty}}=\frac{\lambda}{{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\mathrm{\infty}\right)$ . If p_{∞}>1, then the third regime does not exist, whereas if p_{0} >1, only the first regime exists.
We next turn to set size effects. An interesting property of the model is that $\overline{J}}_{\text{optimal}$ depends only on the probing probability, p_{i}, and on the model parameters – it does not explicitly depend on set size, N. Therefore, the only way in which the model can predict set size effects is through a coupling between N and p_{i}. Such a coupling exists in most studies that use a local task. For example, in delayedestimation tasks, each item is usually equally likely to be probed such that p_{i} = 1/N. For those experiments, the above partitioning of the domain of p_{i} translates to a similar partitioning of the domain of N (Figure 1B). Then, a set size N_{∞}≥0 may exist below which it is optimal to encode items with infinite resource, a region N_{∞}≤N < N_{0} in which it is optimal to encode items with a finite amount of resource, and a region N>N_{0} in which it is optimal to not encode items at all.
Results
Model predictions for delayedestimation tasks
To test the predictions of the model against empirical data, we apply it to the delayedestimation task (Wilken and Ma, 2004; Blake et al., 1997; Prinzmetal et al., 1998), which is currently one of the most widely used paradigms in VWM research. In this task, the observer briefly holds a set of items in memory and then reports their estimate of a randomly probed target item (Figure 2A). Set size effects manifest as a widening of the estimation error distribution as the number of items is increased (Figure 2B), which suggests a decrease in the amount of resource per item (Figure 2C).
To apply our model to this task, we express the expected local behavioral cost as an expected value of the behavioral cost with respect to the error distribution,
where the behavioral cost function c_{behavioral,i}(ε) maps an encoding error ε to a cost and $p\left(\epsilon ;{{\displaystyle \overline{J}}}_{i},\tau \right)$ is the predicted distribution of ε for an item encoded with resource ${{\displaystyle \overline{J}}}_{i}$. We first specify $p\left(\epsilon ;{{\displaystyle \overline{J}}}_{i},\tau \right)$ and then turn to c_{behavioral,i}(ε). As the taskrelevant feature in delayedestimation experiments is usually a circular variable (color or orientation), we make the common assumption that ε follows a Von Mises distribution. We denote this distribution by VM(ε;J), where J is onetoone related to the distribution’s concentration parameter κ (Appendix 1). The distribution of ε for a stimulus encoded with resource ${{\displaystyle \overline{J}}}_{i}$ is found by integrating over J,
Finally, we specify the behavioral cost function ${c}_{\text{behavioral},i}\left(\epsilon \right)$ in Equation 12, which maps an estimation error ε to a behavioral cost. As in most psychophysical experiments, human subjects tend to perform well on delayedestimation tasks even when the reward is independent of their performance. This suggests that the behavioral cost function is strongly determined by internal incentives. A recent paper (Sims, 2015) has attempted to measure this mapping and proposed a twoparameter function. We will test that proposal later, but for the moment we assume a simpler, oneparameter powerlaw function, ${c}_{\text{behavioral},i}\left(\epsilon ;\beta \right)={{\displaystyle \left\epsilon \right}}^{\beta}$, where power β is a free parameter.
To obtain an intuition for the predictions of this model, we plot in Figure 2D for a specific set of parameters the two expected costs per item and their sum, Equation 11, as a function of $\overline{J}$. The expected behavioral cost per item depends on p_{i} and decreases with $\overline{J}$ (colored curves in left panel), while the expected neural cost per item is independent of p_{i} and increases (black line in left panel). The expected total cost per item has a unique minimum (middle panel). The value of $\overline{J}$ corresponding to this minimum, ${{\displaystyle \overline{J}}}_{\text{optimal}}$, increases with p_{i} (Figure 2D, right). Hence, in this example, the optimal amount of resource per item is an increasing function of its probing probability.
We next consider the special case in which each item is equally likely to be probed, that is, p_{i} = 1/N. The values of p_{i} in Figure 2D then correspond to set sizes 1, 2, 4, and 8. When replotting ${{\displaystyle \overline{J}}}_{\text{optimal}}$ as a function of N, we find a set size effect (Figure 2E, right panel) that is qualitatively similar to the empirical result in Figure 2C. An alternative way to understand this predicted set size effect is by considering how the three expected costs across all items, Equation 3, depend on $\overline{J}$. Substituting p_{i} = 1/N in Equation 9, we find that the expected behavioral cost across all items is independent of set size (Figure 2E, left panel, black curve). Moreover, when all items are encoded with the same amount of resource (which is necessarily the optimal solution when p_{i} is identical across items), the expected neural cost across all items equals $N\overline{J}$ and therefore scales linearly with set size (Figure 2E, left panel, colored lines). The sum of these terms has a unique minimum ${{\displaystyle \overline{J}}}_{\text{optimal}}$ (Figure 2E, center panel), which monotonically decreases with set size (Figure 2E, right panel). The costs plotted in Figure 2E can be considered as obtained by multiplying the corresponding costs in Figure 2D by N.
The model thus predicts set size effects in delayedestimation tasks that are fully mediated by individualitem probing probability. The latter notion is consistent with empirical observations. Palmer et al. (1993) reported that "relevant set size" (where irrelevance means p_{i} = 0) acts virtually identically to actual set size. Emrich et al. (2017) independently varied probing probability and set size in their experiment, and found that the former was a better predictor of performance than the latter. Based on this, they hypothesized that set size effects are mediated by probing probability. The predictions of our model are qualitatively consistent with these findings.
Model fits to data from delayedestimation experiments with equal probing probabilities
To examine how well the model accounts for set size effects in empirical data, we fit it to data from six experiments that are part of a previously published benchmark set (E1E6 in Table 1). We use a Bayesian optimization method (Acerbi and Ma, 2017) to estimate the maximumlikelihood parameter values, separately for each individual data set (see Table 2 for a summary of these estimates). The model accounts well for the subjectlevel error distributions (Figure 3A) and the two statistics that summarize these distributions (Figure 3B). The original benchmark set (van den Berg et al., 2014) contained four more data sets, but three of those were published in papers that were later retracted and another one contains data at only two set sizes. Although we decided to leave those four datasets out of our main analyses, the model accounts well for them too (Figure 2—figure supplement 1).
We next compare the goodness of fit of the resourcerational model to that of a descriptive variant in which the amount of resource per item, $\overline{J}$, is assumed to be a powerlaw function of set size (all other aspects of the model are kept the same). This variant is identical to the VPA model in our earlier work, which is one of the most accurate descriptive models currently available (van den Berg et al., 2014). Model comparison based on the Akaike Information Criterion (AIC) (Akaike, 1974) indicates that the data provide similar support for both models, with a small advantage for the resourcerational model (ΔAIC = 5.27 ± 0.70; throughout the paper, X ± Y indicates mean ±s.e.m. across subjects). Hence, the resourcerational model provides a principled explanation of set size effects without sacrificing quality of fit compared to one of the best available descriptive models of VWM. We find that the resourcerational model also fits better than a model in which the total amount of resource is fixed and divided equally across items (ΔAIC = 13.9 ± 1.4).
So far, we have assumed that there is random variability in the actual amount of resource assigned to an item. Next, we test an equalprecision variant of the resourcerational model, by fixing parameter τ to a very small value (10^{−3}). Consistent with the results obtained with the variableprecision model, we find that the rational model has a substantial AIC advantage over a fixedresource model (ΔAIC = 43.0 ± 6.8) and is on equal footing with the powerlaw model (ΔAIC = 2.0 ± 1.7 in favor of the powerlaw model). However, all three equalprecision models (fixed resource, power law, rational) are outperformed by their variableprecision equivalents by over 100 AIC points. Therefore, we will only consider variableprecision models in the remainder of the paper.
To get an indication of the absolute goodness of fit of the resourcerational model, we next examine how much room for improvement there is in the fits. We do this by fitting a nonparametric model variant in which resource $\overline{J}$ is a free parameter at each set size, while keeping all other aspects of the model the same. We find a marginal AIC difference, suggesting that the fits of the rational model cannot be improved much further without overfitting the data (ΔAIC = 3.49 ± 0.93, in favor of the nonparametric model). An examination of the fitted parameter values corroborates this finding: the estimated resource values in the nonparametric model closely match the optimal values in the rational model (Figure 3C).
So far, we have assumed that behavioral cost is a powerlaw function of the absolute estimation error, c_{behavioral}(ε)=ε^{β}. To evaluate the necessity of a free parameter in this function, we also test three parameterfree choices: ε, ε^{2}, and −cos(ε). Model comparison favors the original model with AIC differences of 14.0 ± 2.8, 24.4 ± 4.1, and 19.5 ± 3.5, respectively. While there may be other parameterfree functions that give better fits, we expect that a free parameter is unavoidable here, as the errortocost mapping may differ across experiments (because of differences in external incentives) and also across subjects within an experiment (because of differences in intrinsic motivation). Finally, we also test a twoparameter function that was proposed recently (Equation (5) in Sims [2015]). The main difference with our original choice is that this alternative function allows for saturation effects in the errortocost mapping. However, this extra flexibility does not increase the goodness of fit sufficiently to justify the additional parameter, as the original model outperforms this variant with an AIC difference of 5.3 ± 1.8.
Finally, we use fivefold cross validation to verify the AICbased results reported in this section. We find that they are all consistent (Table 3).
Nonmonotonic relation between total resource and set size
One quantitative feature that sets the resourcerational theory apart from previous theories is its predicted relation between set size and the total amount of invested resource, ${{\displaystyle \overline{J}}}_{\text{total}}={\displaystyle {\sum}_{i=1}^{N}}\text{\hspace{0.17em}}{{\displaystyle \overline{J}}}_{i}$. This quantity is by definition constant in fixedresource models, and in powerlaw models it varies monotonically with set size. By contrast, we find that in the fits to several of the experiments, ${{\displaystyle \overline{J}}}_{\text{total}}$ varies nonmonotonically with set size (Figure 3D, gray curves). To examine whether there is evidence for nonmonotonic trends in the subject data, we next compute an "empirical" estimate ${{\displaystyle \overline{J}}}_{\text{total}}={\displaystyle {\sum}_{i=1}^{N}}\text{\hspace{0.17em}}{{\displaystyle \widehat{\overline{J}}}}_{i}$, where $\hat{\overline{J}}}_{i$ are the bestfitting resource estimates in the nonparametric model. We find that these estimates show evidence of similar nonmonotonic relations in some of the experiments (Figure 3D, black circles). To quantify this evidence, we perform Bayesian paired ttests in which we compare the estimates of ${{\displaystyle \overline{J}}}_{\text{total}}$ at set size 3 with the estimates at set sizes 1 and 6 in the experiments that included these three set sizes (E2 and E4E6). These tests reveal strong evidence that the total amount of resource is higher at set size 3 than at set sizes 1 (BF_{+0}=1.05·10^{7}) and 6 (BF_{+0}=4.02·10^{2}). We next compute for each subject the set size at which ${{\displaystyle \overline{J}}}_{\text{total}}$ is largest, which we denote by N_{peak}, and find a subjectaveraged value of 3.52 ± 0.18. Altogether, these findings suggest that the total amount of resource that subjects spend on item encoding varies nonmonotonically with set size, which is consistent with predictions from the resourcerational model, but not with any of the previously proposed models. To the best of our knowledge, evidence for a possible nonmonotonicity in the relation between set size and total encoding resource has not been reported before.
Predicted effects of probing probability
As we noted before, the model predictions do not explicitly depend on set size, N. Yet, we found that the model accounts well for set size effects in the experiments that we considered so far (E1E6). This happens because in all those experiments, N was directly coupled with probing probability p_{i}, through p_{i} = 1/N. This coupling makes it impossible to determine whether changes in subjects’ encoding precision are the result of changes in N or changes in p_{i}. Therefore, we will next consider experiments in which individual probing probabilities and set size were varied independently of each other (E7E9 in Table 1). According to our model, the effects of N that we found in E1E6 were really effects of p_{i}. Therefore, we should be able to make predictions about effects of p_{i} in E7E9 by recasting the effects of N in E1E6 as effects of p_{i} = 1/N. Given that the amount of resource per item in E1E6 decreases with N, a first prediction is that it should increase as a function of p_{i} in E7E9. A second and particularly interesting prediction is that the estimated total amount of invested resource should vary nonmonotonically with p_{i} and peak at a value p_{peak} that is close to 1/N_{peak} found in E1E6 (see previous section). Based on the values of N_{peak} in experiments E1E6, we find a prediction p_{peak} = 0.358 ± 0.026.
Model fits to data from delayedestimation experiments with unequal probing probabilities
To test the predictions presented in the previous section and, more generally, to evaluate how well our model accounts for effects of p_{i} on encoding precision, we fit it to data from three experiments in which probing probability was varied independently of set size (E7E9 in Table 1).
In the first of these experiments (E7), seven subjects performed a delayedestimation task at set sizes 2, 4, and 8. On each trial, one of the items – indicated with a cue – was three times more likely to be probed than any of the other items. Hence, the probing probabilities for the cued and uncued items were 3/4 and 1/4 at N = 2, respectively, 1/2 and 1/6 at N = 4, and 3/10 and 1/10 at N = 8. The subject data show a clear effect of p_{i}: the higher the probing probability of an item, the more precise the subject responses (Figure 4A, top row, black circles). We find that the resourcerational model, Equation (11), accounts well for this effect (Figure 4A, top row, curves) and does so by increasing the amount of resource as a function of probing probability p_{i} (Figure 4B, left panel, red curves).
In the other two experiments (E8 and E9), the number of cued items and cue validity were varied between conditions, while set size was kept constant at 4 or 6. For example, in one of the conditions of E8, three of the four items were cued with 100% validity, such that p_{i} was 1/3 for each cued item and 0 for the uncued item; in another condition of the same experiment, two of the four items were cued with 66.7% validity, meaning that p_{i} was 1/3 for each cued item and 1/6 for each uncued item. The unique values of p_{i} across all conditions were {0, 1/6, 2/9, 1/4, 1/3, 1/2, 1} in E8 and {0, 1/12, 1/10, 2/15, 1/6, 1/3, 1/2, and 1} in E9. As in E7, responses become more precise with increasing p_{i} and the model accounts well for this (Figure 4A), again by increasing the amount of resource assigned to an item with p_{i} (Figure 4B).
We next examine how our model compares to the models proposed in the papers that originally published these three data sets. In contrast to our model, both Bays (2014) and Emrich et al. (2017) proposed that the total amount of invested resource is fixed. However, while Bays proposed that the distribution of this resource is in accordance with minimization of a behavioral cost function (as in our model), Emrich et al. postulated that the resource is distributed in proportion to each item’s probing probability. Hence, while our model optimizes both the amount of invested resource and its distribution, Bays’ model only optimizes the distribution, and Emrich et al.’s model does not explicitly optimize anything. To examine how the three proposals compare in terms of how well they account for the data, we fit two variants of our model that encapsulate the main assumptions of these two earlier proposals. In the first variant, we compute ${{\displaystyle \overline{J}}}_{\text{optimal}}$ as $\underset{bf\overline{J}}{{\displaystyle \text{argmin}}}\text{\hspace{0.17em}}\left[{\displaystyle \sum _{i=1}^{N}}\text{\hspace{0.17em}}{p}_{i}{{\displaystyle \overline{c}}}_{\text{behavioral}}\left({{\displaystyle \overline{J}}}_{i};\beta ,\tau \right)\right]$ under the constraint $\sum _{i=1}^{N}}\text{\hspace{0.17em}}{{\displaystyle \overline{J}}}_{i}={{\displaystyle \overline{J}}}_{\text{total}$, which is consistent with Bays’ proposal. Hence, in this variant, the neural cost function is removed and parameter λ is replaced by a parameter ${{\displaystyle \overline{J}}}_{\text{total}}$ – otherwise, all aspects of the model are the same as in our main model. In the variant that we use to test Emrich et al.’s proposal, we compute ${{\displaystyle \overline{J}}}_{i}$ for each item as ${p}_{i}{{\displaystyle \overline{J}}}_{\text{total}}$, where p_{i} is the probing probability and ${{\displaystyle \overline{J}}}_{\text{total}}$ is again a free parameter that represents the total amount of resource. Fitting the models to the data from all 47 subjects in E7E9, we find a substantial advantage of our model over the proposal by Emrich et al., with an AIC difference of 18.0 ± 3.9. However, our model cannot reliably be distinguished from the proposal by Bays: either model is preferred in about half of the subjects (our model: 27; Bays: 20) and the subjectaveraged AIC difference is negligible (1.8 ± 2.5 in favor of our model). Hence, the model comparison suggests quite convincingly that subjects distribute their resource nearoptimally across items with unequal probing probabilities, but it is inconclusive regarding the question of whether the total amount of invested resource is fixed or optimized.
As an alternative way to address the question of whether the total amount of resource is fixed, we again fit a nonparametric model to obtain “empirical” estimates of the total amount of invested resource. To this end, we define $\hat{\overline{J}}}_{\text{total}}={\hat{\overline{J}}}_{i}/{p}_{i$, where $\hat{\overline{J}}}_{i$ are the bestfitting values in a nonparametric model, such that $\hat{\overline{J}}}_{\text{total}$ represents the estimated total amount of resource that a subject would invest to encode a display filled with items that all have probing probability p_{i}. We find that these estimates show signs of a nonmonotonicity as a function of p_{i} (Figure 4B, black points), which are captured reasonably well by the resourcerational model (Figure 4B, black curves). Averaged across all subjects in E7E9, the value of p_{i} at which $\hat{\overline{J}}}_{\text{total}$ is largest is 0.384 ± 0.037, which is close to the predicted value of 0.358 ± 0.026 (see previous section). Indeed, a Bayesian independentsamples ttest supports the null hypothesis that there is no difference (BF_{01} = 4.27). Hence, while the model comparison results in the previous paragraph were inconclusive regarding the question of whether the total amount of invested resource is fixed or optimized, the present analysis provides evidence against fixedresource models and confirms a prediction made by our own model.
In summary, the results in this section show that effects of probing probability in E7E9 are well accounted for by the same model as we used to explain effects of set size in E1E6. Regardless of whether total resource is fixed or optimized, this finding provides further support for the suggestion that set size effects are mediated by probing probability (Emrich et al., 2017) or, more generally, by item relevance (Palmer et al., 1993).
Is it ever optimal to not encode an item?
There is an ongoing debate about the question of whether a taskrelevant item is sometimes completely left out of working memory (Adam et al., 2017; Luck and Vogel, 2013; Ma et al., 2014; Rouder et al., 2008). Specifically, slot models predict that this happens when set size exceeds the number of slots (Zhang and Luck, 2008). In resource models, the possibility of complete forgetting has so far been an added ingredient separate from the core of the model (van den Berg et al., 2014). Our normative theory allows for a reinterpretation of this question: are there situations in which it is optimal to assign zero resource to the encoding of an item? We already established that this could happen in delayedestimation tasks: whenever the probing probability is lower than a threshold value ${p}_{0}={\scriptscriptstyle \frac{\lambda}{\left{{\displaystyle \overline{c}}}_{\text{behavioral}}\left(0\right)\right}}$, the optimal amount of resource to invest on encoding the item is zero (see Theory). But what values does p_{0} take in practice? Considering the expected behavioral cost function of a fixedprecision model (a variableprecision model with $\tau \downarrow 0$), we can prove that p_{0} = 0, that is, it is never optimal to invest no resource (Appendix 1). For the expected behavioral cost function of the variableprecision model, however, simulations indicate that p_{0} can be greater than 0 (we were not able to derive this result analytically). We next examine whether this ever happens under parameter values that are representative for human subjects. Using the maximumlikelihood parameters obtained from the data in E7E9, we estimate that p_{0} (expressed as a percentage) equals 8.86 ± 0.54%. Moreover, we find that for 8 of the 47 subjects, p_{0} is larger than the lowest probing probability in the experiment, which suggests that these subjects sometimes entirely ignored one or more of the items. For these subjects, the error distributions on items with p_{i}<p_{0} look uniform (see Figure 4C for an example) and KolmogorovSmirnov tests for uniformity did not reject the null hypothesis in any of these cases (p>0.05 in all tests).
These results suggest that there might be a principled reason why people sometimes leave taskrelevant items out of visual working memory in delayedestimation experiments. However, our model cannot explain all previously reported evidence for this. In particular, when probing probabilities are equal for all items, the model makes an "all or none" prediction: all items are encoded when p_{i}>p_{0} and none are encoded otherwise. Hence, the model cannot explain why subjects in tasks with equal probing probabilities sometimes seem to encode a subset of taskrelevant items. For example, a recent study reported that in a wholereport delayedestimation experiment (p_{i} = 1 for all items), subjects encoded about half of the six presented items on each trial (Adam et al., 2017). Unless additional assumptions are made, our model cannot account for this finding.
Predictions for a global task: wholedisplay change detection
The results so far show that the resourcerational model accounts well for data in a variety of delayedestimation experiments. To examine how its predictions generalize to other tasks, we next consider a change detection task, which is another widely used paradigm in research on VWM. In this task, the observer is sequentially presented with two sets of items and reports if any one of them changed (Figure 5A). In the variant that we consider here, a change is present on exactly half of the trials and is equally likely to occur in any of the items. We construct a model for this task by combining Equations 3, 4, and 8 with an expected behavioral cost function based on the Bayesian decision rule for this task (see Appendix 1), which yields
where $p\left(\text{error}\overline{J}\right)$ is the expected behavioral cost function, which in this case specifies the probability of an error response when a set of items is encoded with resource $\overline{J}$.
In contrast to local tasks, the expected total cost in global tasks cannot be written as a sum of expected costs per item, because the expected behavioral cost – such as $p\left(\text{error}\overline{J}\right)$ in Equation (14) – can only be computed globally, not per item. Consequently, the elements of ${{\displaystyle \overline{J}}}_{\text{optimal}}$ in global tasks cannot be computed separately for each item. This makes resource optimization computationally much more demanding, because it requires solving an Ndimensional minimization problem instead of N onedimensional problems.
We perform a simulation at N = 2 (which is still tractable) to get an intuition of the predictions that follow from Equation (14). For practical convenience, we assume in this simulation that there is no variability in precision, τ↓0, such that λ is the only model parameter. The results (Figure 5B) show that the costminimizing strategy is to encode neither of the items when the amount of reward per correct trial is very low (left panel) and encode them both when reward is high (right panel). However, interestingly, there is also an intermediate regime in which the optimal strategy is to encode one of the two items, but not the other one (Figure 5B, central panel). Hence, just as in the delayedestimation task, there are conditions in which it is optimal to encode only a subset of items. An important difference, however, is that in the delayedestimation task this only happens when items have unequal probing probabilities, while in this change detection task it even happens when all items are equally likely to change.
Simulations at larger set sizes quickly become computationally intractable, because of the reason mentioned above. However, the results at N = 2 suggest that if two items are encoded, the optimal solution is to encode them with the same amount of resource (Figure 5C). Therefore, we conjecture that all nonzero values in ${{\displaystyle \overline{J}}}_{\text{optimal}}$ are identical, which would mean that the entire vector can be summarized by two values: the number of encoded items, which we denote by K_{optimal}, and the amount of resource assigned to each encoded item, which we denote by ${{\displaystyle \overline{J}}}_{\text{optimal}}$. Using this conjecture (which we have not yet been able to prove), we are able to efficiently compute predictions at an arbitrary set size. Simulation results show that the model then predicts that both K_{optimal} and ${{\displaystyle \overline{J}}}_{\text{optimal}}$ depend on set size (Figure 5D, left) and produces response data that are qualitatively similar to human data (Figure 5D, right).
Discussion
Summary
Descriptive models of visual working memory (VWM) have evolved to a point where there is little room for improvement in how well they account for experimental data. Nevertheless, the basic finding that VWM precision depends on set size still lacks a principled explanation. Here, we examined a normative proposal in which expected task performance is traded off against the cost of spending neural resource on encoding. We used this principle to construct a resourcerational model for "local" VWM tasks and found that set size effects in this model are fully mediated by the probing probabilities of the individual items; this is consistent with suggestions from earlier empirical work (Emrich et al., 2017; Palmer et al., 1993). From the perspective of our model, the interpretation is that as more items are added to a task, the relevance of each individual item decreases, which makes it less costefficient to spend resource on its encoding. We also found that in this model it is sometimes optimal to encode only a subset of taskrelevant items, which implies that resource rationality could serve as a principled bridge between resource and slotbased models of VWM. We tested the model on data from nine previous delayedestimation experiments and found that it accounts well for effects of both set size and probing probability, despite having relatively few parameters. Moreover, it accounts for a nonmonotonicity that appears to exist between set size and the total amount of resource that subjects invest in item encoding. The broader implication of our findings is that VWM limitations – and cognitive limitations in general – may be driven by a mechanism that minimizes a cost, instead of by a fixed constraint on available encoding resource.
Limitations
Our theory makes a number of assumptions that need further investigation. First, we have assumed that the expected behavioral cost decreases indefinitely with the amount of invested resource, such that in the limit of infinite resource there is no encoding error and no behavioral cost. However, encoding precision in VWM is fundamentally limited by the precision of the sensory input, which is itself limited by irreducible sources of neural noise – such as Johnson noise and Poisson shot noise (Faisal et al., 2008; Smith, 2015) – and suboptimalities in early sensory processing (Beck et al., 2012). One way to incorporate this limitation is by assuming that there is a resource value ${{\displaystyle \overline{J}}}_{\text{input}}$ beyond which the expected behavioral cost no longer decreases as a function of $\overline{J}$. In this variant, ${{\displaystyle \overline{J}}}_{\text{input}}$ represents the quality of the input and ${{\displaystyle \overline{J}}}_{\text{optimal}}$ will never exceed this value, because any additional resource would increase the expected neural cost without decreasing the expected behavioral cost.
Moreover, our theory assumes that there is no upper limit on the total amount of resource available for encoding: cost is the only factor that matters. However, as the brain is a finite entity, the total amount of resource must obviously have an upper limit. This limit can be incorporated by optimizing ${J}_{\text{optimal}}$ under the constraint ${\sum}_{i=1}^{N}}\text{\hspace{0.17em}}{{\displaystyle \overline{J}}}_{\text{optimal},i}\le {{\displaystyle \overline{J}}}_{max$, where ${{\displaystyle \overline{J}}}_{\text{max}}$ represents the maximum amount of resource that can be invested. While an upper limit certainly exists, it may be much higher than the average amount of resource needed to encode information with the same fidelity as the sensory input. If that is the case, then ${{\displaystyle \overline{J}}}_{\text{input}}$ would be the constraining factor and ${{\displaystyle \overline{J}}}_{\text{max}}$ would have no effect.
Similarly, our theory assumes that there is no lower limit on the amount of resource available for encoding. However, there is evidence that taskirrelevant stimuli are sometimes automatically encoded (Yi et al., 2004; Shin and Ma, 2016), perhaps because in natural environments few stimuli are ever completely irrelevant. This would mean that there is a lower limit to the amount of resource spent on encoding. In contradiction to the predictions of our model, such a lower limit would prevent subjects from sometimes encoding nothing at all. For local tasks, such a lower limit can be incorporated by assuming that probing probability p_{i} is never zero.
We have fitted our model only to data from delayedestimation experiments. However, it applies without modification to other local tasks, such as singleprobe change detection (Luck and Vogel, 1997; Todd and Marois, 2004) and singleprobe change discrimination (Klyszejko et al., 2014). Further work is needed to examine how well the model accounts for empirical data of such tasks. Moreover, it should further examine how the theory generalizes to global tasks. One such task could be wholereport change detection; we presented simulation results for this task but the theory remains to be further worked out and fitted to the data.
A final limitation is that our theory assumes that items are uniformly distributed and uncorrelated. Although this is correct for most experimental settings, items in more naturalistic settings are often correlated and can take nonuniform distributions. In such environments, the expected total cost can probably be further minimized by taking into account statistical regularities (Orhan et al., 2014). Moreover, recent work has suggested that even when items are uncorrelated and uniformly distributed, the expected estimation error can sometimes be reduced by using a "chunking" strategy, that is, encoding similar items as one (Nassar et al., 2018). However, as Nassar et al. assumed a fixed total resource and did not take neural encoding cost into account in their optimization, it remains to be seen whether chunking is also optimal in the kind of model that we proposed. We speculate that this is likely to be the case, because encoding multiple items as one will reduce the expected neural cost (fewer items to encode), while the increase in expected behavioral cost will be negligible if the items are very similar. Hence, it seems worthwhile to examine models that combine resource rationality with chunking.
Variability in resource assignment
Throughout the paper, we have assumed that there is variability in resource assignment. Part of this variability is possibly a result of stochastic factors, but part of it may also be systematic – for example, particular colors and orientations may be encoded with higher precision than others (Bae et al., 2014; Girshick et al., 2011). Whereas the systematic component could have a rational basis (e.g. higher precision for colors and orientations that occur more frequently in natural scenes [Ganguli and Simoncelli, 2010; Wei and Stocker, 2015]), this is unlikely to be true for the random component. Indeed, when we jointly optimize $\overline{J}$and τ in Equation 11, we find estimates of τ that consistently approach 0, meaning that any variability in encoding precision is suboptimal under our proposed cost function. One way to reconcile this apparent suboptimality with the otherwise normative theory is to postulate that maintaining exactly equal resource assignment across cortical regions may itself be a costly process; under such a cost, it could be optimal to allow for some variability in resource assignment. Another possibility is that there are unavoidable imperfections in mental inference (Drugowitsch et al., 2016) that make it impossible to compute ${{\displaystyle \overline{J}}}_{\text{optimal}}$ without error, such that the outcome of the computation will vary from trial to trial even when the stimuli are identical.
Experimental predictions of incentive manipulations
In the present study, we have focused on effects of set size and probing probability on encoding precision. However, our theory also makes predictions about effects of incentive manipulations on encoding precision, because such manipulations affect the expected behavioral cost function.
Incentives can be experimentally manipulated in a variety of ways. One method used in at least two previously published delayedestimation experiments is to make the feedback binary ("correct," "error") and vary the value of the maximum error allowed to receive positive feedback (Zhang and Luck, 2011; Nassar et al., 2018). In both studies, subjects in a "low precision" condition received positive feedback whenever their estimation error was smaller than a threshold value of π/3. Subjects in the "high precision" condition, however, received positive feedback only when the error was smaller than π/12 (Zhang and Luck, 2011) or π/8 (Nassar et al., 2018). Neither of the two studies found evidence for a difference in encoding precision between the low and highprecision conditions. At first, this may seem to be at odds with the predictions of our model, as one may expect that it should assign more resource to items in the highprecision condition. To test whether this is the case, we simulated this experimental manipulation using a behavioral cost function c_{behavioral,i}(ε) that maps values of ε smaller than the feedback threshold to 0 and larger values to 1. The results reveal that the model predictions are not straightforward and that it can actually account for the absence of an effect (Figure 6). In particular, the simulation results suggest that the experimental manipulations in the studies by Zhang and Luck and Nassar et al. may not have been strong enough to measure an effect. Indeed, another study has criticized the study by Zhang and Luck on exactly this point and did find an effect when using an experimental design with stronger incentives (Fougnie et al., 2016).
Another method to manipulate incentives is to vary the amount of potential reward across items within a display. For example, Klyszejko and colleagues performed a local change discrimination experiment in which the monetary reward for a correct response depended on which item was probed (Klyszejko et al., 2014). They found a positive relation between the amount of reward associated with an item and response accuracy, which indicates that subjects spent more resource on encoding items with larger potential reward. This incentive manipulation can be implemented by multiplying the behavioral cost function with an itemdependent factor u_{i}, which modifies Equation (11) to ${\overline{J}}_{\text{optimal}},i\left({r}_{i};\lambda ,\tau \right)=\underset{\overline{J}}{\mathrm{a}\mathrm{r}\mathrm{g}\mathrm{m}\mathrm{i}\mathrm{n}}\text{}\left({u}_{i}{p}_{i}{\overline{c}}_{\mathrm{b}\mathrm{e}\mathrm{h}\mathrm{a}\mathrm{v}\mathrm{i}\mathrm{o}\mathrm{r}\mathrm{a}\mathrm{l}}(\overline{J};\tau )+\lambda \overline{J}\right)$. The coefficients u_{i} and p_{i} can be combined into a single "item relevance" coefficient r_{i} = u_{i}p_{i}, and all theoretical results and predictions that we derived for p_{i} now apply to r_{i}.
A difference between the two discussed methods is that the former varied incentives within a trial and the latter across trials. However, both methods can be applied in both ways. A withintrial variant of the experiments by Zhang and Luck (2011) and Nassar et al. (2018) would be a N = 2 task in which one of the items always has a low positive feedback threshold and the other a high one. Similarly, a betweentrial variant of the experiment by Klyszejko et al. (2014) would be to scale the behavioral cost function of items with a factor that varies across trials or blocks, but is constant within a trial. Our model can be used to derive predictions for these task variants, which to our knowledge have not been previously reported in the published literature.
Neural mechanisms and timescale of optimization
Our results raise the question of what neural mechanism could implement the optimal allocation policy that forms the core of our theory. Some form of divisive normalization (Bays, 2014; Carandini and Heeger, 2012) would be a likely candidate, which is already a key operation in neural models of attention (Reynolds and Heeger, 2009) and visual working memory (Bays, 2014; Wei et al., 2012). The essence of this mechanism is that it lowers the gain when set size is larger, without requiring explicit knowledge of the set size prior to the presentation of the stimuli. Consistent with the predictions of this theory, empirical work has found that the neural activity associated with the encoding of an item decreases with set size, as observed in for example the lateral intraparietal cortex (Churchland et al., 2008; Balan et al., 2008) and superior colliculus (Basso and Wurtz, 1998). Moreover, the work by Bays (2014) has shown that a modified version of divisive normalization can account for the nearoptimal distribution of resources across items with unequal probing probabilities. As set size effects in our model are mediated by probing probability, its predicted set size effects can probably be accounted for by a similar mechanism.
Another question concerns the timescale at which the optimization takes place. In all experimental data that we considered here, the only factors that changed from trial to trial were set size (E1E7) and probing probability (E7E9). When we fitted the model, we assumed that the expected total cost in these experiments was minimized on a trialbytrial basis: whenever set size or probing probability changed from one trial to the next, the computation of ${J}_{\text{optimal}}$ followed this change. This assumption accounted well for the data and, as discussed above, previous work has shown that divisive normalization can accommodate trialbytrial changes in set size and probing probability. However, can the same mechanism also accommodate changes in the optimal resource policy changes driven by other factors, such as the behavioral cost function, c_{behavioral}(ε)? From a computational standpoint, divisive normalization is a mapping from an input vector of neural activities to an output vector, and the shape of this mapping depends on the parameters of the mechanism (such as gain, weighting factors, and a power on the input). As the mapping is quite flexible, we expect that it can accommodate a nearoptimal allocation policy for most experimental conditions. However, topdown control and some form of learning (e.g. reinforcement learning) are likely required to adjust the parameters of the normalization mechanism, which would prohibit instantaneous optimality after a change in the experimental conditions.
Neural prediction
The total amount of resource that subjects spend on item encoding may vary nonmonotonically with set size in our model. At the neural level, this translates to a prediction of a nonmonotonic relation between populationlevel spiking activity and set size. We are not aware of any studies that have specifically addressed this prediction, but it can be tested using neuroimaging experiments similar to previously conducted experiments. For example, Balan et al. used singleneuron recording to estimate neural activity per item for set sizes 2, 4, and 6 in a visual search task (Balan et al., 2008). To test for the existence of the predicted nonmonotonicity, the same recoding techniques can be used in a VWM task with a more finegrained range of set sizes. Even though it is practically impossible to directly measure populationlevel activity, reasonable estimates may be obtained by multiplying singleneuron recordings with set size (under the assumption that an increase in resource translates to an increase in firing rate and not an increase of neurons used to encode an item). A similar method can also assess the relation between an item’s probing probability and the spiking activity related to its neural encoding.
Extensions to other domains
Our theory might apply beyond working memory tasks. In particular, it has been speculated that the selectivity of attention arises from a need to balance performance against the costs associated with spiking (Pestilli and Carrasco, 2005; Lennie, 2003). Our theory provides a normative formalism to test this speculation and may thus explain set size effects in attention tasks (Lindsay et al., 1968; Shaw, 1980; Ma and Huang, 2009).
Furthermore, developmental studies have found that that working memory capacity estimates change with age (Simmering and Perone, 2012; Simmering, 2012). Viewed from the perspective of our proposed theory, this raises the question of why the optimal tradeoff between behavioral and neural cost would change with age. A speculative answer is that a subject's coding efficiency – formalized by the reciprocal of parameter α in Equation 7 – may improve during childhood: an increase in coding efficiency reduces the neural cost per unit of precision, which shifts the optimal amount of resource to use for encoding to larger values. Neuroimaging studies might provide insight into whether and how coding efficiency changes with age, for example by estimating the amount of neural activity required per unit of precision in memory representations.
Broader context
Our work fits into a broader tradition of normative theories in psychology and neuroscience (Table 4). The main motivation for such theories is to reach a deeper level of understanding by analyzing a system in the context of the ecological needs and constraints under which it evolved. Besides work on idealobserver decision rules (Green and Swets, 1966; KordingKörding, 2007; Geisler, 2011; Shen and Ma, 2016) and on resourcelimited approximations to optimal inference (Gershman et al., 2015; Griffiths et al., 2015; Vul and Pashler, 2008; Vul, 2009), normative approaches have also been used at the level of neural coding. For example, properties of receptive fields (Vincent et al., 2005; Liu et al., 2009; Olshausen and Field, 1996), tuning curves (Attneave, 1954; Barlow, 1961; Ganguli and Simoncelli, 2010), neural architecture (Cherniak, 1994; Chklovskii et al., 2002), receptor performance (Laughlin, 2001), and neural network modularity (Clune et al., 2013) have been explained as outcomes of optimization under either a cost or a hard constraint (on total neural firing, sparsity, or wiring length), and are thus mathematically closely related to the theory presented here. However, a difference concerns the timescale at which the optimization takes place: while optimization in the context of neural coding is typically thought to take place at the timescale over which the statistics of the environment change or a developmental timescale, the theory that we presented here could optimize on a trialbytrial basis to follow changes in task properties.
We already mentioned the informationtheory models of working memory developed by Chris R. Sims et al. A very similar framework has been proposed by Chris A. Sims in behavioral economics, who used information theory to formalize his hypothesis of "rational inattention," that is, the hypothesis that consumers make optimal decisions under a fixed budget of attentional resources that can be allocated to process economic data (Sims, 2003). The model presented here differs from these two approaches in two important ways. First, similar to early models of visual working memory limitations, they postulate a fixed total amount of resources (formalized as channel capacity), which is a constraint rather than a cost. Second, even if it had been a cost, it would have been the expected value of a log probability ratio. Unlike neural spike count, a log probability ratio does not obviously map to a biologically meaningful cost on a singletrial level. Nevertheless, recent work has attempted to bridge rational inattention and attention in a psychophysical setting (Caplin et al., 2018).
Materials and methods
Data and code sharing
Request a detailed protocolData from experiments E1E7 (Table 1) and Matlab code for model fitting and simulations are available at http://dx.doi.org/10.5061/dryad.nf5dr6c.
Statistical analyses
Request a detailed protocolBayesian ttests were performed using the JASP software package (JASP Team, 2017) with the scale parameter of the Cauchy prior set to its default value of 0.707.
Model fitting
Request a detailed protocolWe used a Bayesian optimization method (Acerbi and Ma, 2017) to find the parameter vector $\theta =\left\{\beta ,\lambda ,\tau \right\}$ that maximizes the log likelihood function, ${\sum}_{i=1}^{n}}\text{\hspace{0.17em}}logp\left({\epsilon}_{i};{p}_{i},\theta \right)$, where n is the number of trials in the subject’s data set, ε_{i} the estimation error on the i^{th} trial, and p_{i} the probing probability of the probed item on that trial. To reduce the risk of converging into a local maximum, initial parameter estimates were chosen based on a coarse grid search over a large range of parameter values. The predicted estimation error distribution for a given parameter vector θ and probing probability p_{i} was computed as follows. First, $\overline{J}}_{\text{optimal}$ was computed by applying Matlab's fminsearch function to Equation 11. Thereafter, the gamma distribution over J (with mean ${{\displaystyle \overline{J}}}_{\text{optimal}}$ and shape parameter τ) was discretized into 50 equalprobability bins. The predicted (Von Mises) estimation error distribution was then computed under the central value of each bin. Finally, these 50 predicted distributions were averaged. We verified that increasing the number of bins used in the numerical approximation of the integral over J did not substantially affect the results.
Model comparison using crossvalidation
Request a detailed protocolIn the crossvalidation analysis, we fitted the models in the same way as described above, but using only 80% of the data. We did this five times, each time leaving out a different subset of 20% of the data (in the first run we left out trials 1, 6, 11; in the second run we left out trials 2, 7, 12, etc.). At the end of each run, we used the maximumlikelihood parameter estimates to compute the log likelihood of the 20% of trials that were left out. These log likelihood values were then combined across the five runs to give an overall crossvalidated log likelihood value for each model.
Appendix 1
Relation between Fisher information J and concentration parameter κ
As we are only considering stimuli with circular domains, we assume that memory encoding errors follow a Von Mises distribution with a concentration parameter κ,
where I_{0} is the modified Bessel function of the first kind of order 0. We measure encoding precision as Fisher information, J, which measures the performance of the best possible unbiased decoder. Substituting Equation 15 into the definition of Fisher information, we find that J and κ are onetoone related through
Encoding precision J is a monotonically increasing function of κ and therefore invertible. However, the inverse is not analytic, so we use numerical inversion to compute the mapping from J to κ when fitting models.
Mathematical proofs of some properties of the resourcerational model for local tasks
In this section, we prove three properties of the general model that we presented for "local" tasks, that is, tasks in which responses depend on a single item. This model is characterized by Equation 11,
where $\overline{J}\ge 0$, ${p}_{i}\in \left[0,1\right]$, and $\lambda \ge 0$, and we left out the dependence on the parameter τ for notational convenience. We will also use the derivative of the local expected total cost,
where ${{\displaystyle \overline{c}}}_{\text{behavioral}}$ is the derivative of the expected behavioral cost.
We will now prove that the following three claims hold under rather general assumptions about the shape of the expected behavioral cost function in this model:
Claim 1. When neural coding is costly (λ >0), it is optimal to encode items with a finite amount of resource;
Claim 2. It is sometimes optimal not to encode a taskrelevant item;
Claim 3. When each item is equally likely to be probed, p_{i} = 1/N, the optimal amount of resource per item decreases with set size.
Assumptions about the expected behavioral cost
We construct our proofs under two intuitive and general assumptions about the expected behavioral cost function ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\overline{J}\right)$:
Assumption 1. Expected behavioral cost is a monotonically decreasing function of resource: whenever more resource is invested, the expected behavioral cost is lower. This means that ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\overline{J}\right)\le 0$ for all $\overline{J}$.
Assumption 2. A law of diminishing returns: when adding a bit of extra resource, the resulting decrease in ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\overline{J}\right)$ is lower in magnitude when $\overline{J}$ is higher. This means that ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)$ is monotonically increasing, that is, ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}\mathrm{\prime}}\left(\overline{J}\right)$ for all $\overline{J}$. As a consequence, ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)$ takes its lowest value at $\overline{J}=0$ and its largest as $\overline{J}\to \infty $.
Both assumptions are satisfied by the behavioral cost function that we used for fitting human data, namely c_{behavioral}(ε;β)=ε^{β}. Examples of the expected behavioral cost function under this choice and its first and second derivative are presented in Appendix 1—figure 1.
Three scenarios
We now return to the problem of calculating ${{\displaystyle \overline{J}}}_{\text{optimal}}$, Equation 17. We are interested in the value $\overline{J}\in [0,\infty )$ that minimizes the expected total cost, ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)$. We separately consider the following three scenarios: the minimum lies on the left boundary (0), on the right boundary (∞), or in between.
Scenario 1: ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)$ is monotonically decreasing across the domain of $\overline{J}$, so ${{\displaystyle \overline{J}}}_{\text{optimal}}\to \infty $.
When does this happen? The monotonic decrease means that ${\overline{c}}_{\text{total}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\le 0$ for all $\overline{J}$, or equivalently, ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\le \lambda$ for all $\overline{J}$. As we assume ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)$ to be monotonically increasing (Assumption 2), its largest value is attained at $\overline{J}\to \infty $. Therefore, ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\le \lambda$ is equivalent to ${p}_{i}{{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\infty \right)\le \lambda $, or (using Assumption 1) ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\mathrm{\infty}\right)\ge \lambda .$ This means that it is optimal to invest infinite resource when p_{i} exceeds a critical value p_{∞}:
The condition p_{i} ≥p∞ is satisfied when λ = 0. This makes sense: when neural cost plays no role, there is no reason not to invest more. Other than that, the condition will rarely if ever be satisfied, as every expected behavioral cost function that we can think of has the property ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\mathrm{\infty}\right)=0$: as the amount of invested resource approaches infinity, there is no behavioral benefit in investing more resource (note that p_{∞} has a domain [0,∞), not [0,1]). Therefore, unless neural cost plays no role, we do not expect it to be optimal to invest an infinite amount of resource in an item.
In tasks where p_{i} is onetoone related to set size, the above result can be reformulated in terms of set size. In particular, when probing probabilities are equal, ${p}_{i}={\scriptscriptstyle \frac{1}{N}}$, the above result implies that there exists a set size N_{∞} (in general not an integer) below which it is optimal to invest infinite resource in each item:
Scenario 2: ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)$ is monotonically increasing across the domain of $\overline{J}$, so ${{\displaystyle \overline{J}}}_{\text{optimal}}=0$.
The monotonic increase means that ${\overline{c}}_{\text{total}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\ge 0$ for all $\overline{J}$, or equivalently, ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\ge \lambda$ for all $\overline{J}$. As we assume ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)$ to be monotonically increasing (Assumption 2), its smallest value is attained at $\overline{J}\to \infty $. Therefore, ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)\ge \lambda$ is equivalent to ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(0\right)\ge \lambda$, or (using Assumption 1) ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(0\right)\le \lambda$. This means that it is optimal to invest no resource when p_{i} is smaller than or equal to a critical value p_{0}:
A similar condition was derived in our earlier work (de Silva and Ma, 2018) for the case of a fixed total amount of resource (hard constraint).
The condition p_{i}≤ p_{0} is satisfied when p_{i} = 0. This makes sense: when an item never gets probed, one should not invest any resource. More generally, when probing probability is sufficiently low, the behavioral cost function is sufficiently shallow at 0, and neural cost is sufficiently important, it is not worth investing any resource on encoding. The expression for p_{0} also makes clear that the optimal amount of resource is never 0 when the slope of the behavioral cost function at 0 approaches −∞.
In tasks where p_{i} is onetoone related to set size, the above result can be reformulated in terms of set size. In particular, when probing probabilities are equal, ${p}_{i}={\scriptscriptstyle \frac{1}{N}}$, the above result implies that there exists a set size N_{0} (in general not an integer) beyond which it is optimal to not invest any resource in any item:
Intuitively, this means that when set size is too large, the chances of success are too low and one should not even try.
Scenario 3: ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)$ has a stationary point, so ${{\displaystyle \overline{J}}}_{\text{optimal}}$ is finite and nonzero.
We will now consider the remaining scenario, which is the complement of Scenarios 1 and 2; in particular, we can take λ >0 and p_{i} >0. The stationary point of ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)$ will always be a minimum, as the second derivative ${\overline{c}}_{\text{total}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}\mathrm{\prime}}\left(\overline{J}\right)$ is equal to ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}\mathrm{\prime}}\left(\overline{J}\right)$, which is always positive (Assumption 2). At the minimum, we have ${{\displaystyle \overline{c}}}_{\text{total}}\left(\overline{J}\right)=0$, from which it follows that $\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)=\frac{a}{{p}_{i}$ at the minimum. As the lefthand side is monotonically increasing as a function of $\overline{J}$ (Assumption 2), the minimum is either a single point or a single interval, but there cannot be multiple disjoint minima. Graphically, this equation describes the intersection between ${\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)$, which is a monotonically increasing function, and a flat line at a value ${\scriptscriptstyle \frac{\lambda}{{p}_{i}}}$ (Appendix 1—figure 2). The value of at which this intersection occurs necessarily increases with p_{i}.
Three regimes for probing probability
So far, we have assumed a given probing probability p_{i}. Now suppose that for a given ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\overline{J}\right)$ and a given λ, we increase p_{i} from 0 to 1:
The first regime is p_{i}≤ p_{0}. There, Scenario 2 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}=0$: the item does not get encoded at all.
The second regime is p_{0} <p_{i} < p_{∞}; there, Scenario 3 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}$ monotonically increases with p_{i}.
The third regime is p_{i} ≥p_{∞}. There, Scenario 1 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}=\infty $: the item gets encoded with infinite resource.
Even though not all regimes might exist for every parameter combination, the model generally predicts that there is a regime in which ${{\displaystyle \overline{J}}}_{\text{optimal}}$ increases monotonically with p_{i} (Figure 1D).
Three regimes for set size
We can similarly examine the experimentally important special case of equal probing probabilities, ${p}_{i}={\scriptscriptstyle \frac{1}{N}}$:
The first regime is N ≤ N_{∞}. There, Scenario 1 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}=\infty $: all items are encoded with infinite resource.
The second regime is N_{∞}<N <N_{0}. There, Scenario 3 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}$ monotonically decreases with N.
The third regime is N≥ N_{0}. There, Scenario 2 applies and ${{\displaystyle \overline{J}}}_{\text{optimal}}=0$: no items are encoded at all.
Even though not all regimes might exist for every parameter combination, the model generally predicts that there is a regime in which ${{\displaystyle \overline{J}}}_{\text{optimal}}$ decreases monotonically with N (Figure 1E).
Conclusion
In conclusion, given Equation (17) and two additional assumptions, we have proven the following:
Investing infinite resource in an item is only optimal when ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\mathrm{\infty}\right)\ge \lambda$. In practice, this might only happen when neural cost is unimportant (λ = 0). This proves Claim 1.
Investing no resource in an item is optimal when ${p}_{i}{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(0\right)\le \lambda$. This can happen even when the probing probability p_{i} is nonzero. This proves Claim 2.
${{\displaystyle \overline{J}}}_{\text{optimal}}$ is a monotonically increasing function of p_{i}. In particular, if ${p}_{i}={\scriptscriptstyle \frac{1}{N}}$, then ${{\displaystyle \overline{J}}}_{\text{optimal}}$ is a monotonically decreasing function of N. This proves Claim 3.
All three results hold more generally than we have shown here: we can replace the neural cost term $\lambda \overline{J}$ in Equation (17) by any function ${c}_{\text{neural}}\left(\overline{J}\right)$ whose derivative is positive and monotonically increasing. The proofs proceed along the same lines (see below).
Special case: fixedprecision model
For the fixedprecision model (variableprecision model with $\tau \downarrow 0$), Equation (12) in the main text takes the form
We wish to evaluate ${p}_{0}\equiv \frac{\lambda}{\left{{\displaystyle \overline{c}}}_{\text{behavioral}}\left(0\right)\right}$. First, we evaluate the derivative of ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(J\right)$ using the chain rule:
Using Equation (15), the first factor is
where in the last line we used ${I}_{0}\left(\kappa \right)={I}_{1}\left(\kappa \right)$ (see Eq. 9.6.27 in Abramowitz and Stegun [1972]).
We next evaluate the second factor in Equation (19) using Equation (16):
where in the third equality, we used $\frac{d}{d\kappa}(\kappa {I}_{1}(\kappa ))=\kappa {I}_{0}(\kappa )$ (see Eq. 9.6.28 in [Abramowitz & Stegun, 1972]). We now combine Equation (20) and Equation (21) into Equation (19) and the result in turn in the expression for p_{0}. We also realize that the limit $J\downarrow 0$ is, using Equation (16), equivalent to the limit 0. Putting everything together, we find
${p}_{0}=\lambda \underset{\kappa \downarrow 0}{{\displaystyle lim}}\text{\hspace{0.17em}}\left\frac{\kappa \left(1{\scriptscriptstyle \frac{{I}_{1}{{\displaystyle \left(\kappa \right)}}^{2}}{{I}_{0}{{\displaystyle \left(\kappa \right)}}^{2}}}\right)}{{\scriptscriptstyle \frac{1}{2\pi}}{\displaystyle \underset{\pi}{\overset{\pi}{\int}}}\text{\hspace{0.17em}}{c}_{\text{behavioral}}\left(\epsilon \right)\left({\scriptscriptstyle \frac{{I}_{1}\left(\kappa \right)}{{I}_{0}{{\displaystyle \left(\kappa \right)}}^{2}}}{e}^{\kappa cos\epsilon}+{\scriptscriptstyle \frac{cos\epsilon}{{I}_{0}\left(\kappa \right)}}{e}^{\kappa}\right)d\epsilon}\right=0$.
We conclude that in our theory for delayedestimation, assuming the expected behavioral cost function from the fixedprecision model, it is only optimal to invest no resource at all into an item when that item has zero probability of being probed.
Generalization to other neural cost functions
So far, we have assumed that the expected neural cost is linear in resource, Equation (8). Relaxing this assumption, Equation (17) for local tasks becomes
The derivative of the local expected total cost becomes
The three claims above still hold if we modify the two assumptions to
Assumption 1’. $\frac{{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}{{\overline{c}}_{\text{neural}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}\le 0$ for all $\overline{J}$.
Assumption 2’. $\frac{{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}{{\overline{c}}_{\text{neural}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}$ is monotonically increasing for all $\overline{J}$.
The proofs are completely analogous, with ${{\displaystyle \overline{c}}}_{\text{behavioral}}\left(\overline{J}\right)$ replaced by $\frac{{\overline{c}}_{\text{behavioral}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}{{\overline{c}}_{\text{neural}}{\phantom{\rule{thinmathspace}{0ex}}}^{\mathrm{\prime}}\left(\overline{J}\right)}$.
Optimal decision rule for the change detection task
In our simulation of the change detection task, we assume that observers use a Bayesian decision rule. This rule is to report "change" whenever the posterior ratio of change presence over change absence exceeds 1,
where x and y denote the vectors of noisy measurements of the items in the first and second displays, respectively. Under the Von Mises noise assumption, and assuming a flat prior on change presence, this decision rule evaluates to (Keshvari et al., 2013)
where κ_{x,i} and κ_{y,i} denote the concentration parameters of the Von Mises distributions associated with the observations of the items at the i^{th} location in the first and second displays, respectively. The predicted probability of a correct response for a given resource vector, $p\left(\text{error}\overline{J}\right)$, is not analytic, but can easily be computed using Monte Carlo simulations.
Data availability
Data from experiments E1E7 (Table 1) and Matlab code for model fitting and simulations are available at http://dx.doi.org/10.5061/dryad.nf5dr6c.

Data from: A resourcerational theory of set size effects in human visual working memoryAvailable at Dryad Digital Repository under a CC0 Public Domain Dedication.
References

BookHandbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables, Vol. 55New York: Dover Publications.

Clear evidence for item limits in visual working memoryCognitive Psychology 97:79–97.https://doi.org/10.1016/j.cogpsych.2017.07.001

A new look at the statistical model identificationIEEE Transactions on Automatic Control 19:716–723.https://doi.org/10.1109/TAC.1974.1100705

RETRACTED ARTICLE: the plateau in mnemonic resolution across large set sizes indicates discrete resource limits in visual working memoryAttention, Perception, & Psychophysics 74:891–910.https://doi.org/10.3758/s1341401202921

Some informational aspects of visual perceptionPsychological Review 61:183–193.https://doi.org/10.1037/h0054663

An energy budget for signaling in the grey matter of the brainJournal of Cerebral Blood Flow & Metabolism 21:1133–1145.https://doi.org/10.1097/0000464720011000000001

Sensory Communication217–234, Possible principles underlying the transformation of sensory messages, Sensory Communication.

Modulation of neuronal activity in superior colliculus by changes in target probabilityThe Journal of Neuroscience 18:7519–7534.https://doi.org/10.1523/JNEUROSCI.181807519.1998

Noise in neural populations accounts for errors in working memoryJournal of Neuroscience 34:3632–3645.https://doi.org/10.1523/JNEUROSCI.320413.2014

Memory for visual motionJournal of Experimental Psychology: Human Perception and Performance 23:353–369.https://doi.org/10.1037/00961523.23.2.353

Normalization as a canonical neural computationNature Reviews Neuroscience 13:51–62.https://doi.org/10.1038/nrn3136

Component placement optimization in the brainThe Journal of Neuroscience 14:2418–2427.https://doi.org/10.1523/JNEUROSCI.140402418.1994

Cognitive cost as dynamic allocation of energetic resourcesFrontiers in Neuroscience 9:289.https://doi.org/10.3389/fnins.2015.00289

Decisionmaking with multiple alternativesNature Neuroscience 11:693–702.https://doi.org/10.1038/nn.2123

The evolutionary origins of modularityProceedings of the Royal Society B: Biological Sciences 280:20122863.https://doi.org/10.1098/rspb.2012.2863

The same type of visual working memory limitations in humans and monkeysJournal of Vision 15:13–18.https://doi.org/10.1167/15.16.13

Visual shortterm memory compared in rhesus monkeys and humansCurrent Biology 21:975–979.https://doi.org/10.1016/j.cub.2011.04.031

Attention mediates the flexible allocation of visual working memory resourcesJournal of Experimental Psychology: Human Perception and Performance 43:1454–1465.https://doi.org/10.1037/xhp0000398

Strategic tradeoffs between quantity and quality in working memoryJournal of Experimental Psychology: Human Perception and Performance 42:1231–1240.https://doi.org/10.1037/xhp0000211

Variability in the quality of visual working memoryNature Communications 3:1229.https://doi.org/10.1038/ncomms2237

Implicit encoding of prior probabilities in optimal neural populationsAdvances in Neural Information Processing Systems 2010:658–666.

Contributions of ideal observer theory to vision researchVision Research 51:771–781.https://doi.org/10.1016/j.visres.2010.09.027

Rational use of cognitive resources: levels of analysis between the computational and the algorithmicTopics in Cognitive Science 7:217–229.https://doi.org/10.1111/tops.12142

No evidence for an item limit in change detectionPLoS Computational Biology 9:e1002927.https://doi.org/10.1371/journal.pcbi.1002927

Attentional priority determines working memory precisionVision Research 105:70–76.https://doi.org/10.1016/j.visres.2014.09.002

Energy as a constraint on the coding and processing of sensory informationCurrent Opinion in Neurobiology 11:475–480.https://doi.org/10.1016/S09594388(00)002373

The cost of cortical computationCurrent Biology 13:493–497.https://doi.org/10.1016/S09609822(03)001350

Attention and multidimensional discrimination1Perception & Psychophysics 4:113–117.https://doi.org/10.3758/BF03209520

Visual working memory capacity: from psychophysics and neurobiology to individual differencesTrends in Cognitive Sciences 17:391–400.https://doi.org/10.1016/j.tics.2013.06.006

A tutorial on Fisher informationJournal of Mathematical Psychology 80:40–55.https://doi.org/10.1016/j.jmp.2017.05.006

Bayesian inference with probabilistic population codesNature Neuroscience 9:1432–1438.https://doi.org/10.1038/nn1790

Chunking as a rational strategy for lossy data compression in visual working memoryPsychological Review 125:486–511.https://doi.org/10.1037/rev0000101

What limits working memory capacity?Psychological Bulletin 142:758–799.https://doi.org/10.1037/bul0000046

An interference model of visual working memoryPsychological Review 124:21–59.https://doi.org/10.1037/rev0000044

The adaptive nature of visual working memoryCurrent Directions in Psychological Science 23:164–170.https://doi.org/10.1177/0963721414529144

Measuring the effect of attention on simple visual searchJournal of Experimental Psychology: Human Perception and Performance 19:108–130.https://doi.org/10.1037/00961523.19.1.108

Phenomenology of attention: I. color, location, orientation, and spatial frequencyJournal of Experimental Psychology: Human Perception and Performance 24:261–282.https://doi.org/10.1037/00961523.24.1.261

An information capacity limitation of visual shortterm memoryJournal of Experimental Psychology: Human Perception and Performance 40:2214–2242.https://doi.org/10.1037/a0037744

BookIdentifying attentional and decisionmaking components in information processingIn: Nickerson R. S, editors. Attention and Performance VIII. Hillsdale: Erlbaum. pp. 277–296.

A detailed comparison of optimality and simplicity in perceptual decision makingPsychological Review 123:452–480.https://doi.org/10.1037/rev0000028

Toward a rational and mechanistic account of mental effortAnnual Review of Neuroscience 40:99–124.https://doi.org/10.1146/annurevneuro072116031526

Working memory capacity as a dynamic processFrontiers in Psychology 3:567.https://doi.org/10.3389/fpsyg.2012.00567

The development of visual working memory capacity during early childhoodJournal of Experimental Child Psychology 111:695–707.https://doi.org/10.1016/j.jecp.2011.10.007

Implications of rational inattentionJournal of Monetary Economics 50:665–690.https://doi.org/10.1016/S03043932(03)000291

An ideal observer analysis of visual working memoryPsychological Review 119:807–830.https://doi.org/10.1037/a0029856

The Poisson shot noise model of visual shortterm memory and choice response time: normalized coding by neural population sizeJournal of Mathematical Psychology 66:41–52.https://doi.org/10.1016/j.jmp.2015.03.007

Factorial comparison of working memory modelsPsychological Review 121:124–149.https://doi.org/10.1037/a0035234

Variability in encoding precision accounts for visual shortterm memory limitationsProceedings of the National Academy of Sciences 109:8780–8785.https://doi.org/10.1073/pnas.1117465109

Fechner's law in metacognition: A quantitative model of visual working memory confidencePsychological Review 124:197–214.https://doi.org/10.1037/rev0000060

Is the early visual system optimised to be energy efficient?Network: Computation in Neural Systems 16:175–190.https://doi.org/10.1080/09548980500290047

Measuring the crowd withinPsychological Science 19:645–647.https://doi.org/10.1111/j.14679280.2008.02136.x

Explaining human multiple object tracking as resourceconstrained approximate inference in a dynamic probabilistic modelAdvances in Neural Information Processing Systems 22:1–9.

A bayesian observer model constrained by efficient coding can explain 'antiBayesian' perceptsNature Neuroscience 18:1509–1517.https://doi.org/10.1038/nn.4105

A detection theory account of change detectionJournal of Vision 4:11–35.https://doi.org/10.1167/4.12.11

The number and quality of representations in working memoryPsychological Science 22:1434–1441.https://doi.org/10.1177/0956797611417006
Decision letter

Stephanie PalmerReviewing Editor; University of Chicago, United States

Timothy E BehrensSenior Editor; University of Oxford, United Kingdom
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Ecological rationality in human working memory and attention" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Timothy Behrens as the Senior Editor. The following individual involved in review of your submission has agreed to reveal her identity: Jacqueline Gottlieb (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
In this manuscript, van den Berg and Ma propose a normative theory of a phenomenon called setsize effect – the fact that attention and working memory performance degrades as a function of the number of items. The setsize effect is fundamental in attention and working memory research, and many previous models have successfully described it by using ad hoc assumptions about the relation between precision and setsize. The authors' main innovation is to explain this relation in a normative framework. Specifically, they propose that the degradation in precision as a function of set size reflects a tradeoff between the benefits of precise encoding and the neural costs that this encoding demands. The authors develop the model, use it to fit several existing data sets, and offer an extensive discussion of the limitations of their model and its relation to the previous theoretical and empirical literatures.
The paper incorporates a genuinely new idea and is clearly written and quite thorough in its analysis and discussion. The Results section details a substantive amount of work, in which the model was fit to data from several experiments where working memory precision was measured in the context of local delayedestimation, global estimation, change detection, change localization and present/absent visual search, and was quantitatively evaluated against other (often equally good) prediction schemes.
Some major revisions are requested to clarify situations in which not all items are stored, to answer questions about the model parameter τ, and to expand the discussion about how the flexibility of resource allocation in the model could be mechanistically realized in the brain.
Essential revisions:
1) The number of encoded items: The paper mentions a possible hard constraint on the number of items encoded. Depending on set size and cost functions, when (if ever) is it "optimal" not to encode some of the items? Some new modeling results should be shown here to shore this up and a longer discussion of this point should be added to the Results and Discussion.
Furthermore, please expand the discussion of the costs of encoding more items versus fewer. Intuitively, it is obvious that it is more costly to encode 8 items than 2, but there are numerous reasons why this could be the case. Are there any experimental data in support of this assumption? In particular, are there any experimental data showing that encoding more items results in higher firing rates (at the population level)?
Michael Frank and his group have recently made a very principled attempt to characterize the nature of the cost in these tasks, in terms of how participants might group or chunk the items. It may be beyond the scope of the current paper to attempt to outline a theory that explains these results, but perhaps readers would enjoy some more elaborate discussion of this issue.
2) τ: Please provide additional model results that show how the fits look when τ = 0. In particular, please show the goodness of fit of the rational model (with τ = 0) as compared to the fit of a model with a "hard" constraint on resources? In general, it was confusing that the theory is initially described in absence of τ, while τ is used for actual fits to data. It would have been easier to understand if the theory had been evaluated in presence of τ, and its effect studied within that theoretical framework.
3) The speed of policy update: An assumption of the model that may be problematic is that people must almost instantaneously optimize their encoding precision when setsizes change unpredictably from trial to trial. The theory predicts that, when a trial contains 2 stimuli and the next trial contains 4 stimuli, the participant instantaneously lowers the encoding precision to the new (near) optimal level. This process sounds pretty demanding itself, especially considering that humans may also play an active role in determining their intrinsic motivation or deciding which items to memorize, which may further slow down the adjustment process. The flexibility in allocating resources that is implied by this model seems to be at odds with the slowness of cognitive control, welldocumented by task switching costs. The authors touch on this point in the very last line of the paper, where they note that divisive normalization can provide a rapid adjustment mechanism. Even though the discussion is already long, it would be good to hear more about this point, and a comparison between a hardwired allocation mechanism and slower but more flexible cognitive control strategies.
Furthermore, in the Frank paper described above, they use a task in which there is binary feedback that depends on the liberal vs. conservative error criterion (and they don't report major differences in performance as a function of this). One might suspect that participants would fail to adapt their policy even in an incentivecompatible version of this task which systematically varied these behavioral costs, and this would present a challenge to the authors' theory as described here. Please add some text to the Discussion addressing this point.
4) Total precision vs. setsize: A novel and interesting prediction of the model is the nonmonotonic relation between setsize and total precision (Figure 3B). Although the authors state that this point requires more empirical documentation, are the model results consistent with a nonmonotonic encoding of target location that was reported by Balan et al. (2008) in monkey area LIP? That study found that, in a covert visual search task with different set sizes, the fidelity of target location encoding by area LIP was higher at set size 4 than at set size 2 or set size 6 (see Figure 5) – a nonlinearity that was puzzling at the time but may gain new significance in light of this paper. Please add some discussion of this point to the manuscript.
https://doi.org/10.7554/eLife.34963.022Author response
Essential revisions:
1) The number of encoded items: The paper mentions a possible hard constraint on the number of items encoded. Depending on set size and cost functions, when (if ever) is it "optimal" not to encode some of the items? Some new modeling results should be shown here to shore this up and a longer discussion of this point should be added to the Results and Discussion.
This is an interesting question, because its answer can possibly provide a principled bridge between slotbased and resource models of VWM. We now address this question in three different places. First, a mathematical analysis of the general conditions under which it is optimal to not encode an item is provided in Appendix 1. Second, the question is addressed in the context of delayedestimation tasks in the new Results section “Is it ever optimal to not encode an item?”.Third, for the changedetection task, the question is addressed in the new section “Predictions for a global task: wholedisplay change detection”.
Furthermore, please expand the discussion of the costs of encoding more items versus fewer. Intuitively, it is obvious that it is more costly to encode 8 items than 2, but there are numerous reasons why this could be the case. Are there any experimental data in support of this assumption? In particular, are there any experimental data showing that encoding more items results in higher firing rates (at the population level)?
For many choices of spike variability, the total precision of a set of stimuli encoded in a neural population is proportional to the trialaveraged neural spiking rate (e.g., Paradiso, 1988; Seung and Sompolinsky, 1993; Ma et al., 2006). Based on this theoretical argument, it is expected that it is more costly (in terms of neural spiking) to encode 8 items compared to 2, if they are encoded with the same precision.
However, it is important to keep in mind that our model does not predict that the total spiking rate will increase with set size, because it generally predicts the precision per item (i.e., spike rate per item) to decrease with set size, which is consistent with physiological evidence (e.g., Churchland et al., 2008; Balan et al., 2008; Basso and Wurtz, 1998). The maximumlikelihood fits suggest that the total amount of invested resource varies nonmonotonically with set size, which predicts that the populationlevel spiking activity also varies nonmonotonically with set size. We are not aware of any work that strongly supports or rejects this prediction (see also our response below to the point about the Balan et al. paper). We address this point as follows in a new discussion section “Neural prediction”.
Michael Frank and his group have recently made a very principled attempt to characterize the nature of the cost in these tasks, in terms of how participants might group or chunk the items. It may be beyond the scope of the current paper to attempt to outline a theory that explains these results, but perhaps readers would enjoy some more elaborate discussion of this issue.
We assume that this comment refers to the recent paper by Nassar, Helmers, and Frank. If we understand correctly, this paper is currently in press, so we base our response to this comment on the preprint that is available on bioRxiv.
We agree that this paper has several connections with our own study, and we now refer to it at two different places. First, in the Introduction:
“Finally, Nassar and colleagues have proposed a normative model in which a strategic tradeoff is made between the number of encoded items and their precision: when two items are very similar, they are encoded as a single item, such that there is more resource available per encoded item (Nassar et al., 2018). […] However, just as in much of the work discussed above, this theory assumes a fixed resource budget for item encoding, which is not necessarily optimal when resource usage is costly.”
And then again in the “Limitations” section in the Discussion:
“A final limitation is that our theory assumes that items are uniformly distributed and uncorrelated. […] Hence, it seems worthwhile to examine models that combine resource rationality with chunking.”
2) τ: Please provide additional model results that show how the fits look when τ = 0. In particular, please show the goodness of fit of the rational model (with τ = 0) as compared to the fit of a model with a "hard" constraint on resources?
We have added this analysis:
“So far, we have assumed that there is random variability in the actual amount of resource assigned to an item. […] Therefore, we will only consider variableprecision models in the remainder of the paper.”
As in the variableprecision model, the optimal amount of resource per item decreases with set size in the equalprecision variant of the rational model:
Since we think that focusing too much on equalprecision results distracts a bit from the main story, we decided not to include this plot in the paper. As we explain in response to a later comment, the difference between the equalprecision and variableprecision models is mainly in the predicted kurtosis (“peakiness”) of the error distribution, not in the variance of these distributions (let alone in how the variance changes with set size). Hence, the equalprecision vs.variableprecision question is orthogonal to our main question.
In general, it was confusing that the theory is initially described in absence of τ, while τ is used for actual fits to data. It would have been easier to understand if the theory had been evaluated in presence of τ, and its effect studied within that theoretical framework.
Sorry, this was indeed confusing. In the rewritten “Theory” section, we explicitly indicate which equations depend on τ, by including it in the function arguments.
3) The speed of policy update: An assumption of the model that may be problematic is that people must almost instantaneously optimize their encoding precision when setsizes change unpredictably from trial to trial. The theory predicts that, when a trial contains 2 stimuli and the next trial contains 4 stimuli, the participant instantaneously lowers the encoding precision to the new (near) optimal level. This process sounds pretty demanding itself, especially considering that humans may also play an active role in determining their intrinsic motivation or deciding which items to memorize, which may further slow down the adjustment process. The flexibility in allocating resources that is implied by this model seems to be at odds with the slowness of cognitive control, welldocumented by task switching costs. The authors touch on this point in the very last line of the paper, where they note that divisive normalization can provide a rapid adjustment mechanism. Even though the discussion is already long, it would be good to hear more about this point, and a comparison between a hardwired allocation mechanism and slower but more flexible cognitive control strategies.
This is an important issue, which we now discuss in the new Discussion section “Neural mechanisms and timescale of optimization”.
Furthermore, in the Frank paper described above, they use a task in which there is binary feedback that depends on the liberal vs. conservative error criterion (and they don't report major differences in performance as a function of this). One might suspect that participants would fail to adapt their policy even in an incentivecompatible version of this task which systematically varied these behavioral costs, and this would present a challenge to the authors' theory as described here. Please add some text to the Discussion addressing this point.
The experiment by Frank et al. used feedback threshold of π/3 (“low precision” condition) and π/8 (“high precision” condition) and found no difference in absolute estimation error between these two conditions. This would be at odds with any model that predicts that encoding precision is higher in the “high precision” condition, which is what one may expect to happen in our model. However, it turns out that the predictions for this experiment are not that straightforward and that the model can actually account for the lack of an effect. The short explanation is that there is a threshold region in which the prediction barely changes as a function of threshold, due to the performance benefit of adding extra resource is almost exactly outdone by the added neural cost. For a more detailed explanation, we refer to the new Figure 6.
This point is now also discussed in a new Discussion section “Experimental predictions of incentive manipulations”.
4) Total precision vs. setsize: A novel and interesting prediction of the model is the nonmonotonic relation between setsize and total precision (Figure 3B). Although the authors state that this point requires more empirical documentation, are the model results consistent with a nonmonotonic encoding of target location that was reported by Balan et al. (2008) in monkey area LIP? That study found that, in a covert visual search task with different set sizes, the fidelity of target location encoding by area LIP was higher at set size 4 than at set size 2 or set size 6 (see Figure 5) – a nonlinearity that was puzzling at the time but may gain new significance in light of this paper. Please add some discussion of this point to the manuscript.
We thank the reviewer for the reference, as we were not aware of that paper. However, after a careful study of the results reported in that paper, we don’t see how the nonmonotonic trend in Figure 3 can be linked to the predicted nonmonotonicity in the total amount of invested resource. The nonmonotonicity in the Balan paper shows that the stimulus identity (target/distractor) can be decoded more accurately from neural data in N=4 trials compared to N=2 and N=6 trials. However, we do not see how decoding accuracy of a single item relates to the total amount of resource invested in all items. Although it would have been nice if the Balan paper backs up the nonmonotonicity prediction, we believe that linking our prediction to their result would be a bit misleading, so we decided to not include this point. (However, if we misunderstood the reviewer’s suggestion, we would of course be happy to have another look at it after some clarification).
Nevertheless, the Balan paper is relevant to our work for other reasons and we now cite it at two different places in the Discussion. First, in the Discussion section about experimental predictions and, second, in the rewritten part about Neural mechanisms (see responses to previous comments).
https://doi.org/10.7554/eLife.34963.023Article and author information
Author details
Funding
Vetenskapsrådet (201500371)
 Ronald van den Berg
Marie SkłodowskaCurie Actions, Cofund (INCA 600398)
 Ronald van den Berg
National Institutes of Health (R01EY020958)
 Wei Ji Ma
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by grant R01EY020958 from the National Institutes of Health, grant 2015–00371 by the Swedish Resarch Council, and grant INCA 600398 by Marie Sklodowska Curie Actions. We thank all authors of the papers listed in Table 1 for making their data available.
Senior Editor
 Timothy E Behrens, University of Oxford, United Kingdom
Reviewing Editor
 Stephanie Palmer, University of Chicago, United States
Publication history
 Received: January 10, 2018
 Accepted: July 28, 2018
 Accepted Manuscript published: August 7, 2018 (version 1)
 Version of Record published: August 27, 2018 (version 2)
Copyright
© 2018, van den Berg et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 2,681
 Page views

 421
 Downloads

 8
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.