Budgeted and Non-budgeted Causal Bandits

Vineet Nair, Vishakha Patil, Gaurav Sinha

Introduction

Causal Bayesian Networks (CBN) [Pea09] have become the popular choice to model causal relationships in many real-world systems such as online advertising, gene interaction networks, brain functional connectivity, etc. The underlying directed acyclic graph (DAG) of a CBN is called its causal graph. The nodes of this graph are labeled by random variables The joint distribution of these random variables factorizes over the graph. representing the underlying system, and edges between these variables capture direct causal relationships. Once the causal graph is known, any external manipulations on the system that forcibly fixes some target variables can be modeled via an operation called intervention. An intervention simulates the effect of such a manipulation of the target variables on other system variables by disconnecting the target variables from their parents A process known as causal surgery. and setting them to the desired value.

Two key questions in causal learning are: 1) learning the causal graph itself, and 2) finding the intervention that optimizes some variable of interest (often called reward variable) assuming that the causal graph is known. In this work, we focus on the second question by modeling the causal learning problem as an extension of the Stochastic Multi-armed Bandit Problem (mab) [Rob52]. The Mab problem is a popular model used to capture decision-making in uncertain environments where a decision-maker is faced with kk choices (called arms) and at each time step the decision-maker has to choose one out of the kk arms (pull an arm). The arm that is pulled gives a reward drawn from an underlying distribution which is unknown to the decision-maker beforehand.

We study the Mab problem with dependencies between the arms modelled via a causal graph. This model, called causal bandits, was studied in the recent works of [BFP15, LLR16, SSDS17, SSK+17, LB18, YHS+18, LB19, LMTY20], where the interventions are modelled as the arms of the bandit and the influence of the arms on the reward is assumed to conform to a known causal graph. In addition to the possible interventions allowed, the set of arms also contains the empty intervention called the observational arm, where the algorithm does not perform any intervention on the causal graph. The goal of a causal bandit algorithm is to learn the intervention that maximizes the reward.

We explain the causal bandit problem with a motivating example from the marketing domain for which a simple causal graph is shown in Figure 1. An e-commerce company sells a product online and makes a profit whenever a customer purchases their product. This corresponds to the green (reward) node labeled Purchase in the causal graph. On every new customer visit, the product webpage is rendered using some values of the blue (intervenable) nodes, chosen from an underlying distribution. For example, a customer might see a large image, small description, no promotions, and comparison with a competing brand. The red nodes capture actions taken by the customer before they make any purchase decision. For example, based on the rendered product page, a customer might want to get more information from the reviews before deciding to add the product to their shopping cart. Note that, while the blue nodes are actionable and can be manipulated to increase the chances of purchase, the red nodes are not directly manipulable and can only be passively observed for any given values of the blue nodes. For example, the company might take an action by always offering a promotion, seeing which customers might decide to skip going through the reviews and directly add the product to their shopping cart. The objective here is to learn the intervention (on blue nodes) that maximizes the chances of the product being purchased.

However, in many situations, interventions are costly [KDV17, KSB17, LKDV18, AKMM20]. Consider the marketing example above where observational data from this graph can be collected via independent customer visits whereas to get an interventional sample, one needs to render a specific page configuration that would require additional expenditure. But recent works suggest that in many scenarios the effect of interventions can be efficiently estimated using observational samples [TP02, BGK+20, Pea09]. Hence, in the causal bandit framework, for a fixed budget, there is a trade-off between the more economical observational arm and the high-cost interventional arm. This is because the observational arm, though less rewarding, aids in the exploration of the possibly high rewarding interventional arms. This motivates the study of observation/intervention trade-off in the budgeted bandit setting.

We study the problem of finding the best intervention in a causal graph in two settings: with and without budget constraints. Further, we study these problems with two objectives that are common in the Mab literature: simple regret minimization and cumulative regret minimization.

Budgeted Setting: In Sections 3 and 4 we consider a class of causal graphs that we call no-backdoor graphs (see Section 3 for the definition). A special instance of the no-backdoor graph class is the parallel graph model defined in [LLR16]: G\mathcal{G} consists of M+1M+1 nodes, X={Y,X1,…,XM}\mathcal{X}=\{Y,X_{1},\ldots,X_{M}\}, and the only edges in G\mathcal{G} are from each XiX_{i} to YY. For this, [LLR16] propose an algorithm called the parallel bandit algorithm (PB-ALG henceforth). We observe that PB-ALG in fact works for the more general class of no-backdoor graphs.

We study the causal bandit problem for no-backdoor graphs in the budgeted bandit setting [TTCRJ12], where a budget BB is specified and the ratio of the cost of the intervention to the cost of observation is γ≥1\gamma\geq 1. The goal of an algorithm is to find the best intervention such that the total cost of arm pulls does not exceed BB. In Section 3, we first study this problem with the goal of minimizing simple regret. Note that PB-ALG does not take into account the cost of interventions and is only optimal in the non-budgeted setting. We show that when γ\gamma is higher than a threshold (unknown to the algorithm), the simple algorithm OBS-ALG that plays the observational arm every time achieves better simple regret in terms of BB than PB-ALG. Next, we propose γ\gamma-NB-ALG (Algorithm 1) which determines this unknown threshold online and successfully manages to trade-off interventions with observations for a specified budget.

In Section 4, we study the cumulative regret minimization (CRM) problem in the above setting and give the CRM-NB-ALG algorithm. CRM-NB-ALG is based on the Fractional−KUBE\mathtt{Fractional-KUBE} algorithm (F−KUBE\mathtt{F-KUBE} henceforth) given in [TTCRJ12] for budgeted bandits with no side-information. CRM-NB-ALG achieves constant regret if the observational arm is the optimal arm and otherwise achieves logarithmic regret which is better than that of F−KUBE\mathtt{F-KUBE} in terms of instance-specific constants.

Non-Budgeted Setting: In Section 5, we study the problem of minimizing the cumulative regret for general causal graphs in the non-budgeted setting. We assume that the distribution of parents of the reward variable for each intervention is known to the algorithm. This assumption though limiting in the practical setting is also made in the recent work of [LMTY20] (which studies the same problem) as well as in the work of [LLR16]. [LMTY20] proposed an algorithm called C-UCB which has a worst-case regret guarantee of O(knT)O(\sqrt{k^{n}T)} where kk is the number of distinct values that each of the nn parents of the reward variable can take. For the same problem, we propose C-UCB-2 (Algorithm 3) and show it has constant expected cumulative regret in terms of instance parameters which is a significant improvement.

Model and Notations

An algorithm for this problem is a sequential decision-making process that at each time tt performs an intervention at∈Aa_{t}\in\mathcal{A} and observes reward Yt∈{0,1}Y_{t}\in\{0,1\}. For each intervention a∈Aa\in\mathcal{A}, where a=do(X=x)a=do(\mathbf{X}=\mathbf{x}), the expected reward of aa is denoted μa=E[Y∣do(X=x)]\mu_{a}=E[Y\mid do(\mathbf{X}=\mathbf{x})]. We study a budgeted as well as a non-budgeted variant of this problem.

Simple regret: Let ALG\mathtt{ALG} be an algorithm for the above problem that outputs arm aBa_{B} when the budget given is BB. Then the simple regret of ALG\mathtt{ALG} with budget BB, denoted r(B)r(B), is

