Theoretical Paper

Journal of the Operational Research Society (2003) 54, 89–100. doi:10.1057/palgrave.jors.2601484

On the design of lottery games

R Hartley1 and G Lanot1

1Economics Department, Keele University, Keele, UK

Correspondence: R Hartley, Roger Hartley, Economics Department, Keele University, Keele ST5 5BG, UK. E-mail: eca26@keele.ac.uk

Received April 2000; Accepted May 2002.

Top

Abstract

We describe a model of participation in lottery games designed to address the optimisation of tax revenue in state-sponsored lotteries. The model treats participants dynamically and examines a long-run equilibrium. A novel high frequency approximation is used to turn the problem into a static, state-contingent deterministic programming problem. We demonstrate that the solution of this problem has qualitatively plausible properties and then calibrate the model against the United Kingdom National Lottery (UKNL). The results suggest that the current design of the UKNL may not be maximising tax revenue.

Keywords:

lottery, dynamic optimisation, simulations

Top

Introduction

State-sponsored lotteries are extensively used to raise revenue for good causes such as education and the arts. Worldwide the annual value of tickets bought in 1998 exceeded $51bn (see La Fleur and La Fleur1 for references). States often contract out the running of such games to an operating company with a proportion of the cost of each ticket going to good causes or to general tax revenue and the contractor taking a further proportion as operating costs; most of the remainder is returned to participants as prizes. For example, in 1998 out of every pound spent on the UK National Lottery (UKNL), about 50 pence went to prizes, 28 pence to good causes, 13 pence to Customs and Excise as excise duties, 5 pence to the retailer as commission, 3 pence to cover the operating company (Camelot) operating costs, and Camelot kept 1 pence as profits. The tax rate and the level of participation will affect the total good-causes revenue. Participation depends on the tax rate and the design of the lottery, in particular, the probability of winning a jackpot prize. In order to study the effects of changing such parameters on the good-causes revenue it is necessary to model their effect on participation. In this paper we present such a model and use it to address the issue of whether the current tax rate and jackpot probability are the appropriate ones to maximize good-causes revenue in the UKNL.

A common feature of lotteries is that, if there are no winners in a given draw, the jackpot prize pool from that draw is added to the pool for the next draw: a rollover. The increased pool typically induces a higher level of participation. This suggests good-causes revenue may be increased by making rollovers more likely and this can be achieved by reducing the probability of winning a jackpot. However, doing so may also have a disincentive effect on participation in draws that do not feature a rollover and thus offset any increase in revenue. To capture such a trade-off effectively requires a dynamic model of individual participation. One possible approach would be to fit a 'black box' model (eg, regression-based) using standard statistical methods. However, the similarity and stability of lottery designs means that the available data offers little natural variation and poor prospects of a good fit, particularly if one controls for differences in country-specific propensities to gamble. Our approach is to work from first principles and build a model of individual participation based on standard inter-temporal models of individual behaviour under uncertainty. Lack of variation in the data forces us to adopt the simplest model parameterization. Nonetheless, this is sufficient to capture the principal features of observed behaviour. Since quantitative verification of the model is not possible, we check that the qualitative predictions of the model are plausible and consistent with casual observation. Our final step is to aggregate the behaviour of individual agents into a full model of demand for lottery tickets, which enables us to examine the effects of changing lottery design and tax rates.

There is little literature on the issues discussed here and (to our knowledge) none which bases the optimal design problem on a dynamic model of individual behaviour. Farrell et al2 equate the effective price of a ticket in the UKNL to the ticket price less the expected value of the prizes and use the variation in the latter in rollover weeks to estimate the elasticity of demand. Farrell et al3 analyse the same data using use time-series methods. Although both these papers find values for the short- and long-run elasticities roughly consistent with revenue maximisation, they implicitly assume risk neutrality and neither addresses inter-temporal substitution by participants. Furthermore, they do not disentangle the effects of the tax rate and the lottery design and, using point estimates, they cannot investigate non-local changes of the parameters.

Scoggins4 offers rank-dependent utility as an explanation of the complex prize structure. An interesting paper by Simon5 attempts to construct a model that fits the first couple of years of participation in the UKNL. Although it is too tuned to UK data to be useful as a general tool for addressing design problems, perhaps its most interesting contribution is an attempt to model how observed attrition in participation could be counteracted by the high probability of winning small fixed prizes. This is usually offered as justification for the presence of small fixed prizes. We have not explicitly included this (small) effect in our model as it will not show up in the high-frequency limit we employ.

We conclude the Introduction by giving an informal description of our model. Full details of each component of the model are then set out in subsequent sections.

