2.1 Discrete Choice Experiment design
Discrete Choice Experiments (DCEs) present participants with a series of choices between options that have multiple attributes and then preferences over the attributes are extracted using a technique based on Thurstone’s Random Utility Theory (RUT), which was extended by McFadden in the late 20th century [15]. One of RUT’s core assumptions is that choice is underpinned by a ‘utility’, which has a systematic and random component. Systematic components comprise latent values attached to attributes (and levels of each attribute) in choice alternatives, as well as covariates that are influential in selection. In contrast, random components comprise unknown factors that may impact decisions [16].
DCEs have been increasingly applied in public health research to understand citizens’ preferences with respect to interventions or programmes [17], and empirical research has demonstrated their robust external validity: for example, they have been used to accurately predict medical treatment [18] and vaccination behaviour [19]. During the pandemic, DCEs have been used to understand which vaccine characteristics are most influential in decision-making [20]; health and economic trade-offs in lockdowns [21]; and exit strategies [22].
Given its application across domains in public health research and strong external validity, a DCE was chosen as the core component within this study.
Attributes, experimental design and operationalisation
In contrast to many other public health DCEs – which often rely on primary qualitative research to inform design (see [17]) – this study’s DCE was envisaged and designed in collaboration with individuals responsible for the strategic planning of the UK’s mass vaccination programme. The design of the DCE was informed by three sessions in the weeks commencing 8th and 15th March 2021; representatives from Kantar Public UK, PHE, NHSE and DHSC were present at each session.
The final design comprised four attributes (1 x 5 levels, 3 x 3 levels), each of which can be seen in Table 1. The rationale for the definition of attributes and levels was as follows:
Mode of delivery. Mechanisms for vaccination through Primary Care Networks (PCNs) and Vaccination Centres (VCs) were already well established as part of the first phase of the programme, and it was assumed that these would continue to play an important role in providing access for young people. There was also recognition of the likely need for more targeted, localised services: recruitment of local pharmacies had recently begun, and there was interest in the potential utility of mobile services (which could be deployed in convenient locations) as well as drive-through options. Therefore, the rationale for the inclusion of this attribute was to understand whether these options should be employed – or scaled up or down – according to young people’s preferences.
Appointment time. Existing research shows that greater convenience plays a role in reducing hesitancy [23]; but evidence on whether extended opening hours are advantageous for young people specifically is lacking. Additionally, as noted, COVID-19 presents a novel situation in which past experiences might not apply. Services had been commissioned to operate seven days a week with extended opening hours to heighten convenience. The rationale for the inclusion of this attribute was therefore to identify any strong preferences among young people for appointments outside of normal working hours, further to inform resource planning and allocation.
Proximity. Evidence for a distance decay effect, whereby people who live further away from healthcare facilities have lower levels of usage after adjustment for need [24], suggests that more proximate vaccination locations will result in higher levels of uptake. Vaccination Centres had been specifically situated within 45 minutes of 99% of the population in England. As such, the rationale for the inclusion of this attribute was to identify any impacts of travel times ranging from 0 to 45 minutes on young people’s propensity to attend a vaccination appointment.
SMS invitation sender. Unpublished observations had indicated that SMS text messages from friends and family may be effective in motivating attendance at appointments among younger people; General Practice (GP) text messaging systems were already in operation, and a national SMS booking system was about to go live. Consequently, the rationale for the inclusion of this attribute was to provide evidence of any difference in motivational impact due to these SMS senders, to inform planning.
Table 1: DCE attributes and levels
Attribute
|
Level 1
|
Level 2
|
Level 3
|
Level 4
|
Level 5
|
Mode of delivery
|
Local vaccination centre
|
Nearby GP surgery (Primary Care Network)
|
Nearby pharmacy
|
Drive-thru
|
Mobile/ pop-up
|
Appointment time
|
Monday to Friday, 9am-5pm
|
Monday to Friday, after hours (before 9am or after 5pm)
|
Weekends
|
|
Proximity from one’s home
|
Less than 15 minutes
|
Between 15 and 30 minutes
|
Between 30 and 45 minutes
|
|
|
SMS invitation sender
|
NHS
|
Your GP
|
Best friend
|
|
|
Table 1 legend: the five attributes (and their levels) included in the DCE design
A full factorial design with these attributes/levels would have included n=135 profiles (51 * 33), the presentation of which would have been infeasible in this study. Therefore, a d-optimal fractional factorial design was generated using the choiceDes package in R Statistical Software.
The final design was unlabelled and comprised n=6 paired choice sets, each of which contained an opt-out to maximise external validity. Following the finalisation of the design, each of the choices was translated into an image of an invitation SMS message for presentation in the DCE, an example of which can be seen in Figure 1 below.
The use of pictorial choice options arguably delivers a more natural and relatable set of stimuli for participants than the more usual tabular format for DCEs. As Kahneman and Tversky have argued, “the method of hypothetical choices… relies on the assumption that people often know how they would behave in actual situations of choice, and on the further assumption that the subjects have no special reason to disguise their true preferences” [25]. The method used in this study brings participants closer to the “actual situation of choice” than is often the case with DCEs, enhancing its external validity.
Before starting the DCE, participants were provided with an overview of the scenario (detailing the vaccination programme) and their choice task, which involved selecting the vaccination appointment that they would be most likely to book based on its characteristics, or selecting ‘neither appointment’. The introduction to the DCE and the complete choice set can be seen in the Appendix, and an example of a paired choice set can be seen in Figure 2.
In the experiment, the order in which the pairs were shown to participants was randomised to minimise the influence of order effects [26].
2.2 Participants
Sample size requirement
There is no scientific consensus on the sample size required for a sufficiently powered DCE. However, rules of thumb have been proposed in the literature, the most common of which is that from Johnson and Orme [27]. According to the authors, the sample size required for a main effects DCE model can be calculated using the following equation:

Where: a represents the number of alternatives (2); t represents the number of choice tasks (6), and c represents the number of levels in the largest attribute (5). According to this rule of thumb, we required a minimum of 209 participants, a total which we exceeded in our final sample.
Sample profile
This study was conducted online from 25 March to 2 April 2021, with sample sourced from LifePoints (Kantar’s online access panel).
The sample for the experiment comprised n=2,012 adults aged 18-29 years who were currently living in the UK and had not been vaccinated at the time of interview. To ensure that the sample was nationally representative of this age group in terms of key demographic characteristics, we enforced flexible parallel quotas on age and ethnicity. These quotas were based on mid-year population statistics from the ONS [28]. Quota targets and achieved sample can be seen in Table 1.
Table 2: Quota targets
Age
|
Target %
|
Achieved %
|
Ethnicity
|
Target %
|
Achieved %
|
18-21
|
30%
|
32%
|
White
|
81%
|
81%
|
22-25
|
34%
|
33%
|
Asian, Black Mixed/Other
|
19%
|
19%
|
26-29
|
36%
|
35%
|
|
|
|
Table 2 legend: Age and ethnicity parallel quota targets, and the profile of the final achieved sample
As an incentive for participation in the study, all participants were provided with LifePoints’ reward points.
2.3 Statistical methods
Participants’ choices were analysed using mixed logit models (alternatively termed a random parameters model), adjusted for their panel nature. Attributes were set as random parameters – each with a normal distribution – to allow for preference heterogeneity across participants [29]. A likelihood ratio test was conducted to test for a difference in fit between a model which allowed for correlations between random parameters using Choleski decomposition; however, this test did not indicate significant improvement (χ2(55) = 31.149, p = 0.996), so the non-correlated model was selected for use. This model was estimated in R statistical software using the mlogit package [30], and can be written as:

Where:
α denotes the model alternative specific constant (ASC; the systematic preference for ‘opting in’ to appointment options).
β1 – β5 denote individual-specific coefficients representing the effect of vaccination delivery modes (w1 represents ‘Vaccination centre’, w2 represents ‘GP surgery’, w3 represents ‘Nearby pharmacy’, w4 represents ‘Drive-thru’, w5 represents ‘Mobile/pop-up’) on selection.
β6 – β8 denote individual-specific coefficients representing the effect of appointment times (x1 represents ‘Monday to Friday, 9am-5pm’; x2 represents ‘Monday to Friday, after hours (before 9am or after 5pm)’, x3 represents ‘Weekends’) on selection.
β9 – β11 denote individual-specific coefficients representing the effect of venue proximity (y1 represents ‘Less than 15m’, y2 represents ‘Between 15 and 30m’, y3 represents ‘Between 30 and 45m’) on selection.
β12 – β14 denote individual-specific coefficients representing the effect of invitation sources (z1 represents ‘NHS’, z2 represents ‘YourGP’, z3 represents ‘Best friend’) on selection.
ε is the random error term, representing the non-systematic component in selection.
All attributes were dummy coded, such that 1 represented their presence in each choice card, while 0 represented their absence. Coefficients’ signs reflect whether a level has a positive or a negative effect on utility compared to the reference category; further, their absolute values indicate their relative importance in selection, again compared to the reference category. To facilitate ease of interpretation, coefficients were exponentiated to generate odds ratios.