An algorithm whose objective is to minimize the simple regret is a pure-exploration algorithm and its goal is to identify the best arm without having to restrict the number of times a sub-optimal arm may be played using the budget BB. In many applications, we may require that a sub-optimal arm should not be pulled too many times right from the start. This motivates the definition of cumulative regret.

Cumulative regret: Let GALG(B)G_{\mathtt{ALG}}(B) be the expected reward accumulated by algorithm ALG\mathtt{ALG} with budget BB, and let GB=max⁡ALGGALG(B)G_{B}=\max_{\mathtt{ALG}}G_{\mathtt{ALG}}(B). Then, the cumulative regret of an algorithm ALG\mathtt{ALG} with budget BB, denoted RALG(B)R_{\mathtt{ALG}}(B), is

An algorithm for cumulative regret minimization has to carefully trade-off between exploration vs. exploitation. Hence, an algorithm with good simple regret guarantees may not have good cumulative regret guarantees and vice-versa.

2 Non-budgeted Causal Bandits

In the non-budgeted variant of the problem, the cost associated with every intervention is the same, i.e., we can assume γ=1\gamma=1 for all a∈Aa\in\mathcal{A}. In Section 5, we study this problem with the objective of minimizing the expected cumulative regret when the time horizon TT is unknown but finite. The regret notion is defined as in Equation 2 with B=TB=T. Observe that GT=T⋅max⁡a∈AμaG_{T}=T\cdot\max_{a\in\mathcal{A}}\mu_{a}. The cumulative regret of an algorithm ALG\mathtt{ALG} after TT rounds, denoted RALG(T)R_{\mathtt{ALG}}(T), is then defined as

The goal of any algorithm for such a setting is to minimize the expected cumulative regret E[RALG(T)]E[R_{\mathtt{ALG}}(T)] where the expectation is taken over the randomness in the rewards as well as in the algorithm.

Budgeted mab: No-backdoor Graphs

In this section and Section 4, we assume that the interventions are of size 11. Formally, let {X1,…,\{X_{1},\ldots, XM}X_{M}\} ⊆X\subseteq\mathcal{X} be the set of intervenable nodes such that Xi∈{0,1}X_{i}\in\{0,1\} for all i∈[M]i\in[M]. An intervention in this setting is defined as explicitly setting the value of a single node XiX_{i} as either 00 or 11. When not intervened upon, Xi∼Bernoulli(pi)X_{i}\sim\text{Bernoulli}(p_{i}). Hence, we have 2M+12M+1 interventions in total: 2M2M interventions correspond to setting each of the MM variables XiX_{i} to either 00 or 11, denoted do(Xi=0)do(X_{i}=0) and do(Xi=1)do(X_{i}=1) respectively, and the last intervention corresponds to the empty intervention, do()do(). Moreover, we assume that there are no backdoor paths from XiX_{i} to the reward variable. This implies E[Y∣do(Xi=x)]=E[Y∣Xi=x]=μi,xE[Y\mid do(X_{i}=x)]=E[Y\mid X_{i}=x]=\mu_{i,x} (see Section 3.3.1 [Pea00]). We call a causal graph G\mathcal{G} satisfying this property as a no-backdoor graph ( NB\mathtt{NB}-graph).

For ease of notation, we denote the intervention do(Xi=x)do(X_{i}=x) by ai,xa_{i,x} where i∈[M]i\in[M] and x∈{0,1}x\in\{0,1\}, and the empty intervention as a0a_{0}. The set of interventions is then A={ai,x∣i∈[M],x∈{0,1}}⊎{a0}\mathcal{A}=\{a_{i,x}\mid i\in[M],x\in\{0,1\}\}\uplus\{a_{0}\}. The expected reward for the intervention ai,xa_{i,x} and a0a_{0} are μi,x=E[Y∣Xi=x]\mu_{i,x}=E[Y\mid X_{i}=x] and μ0=E[Y]\mu_{0}=E[Y] respectively. Throughout Sections 3 and 4, i∈[M]i\in[M] and x∈{0,1}x\in\{0,1\}. Also, we use aa to denote an intervention in A\mathcal{A} when we do not differentiate between ai,xa_{i,x} and a0a_{0}. We study the budgeted causal bandit problem for no-backdoor graphs. As stated in Section 2, an algorithm for this problem is given as input the graph G\mathcal{G}, the set of intervenable nodes {X1,…,XM}\{X_{1},\ldots,X_{M}\}, a budget BB, and γ\gamma which is the cost for pulling an arm ai,xa_{i,x}. The algorithm does not know pip_{i} for any ii. Note that if γ≥B\gamma\geq B then trivially the algorithm can only make observations.

In Section 3.1, we show that for γ=Ω(1p⋅m(p))\gamma=\Omega(\frac{1}{p\cdot m(\mathbf{p})}) the simple algorithm that plays the observation arm for BB rounds achieves better expected simple regret than PB-ALG. Since pip_{i} for all ii is unknown, the threshold 1p⋅m(p)\frac{1}{p\cdot m(\mathbf{p})} is a priori unknown to an algorithm. Hence, in Section 3.2 we propose an algorithm that estimates this threshold online, and trades-off between interventions and observations dependent on γ\gamma and the threshold to minimize the expected simple regret.

Here, we analyze the simple-regret of the observational algorithm (OBS-ALG) which plays the arm a0a_{0} for all the rounds, and at the end of BB rounds outputs the arm a∈Aa\in\mathcal{A} with the highest empirical mean estimate. The empirical estimate of μi,x\mu_{i,x} is computed as the average of the rewards accrued in those rounds where XiX_{i} was sampled as xx. Theorem 1 shows the dependence of the expected simple regret of OBS-ALG on pp.

The expected simple regret of OBS-ALG with budget BB is O(1pBlog⁡(pMB))O\left(\sqrt{\frac{1}{pB}\log(pMB)}\right).

The proof of Theorem 1 is in Section 8.2. Theorem 1 is proved by crucially leveraging the fact that the arm a0a_{0} aides in the exploration of all the other 2M2M arms, which is the side-information available in NB\mathtt{NB}-graphs. Observe that the guarantee of observational algorithm is better than that of PB-ALG in [LLR16] if γ=Ω(1p⋅m(p))\gamma=\Omega(\frac{1}{p\cdot m(\mathbf{p})}).

2 Observation-Intervention Trade-off

γ\gamma-NB-ALG (Algorithm 1) trades-off between observations and interventions depending on the value of γ\gamma to minimize the expected simple regret. The idea behind γ\gamma-NB-ALG is that if γ\gamma is larger than the threshold 1p⋅m(p)\frac{1}{p\cdot m(\mathbf{p})} then performing only observations gives a better regret (as stated at the end of Section 3.1), whereas if γ\gamma is less than this threshold then the algorithm follows the strategy of PB-ALG by playing the interventions in set AA (see step 11 of γ\gamma-NB-ALG ) for an equal number of times in the remaining rounds. At steps 13-14, the empirical estimates of only the arms in AA are updated. Since p\mathbf{p} and pp are not known a priory, the algorithm has to estimate the threshold online as done in Step 6 of γ\gamma-NB-ALG. Note that at step 6, p^=(p^1…p^M)\widehat{\mathbf{p}}=(\widehat{p}_{1}\ldots\widehat{p}_{M}), where p^i=p^i,1\widehat{p}_{i}=\widehat{p}_{i,1}, and m(p^)m(\widehat{\mathbf{p}}) is defined similar to m(p)m(\mathbf{p}).