Lottery draws take place at regular intervals and, at each draw, consumers decide how many tickets to purchase. The cost is measured in foregone expenditure on other goods. Each ticket displays a set of positive integers (1 to 49 for the UKNL) from which participants must choose a 'combination' ie, a subset of the available integers (of size 6 for the UKNL). After all purchases of tickets have been made, the Lottery organisers randomly draw a combination according to a uniform distribution (see Haigh6,7 and references therein). A ticket holder whose chosen combination matches that drawn is a jackpot winner and his/her prize is determined by dividing the jackpot prize pool equally amongst all the jackpot winners. The jackpot prize pool is a fixed proportion of the prize pool augmented by any amount rolled over from previous draws. Such rollovers occur if there are no jackpot winners in a given draw, in which case the jackpot prize pool is held over and added to the jackpot prize poo1 in the next draw. A run of draws with no winners can lead to such a large jackpot prize pool that behaviour regarded as undesirable (such as attempting to buy all combinations) can occur. To avoid this, a limit on the number of consecutive draws that can be rolled over is imposed. If this limit is reached, we assume, for expositional convenience, that the money in the accumulated jackpot prize pool is simply lost from the system. The remainder of the prize pool is allocated by looking at partial matches of the subsets chosen with that drawn. In our model, we concentrate on jackpot prizes where winning is rare (1 in 14 million chance for the UKNL) with values large enough to change the life of winners and other prizes that are low in value and common (more than 1 in 60 chance for fixed-value prizes in the UKNL).

Knowing the prize structure, we model participants' decisions of how many tickets to purchase as a dynamic optimisation problem. Since the number of tickets bought is very large (over 60 million per draw in the UKNL), we assume that consumers treat rollovers as exogenous, ignoring the minuscule effect their individual actions have on their occurrence. Nevertheless, aggregate levels of participation determine (stochastically) the rollover process. This results in a circular model in which the rollover process affects participation and the latter affects the rollover process. Rather than specify the dynamics of the adjustment process, we analyse the equilibrium in which we have rollover-contingent participation levels that are optimal given the rollover process determined by those participation levels. Once the full network of terminals was installed, participation levels rapidly stabilised in the UKNL, suggesting that equilibrium had been reached.3

In the next section, we develop our model of consumer choice as a stochastic dynamic optimisation problem. In the third section and the Appendix, we analyse the limit as the time interval between draws approaches zero. By regarding this limit as a first order approximation to the dynamic problem, we are able to recast the consumer's optimisation problem as a static, deterministic, non-linear programming problem. In the fourth section we study the Kuhn–Tucker conditions for this problem and show that the unique optimal solution captures several qualitative properties of the behaviour we are seeking to model. The next section closes the model by relating rollovers to participation levels. This necessitates the construction of a simple sub-model of how combinations are chosen. Finally, we apply the model to the UKNL and draw some tentative conclusions on the current values of the parameters.

Top

Individual behaviour

We will suppose participants live indefinitely (or participants leaving the system are replaced with others having similar tastes) and divide their expenditure in period t between yt on lottery tickets and xt on other goods so as to maximise the expected value of the lifetime utility:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where delta=e- betaDeltat is the discount factor, beta is the instantaneous rate of discount, Deltat is the time interval between draws and u is the instantaneous, von Neumann–Morgenstern rate-of-utility function in period t.

We have included lottery expenditure as a continuous variable, which will permit the use of calculus-based methods in the sequel. As tickets are indivisible, this raises the question of the interpretation of non-integral y. If restricted to purchasing whole numbers of tickets, it is not hard to see that, even in the absence of rollovers, a participant's optimal number of tickets will typically vary from week to week. However, if purchases are made more frequently (as we shall assume), the expected utility will be nearly the same as that of a stationary policy in which the (non-integral) long-run average number of tickets is purchased in each draw.

Friedman and Savage8 explain how gambling can be reconciled with risk averse behaviour in a static expected utility framework by including convex sections in the utility function. Dynamic versions of this approach based on (1) are easy to write down but Farrell and Hartley9 show that repeated purchase of lottery tickets cannot be explained using expected utility functions unless consumers derive direct value ('fun') from their purchases of tickets (or capital markets are imperfect). We shall also assume an absence of income effects in lottery purchases since this allows us to simplify the utility function and is broadly consistent with observation. The most general form of utility function consistent with these observations is

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where v captures the consumer's attitude to risk and will be assumed increasing and concave, whilst the 'fun function' f expresses the direct consumption value of lottery purchases in terms of expenditure on other goods. We also permit the consumption value to depend on external factors, ae, such as a rollover draw, as well as past history, ah, including the consumer's past winnings, if any. We will specify these variables more fully in the sequel. Under such a specification, participation is induced both by the direct enjoyment gained from purchasing tickets (which is independent of whether a prize is won) as well as the extra wealth in the event of winning.

Rollover technology

We allow expenditure choices to be contingent on the number of weeks that the jackpot has been rolled over, to reflect the observed increase in ticket sales in rollover weeks.2 We choose this way of allowing for rollovers to simplify the model. In fact making purchases contingent on the actual amount of money rolled over makes no difference in the high frequency limit used in the sequel. Furthermore, measured levels of participation contingent on weeks rolled over is subject to rather little variation.2

Let the (non-negative) integer-valued random variable Rt denote the number of weeks the jackpot has been rolled over into week t. More precisely, Rt=k if and if only the jackpot was won in week t- k- 1 but not in weeks t- k,...,t- 1 (or, if the jackpot has never been won and t=k+1). It follows that the transition probability P(Rt+1=j|Rt=k) is zero if jnot in{0,k+1}. Thus, if K is the maximum permitted value of Rt, the latter follows a Markov chain with Transition Probability Matrix (TPM):

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where pik=P(Rt+1=k+1|Rt=k)>0 for k=0,..., K- 1, and the rows and columns are entered in the order 0,1,2,...,K- 1,K. Since Q is uni-chained and acyclic, there is a unique limiting probability vector p=(p0,p1,...,pK)greater than or equal to0 which satisfies QTp=p and sumj=0Kpj=1. It is readily verified that p0=sigma- 1 and pk=sigma- 1pi0...pik- 1 for kgreater than or equal to1, where sigma=1+sumj=0K+1pi0...pij.

Our assumptions on how ae and ah affect f are designed to reflect certain observations. Participation is typically higher in rollover draws than in normal draws. Higher potential winnings are an obvious explanation but it is also possible that the increased publicity and general excitement of participation associated with rollover draws increase the direct consumption value. To maintain the simplicity of the model, we ignore other environmental factors so that ae depends only on Rt and write fk(y) for the value of expenditure y on tickets when Rt=k. We also assume that the fun associated with purchasing tickets derives predominantly from contemplation of winning a prize large enough to change the winner's lifestyle substantially. We reflect this by including only a past jackpot win in the history ah and, in particular, setting fequivalent to0 for a consumer who has won the jackpot in the past. This will mean that such a consumer no longer participates.

Since participation and other consumption will typically depend on the rollover state and past winnings, if any, they are random variables. We write Yt[Xt] for expenditure on lottery tickets [other goods] in period t where the process {Xt,Yt}t=0infinity is adapted to the rollover process {Rt}t=0infinity and the individual winning process.

Thus, consumers maximise:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where the stopping time T denotes the period in which the consumer wins the jackpot. Note that T is a random variable that depends (probabilistically) on {Yt}t=0infinity, and takes the value +infinity for a non-participant.

Constraint

A consumer's expenditure is subject to the lifetime budget constraint:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where r=exp(betaDeltat)- 1 is the rate of interest, m is income per period (assumed constant) and Wt* (y) is a non-negative random variable representing winnings in period t (taking the value {0} if the consumer does not win in period t) if y is spent on lottery tickets. We interpret the constraint as holding with probability one. This prevents borrowing on the strength of future lottery wins since for any alt epsilon>0 and pattern of participation, there is a positive probability that the net present value at t=1 of future winnings is less than alt epsilon.

For t>T, lottery winnings cease so that (3a) can be written:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where Gt(y) is the random variable representing the value of a non-jackpot prize in period t when y is spent on lottery tickets and takes the value 0 if no such prize is won. Similarly, ignoring the fact that a ticket cannot win more than one type of prize, Wk is a random variable representing the value of a jackpot prize in state k.

To solve the consumer's problem, we first look at the situation faced by a jackpot winner. This involves examining the optimisation problem for a consumer conditional on T=t^ and Wk(t^)=w. Discounted to the first period, this yields optimal utility of Deltatdeltat^- 1psi (w;Deltat) where:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

We have assumed that winners continue to receive per-period income m (or its equivalent from quitting a job in monetary value). Since v (dot) is concave, we may use the fact that {r/(1+r)t}t=1infinity is a probability distribution, and that the constraint says that the mean value of {xt}t=1infinity is x¯=m+wr/Deltat(1+r), to apply Jensen's inequality and deduce that the optimal solution is xt=x¯ for every t. Hence

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

We shall focus on the case when Delta t is small and note that, as Deltatright arrow0+, Delta tpsi(w;Deltat)right arrow(1/beta)v(m+betaw).