In Theorem 2, we bound the expected simple regret of γ\gamma-NB-ALG which depends upon γ\gamma and the value of the threshold.

If γ≥1p⋅m(p)\gamma\geq\frac{1}{p\cdot m(\mathbf{p})} then the expected simple regret of γ\gamma-NB-ALG is O(1pBlog⁡(pMB))O\left(\sqrt{\frac{1}{pB}\log(pMB)}\right), and if γ≤1p⋅m(p)\gamma\leq\frac{1}{p\cdot m(\mathbf{p})} then it is O(γ⋅m(p)Blog⁡MBγ⋅m(p))O\left(\sqrt{\frac{\gamma\cdot m(\mathbf{p})}{B}\log\frac{MB}{\gamma\cdot m(\mathbf{p})}}\right).

The proof of Theorem 2 is in Section 8.3. Observe that the expected simple regret of γ\gamma-NB-ALG is equal to that of PB-ALG if γ≤1p⋅m(p)\gamma\leq\frac{1}{p\cdot m(\mathbf{p})}, and is equal to the that of OBS-ALG if γ>1p⋅m(p)\gamma>\frac{1}{p\cdot m(\mathbf{p})}. For γ=O(1p⋅m(p))\gamma=O(\frac{1}{p\cdot m(\mathbf{p})}) the optimality of the regret up to log factors follows from Theorem 2 in [LLR16] where they show a Ω(γmB)\Omega\left(\sqrt{\frac{\gamma m}{B}}\right) lower bound on the expected simple regret. The lower bound is shown in non-budgeted setting, which translates to Ω(γmB)\Omega\left(\sqrt{\frac{\gamma m}{B}}\right) lower bound in our setting if γ=O(1p⋅m(p))\gamma=O(\frac{1}{p\cdot m(\mathbf{p})}). The experiment 2 in Section 6 shows that the performance of γ\gamma-NB-ALG matches or is better than the performance of PB-ALG for all values of γ\gamma, which validates our theoretical claim.

Cumulative-Regret in No-backdoor Graphs with Budget

At the end of tt rounds CRM-NB-ALG computes μ^i,x(t)\widehat{\mu}_{i,x}(t) and μ^0(t)\widehat{\mu}_{0}(t) which are empirical estimates of μi,x\mu_{i,x} and μ0\mu_{0} respectively, as follows:

Based on this estimates the CRM-NB-ALG computes the weighted UCB estimate μ‾i,x(t)\overline{\mu}_{i,x}(t) and μ‾0(t)\overline{\mu}_{0}(t) for the arms ai,xa_{i,x} and a0a_{0} as follows:

In each round CRM-NB-ALG first ensures arm a0a_{0} is pulled at least β2log⁡T\beta^{2}\log T times (steps 4-5), where β\beta is set as in steps 11-14 and otherwise pulls the arm with the highest weighted UCB estimate (steps 6-8). Ensuring the arm a0a_{0} is pulled at least β2log⁡T\beta^{2}\log T times at the end of TT rounds delicately balances the exploration-exploit trade-off: the causal side-information by pulling the arm a0a_{0} ensuring free exploration of the other 2M2M interventions and the loss experienced in pulling the arm a0a_{0} (if a0a_{0} is the sub-optimal arm). The reason for setting β\beta as in steps 11-14 is explained after Theorem 3, which bounds the expected cumulative regret of CRM-NB-ALG. Observe that CRM-NB-ALG halts once it has exhausted its entire budget BB. Crucially though, the decisions of CRM-NB-ALG do not depend on the budget, i.e., it is budget oblivious. But note that the decisions of the algorithm do take into account the cost of an intervention, i.e., the algorithm is not cost-oblivious.

Before stating Theorem 3, we introduce a few more notations which are used in the theorem. Let vi,x=μi,xγv_{i,x}=\frac{\mu_{i,x}}{\gamma}, and v0=μ0v_{0}=\mu_{0}, and a∗=arg⁡max⁡a∈A{va}a^{*}=\arg\max_{a\in\mathcal{A}}\{v_{a}\}. Further, let Δa=μa∗−μa\Delta_{a}=\mu_{a^{*}}-\mu_{a} and da=va∗−vad_{a}=v_{a^{*}}-v_{a} for each a∈Aa\in\mathcal{A}. Note that there could be a∈Aa\in\mathcal{A} such that Δa<0\Delta_{a}<0.

If a∗=a0a^{*}=a_{0} then the expected cumulative regret of the algorithm is O(1)O(1) and otherwise the expected cumulative regret of the algorithm is of order ∑Δi,x>0Δi,x(max(0,1+8ln⁡B(1di,x2−pi,x3d02))+π23)+Δ0(50ln⁡Bd02+1+π23)\sum_{\Delta_{i,x}>0}\Delta_{i,x}\left(\text{max}\left(0,1+8\ln B\left(\frac{1}{d_{i,x}^{2}}-\frac{p_{i,x}}{3d_{0}^{2}}\right)\right)+\frac{\pi^{2}}{3}\right)+\Delta_{0}\left(\frac{50\ln B}{d_{0}^{2}}+1+\frac{\pi^{2}}{3}\right).

The optimal arm a∗a^{*} is equal to a0a_{0} if the ratio of the expected reward of any intervention to expected reward of a0a_{0} is at most γ\gamma. In particular, if maxi,xμi,xμ0≤γ\frac{\text{max}_{i,x}\mu_{i,x}}{\mu_{0}}\leq\gamma then a∗=a0a^{*}=a_{0} and in this case the expected cumulative regret of CRM-NB-ALG is bounded by a constant. The proof of Theorem 3 is given in Section 8.4. For the value of β\beta set as in steps 11-14, we show that if a∗≠a0a^{*}\neq a_{0} then 89d02≤E[β2]≤50d02\frac{8}{9d_{0}^{2}}\leq E[\beta^{2}]\leq\frac{50}{d_{0}^{2}} (see Lemma 8.6 in Section 8.4). This in particular ensures that if a∗≠a0a^{*}\neq a_{0} then the expected number of pulls of a sub-optimal arm ai,xa_{i,x} is at most max(0,1+8ln⁡B(1di,x2−pi,x3d02))+π23\text{max}\left(0,1+8\ln B\left(\frac{1}{d_{i,x}^{2}}-\frac{p_{i,x}}{3d_{0}^{2}}\right)\right)+\frac{\pi^{2}}{3}. Hence note that, if 1di,x2≥pi,x3d02\frac{1}{d_{i,x}^{2}}\geq\frac{p_{i,x}}{3d_{0}^{2}} then this sub-optimal arm is pulled at most a constant number of times. In Section 6, we show via simulation that CRM-NB-ALG performs much better than F−KUBE\mathtt{F-KUBE} even for small values of γ\gamma. Note that F−KUBE\mathtt{F-KUBE} does not take side information into account.

Cumulative Regret in General Graphs

In this section, we study the non-budgeted version of the causal bandit problem for general graphs (see Section 2.2) with the goal of minimizing the expected cumulative regret. This problem was studied in the recent work of [LMTY20] who gave a UCB based algorithm, called C-UCB, which has a worst-case regret bound of knT\sqrt{k^{n}T} when each of the nn parent nodes of YY (the reward variable) in the graph can take one of kk values. For the same problem, we propose an algorithm called C-UCB-2 , which has constant regret in terms of instance-parameters. Additionally, C-UCB in [LMTY20] takes the time horizon TT as input, but our algorithm C-UCB-2 works for any unknown (but finite) time horizon.