The consumer's problem can therefore be written as:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Top

An approximation to the consumer's problem

Since lottery draws occur at least weekly, the discount factor delta will be close to one. In a Markov reward process with this property, a good approximation to the expected discounted reward per period is the long-run average reward obtained by taking the expectation of the reward per period with respect to the long-run stationary probabilities. (See, for example, Whittle,10 p. 123). However, the objective function and constraint in CP are complicated by the presence of additional terms representing the possibility of exiting the reward process upon winning the jackpot. Since the probability of winning the jackpot is so low, for a consumer who buys the same small number of tickets each week, wins can be regarded as events in a (discrete approximation to a) Poisson process with expected inter-event times of the order of 107 periods. When considering such rare events, it is quite inappropriate to approximate the expected discounted value of the jackpot by the expected average reward as the discount rate beta, though small, is much larger than the quit rate rho. To cope with these two distinct components of the reward, we adopt a hybrid approach in which consumers receive the long-run average reward per period until they win the jackpot, the value of which is discounted to the current period.

The consumer's problem CP is actually slightly more complicated than the preceding description suggests, for the probability of winning the jackpot will depend on the number of tickets purchased and this is contingent on the rollover state. Hence, quitting probabilities are state-dependent and the corresponding process is not Poisson. Nevertheless, the expected number of periods until a jackpot win will still be very large and, if we assume that the pre-quit process can be regarded as in the steady state, it is plausible that exiting can still be well-approximated by a Poisson process. In this process, we take the probability of a jackpot win leading to a quit as the expected probability of winning, taken with respect to the steady state probabilities. A formal statement of these results, with proofs, is offered in the Appendix.

Objective function

For consumers who have not yet won the jackpot we concentrate on stationary decision rules in which consumers' decisions depend on the rollover state k but not on t and write yk[xk] for the rate of expenditure on lottery tickets [other goods] when Rt=k. To allow for winning the jackpot, the environmenzt of any particular consumer can be summarised by the extended state space:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where, at period t, the extended state is (k,0) if the consumer has not yet won the jackpot and Rt=k is (k,1) if the consumer won the jackpot in period t- 1 and Rt- 1=k and is omega macr if the consumer won the jackpot in a period before t- 1. Writing D for the order-(K+1) diagonal matrix with dkk=yk and rho for the probability that a single ticket wins the jackpot, the transition matrix for the extended state space becomes:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where every component of e is 1 and the 0 s represent appropriately dimensioned zero matrices or vectors. (Note that since Q is a transition matrix, the elements in each row sum to one). We have used the fact that the probability of winning the jackpot when Rt=k is rhoDeltat yk to a very good approximation. Define the (K+1)-vectors:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Then the expected value of the objective function of the consumer's problem CP is the (0,0)th (in the extended state space) component of:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where A=sumt=1infinityDeltatdeltat- 1Pt- 1.

In the appendix, we demonstrate that this component is well approximated by:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

This result demonstrates our earlier discussion. The first term is the net present value of receiving the long run expected reward every period for a long time. The second term is the expected utility of winning the jackpot when the rollover variable is in the steady state, discounted to the initial time.

Constraint

We shall treat the budget constraint in (CP) in a similar way. We assume that, conditional on Rt=k, Gt(yk)>0 with probability Deltatphi(yk) where Deltatphi(yk) is the probability of winning a non-jackpot prize. Similarly, conditional on Gt(yk)>0 and the rollover state being k, we assume that Gt(yk) is distributed as Gk: the non-jackpot prize.

We write gk for E[Gk]. In the Appendix, we demonstrate that the left-hand side of the constraint in CP is well approximated by:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

which says that the budget constraint must be satisfied on average using the limiting probabilities.

Since v is increasing in xk, the consumer's optimisation problem becomes:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Top

Properties of the optimal solution

An important test of our model of individual participation is that the optimal solution of CP* exhibits plausible qualitative properties if reasonable assumptions are made about the preferences. We now turn to such a set of assumptions. These are stated in terms of derivatives so we suppose throughout that v,f0,f1,...,fK and phi are twice differentiable for positive arguments.

Assumptions:

  1. For all x>0, we have v' (x)>0 and v"(x)<0.
  2. For k=0,...,K, we have fk(0)=0 and 0<f'k(y)<1, f"k(y)<0 for all y>0.
  3. For Kgreater than or equal tok>jgreater than or equal to0 and all y>0, we have f'k(y)greater than or equal tof'j(y).
  4. We have phi<(0)=0, and phi'(y)>0, phi"(y)less than or equal to0 for all y>0.
  5. For Kgreater than or equal tok>jgreater than or equal to0, we have gkgreater than or equal togj, and Epsi(Wk)greater than or equal toEpsiWj).

Assumption 1 says that 'other goods' are desirable and participants are risk averse. The sign constraints on the derivative in Assumption 2 state that lottery tickets are desirable but that the additional fun derived from the purchase of an extra ticket falls as the number of tickets purchased rises. The requirement f'k(y)<1 implies that the 'fun' generated by purchasing tickets does not exceed the foregone consumption of other goods. Hence, both fun and the prospect of winning prizes are required to induce lottery participation; neither alone is sufficient. Assumption 3 specifies that marginal fun is non-decreasing in the size of the rollover, and Assumption 4 says that phi is concave. This will hold if participants choose the numbers on their tickets optimally. When the number of tickets purchased is small, we expect phi to be (nearly) linear. Assumption 5 requires the expected size of prizes to be non-decreasing in the rollover state variable. The total prize pool always increases but the number of participants and therefore, on average, the number of winners sharing this pool will also increase. Our assumption, which is supported by empirical evidence, is that this effect is not large enough to fully offset the increase in the prize pool.

Assumptions 1 and 2 mean that v(x+fk(y)) is a strictly concave function of (x,y) for k=0,...,K. Hence, this is also true for the objective function in CP*. Assumption 4 means that the left-hand side of the equality constraint is a convex function of (x,y). Since the objective function is increasing in xk the constraint may be replaced by a weak inequality. Furthermore, this inequality can be satisfied strictly (by setting x and y small and assuming m>0). It follows that the first-order Kuhn–Tucker conditions are necessary and sufficient and have a unique solution.11 For k=0,...,K, these conditions can be written in the form:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where lambda is a Kuhn–Tucker multiplier for the first constraint. We have cancelled pk from both sides of these inequalities, using the fact that pk>0 for all k.

We shall consider only solutions that are interior for other goods. Thus, all xk>0 and, since v' is strictly decreasing by Assumption 1, we derive:

Optimality condition 1

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

It follows from Assumption I that lambda>0 and, hence, that the second part of the optimality condition can be written:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Assumptions 2, 3 and 5 imply that:

  1. xik(y) is strictly decreasing in y for a given k,
  2. xik(y) is strictly increasing in k for a given y.

From (a) and (8), we deduce that, if xik(0)>1, then yk>0 and xik(yk)=1; otherwise, yk=0. When xik(0)>1 for some k=0,...,K, we can define J to be the smallest k such that xik(0)>1. Then, (a) and (b) lead to:

Optimality condition 2

If xiK(0)less than or equal to1, then y0=dotdotdot=yK=0, (Case NP) otherwise, yk=0 for k<J and yk>0 with xik(yk)=1 for kgreater than or equal toJ. (Case P)

Case P corresponds to lottery participation (when the rollover state variable is at least J) and NP to non-participation for any rollover state.

Properties (a) and (b) further imply that yk>yj when k>jgreater than or equal toJ. Now observe that, by integrating the inequality in Assumption 3 from 0 to y, and using Assumption 2, we obtain fk(y)greater than or equal tofj(y) for all ygreater than or equal to0 when k>j. Hence,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

using Assumption 2. In view of Optimality condition 1, we can conclude that xk<xj.

If k>j, then it follows from Assumption 2 that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Suppose, further, that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

ie, additional participation following an increase in the rollover state is induced solely by the enlarged prize pool rather than by any additional fun. Then, using Optimality condition 1,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Theorem 1 The optimal solution of CP* has one of the following forms:

Case NP: For all k=0,..., K, we have yk=0 and xk=m.

Case P: There exits a Jset symbol{0,..., K} such that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

For a potential participant the theorem says that positive participation will take place if the number of rollovers is at least J, as the rollover state increases, more is spent on lottery tickets, less on other goods. Nevertheless total expenditure increases if the fun function is independent of the rollover state.

Top

Equilibrium/consistency

In this section we recognise that the collective decisions of all consumers will determine aggregate lottery expenditure and this influences both the transition matrix Q and the (random) jackpot prize W=(W0,...,WK). Current lack of suitable data on heterogeneous purchases forces us to consider a homogenous population of identical individuals who engage in the gamble whatever the rollover state. Nonetheless, the arguments and formulae we develop below are readily extended to a population of heterogeneous agents.