if Ny,t≥1N_{\mathbf{y},t}\geq 1 and otherwise μ^y(t)=0\widehat{\mu}_{\mathbf{y}}(t)=0. The algorithm also computes the empirical estimate μ^a(t)\widehat{\mu}_{a}(t) and the UCB estimate μ‾a(t)\overline{\mu}_{a}(t) for all aa using μ^y(t)\widehat{\mu}_{\mathbf{y}}(t) at the end of every round as follows:

The quantity log⁡(knt2/2)tζa\sqrt{\frac{\log(k^{n}t^{2}/2)}{t}}\zeta_{a} is called the upper confidence radius around the empirical estimate μ^a(t)\widehat{\mu}_{a}(t) at the end of tt rounds. We remark here that the difference between our algorithm C-UCB-2 and C-UCB by [LMTY20] is that C-UCB maintains a UCB estimate for each parent value tuple y\mathbf{y}, whereas C-UCB-2 maintains a UCB estimate for each intervention a∈Aa\in\mathcal{A}.

Theorem 4 bounds the expected cumulative regret of C-UCB-2 . In Theorem 4, Δa=max⁡b∈Aμb−μa\Delta_{a}=\max_{b\in\mathcal{A}}\mu_{b}-\mu_{a}.

Observe that in Theorem 4, LaL_{a} is a constant based on problem instance parameters, and hence Theorem 4 proves that C-UCB-2 achieves instance dependent constant regret. Theorem 4 is proved by showing that the expected number of pulls of a sub-optimal arm a∈Aa\in\mathcal{A} after time LaL_{a} is at most 2π23\frac{2\pi^{2}}{3} (proof in Section 8.5). In Section 6, we show via simulations that the expected cumulative regret of C-UCB-2 is better than that of C-UCB, and the experiment also validates that the regret of C-UCB-2 is a constant.

Algorithm Simulations

Experiment 11 (OBS-ALG vs. PB-ALG ): This experiment compares the performance of OBS-ALG with PB-ALG for a fixed budget on a parallel graph with M=50M=50, i.e the reward variable YY has 5050 parents X1,…,X50X_{1},\ldots,X_{50}: Xi∼Bernoulli(pi)X_{i}\sim Bernoulli(p_{i}) for i∈i\in. The rewards variable YY depends on XiX_{i}’s as follows (unknown to both the algorithms): if X1=1X_{1}=1 then Y∼Bernoulli(0.5+ϵ)Y\sim Bernoulli(0.5+\epsilon) and otherwise Y∼Bernoulli(0.5−ϵ′)Y\sim Bernoulli(0.5-\epsilon^{\prime}), where ϵ=0.3\epsilon=0.3, and ϵ′=p1ϵ1−p1∼0.006\epsilon^{\prime}=\frac{p_{1}\epsilon}{1-p_{1}}\sim 0.006. The chosen causal graph structure is the same as in the experiments of [LLR16]. Throughout the experiment pi=0.5p_{i}=0.5 for i∈i\in, and p1=p2p_{1}=p_{2} is the minimum probability pp. The value of p1p_{1} and p2p_{2}, i.e. the minimum probability pp is increased from 0.020.02 to 0.30.3. Figure 3 plots the simple regret of these algorithms with respect to minimum probability. The regret is computed by averaging it over 10001000 independent runs. The budget BB is fixed to a moderate value of 100100 and the cost of intervention γ\gamma to 11. The plot in Figure 3 shows an inverse relationship between simple regret of OBS-ALG and pp as proved in Theorem 1, whereas the simple regret of PB-ALG does not depend on pp. Recall that the expected simple regret of PB-ALG depends on m(p)m(\mathbf{p}) (and not on pp) and for the pip_{i}’s as stated before, the quantity m(p)=2m({\bf p})=2, does not change. Note that the performance of γ\gamma-NB-ALG is best for γ=1\gamma=1 and m(p)=2m({\bf p})=2. Finally, also observe that after a threshold value of pp, OBS-ALG starts performing much better than PB-ALG as can be seen from the plot.

Experiment 22 (γ\gamma-NB-ALG vs. PB-ALG): This experiment compares the performance of γ\gamma-NB-ALG and PB-ALG on a parallel graph with M=50M=50, i.e the reward variable YY has 5050 parents X1,…,X50X_{1},\ldots,X_{50}: Xi∼Bernoulli(pi)X_{i}\sim Bernoulli(p_{i}) for i∈i\in, p1=p2=0.02p_{1}=p_{2}=0.02, and pi=0.5p_{i}=0.5 for i∈i\in. For this choice of pip_{i}’s the PB-ALG algorithm asymptotically achieves its best regret. The rewards variable YY depends on XiX_{i}’s as follows (unknown to both the algorithms): if X1=1X_{1}=1 then Y∼Bernoulli(0.5+ϵ)Y\sim Bernoulli(0.5+\epsilon) and otherwise Y∼Bernoulli(0.5−ϵ′)Y\sim Bernoulli(0.5-\epsilon^{\prime}), where ϵ=0.3\epsilon=0.3, and ϵ′=p1ϵ1−p1∼0.006\epsilon^{\prime}=\frac{p_{1}\epsilon}{1-p_{1}}\sim 0.006. Under these settings, [LLR16] demonstrated a faster exponential decay of simple regret compared to the non-causal algorithms. Since in this experiment we compare γ\gamma-NB-ALG to PB-ALG (adapted to the budgeted version), we choose the same causal graph and distribution. The part a of this experiment in Figure 3 compares the simple regret of the two algorithms when γ=60\gamma=60 and the budget is increased to 30003000. The regret is computed by averaging it over 10001000 independent runs. The part b of this experiment in Figure 5 illustrates the effect on the simple regret of the algorithms as γ\gamma increases from 11 to 7575. In Figure 5, observe that till a threshold value of γ\gamma both algorithms have very close simple regret and post the threshold, γ\gamma-NB-ALG trades off between observations and interventions to yield a much better simple regret.

Experiment 33 (F−KUBE\mathtt{F-KUBE} vs. CRM-NB-ALG): This experiment compares the performance of F−KUBE\mathtt{F-KUBE} and CRM-NB-ALG . The model is as in Experiment 2, except ϵ=0.5\epsilon=0.5, i.e. the best arm has reward 11. If the reward distribution is the same as in experiment 11 the cumulative regret of CRM-NB-ALG even with γ=1.1\gamma=1.1 converges very quickly to a small constant. This is attributed to the fact that the observation arm is closer to being optimal (i.e. d0d_{0} is smaller). Even though this validates the better performance of our algorithm, for a better visual description we set the expected reward of the best arm to 11. Even with this reward distribution, the performance of CRM-NB-ALG is much better than F−KUBE\mathtt{F-KUBE}. Figures 5, 7, and 7 illustrate the cumulative regrets of both the algorithms for γ\gamma equal to 1,1.11,1.1 and 1.51.5 respectively as the budget is increased. The regret is computed by averaging over 5050 independent runs. Notice that CRM-NB-ALG yields a much better regret in all three cases and its regret is constant for γ=1.5\gamma=1.5.