We assume that there are n individuals sharing the same preferences, reflected in the utility function v and the fun functions fk, and the same income m. In rollover state k, the number of tickets sold is nk=nyk and this affects both Q and W through the probability distribution of the number of prizewinners; if this number is 0 the jackpot prize pool is rolled over. To obtain this distribution, the simplest assumption to make would be that each combination of numbers (6 out of 49 for the UKNL) is chosen independently and with equal probability (ie, =rho). In this case, the number of winners would follow the binomial distribution B(rho,nk) or, to a very good approximation, the Poisson distribution P(rhonk). However, the equal probability assumption conflicts with observation in that it significantly under-predicts the probability that there are no winners2,7 and a consequent rollover.

We will therefore assume that any ticket design induces a probability distribution chi over the set I of permitted combinations, such that the probability of the combination iset symbolI being chosen on a randomly selected ticket is chi(i), independently of all other choices. Since all combinations are drawn with probability |I|- 1, the distribution Zk(chi) of the number of winners when the rollover state variable is k is an equally weighted combination of Poisson distributions. In particular,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

The uniform case corresponds to chi(i)=rho=|I|- 1 for all iset symbolI; a simple, two-parameter generalisation is discussed in the next section. The probability of a rollover is given by

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

and the expression for the stationary probabilities implies that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

We now turn to the distribution of Wk, the jackpot prize from the perspective of a participant who will use chi(dot) to choose their combination. The total jackpot prize pool in rollover state k is a proportion alpha of the expenditure on tickets net of taxes and operating costs augmented, if kgreater than or equal to1 by the jackpots rolled over from the previous k- 1 draws. Hence the total jackpot prize in this state is alpha(1- tau- kappa)n¯k, where n¯k=suml=0k- 1nl, and kappa is the proportion of revenue attributed to operating cost and profit. For any ticket, conditional on the combination iset symbolI being drawn, the number of winning tickets excluding that ticket is a random variable Zik where Zikapprox P{chi(i)(nk- 1)}. It follows that a winning ticket will receive a jackpot prize of alpha(1- tau- kappa)n¯k/(Zik+1). Since the probability that i is chosen, given that a ticket wins, is chi(i), we deduce that the probability distribution of the size of jackpot prize is:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

(We have ignored the highly unlikely possibility that any participant wins two or more jackpot prizes in a given draw.) Using (12), we can write, by rearranging the series, for k=0,...,K:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Since non-jackpot prizes are not subject to rollovers, a proportion (1- alpha)(1- tau- kappa) of the total expenditure on tickets is allocated to non-jackpot prizes and the average prize is worth gk in state k. As all tickets are equally likely to win, the probability of winning must be (1- alpha)(1- tau- kappa)/gk. We shall further assume that phi is linear (since y is small):

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

in which case the term phi (yk)gk in the constraint of CP* can be replaced with (1- alpha)(1- tau- kappa)yk. We can therefore write the (necessary and sufficient) first-order conditions for households in state k as:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Note that our assumption of homogeneity forces yk>0 for all k so that the inequality in (15) is satisfied as an equation. The equilibrium solutions (x,y) are the solutions of equations (10), (11) and (13) to (16).

Top

Simulations

In this section, we apply the model to the UKNL and so we need to select specific functional forms for the utility function, fun functions, and the distribution chi. The limited variability in the data dictates the use of parsimoniously parameterized functions and, for the utility and fun functions, we choose functions with a single parameter. Individual preferences are assumed to exhibit constant relative risk aversion (CRRA):

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

The parameter theta characterizes relative risk aversion (defined as - xv"(x)/v' (x)) and is assumed constant. The literature tends to favor CRRA as a parsimonious and accurate description of individual behaviour.12,13

A convenient model for the fun function is a power function.14 However, in order to satisfy Assumptions 2 and 3, we need to transform it into the form

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where 0<etajless than or equal toetak<1 for Kgreater than or equal tok>jgreater than or equal to0. For simplicity, in the simulations we assume that etak=eta, forallk. These are the conditions in the last part of Theorem 1, which allow us to deduce that expenditure increases with the rollover.

The distribution chi needs to be parsimoniously parameterised and not tied too closely to a particular ticket design (since changing rho entails changing the design). A simple way to do this is to assume that a proportion 1- zeta, 0<zetaless than or equal to1, of the |I| possible combinations are never chosen and that the remaining combinations are equally likely to be chosen, with probability rho/zeta. A slightly more sophisticated version assumes that some participants use this distribution and others choose uniformly. Indeed, in many lotto games, computerized pseudo-random generation of combinations is offered to participants (called 'Lucky Dip' in the UKNL). Let gammaset symbol(0,1) be the proportion of tickets on which combinations are chosen in this way. Then chi(i)=gammarho(=xi1) if i is a forbidden combination, and chi(i)=(1- gamma+zetagamma)(rho/zeta)(=xi2), otherwise. It follows from (10) that the probability of a rollover in state k is:

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

and (13) can be written

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where, for i=1,2,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

We take alpha and kappa, the jackpot proportion and operating costs, from the current design and consider values of tau, the tax rate, and rho in a range of 4% either side of the current values, while we put income equal to the average UK weekly household expenditure and N to the number of UK households. We assume an annual interest rate of 10% and set the maximum number of rollovers K to 2 reflecting the maximum observed value. This leaves setting the values of the parameters: zeta, the proportion of combinations available for selection, theta, the coefficient of relative risk aversion, and eta, the fun function parameter, to complete the calibration. The values are chosen such that (i) sales in normal (non-rollover) weeks (ii) sales in single rollover weeks, and (iii) the long run proportion of normal weeks agree with observed values. We note that the fitted value of theta falls within the commonly observed range of values (ie, between 0.5 and 3).

In Figure 1 we present contour plots in the (tau+kappa,rho) plane for the parameter values in Table 1. The panels illustrate how local changes in the design affect the expected tax revenue (top left panel), the expenditure per household on the lottery during a non-rollover week (top right panel) and during a single rollover week (middle right panel), the expected number of jackpot winners during a non-rollover week (middle left panel), the expected size of a jackpot in a non-rollover week (bottom right panel), and the long run proportion of non-rollover weeks (bottom left panel).

Figure 1.
Figure 1 - Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Simulation of design changes, contour plots.

Full figure and legend (61K)


With the parameters set at the calibrated values, we see that increasing rho, and more surprisingly, decreasing tau, leads to an increase in expected tax revenue. In the former case, although the number of rollovers and therefore of potential additional revenue associated with more rollovers falls, this is more than offset by the increased participation induced by the increased chance of a jackpot win; the number of winners also increases, indeed sufficiently so that expected (normal week) jackpots fall in spite of the enlarged jackpot prize pool. Decreasing tau also increases participation sufficiently to offset the reduced tax take per ticket. Normal and rollover week expenditure rise (in contrast to increasing rho), the number of rollovers falls and the expected number of winners and the expected jackpot both increase. Changing the parameters from the current values in a direction of steepest ascent of the expected tax revenue function has the same effect on all the graphed functions as increasing rho. (The expected jackpot decreases slightly).

Before drawing policy conclusions from these results, it is appropriate to examine the sensitivity of the results to the parameter values and functional forms. Indeed, it is precisely because they will be particularly sensitive to such model details that we have refrained from suggesting optimal values. To investigate these issues we have experimented by varying the parameter values (re-calibrating where appropriate) and with different functional forms such as a constant absolute risk aversion utility function. The results of these experiments are too numerous to include in the paper but most of the qualitative features of the functions plotted in Figure 1 remain robust to these changes (details are available from the authors). In particular, the expected tax revenue function always appears to be upward sloping at current parameter values in a west to northwest direction. We suggest that this evidence is sufficient to prompt further investigation into whether the tax rate in the UKNL may be too high and whether the probability of winning the jackpot may be too low.

Top

References

  1. La Fleur T, La Fleur B (1999). La Fleur's 1999 World Lottery Almanac. TLF Publications Incorporated: Maryland, USA.
  2. Farrell L, Lanot G, Hartley R, Walker I (1999). The demand for lotto: the role of conscious selection. J Bus Econ Statist 18: 228–241.
  3. Farrell L, Morgenroth E, Walker I (1999). A time series analysis of UK lottery sales: the long run price elasticity. Oxf Bull Econ Statist 61: 513–526. | Article |
  4. Scoggins JF (1995). The lotto and expected net revenue. Nat Tax J 5: 61–70.
  5. Simon J (1998). Dreams and disillusionment: A dynamic model of lottery demand. Four essays and a note on the demand for lottety tickets and how lotto players choose their numbers, unpublished Ph.D., Department of Economics, European University Institute (Florence, Italy).
  6. Haigh J (1996). Lottery—the first 57 draws. Roy Statist Soc News 23: 1–2.
  7. Haigh J (1997). The statistics of the national lottery. J Roy Statist Soc (Series A) 160: 187–206. | Article |
  8. Friedman M, Savage LJ (1948). Utility analysis of choices involving risk. J Polit Econ 56: 279–304. | Article |
  9. Farrell L, Hartley R (2002). Can Friedman-Savage utility functions explain gambling? American Economic Review 92: 613–624. | Article |
  10. Whittle P (1982). Optimization Over Time. John Wiley and Sons: Chichester.
  11. Rockafellar RT (1970). Convex Analysis. Princeton University Press: Princeton.
  12. Deaton A, Muellbauer J (1980). Economics and Consumer Behavior. Cambridge University Press: Cambridge.
  13. Hirshleifer J, Riley JG (1992). The Analytics of Uncertainty and Information. Cambridge University Press: Cambridge.
  14. Johnson JEV, Shin HS (1995). A violation of dominance and the consumption value of gambling. Department of Economics. University of Southampton, WP 9525.
  15. Cox DR, Miller HD (1965). The Theory of Stochastic Processes. Chapman and Hall: London.
Top

Appendices

Appendix

We will use the following ergodic results for which we offer an elementary proof.

Lemma Suppose {at}t=1infinity is a sequence of independent and identically distributed random (K+1)-vectors and Ea1 exists. Let Q* be the transition matrix in which every row is p (the stationary distribution of Q). Then, as Deltatright arrow0+

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

with probability one.

Proof. It is well known15 that Qtright arrowQ* as tright arrowinfinity. Hence, QQ*=Q*Q=Q*Q*=Q*. It follows inductively that, if S=Q- Q*, then Qt=St+Q* for tgreater than or equal to1. Consequently, Stright arrow0 as tright arrowinfinity, which means that I- S is non-singular. Thus, using delta=e- betaDeltat,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Write

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

The latter inverse exists since delta(I- rhoDeltat D)Q is a sub-stochastic matrix. Note that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

from which it follows that

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

For ngreater than or equal to1, write sn=(1/n)sumt=1nat. By the strong law of large numbers, for almost all sample paths omega, sn(omega)right arrowEa1. For such an omega, we shall write

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

and prove that S(Deltat)right arrowAinfinityEa1 as Deltatright arrow0+. Note that by rearranging the sums,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

A similar argument gives

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

Taking the difference of these results gives

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

and the proof is completed by showing that the term in braces approaches 0 as Deltatright arrow0+.

To achieve this, we use the vector and matrix norms parallelxparallel=maxk|xk|, parallelAparallel=maxk'sumk|ak'k|, and the fact that parallelAxparallelless than or equal toparallelAparallelparallelxparallel together with the standard results on matrix norms. We have parallelQparallel=1 and parallelI- rhoDeltat Dparallel=1- rhoDeltaty, where y=minkyk. Since

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

there exists a mu>0 and delta1>0 such that, if 0<Deltat<delta1, then 1- delta(1- rhoDeltat y)>muDeltat.

Let alt epsilon>0. Since st(omega)right arrowEa1 as tright arrowinfinity there is an N such that parallelst(omega)- Ea1parallel<mu2(alt epsilon/2) for all t>N. Define

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

If Deltat<min{delta1,delta2}, then

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

which concludes the proof.

For igreater than or equal to1,

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where the missing term is chosen to make the rows sum to 1. It follows that the (0,0)th component (in the extended state space) of (5) is the 0th (in the original state space) component of

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

using the lemma. (We can drop the 'almost surely' qualification in this case.) For the UKNL rho- 1 is close to 14times106 and if interest rates are less than 10% then beta does not exceed 1/500. Hence for purchases of a few tickets a week, the inverse matrix in the expression above is close to I. Furthermore, since v is concave and increasing, it is bounded above by an affine function and below by a constant, and this allows us to apply the dominated convergence theorem when taking the limit of DeltatEbulletpsi(Wk;Deltat) as Deltatright arrow0+. Expression (6) is a rewriting of the limit in (A2) and in using it in our model, we are assuming that draws are sufficiently frequent for us to apply this limit. One observation that supports this assumption is that winning small prizes has little effect on expenditure.

We can use the same approach for the constraint. In time period t, if the rollover state is k, the probability of winning a non-jackpot prize is Delta t phi(yk) and if won, the prize is gtk, an independent copy of g1k. Thus the left-hand side of the constraint of CP* is the (0,0)th component of

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

,

where bold utildet=(x0+y0- m- phi(y0)gt0,...,xK+yK- m- phi(yK)gtK).

We can use the expression (A1) for Pt above to rewrite this as the 0th component of

Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact help@nature.com or the author

where D~ is the diagonal matrix with d~kk=phi(yk). By the same argument as above, we can replace the inverse matrix in this expression with I, thus justifying (7).

Top

Acknowledgements

Financial support from the ESRC under research grant R000236821 is acknowledged. ESRC Research Fellowship Award H53627501695 supported the second author. We would like to thank Lisa Farrell, Richard Cornes and Ian Walker for discussions and comments, as well as seminar participants at Edinburgh, Keele, Oxford and Toulouse.

Extra navigation

.

Society resources

ADVERTISEMENT
JORS-Link to full archive