Experiment 44 (C-UCB vs. C-UCB-2): This experiment compares the performance of C-UCB and C-UCB-2 . The causal graph used in this experiment is as shown in figure 9. Notice that this graph has a backdoor path from X2X_{2} to YY and therefore algorithms such as PB-ALG , γ\gamma-NB-ALG and CRM-NB-ALG , which are for no-backdoor graphs, cannot be used. Our conditional probabilities for nodes (P(node∣Pa(node))P(\text{node}|Pa(\text{node}))) are given in Table 1.

The conditional distribution of the reward variable YY was chosen as Y∣w1,w2=θ1X1+θ2X2+ϵY|w_{1},w_{2}=\theta_{1}X_{1}+\theta_{2}X_{2}+\epsilon, where θ1\theta_{1} and θ2\theta_{2} are fixed to 0.250.25 (similar to that in [LMTY20]) and ϵ\epsilon is distributed as N(0,0.01)\mathcal{N}(0,0.01). Here N(0,0.01)\mathcal{N}(0,0.01) denotes the normal distribution with mean 00 and standard deviation 0.010.01. The conditional probabilities in the above table are chosen to be close to each other in order to ensure that the expected rewards for all the arms are competitive and the algorithm takes longer to distinguish between them. The expected reward of the four arms do(Xi=x),i∈,x∈{0,1}do(X_{i}=x),i\in,x\in\{0,1\} are given in the Table 2.

Figure 9 shows a comparison between the cumulative regret incurred by both algorithms for values of TT in the range $.Theregretiscomputedbyaveragingover. The regret is computed by averaging over500$ independent runs. Notice that the regret of C-UCB is much higher than C-UCB-2 and also grows with time. Moreover, the regret of C-UCB-2 grows a little initially and then becomes constant as proved in Theorem 4.

Discussion and Future Work

The Mab problem can be used to model several real-world scenarios where additional information besides the reward of the pulled arms is available and hence the study of the Mab problem with side-information has been an area of significant interest in the research community. One of the most prominent models with side-information is the contextual Mab problem where the algorithm receives extra information (called context) before each arm pull [LPP10]. A class of bandit problems where the side-information obtained conforms to a feedback graph has also been studied in the literature [ACBDK15]. The special case of parallel causal graphs studied in [LLR16] and in this work is in fact captured by such a model, but as shown by [LLR16] their regret bounds are not optimal in this setting.

In this work, we study the the causal bandit problem for no-backdoor graphs in the budgeted bandit framework. In this setting, observations are cheaper compared to interventions, which is practically well-motivated. In Sections 3 and 4 we provided two algorithms, γ\gamma-NB-ALG and CRM-NB-ALG , that minimized the expected simple regret and expected cumulative regret respectively. [SSDS17] also studies the best intervention identification problem via importance sampling under budget constraint. But in contrast to our work, they consider soft interventions on a single node VV, and also assume that the interventional distributions and the marginals of the parent distribution of the node VV are known. This is incomparable with hard interventions on no-backdoor graphs, where interventions can be performed on different variables and the parent distributions of the intervened nodes are not known. Also their setting is parameterized by B′B^{\prime} and TT, where B′B^{\prime} is the upper bound on the average cost of sampling and TT is the total number of samples that the algorithm draws. This can be mapped to our setting by setting the budget to be B′TB^{\prime}T. In our budgeted setting TT is not given as input to the algorithm, and this is important for the trade-off between observations and interventions.

In the non-budgeted setting, we showed that our algorithm C-UCB-2 has constant expected cumulative regret in terms of instance-parameters. We conjecture that the worst-case regret bound of our algorithm matches that in [LMTY20], and resolving that remains open. The work by [SB17] studies a similar problem as that in our work and experimentally show the effectiveness of Thompson Sampling but do not provide any theoretical guarantees.

Finally, many of the works in the literature such as those of [LMTY20] and [LLR16] assume that the parent distribution for each intervention is known to the algorithm. We only make this assumption in Section 5. This assumption is limiting in practice and showing a non-trivial regret guarantee for settings without this assumption remains an important open direction.

Proofs of Theorems

We require the following two versions of the Chernoff-Hoeffeding inequality in our proof.

Suppose X1,…,XTX_{1},\ldots,X_{T} are independent random variables taking values in the interval $,andlet, and letX=\sum_{t\in[T]}X_{t}andand\overline{X}=\frac{\sum_{t\in[T]}X_{t}}{T}.Thenforany. Then for any\varepsilon\geq 0$ the following holds:

2 Proof of Theorem 1

where YtY_{t} is value of YY sampled in round tt. Notice that μ^i,x\widehat{\mu}_{i,x} is the empirical estimate of μi,x\mu_{i,x} computed by OBS-ALG at the end of BB rounds. Similarly the empirical estimate of μ0\mu_{0}, denoted μ^0\widehat{\mu}_{0}, is computed by OBS-ALG at the end of BB rounds as follows:

Finally, also let p^i=p^i,1\widehat{p}_{i}=\widehat{p}_{i,1}. The proof of the theorem is completed using the following lemma.

At the end of BB rounds played by OBS-ALG the following hold:

1) Part 1 directly follows from Lemma 8.1.

2) Observe that E[p^i]=piE[\widehat{p}_{i}]=p_{i}, and hence from Lemma 8.1, for an i∈[M]i\in[M] at the end of BB rounds we have

Since ε≤p/2\varepsilon\leq\sqrt{p/2} (from Equation 4), εBp2≤pB2\varepsilon B\sqrt{\frac{p}{2}}\leq\frac{pB}{2}. This implies

Hence from Equations 5 and 6, for a fixed (i,x)(i,x) the following holds:

3) Notice that p^i,xB\widehat{p}_{i,x}B is the number of times XiX_{i} was sampled as xx in BB rounds. In particular, part 2 of Lemma 8.2 bounds the probability that the number of times XiX_{i} was sampled as xx is small. We use this to prove part 3. First observe that from Lemma 8.1 we have

In particular, Equation 7 bounds the error probability of estimating μ^i,x\widehat{\mu}_{i,x} conditioned on the event that XiX_{i} has been sampled as xx sufficiently many times. Next by law of total probability, for any fixed (i,x)(i,x),

Hence, from Equation 7 and part 2 of Lemma 8.2 we have

Let U0U_{0} be the event that ∣μ^0−μ0∣≤ε|\widehat{\mu}_{0}-\mu_{0}|\leq\varepsilon, and for any i,xi,x let Ui,xU_{i,x} be the event ∣μ^i,x−μi,x∣≤ε|\widehat{\mu}_{i,x}-\mu_{i,x}|\leq\varepsilon. Also let U=(∩i,xUi,x)∩U0U=(\cap_{i,x}U_{i,x})\cap U_{0}, U‾\overline{U} denote the compliment of UU. Then applying union bound on the events in part 1 and 3 in Lemma 8.2, we have that

Let a∗=arg⁡max⁡a∈A(μa)a^{*}=\arg\max_{a\in\mathcal{A}}(\mu_{a}). Note that if event U‾\overline{U} holds then the simple regret of OBS-ALG , rOBS-ALG (B)≤1r_{\text{{OBS-ALG} }}(B)\leq 1. On the other hand, if the event UU holds, and aBa_{B} is the arm output by the algorithm, then rOBS-ALG (B)=μa∗−μaB≤2εr_{\text{{OBS-ALG} }}(B)=\mu_{a^{*}}-\mu_{a_{B}}\leq 2\varepsilon. Setting δ=16Me−ε2pB\delta=16Me^{-\varepsilon^{2}pB}, and substituting the value of ε\varepsilon, we have δ=116Mp2B2\delta=\frac{1}{16Mp^{2}B^{2}}. Hence, the expected simple regret is at most:

3 Proof of Theorem 2

For convenience, we denote m(p)m(\mathbf{p}) and m(p^)m(\widehat{\mathbf{p}}) as mm and m^\widehat{m} respectively. Throughout the proof we assume that BB is such that: a) B≥max⁡(γm,pM)B\geq\max(\gamma m,pM) and b) B≥max⁡(16p2log⁡2MBγm,16p2log⁡2pMB)B\geq\max(\frac{16}{p^{2}}\log\frac{2MB}{\gamma m},\frac{16}{p^{2}}\log 2pMB). Note that the two constraints hold for sufficiently large BB. To begin with observe that if γ=θ(1p⋅m(p))\gamma=\theta(\frac{1}{p\cdot m(\mathbf{p})}) then O(1pBlog⁡(pMB))=O(γmBlog⁡MBγm)O\left(\sqrt{\frac{1}{pB}\log(pMB)}\right)=O\left(\sqrt{\frac{\gamma m}{B}\log\frac{MB}{\gamma m}}\right). Hence, it is sufficient to show that if γ≤15p⋅m(p)\gamma\leq\frac{1}{5p\cdot m(\mathbf{p})} then the expected simple regret of γ\gamma-NB-ALG is O(γmBlog⁡MBγm)O\left(\sqrt{\frac{\gamma m}{B}\log\frac{MB}{\gamma m}}\right) and if γ≥5p⋅m(p)\gamma\geq\frac{5}{p\cdot m(\mathbf{p})} then the expected simple regret of γ\gamma-NB-ALG is O(1pBlog⁡(pMB))O\left(\sqrt{\frac{1}{pB}\log(pMB)}\right). Theorem 2 is proved using Lemmas 8.3 and 8.4.

The following lemma is similar to Lemma 8 in [LLR16].

From Lemmas 8.3 and 8.4 it follows that if F=0F=0 then at the end of B/2B/2 rounds the following holds:

This implies that if F=0F=0 then at then end of B/2B/2 rounds the following holds:

Case a (γ<15p⋅m\gamma<\frac{1}{5p\cdot m}): We condition on F=0F=0. Hence, from the argument above it follows that Equation 9 holds. Hence, γ<15p⋅m≤1p^⋅m^\gamma<\frac{1}{5p\cdot m}\leq\frac{1}{\widehat{p}\cdot\widehat{m}}. This implies at step 6 in γ\gamma-NB-ALG , p^⋅m^<1γ\widehat{p}\cdot\widehat{m}<\frac{1}{\gamma}, and γ\gamma-NB-ALG executes steps 11-14. That is γ\gamma-NB-ALG makes B4γ\frac{B}{4\gamma} interventions in the remaining rounds. The algorithm constructs set A={ai,x∣p^i,x≤1m^}A=\{a_{i,x}\mid\widehat{p}_{i,x}\leq\frac{1}{\widehat{m}}\}. Now for arms in AA, μ^i,x\widehat{\mu}_{i,x} is computed as in step 14 of γ\gamma-NB-ALG, i.e for ai,x∈Aa_{i,x}\in A

Notice that ∣A∣≤m^|A|\leq\widehat{m} (from the definition of m(p^)m(\widehat{\mathbf{p}})). Hence

Thus from Lemma 8.1 for each arm ai,x∈Aa_{i,x}\in A and any ε>0\varepsilon>0

Also for arms not in AA, μ^i,x\widehat{\mu}_{i,x} is computed as in step 3 of γ\gamma-NB-ALG, i.e. for ai,x∉Aa_{i,x}\notin A

The last inequality holds since γ≥1\gamma\geq 1. Using Equations 10 and 11 we have for any arm a∈Aa\in\mathcal{A},

Substituting ε=8γmBlog⁡MBγm\varepsilon=\sqrt{\frac{8\gamma m}{B}\log\frac{MB}{\gamma m}} we have

To get the last inequality, we use that 8M3(γmB)4≤8γmBlog⁡MBγm\frac{8}{M^{3}}\left(\frac{\gamma m}{B}\right)^{4}\leq\sqrt{\frac{8\gamma m}{B}\log\frac{MB}{\gamma m}}, as M≥1M\geq 1 and B≥γmB\geq\gamma m. Finally, we use Equation 12 and Lemma 8.3 to bound the expected simple regret of γ\gamma-NB-ALG in this case as follows:

In last but one line of the above equation, we use that BB satisfies B≥4p2log⁡2MBγmB\geq\frac{4}{p^{2}}\log\frac{2MB}{\gamma m} and B≥γmB\geq\gamma m implying 2Me−p216B2Me^{-\frac{p^{2}}{16}B} is at most 32γmBlog⁡MBγm\sqrt{\frac{32\gamma m}{B}\log\frac{MB}{\gamma m}}.

Case b (γ≥5p⋅m(p)\gamma\geq\frac{5}{p\cdot m(\mathbf{p})}): Again we condition on F=0F=0, and hence Equation 9 holds. Hence, γ≥5p⋅m(p)≥1p^⋅m(p^)\gamma\geq\frac{5}{p\cdot m(\mathbf{p})}\geq\frac{1}{\widehat{p}\cdot m(\widehat{\mathbf{p}})}. This implies at step 6 in γ\gamma-NB-ALG , p^⋅m(p^)≥1γ\widehat{p}\cdot m(\widehat{\mathbf{p}})\geq\frac{1}{\gamma}, and γ\gamma-NB-ALG executes steps 7-9. That is it plays the arm a0a_{0} for BB rounds. Thus, from the analysis of Theorem 1 we have that (see Equation 8)

We use Equation 13 and Lemma 8.3 to bound the expected simple regret of γ\gamma-NB-ALG in this case as follows:

Again in the last but one line of the above equation, we use that 4log⁡MBp2B≤1\frac{4\log MB}{p^{2}B}\leq 1 and hence 2Me−p216B2Me^{-\frac{p^{2}}{16}B} is at most 8pBlog⁡(16pMB)\sqrt{\frac{8}{pB}\log(16pMB)}.

4 Proof of Theorem 3

The proof of Theorem 3 requires the the following lemmas.

1. Since β≥1\beta\geq 1, at the end of TT rounds arm a0a_{0} is pulled by CRM-NB-ALG at least (ln⁡T)2(\ln T)^{2} times. Hence, NT0≥(ln⁡T)N^{0}_{T}\geq(\ln T), and from Lemma 8.1 we have

3. Recall that the effective number of arm pulls of arm ai,xa_{i,x} at the end of TT rounds is

Hence, ETi,x=NTi,x+p^i,xNT0E^{i,x}_{T}=N^{i,x}_{T}+\widehat{p}_{i,x}N^{0}_{T}, where p^i,x\widehat{p}_{i,x} is as defined in part two of this lemma. Hence for any i,xi,x at the end of TT rounds if p^i,x≥p2\widehat{p}_{i,x}\geq\frac{p}{2} then ETi,x≥pNT02E^{i,x}_{T}\geq\frac{pN^{0}_{T}}{2}. Further, as NT0≥ln⁡TN^{0}_{T}\geq\ln T, it follows that at the end of TT rounds if p^i,x≥p2\widehat{p}_{i,x}\geq\frac{p}{2} then ETi,x≥pln⁡T2E^{i,x}_{T}\geq\frac{p\ln T}{2}. Hence, from the definition of μ^i,x(T)\widehat{\mu}_{i,x}(T) and Lemma 8.1, at the end of TT rounds we have for any fixed i,xi,x:

The last line in the above inequality follows from Equation 15 and part 2 of this lemma. ∎

Recall that β\beta is set as in steps 11-14 in CRM-NB-ALG . We begin by making the following easy to see observations.

If a∗≠a0a^{*}\neq a_{0} then d0=μa∗γ−μ0d_{0}=\frac{\mu_{a^{*}}}{\gamma}-\mu_{0}.

Let μ^∗=max⁡i,x(μ^i,x(T))\widehat{\mu}^{*}=\max_{i,x}(\widehat{\mu}_{i,x}(T)) (as computed in step 11 of CRM-NB-ALG ). If ∣μ^0(T)−μ0∣≤d04|\widehat{\mu}_{0}(T)-\mu_{0}|\leq\frac{d_{0}}{4} and ∣μ^i,x(T)γ−μi,xγ∣≤d04|\frac{\widehat{\mu}_{i,x}(T)}{\gamma}-\frac{\mu_{i,x}}{\gamma}|\leq\frac{d_{0}}{4} for all (i,x)(i,x) then d02≤μ^∗γ−μ^0(T)≤3d02\frac{d_{0}}{2}\leq\frac{\widehat{\mu}^{*}}{\gamma}-\widehat{\mu}_{0}(T)\leq\frac{3d_{0}}{2}, and 329d02≤β2≤32d02\frac{32}{9d_{0}^{2}}\leq\beta^{2}\leq\frac{32}{d_{0}^{2}}. Notice that since T≥e50d02T\geq e^{\frac{50}{d_{0}^{2}}}, 32d02≤ln⁡T\frac{32}{d_{0}^{2}}\leq\ln T.

Let U0U_{0} be the event that ∣μ^0−μ0∣≤d04|\widehat{\mu}_{0}-\mu_{0}|\leq\frac{d_{0}}{4}, and for any i,xi,x let Ui,xU_{i,x} be the event ∣μ^i,xγ−μi,xγ∣≤d04|\frac{\widehat{\mu}_{i,x}}{\gamma}-\frac{\mu_{i,x}}{\gamma}|\leq\frac{d_{0}}{4}. Also let U=(∩i,xUi,x)∩U0U=(\cap_{i,x}U_{i,x})\cap U_{0}, and let U‾0\overline{U}_{0}, U‾i,x\overline{U}_{i,x}, and U‾\overline{U} denote the compliment of the events U0,Ui,xU_{0},U_{i,x}, and U‾\overline{U} respectively. From parts 1 and 3 of Lemma 8.5, we have

Since TT satisfies Tp2d0216ln⁡T≥15M\frac{T^{\frac{p^{2}d_{0}^{2}}{16}}}{\ln T}\geq 15M, this implies 32δ9d02≤249d02\frac{32\delta}{9d_{0}^{2}}\leq\frac{24}{9d_{0}^{2}}, and hence E[β2]≥89d02E[\beta^{2}]\geq\frac{8}{9d_{0}^{2}}. Similarly, from part 2 of Observation 8.1 we have that the event UU implies β2≤32d02\beta^{2}\leq\frac{32}{d_{0}^{2}}. Here, we use that if UU does not hold then β2≤ln⁡T\beta^{2}\leq\ln T. Hence

Since TT satisfies Tp2d0216ln⁡T≥15M\frac{T^{\frac{p^{2}d_{0}^{2}}{16}}}{\ln T}\geq 15M, we have δln⁡T≤18d02\delta\ln T\leq\frac{18}{d_{0}^{2}}, and hence E[β2]≤50d02E[\beta^{2}]\leq\frac{50}{d_{0}^{2}}. ∎

Suppose the algorithm pulls the arms for TT rounds and if a∗≠ai,xa^{*}\neq a_{i,x}. Then

For ease of notation we denote E[NTi,x∣T]E[N^{i,x}_{T}|T] as E[NTi,x]E[N^{i,x}_{T}]. Observe that

We require the following observation which is easy to prove.

We continue by taking expectation on both sides of Equation 17 and use Observation 8.2,

If μ^a∗(s,t)γ+2ln⁡tγ2s≤μ^i,x(sj,t)γ+2ln⁡tγ2sj\frac{\widehat{\mu}_{a^{*}}(s,t)}{\gamma}+\sqrt{\frac{2\ln t}{\gamma^{2}s}}\leq\frac{\widehat{\mu}_{i,x}(s_{j},t)}{\gamma}+\sqrt{\frac{2\ln t}{\gamma^{2}s_{j}}} is true then at least one of the following events is true

The probability of the events in Equations 19a and 19b can be bounded using Lemma 8.1,

If a∗=a0a^{*}=a_{0} then using the exact arguments as above we can show that Equation 20 still holds. Hence, using Equations 18 and 20 we have if a∗≠ai,xa^{*}\neq a_{i,x} then

The arguments used to bound E[NT0∣T]E[N^{0}_{T}|T] (denoted E[NT0]E[N^{0}_{T}] for convenience), when a∗≠a0a^{*}\neq a_{0} is similar. In this case the equation corresponding to Equation 18 is

Finally using Equations 21 and 22, we have

If a∗=a0a^{*}=a_{0} and suppose the algorithm pulls the arms for TT rounds then

For convenience, we denote E[NTi,x∣T]E[N^{i,x}_{T}|T] and E[NT0∣T]E[N^{0}_{T}|T] as E[NTi,x]E[N^{i,x}_{T}] and E[NT0]E[N^{0}_{T}] respectively. At the end of TT rounds we have

Taking expectation on both sides of the above equation and rearranging the terms we have,

Before we bound the regret of the algorithm we make the following observation regarding TT, which is the number of rounds CRM-NB-ALG pulls the arms before exhausting the budget BB:

Now are ready to bound the expected cumulative regret of CRM-NB-ALG for the two cases:

Case a (a∗=a0a^{*}=a_{0}): In this case we bound the expected cumulative regret of CRM-NB-ALG for BB satisfying

Observe that the constraint on BB in Equation 24 is satisfied for any large BB. We begin by making the following observation which shows that in this case the expected number of pulls of a sub-optimal arm is bounded by a constant for any large BB. Observe that the constraint on BB in Observation 8.3 is satisfied for any large BB.

Let a∗=a0a^{*}=a_{0}, and TT be the number of rounds CRM-NB-ALG pulls the arms before the budget BB is exhausted, where BB satisfies the constraint in Equation 24. Then ET[NTi,x]≤π23E_{T}[N^{i,x}_{T}]\leq\frac{\pi^{2}}{3}.

From Lemmas 8.7 and 8.8 for any TT satisfying

we have E[NTi,x∣T]≤π23E[N^{i,x}_{T}|T]\leq\frac{\pi^{2}}{3}. Notice that the constraint on TT in Equation 25 is the same as the constraint on Bγ\frac{B}{\gamma} in Equation 24. Moreover, observe that if Bγ\frac{B}{\gamma} satisfies the constraint in Equation 24 then T≥BγT\geq\frac{B}{\gamma} satisfies Equation 25 with probability 11. Hence, ET[NTi,x∣T]≤π23E_{T}[N^{i,x}_{T}|T]\leq\frac{\pi^{2}}{3}. ∎

Next observe that in this case GBG_{B} (see Equation 2) is Bμ0B\mu_{0}, i.e the optimal solution is to play arm a0a_{0} in all the rounds. We require the following observation which lower bounds ET[T]E_{T}[T] in terms of BB, which is the total number of rounds played by the optimal solution.

Let a∗=a0a^{*}=a_{0}, and TT be the number of rounds CRM-NB-ALG pulls the arms before the budget BB is exhausted, where BB satisfies the constraint in Equation 24. Then ET[T]≥B−1−2Mπ2(γ−1)3E_{T}[T]\geq B-1-\frac{2M\pi^{2}(\gamma-1)}{3}.

Let catc_{a_{t}} denote the cost of arm ata_{t} pulled at time t≤Tt\leq T. That is cat=γc_{a_{t}}=\gamma if at=ai,xa_{t}=a_{i,x} and cat=1c_{a_{t}}=1 if at=a0a_{t}=a_{0}. Then the following is always true, as CRM-NB-ALG pulls arms till the budget is the exhausted:

Taking expectation over TT and the sequence of arm pulls {at}\{a_{t}\} made by CRM-NB-ALG , on both sides of the above equation, we have

Finally we bound the expected cumulative regret of CRM-NB-ALG when a∗=a0a^{*}=a_{0} as follows:

Thus, from Observations 8.3 and 8.4, we have

Observe that the expected regret of CRM-NB-ALG is bounded by a constant for large BB and hence O(1)O(1).

Case b (a∗≠a0a^{*}\neq a_{0}): In this case we bound the expected cumulative regret of CRM-NB-ALG for BB satisfying B≥max⁡(L,e50d02)B\geq\max(L,e^{\frac{50}{d_{0}^{2}}}), where LL is as in Lemma 8.6. Observe that the constraint is satisfied for any large BB. Let TT be the number of rounds CRM-NB-ALG pulls the arms before exhausting the budget BB. Then from Equation 25, we have T≥max⁡(L,e50d02)T\geq\max(L,e^{\frac{50}{d_{0}^{2}}}). Hence, from Lemmas 8.6 and 8.7, and as T≤BT\leq B (from Equation 23), we have for a∗≠ai,xa^{*}\neq a_{i,x}

Also observe that in this case GBG_{B} is at most Bμa∗γ\frac{B\mu_{a^{*}}}{\gamma}. Below we bound the expected cumulative regret of CRM-NB-ALG when a∗≠a0a^{*}\neq a_{0}

Hence, we have that the expected cumulative regret of CRM-NB-ALG is:

5 Proof of Theorem 4

Throughout this proof aa and y\mathbf{y} indexes the sets A\mathcal{A} and SnS^{n} respectively. Let δ\delta, L1L_{1}, L2,aL_{2,a} and LaL_{a} for all aa, be as in the theorem statement. Let a∗=arg⁡max⁡a(μa)a^{*}=\arg\max_{a}(\mu_{a}). As is standard in MAB literature, we assume without loss of generality that a∗a^{*} is unique. Further, let Δa=μa∗−μa\Delta_{a}=\mu_{a^{*}}-\mu_{a}. The regret upper bound is proved using Lemmas 8.9 and 8.10.

Let TT be the number of rounds C-UCB 2 has pulled the arms. Then for T≥L1T\geq L_{1} the following holds:

1. Part 1 of the lemma follows from Lemma 8.1.

2. Using Lemma 8.1 again, it follows that for all y{\mathbf{y}} such that cy>0c_{\mathbf{y}}>0, and for all εy≥0\varepsilon_{\mathbf{y}}\geq 0,

Hence, for all y{\mathbf{y}} such that cy>0c_{\mathbf{y}}>0, using the law of total probability we have

implies there is a y\mathbf{y} such that cy>0c_{\mathbf{y}}>0 and {∣μ^y(T)−μy∣≥εy}\{|\widehat{\mu}_{\mathbf{y}}(T)-\mu_{\mathbf{y}}|\geq\varepsilon_{\mathbf{y}}\}. Hence, using part 2 of this lemma and applying union bound over all y\mathbf{y} such that cy>0c_{\mathbf{y}}>0, we have for every aa

where δ=min⁡cy>0cy\delta=\min_{c_{\mathbf{y}}>0}c_{\mathbf{y}}. Since T≥L1T\geq L_{1}, T≥2log⁡(knT2)δ2T\geq\frac{2\log(k^{n}T^{2})}{\delta^{2}}. This implies kne−δ2T/2≤1T2k^{n}e^{-\delta^{2}T/2}\leq\frac{1}{T^{2}}, and

Let a∈Aa\in A be a sub-optimal intervention. Then the expected number of times intervention aa is made after La=max⁡{L1,L2,a}L_{a}=\max\{L_{1},L_{2,a}\} rounds is at most 2π23\frac{2\pi^{2}}{3}.

For ease of notation, we denote log⁡(knt2/2)tζa\sqrt{\frac{\log(k^{n}t^{2}/2)}{t}}\zeta_{a} as ca,tc_{a,t}. Note that ca,tc_{a,t} is the confidence radius of intervention aa C-UCB-2 maintains at the end of tt rounds. Further, let Na,T′N^{\prime}_{a,T} denote the number of times the algorithm performs intervention aa from time La+1L_{a}+1 to time T≥LaT\geq L_{a}, and also let ata_{t} denote the intervention performed at time tt. Hence,

Note that at=aa_{t}=a implies μˉa∗(t−1)≤μˉa(t−1)\bar{\mu}_{a^{*}}(t-1)\leq\bar{\mu}_{a}(t-1) i.e. μ^a∗(t−1)+ca∗,t−1≤μ^a(t−1)+ca,t−1\widehat{\mu}_{a^{*}}(t-1)+c_{a^{*},t-1}\leq\widehat{\mu}_{a}(t-1)+c_{a,t-1} . Hence from Equation 30, we have

The event μ^a∗,t+ca∗,t≤μ^a,t+ca,t\widehat{\mu}_{a^{*},t}+c_{a^{*},t}\leq\widehat{\mu}_{a,t}+c_{a,t} implies that at least one of the following events is true

Since t≥La≥L1t\geq L_{a}\geq L_{1}, using Lemma 8.9 the probability of the events in Equations 31 and 32 can be bounded as:

The event in equation 33 {μa∗<μa+2ca,t}\big\{\mu_{a^{*}}<\mu_{a}+2c_{a,t}\big\} can be written as {μa∗−μa−2log⁡(knt2/2)tζa<0}\Big\{\mu_{a^{*}}-\mu_{a}-2\sqrt{\frac{\log(k^{n}t^{2}/2)}{t}}\zeta_{a}<0\Big\}. Substituting Δa=μa∗−μa\Delta_{a}=\mu_{a^{*}}-\mu_{a} and since t≥La≥L2,at\geq L_{a}\geq L_{2,a}, we have

Now we bound the expected cumulative regret of C-UCB-2. From Equation 3 in Section 2, we have at the end of TT rounds

The inequality in the last line of the above equation follows from Lemma 8.10.

Acknowledgements

Vineet Nair is thankful to be supported by the European Union’s Horizon 2020 research and innovation program under grant agreement No 682203 -ERC-[ Inf-Speed-Tradeoff]. Vishakha Patil gratefully acknowledges the support of a Google PhD Fellowship.

References