Hostname: page-component-7c8c6479df-7qhmt Total loading time: 0 Render date: 2024-03-19T05:34:06.226Z Has data issue: false hasContentIssue false

Injection fears and COVID-19 vaccine hesitancy

Published online by Cambridge University Press:  11 June 2021

Daniel Freeman*
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK NIHR Oxford Health Biomedical Research Centre (BRC), Oxford, UK
Sinéad Lambe
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK
Ly-Mee Yu
Affiliation:
Nuffield Department of Primary Care, University of Oxford, Oxford, UK
Jason Freeman
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK
Andrew Chadwick
Affiliation:
Department of Communication and Media, Online Civic Culture Centre, Loughborough University, Loughborough, UK
Cristian Vaccari
Affiliation:
Department of Communication and Media, Online Civic Culture Centre, Loughborough University, Loughborough, UK
Felicity Waite
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK NIHR Oxford Health Biomedical Research Centre (BRC), Oxford, UK
Laina Rosebrock
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK
Ariane Petit
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK
Samantha Vanderslott
Affiliation:
Oxford Vaccine Group, Department of Paediatrics, University of Oxford, Oxford, UK
Stephan Lewandowsky
Affiliation:
School of Psychological Science, University of Bristol, Bristol, UK
Michael Larkin
Affiliation:
Department of Psychology, Life and Health Sciences, Aston University, Birmingham, UK
Stefania Innocenti
Affiliation:
Smith School of Enterprise and the Environment, University of Oxford, Oxford, UK
Helen McShane
Affiliation:
The Jenner Institute, Nuffield Department of Medicine, University of Oxford, Oxford, UK NIHR Oxford Biomedical Research Centre (BRC), Oxford, UK
Andrew J. Pollard
Affiliation:
Oxford Vaccine Group, Department of Paediatrics, University of Oxford, Oxford, UK NIHR Oxford Biomedical Research Centre (BRC), Oxford, UK
Bao Sheng Loe
Affiliation:
The Psychometrics Centre, University of Cambridge, Cambridge, UK
*
Author for correspondence: Professor Daniel Freeman, Department of Psychiatry, University of Oxford, Warneford Hospital, Oxford, UK, OX3 7JX. E-mail: daniel.freeman@psych.ox.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Background

When vaccination depends on injection, it is plausible that the blood-injection-injury cluster of fears may contribute to hesitancy. Our primary aim was to estimate in the UK adult population the proportion of COVID-19 vaccine hesitancy explained by blood-injection-injury fears.

Methods

In total, 15 014 UK adults, quota sampled to match the population for age, gender, ethnicity, income and region, took part (19 January–5 February 2021) in a non-probability online survey. The Oxford COVID-19 Vaccine Hesitancy Scale assessed intent to be vaccinated. Two scales (Specific Phobia Scale-blood-injection-injury phobia and Medical Fear Survey–injections and blood subscale) assessed blood-injection-injury fears. Four items from these scales were used to create a factor score specifically for injection fears.

Results

In total, 3927 (26.2%) screened positive for blood-injection-injury phobia. Individuals screening positive (22.0%) were more likely to report COVID-19 vaccine hesitancy compared to individuals screening negative (11.5%), odds ratio = 2.18, 95% confidence interval (CI) 1.97–2.40, p < 0.001. The population attributable fraction (PAF) indicated that if blood-injection-injury phobia were absent then this may prevent 11.5% of all instances of vaccine hesitancy, AF = 0.11; 95% CI 0.09–0.14, p < 0.001. COVID-19 vaccine hesitancy was associated with higher scores on the Specific Phobia Scale, r = 0.22, p < 0.001, Medical Fear Survey, r = 0.23, p = <0.001 and injection fears, r = 0.25, p < 0.001. Injection fears were higher in youth and in Black and Asian ethnic groups, and explained a small degree of why vaccine hesitancy is higher in these groups.

Conclusions

Across the adult population, blood-injection-injury fears may explain approximately 10% of cases of COVID-19 vaccine hesitancy. Addressing such fears will likely improve the effectiveness of vaccination programmes.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

Vaccination is now a key global public health intervention to combat the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. A variable proportion of the population in each country will seek to delay or avoid being vaccinated, thereby limiting the success of vaccination programmes. In the Oxford Coronavirus Explanations, Attitudes, and Narratives Survey II (OCEANS-II), key factors associated with vaccine hesitancy in the UK population were lower perception of risk from the virus, less awareness of the collective benefits of vaccination, doubts about the efficacy of vaccination and worry about potential side effects, particularly in the context of the rapid development and testing of the vaccines (Freeman et al., Reference Freeman, Loe, Chadwick, Vaccari, Waite, Rosebrock and Lambe2020). In subsequent qualitative interviews, participants raised an issue that had not been assessed in the survey: anxiety about injections. Fears about injection in the short term might prevent the long-term benefits of vaccination. In children, adolescents and adults, fear of injection falls within the unitary specific phobia subtype of blood-injection-injury fears (Kendler, Aggen, Werner, & Fried, Reference Kendler, Aggen, Werner and Fried2020; Loken, Hettema, Aggen, & Kendler, Reference Loken, Hettema, Aggen and Kendler2014; Muris, Schmidt, & Merckelbach, Reference Muris, Schmidt and Merckelbach1999; Wenzel & Holt, Reference Wenzel and Holt2003). This is persistent excessive fear of blood, needles or invasive medical procedures, leading to avoidance or endurance with intense anxiety. Uniquely among anxiety disorders, part of the typical response pattern for this phobia is a drop in heart rate and blood pressure leading to fainting. A fear of fainting may cause great reluctance to join a line of people waiting for a vaccine injection. In our latest UK survey (Freeman et al., Reference Freeman, Loe, Yu, Freeman, Chadwick, Vaccari and Lambe2021), we therefore included assessment of blood-injection-injury fears in order to determine the degree to which they may be a factor in hesitancy about COVID-19 vaccination, which currently depends on injection.

There is a spectrum of severity of blood-injection-injury fears in the general population, with heritability estimated at approximately 40% (Van Houtem et al., Reference Van Houtem, Laine, Boomsma, Ligthart, van Wijk and De Jongh2013). Determined by diagnostic interview, lifetime prevalence of blood-injection-injury phobia in the adult population is approximately 3–4% (Stinson et al., Reference Stinson, Dawson, Chou, Smith, Goldstein, Ruan and Grant2007; Wardenaar et al., Reference Wardenaar, Lim, Al-Hamzawi, Alonso, Andrade, Benjet and De Jonge2017). Rates are higher in women and in those of younger age. Self-identified fear of needles is much more common. One convenience sample study found that over 20% of parents and over 60% of children reported a fear of needles. Such fears were identified as the primary reason for immunisation non-compliance for 7% of parents and 8% of children (Taddio et al., Reference Taddio, Ipp, Thivakaran, Jamal, Parikh, Smart, Sovran and Katz2012). Fear of injection has been cited by both the general public and health professionals as a reason for vaccination refusal, including for influenza, tetanus and pneumococcal infections (Johnson, Nichol, & Lipczynski, Reference Johnson, Nichol and Lipczynski2008; Yaqub, Castle-Clarke, Sevdalis, & Chataway, Reference Yaqub, Castle-Clarke, Sevdalis and Chataway2014). In a US survey conducted in June 2020, 11.8% of those who were hesitant about COVID-19 vaccination reported that a reason was a dislike of needles and injections (Ruiz & Bell, Reference Ruiz and Bell2021). A higher proportion (43.8%) identified fear of dangerous side effects as the reason. In a US study of 9000 older adults conducted in November 2020, 1.7% were concerned about receiving a COVID-19 vaccine because of a fear of needles (Nikolovski et al., Reference Nikolovski, Koldijk, Weverling, Spertus, Turakhia, Saxon and Navar2021). These studies suggest that, although fear of injection is not the dominant reason for vaccine hesitancy, it may be a contributory factor.

Our primary objective in this study was to determine the degree to which blood-injection-injury fears may explain COVID-19 vaccine hesitancy in the UK adult population. We also set out to estimate the prevalence of significant blood-injection-injury fears in the UK adult population. Furthermore, we wanted to test the extent to which injection fears may account for socio-demographic (age, gender, ethnicity and income) differences in COVID-19 vaccine hesitancy. Rather than merely asking a single question about fear of needles, we used established assessments. Our focus was on self-report, rather than diagnostic interview, since it is likely to be subjective thoughts and feelings across the spectrum of fear of injection that may have an impact regardless of whether they match diagnostic criteria. In a large survey of the adult general population, we included a self-report assessment of blood-injection-injury fears that can be used as a screening tool for the potential presence of a phobia. To validate results, we included an additional blood-injection fears questionnaire. We expected blood-injection-injury fears to explain a small amount of variance in vaccine hesitancy across the population. However, individuals with pronounced fears would be substantially more likely to be hesitant about taking a vaccine. We also hypothesised that injection fears may help explain (i.e. partially mediate) why vaccine hesitancy is somewhat higher in particular demographic groups (younger age, females, ethnic minority groups and lower income).

Methods

Participants

Oxford Coronavirus Explanations, Attitudes, and Narratives Survey (OCEANS) III (Freeman et al., Reference Freeman, Loe, Yu, Freeman, Chadwick, Vaccari and Lambe2021) is an online survey with a quota sampled UK participant group of 15 014 adults (18+ years old), conducted from 19 January 2021 to the 5 February 2021 via a market research company, Lucid. The quotas were based on UK Office for National Statistics population estimate data for gender, age, ethnicity, income and region. We did not wish to focus in this study only on those who had not been vaccinated by a particular point in time, since such results would be less generalisable, COVID-19 vaccination is highly unlikely to be a ‘one off’ event (there are two doses and it is likely to need renewal), and people who have been vaccinated still influence the rest of the population. Invited respondents did not know the topic of the survey before provisional agreement to complete it. They were simply told that there was a new survey and informed of the time period for it to be completed. Only after agreeing to participate did they see the online introduction. There was no mention of phobias in the rationale for the study provided, which was ‘There are now approved vaccines for COVID-19 that will be rolled out in the UK over the coming months. We want to learn about people's views about vaccination for COVID-19. In particular, we want to find out how many people would or would not wish to be vaccinated and the reasons behind their decision’. OCEANS-III was approved by the University of Oxford Central University Research Ethics Committee, and all participants provided informed consent online.

Assessments

Initial vaccine hesitancy question

After agreeing to take part, participants completed a single question: ‘If the vaccine was available at my GP surgery I would: 1. Get it as soon as possible/2. Get it when I have time/3. Delay getting it/4. Avoid getting it for as long as possible/5. Never get it/6. Don't know’. This question was found in our OCEANS-II study to have the highest loading (0.95) on the COVID-19 vaccine hesitancy latent factor but for that precise reason was not included in the primary outcome measure (see Freeman et al., Reference Freeman, Loe, Chadwick, Vaccari, Waite, Rosebrock and Lambe2020 supplementary materials).

Oxford COVID-19 Vaccine Hesitancy Scale (Freeman et al., Reference Freeman, Loe, Chadwick, Vaccari, Waite, Rosebrock and Lambe2020)

Seven items [e.g. ‘Would you take a COVID-19 vaccine (approved for use in the UK) if offered?’] are rated on item specific response options (Saris, Krosnick, Revilla, & Shae, Reference Saris, Krosnick, Revilla and Shae2010), coded from 1 to 5. A ‘Don't know’ option is also provided, which is excluded from scoring. Higher scores indicate a higher level of vaccine hesitancy. In the current study the Cronbach's alpha of the scale was 0.97 (N = 15 014).

Specific Phobia Scale – blood-injection-injury (Ovanessian et al., Reference Ovanessian, Fairbrother, Vorstenbosch, McCabe, Rowa and Antony2019)

This is a self-report screening scale for different types of specific phobia based on the DSM-V criteria (APA, 2013). We used the 14 items assessing blood-injection-injury phobia. Each item (e.g. ‘Watching someone else get an injection’, ‘Getting minor surgery’, ‘Giving blood’ and ‘Receiving an injection’) is rated on a 5 point rating scale for the degree of fear triggered [0 (no fear) to 4 (extreme fear)]. Scores on the subscale can range from 0 to 56, with higher scores indicating a higher level of fear. A score of 20 or higher indicates the possible presence of a blood-injection-injury specific phobia (Ovanessian et al., Reference Ovanessian, Fairbrother, Vorstenbosch, McCabe, Rowa and Antony2019). In the current study, the Cronbach's alpha of the scale was 0.94 (N = 15 014).

Medical fear survey – short version – injections and blood (Olatunji et al., Reference Olatunji, Ebesutani, Sawchuk, McKay, Lohr and Kleinknecht2012)

Each of the four items [e.g. ‘Receiving a hypodermic (i.e. into the skin) injection in the arm’, ‘Having blood drawn from your arm’] is rated on a 4-point rating scale (0 = no fear to 3 = intense fear) for the degree of fear provoked. Scores can range from 0 to 12, with higher scores indicating greater fear. In the current study, the Cronbach's alpha of the scale was 0.90 (N = 15 014).

Injection fears

We also created from the two scales a factor score specifically for injection fears, using the four items about injection [‘Watching someone else get an injection’, ‘Receiving an injection’, ‘Receiving a hypodermic (i.e. into the skin) injection in the arm’, ‘Seeing someone receiving an injection in the arm’]. Higher scores indicate greater fear of injection.

Statistical analysis plan

Statistical analyses were carried out in the R software environment for statistical computing and graphics, Version 4.0.3.

Structural equation modelling (SEM) was used to evaluate the factorial structure of the Oxford COVID-19 Vaccine Hesitancy Scale, Specific Phobia Scale, Medical Fear Survey and the combined four items measuring fear of injection, and to assess the potential mediation effect of fear of injection for associations between vaccine hesitancy and demographic factors. SEM is a collection of statistical techniques composed of two parts. The first part is the confirmatory measurement model, also known as the confirmatory factor analysis (CFA), which estimates the relations among latent constructs and their observed indicators. We assessed the measurement models of the existing and new injection fear scale based on factor loadings, inter-correlation between factors and several goodness-of-fit indices to ensure appropriate specification. The second part is the structural model, which estimates the relations among constructs (Kline, Reference Kline2015). This combination allows the elimination of unreliability of measurement in the models and can be used to make assumptions about how latent constructs derived from observed variables are associated with each other. The mediation analysis was also conducted under the SEM framework and the advantages of such an approach are discussed in Gunzler, Chen, Wu, and Zhang (Reference Gunzler, Chen, Wu and Zhang2013).

With non-normal distributions and missing values within the data, which in this case is considered either missing completely at random or missing at random due to test design (Little & Rubin, Reference Little and Rubin2019; Rubin, Reference Rubin1976), we utilised the full information maximum likelihood estimation procedure with robust (Huber–White) standard errors and a scaled test statistic to estimate the SEMs, which has been described as a more efficient and less biased approach than other ad hoc missing data techniques (Enders, Reference Enders2001; Enders & Bandalos, Reference Enders and Bandalos2001). The SEM supports the data when the root mean square error of approximation (RMSEA) and standardised root mean square residual (SRMR) are <0.08 (Kline, Reference Kline2005). Additional fit indices included the comparative fit index (CFI) and Tucker–Lewis Index (TLI), all of which should exceed 0.90 (Kline, Reference Kline2015), with RMSEA < 0.06, and CFI and TLI > 0.95 indicating good model data fit (Hu & Bentler, Reference Hu and Bentler1999). The confidence intervals (CIs) of the indirect effects in the mediational model using SEM were tested using a Monte Carlo method with 20 000 replications as described by MacKinnon, Lockwood, and Williams (Reference MacKinnon, Lockwood and Williams2004). The Holm method was used to adjust the p values in the mediational model using SEM to control the type I error rate due to multiple comparisons (Holm, Reference Holm1979). The lavaan and the semTools R packages were used to conduct the SEM analyses (Jorgensen, Pornprasertmanit, Schoemann, & Rosseel, Reference Jorgensen, Pornprasertmanit, Schoemann and Rosseel2020; Rosseel, Reference Rosseel2012).

We used the population attributable fraction (PAF) to evaluate vaccine hesitancy due to blood-injection-injury phobia. A historical review of the PAF is described by Poole (Reference Poole2015). The general idea of PAF is to understand the proportion of unfavourable outcomes that would have been prevented if the exposure of interest was eliminated from the population (Levin, Reference Levin1953). Analysis of PAF was conducted using the AF R package (Dahlqwist, Zetterqvist, Pawitan, & Sjölander, Reference Dahlqwist, Zetterqvist, Pawitan and Sjölander2016). We dichotomised (1, 2 = no hesitancy, 3, 4, 5 = hesitancy) the single vaccine hesitancy item that opened the survey (‘If the vaccine was available at my GP surgery I would…’). Participants who answered 6 (‘don't know’) were defined as missing and excluded from the analysis. The presence or absence of having a blood-injection-injury phobia was calculated based on the total score of the Specific Phobia Scale. Blood-injection-injury phobia was used as the exposure variable and was dichotomised according to the recommended cut-offs by Ovanessian et al. (Reference Ovanessian, Fairbrother, Vorstenbosch, McCabe, Rowa and Antony2019). We used a logistic regression model to first estimate the log odds of vaccine hesitancy levels based on the presence of a phobia and included demographic associations (age, ethnicity and gender) as confounding controls. These demographic characteristics were known to be significant predictors as reported in OCEANS-II (Freeman et al., Reference Freeman, Loe, Chadwick, Vaccari, Waite, Rosebrock and Lambe2020). The beta coefficient estimates from the logistic regression model were then used to estimate the adjusted PAF with the formula described by Dahlqwist et al. (Reference Dahlqwist, Zetterqvist, Pawitan and Sjölander2016). The adjusted PAF was calculated as the proportion of vaccine hesitancy in the population that is attributable to blood-injection-injury phobia (i.e. the proportion of cases that might be prevented in the population if the exposure variable was eliminated), after accounting for demographic associations. The standard errors were estimated using the sandwich formula with the delta method (Greenland & Drescher, Reference Greenland and Drescher1993; Sjölander & Vansteelandt, Reference Sjölander and Vansteelandt2011).

Results

A summary of the socio-demographic characteristics of the participants is provided in Table 1. For the initial vaccine hesitancy question – ‘If the vaccine was available at my GP surgery’ – 11 012 (73.3%) individuals reported that they would get it as soon as possible, 1451 (9.7%) reported that that they would get it when they have time, 672 (4.5%) reported that they would delay getting it, 771 (5.1%) reported that they would avoid getting it for as long as possible, and 632 (4.2%) said they would never get it. In total, 476 (3.2%) said that they did not know. Hence, 2075 (13.8%) participants were classed as hesitant (scoring 3, 4 or 5).

Table 1. Socio-demographic information (N = 15 014)

Prevalence of blood-injection-injury fears

There was a high correlation between fears as assessed with the Specific Phobia Scale and Medical Fears Survey, r = 0.81, p < 0.001, n = 15 014. The injection factor score correlated highly with the Specific Phobia Scale score, r = 0.89, p < 0.001, n = 15 014 and Medical Fears Survey score, r = 0.92, p < 0.001, n = 15 014. The mean score of the participant group on the specific phobia scale was 13.19, s.d. = 11.38, n = 15 014. In total, 3927 (26.2%) screened positive for a specific blood-injection-injury phobia and 11 097 (73.8%) screened negative.

Creation of factor scores

Participants who had four or more missing items (i.e. answered don't know to over half the items) on the Oxford COVID-19 Vaccine Hesitancy Scale were removed from further analysis. This resulted in a sample size of 14 820 individuals, which was used for the CFA models. The CFA model for the Oxford COVID-19 Vaccine Hesitancy Scale indicated an excellent model fit (χ2 = 1112.84, df = 14, p < 0.001, CFI = 0.98, TLI = 0.97, RMSEA = 0.07, SRMR = 0.01). The initial CFA model indicated that the 14-item Specific Phobia Scale was unsatisfactory (χ2 = 11 328.93, df = 77, p < 0.001, RMSEA = 0.10, SRMR = 0.05, CFI = 0.86, TLI = 0.83). Therefore, in a post-hoc analysis, we evaluated the model adequacy based on the modification index (MI). MI is an estimate of the amount of chi-square which would be reduced if a parameter was added or the restriction of a specific parameter was removed in the CFA model (Sörbom, Reference Sörbom1989). Specifically, we investigated MIs of correlated measurement errors between given pairs of items. Freely estimating covariances of correlated measurement errors with high MI results in a substantial decrease in the chi-square, thereby increasing model fit. However, the use of correlated measurement errors comes at an expense of theoretical accuracy since the covariance of the error term is from at least one unknown common source (Fornell, Reference Fornell1983; Gerbing & Anderson, Reference Gerbing and Anderson1984). It indicates that the specific pair of items has unexplained common sources of variance that are not captured by the latent variable. Given that the addition of the correlated measurement errors does not improve the interpretation of the model substantially, we instead used an alternative strategy, which was to remove one of the paired items that displayed high MI. We re-ran the analysis after the removal of each item to inspect the model's goodness-of-fit, which resulted in four items being removed (item 3, 8, 10, 14). The final model indicated acceptable goodness of fit results (χ2 = 2831.13, df = 35, p < 0.001, CFI = 0.93, TLI = 0.92, RMSEA = 0.07, SRMR = 0.04). The CFA model for the Medical Fear Survey indicated an excellent model fit (χ2 = 187.50, df = 2, p < 0.001, CFI = 0.99, TLI = 0.96, RMSEA = 0.08, SRMR = 0.02). The CFA model for the four-item injection fears subscale indicated an excellent model fit (χ2 = 174.98, df = 2, p < 0.001, CFI = 0.99, TLI = 0.96, RMSEA = 0.08, SRMR = 0.02).

Associations of blood-injection-injury fears and vaccine hesitancy

The Oxford COVID-19 Vaccine Hesitancy Scale factor score correlated positively with the Specific Phobia Scale factor score, 0.22, p < 0.001, n = 14 820; the Medical Fear Survey factor score, r = 0.23, p < 0.001, n = 14 820; and the injection fears factor score, r = 0.25, p < 0.001, n = 14 820. (In a sensitivity analysis, the size of the associations was unchanged when controlling for condition in the OCEANS-III trial.) The individuals screening positive for a specific phobia were more likely to be hesitant (830/3754; 22.1%) on the survey opening hesitancy item (scoring 3, 4 or 5) than those screening negative for a specific phobia (1245/10 784; 11.5%), odds ratio = 2.18, 95% CI 1.97–2.40, p < 0.001.

Explaining vaccine hesitancy and demographic associations

Mediation analysis using SEM was used to investigate the hypothesis that fear of injection partially mediates the effect of demographic associations on vaccine hesitancy. Hence, we tested a mediational model under the SEM framework consisting of age, gender, ethnicity and income. The sample size used in the analysis was n = 14 149. Participants who had missing data from the demographic variables and those who selected ‘other’ for gender (n = 32) were excluded. We also combined groups within ethnicity and income to increase the sample size at a group level. A breakdown of the demographic information is presented in Table 2.

Table 2. Socio-demographic information for mediation analysis (N = 14 149)

The mediation results are reported in Table 3. According to the mediation model, while holding demographic information constant, participants with a fear of injection had higher levels of vaccine hesitancy (B = 0.20, adjusted p < 0.01). The mediation model indicated that gender was not a significant predictor of fear of injection and vaccine hesitancy after controlling for age, ethnicity and income. Thus, further mediation analysis for gender was not conducted.

Table 3. Mediation results using SEM

CI, confidence intervals; s.e., standard error.

a Monte Carlo CI. Reference group for ethnicity: White. Reference group for income: less than 20k. Adjusted p values using the Holm method.

Age was a significant predictor of fear of injection (B = −0.015, adjusted p < 0.01) and the indirect effect of age and vaccine hesitancy via fear of injection was also found to be significant (B = −0.003, adjusted p < 0.01). Age remained a significant predictor of vaccine hesitancy after controlling for fear of injection (mediator), indicating a partial mediation (B = −0.013, adjusted p < 0.01).

The mediation for ethnicity indicated that relative to White ethnicity, only Black and Asian ethnicities were significant predictors of fear of injection (B = 0.208, adjusted p < 0.01; B = 0.168, adjusted p < 0.01). The relative indirect effects between ethnicity and vaccine hesitancy via fear of injection were also found to be significant only for Black and Asian ethnicities (B = 0.042, adjusted p < 0.01; B = 0.034, adjusted p < 0.01). The relative direct effects between ethnicity and vaccine hesitancy remained significant for Black and Asian ethnicities (B = 0.70, adjusted p < 0.01; B = 0.14, adjusted p < 0.01) after controlling for the fear of injection, indicating a partial mediation for these ethnic groups. Relative to White ethnicity, White other, Other mixed background and Other ethnic groups were also associated with higher levels of vaccine hesitancy (B = 0.22, adjusted p < 0.01; B = 0.36, adjusted p < 0.05; B = 0.25, adjusted p < 0.01), but the relationships between these ethnic groups and vaccine hesitancy were not (partially) mediated by fear of injection.

Relative to the income group with less than £20k, only those in the greater than £60k income group was a significant predictor of fear of injection (B = −0.096, adjusted p < 0.01). The relative indirect effect between those in the greater than £60k income group and vaccine hesitancy via fear of injection was also significant (B = −0.02, adjusted p < 0.01). The relative direct effect for the greater than £60k income group remained significant (B = −0.30, adjusted p < 0.01) after controlling for the fear of injection, indicating a partial mediation. Relative to the less than £20k income group, those in the £20k to £40k and £41k to £60k income groups were associated with lower levels of vaccine hesitancy (B = −0.11, adjusted p < 0.01; B = −0.25, adjusted p < 0.01), but the relationships were not (partially) mediated by fear of injection. The relative direct effects for those in the Prefer not to say income group were not significant for both vaccine hesitancy and fear of injection. Thus, further mediation analysis for the Prefer not to say income group was not necessary. The total effect combining all (relative) direct and indirect effects was significant (B = 1.011, adjusted p < 0.01). Approximately 16.1% (R 2 = 0.161) of the variance in vaccine hesitancy was accounted for by the demographic variables in the mediation model. The model showed an excellent fit (χ2 = 2701.92, df = 142, p < 0.001, RMSEA = 0.04, SRMR = 0.01, CFI = 0.98, TLI = 0.97). Figure 1 presents the mediational model using SEM.

Fig. 1. Mediational model using SEM for demographic associations with vaccine hesitancy. FOI, fear of injection; VAC.HES, vaccine hesitancy; ETH.B, Black ethnicity; ETH.WO, White other ethnicity; ETH.OM, other mixed background ethnicity; ETH.A, Asian ethnicity; ETH.O, other ethnicity; INC.2040, income group between £20k and £40k; INC.4060, income group between £40k and £60k; INC.G60, income group greater than £60k; INC.PNTS, prefer not to say income group. Reference group for ethnicity is White. Reference group for gender is male.

Population attributable fraction

The PAF analysis was conducted based on a sample size of n = 13 875 after removing participants who answered ‘don't know’ to the initial vaccine hesitancy item and those with missing responses to the demographics information, including those who answered ‘other’ to gender. According to the logistic model, the log of the odds of people being vaccine hesitant was positively related to the presence of blood-injury-injection phobia (B = 0.44, s.e. = 0.06, 95% CI 0.33–0.55, p < 0.001). The odds for participants with a blood-injection-injury phobia being vaccine hesitant were 1.55 ( = e 0.4409) times greater than the odds for participants with no blood-injection-injury phobia. The estimated PAF indicated that approximately 11.5% of all participants who are vaccine hesitant might have been prevented if fear of injection was absent (AF = 0.1149; s.e. = 0.02, 95% CI 0.09–0.14, p < 0.001). The odds ratios for different factors associated with vaccine hesitancy are summarised in a forest plot generated from a single model that contained all factors (see Fig. 2).

Fig. 2. Forest plot of associations with vaccine hesitancy (all factors included in the same model). Reference group for gender is male and for ethnicity is White.

Discussion

Blood-injection-injury fears are remarkably common in the UK adult general population. Our findings clearly show that these fears have a small to medium association with COVID-19 vaccine hesitancy. Approximately 5% of variance in COVID-19 hesitancy scores across the adult population may be accounted for by injection fears. At an individual level, screening positive for a blood-injection-injury phobia doubles the odds of also reporting COVID-19 vaccine hesitancy. Around a quarter of the population screen positive for potentially having such a phobia. A PAF calculation indicated that if blood-injection-injury phobia were eliminated, then just over 10% of instances of people who are COVID-19 vaccine hesitant might also be removed. It was also found that two demographic factors repeatedly found to be associated with vaccine hesitancy (age and ethnicity) were explained to a small degree by higher levels of injection fears in these groups. Overall, the results indicate that fear of injection may have a small role to play in the occurrence of vaccine hesitancy in adults. It is a highly plausible causal connection, and one which we find is voiced by individuals who are vaccine hesitant. However, the key caveat of this study is that all the analyses were correlational, and causation cannot be determined. Only future studies that lessen fears of injection and test effects on vaccine hesitancy can determine causation with certainty.

Our survey used a non-probability online quota sampling method, which will have introduced bias into who was approached to take part, and hence caution is required particularly for prevalence estimates. Injection fears will also be higher in children. Given the plausible causal connection, however, in our view the study findings indicate that injection fears should be addressed in vaccination programmes. The fewer the barriers to vaccine acceptance, the more successful the programme. Fear of injection is likely to be a barrier for some individuals but it can be addressed in a number of ways. First, alternative needle-free delivery routes of COVID-19 vaccination could be made available. For example, there has been a nasal spray developed for childhood flu and oral vaccines developed for rotavirus and for polio. Second, administration of the vaccine by injection could be made more tolerable, producing positive narratives in the general public of the process. For example, looking away from the injection may reduce fear (Mithal et al., Reference Mithal, Simmons, Cornelissen, Wong, Riddell, McMurtry and Taddio2018). Pain mitigation strategies (e.g. acute exercise) may also be helpful (Edwards & Booy, Reference Edwards and Booy2019). Third, there could be wider provision of the brief psychological treatments (exposure therapy and applied tension) for blood-injection-injury fears that show large treatment effects (Ayala, Meuret, & Ritz, Reference Ayala, Meuret and Ritz2009; Hiermeier & Mofrad, Reference Hiermeier and Mofrad2020; McMurtry et al., Reference McMurtry, Taddio, Noel, Antony, Chambers, Asmundson and Scott2016). The success of vital COVID-19 vaccination programmes is dependent on uptake. As such, we should not underestimate the importance of making the thought of the jab less anxiety-provoking for the millions of people who are fearful of injections.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721002609.

Acknowledgements

The study was jointly funded by the NIHR Oxford Biomedical Research Centre (IS-BRC-1215-20008) and NIHR Oxford Health Biomedical Research Centre (BRC-1215-20005). Daniel Freeman is a National Institute for Health Research (NIHR) Senior Investigator. The views expressed are those of the authors and not necessarily those of the National Health Service, NIHR or the Department of Health and Social Care. Stephan Lewandowsky received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 964728 (JITSUVAX). Felicity Waite is funded by a Wellcome Trust Clinical Doctoral Fellowship (102176/B/13/Z).

Conflict of interest

AJP is Chair of UK Dept. Health and Social Care's (DHSC) Joint Committee on Vaccination & Immunisation (JCVI), but does not participate in discussions on COVID-19 vaccines, and is a member of the WHO's SAGE. The views expressed in this paper do not necessarily represent the views of DHSC, JCVI or WHO. Oxford University has entered into a partnership with Astra Zeneca for the development of a coronavirus vaccine.

References

American Psychiatric Association (2013). Diagnostic and statistical manual of mental disorders (5th ed.). Washington, DC: American Psychiatric Publishing.Google Scholar
Ayala, E. S., Meuret, A. E., & Ritz, T. (2009). Treatments for blood-injury-injection phobia. Journal of Psychiatric Research, 43, 12351242.CrossRefGoogle ScholarPubMed
Dahlqwist, E., Zetterqvist, J., Pawitan, Y., & Sjölander, A. (2016). Model-based estimation of the attributable fraction for cross-sectional, case–control and cohort studies using the R package AF. European Journal of Epidemiology, 31, 575582.CrossRefGoogle Scholar
Edwards, K. M., & Booy, R. (2019). Pain control at the vaccine injection site: New insights. International Journal of Travel Medicine, 26(6), taz053. https://doi.org/10.1093/jtm/taz053.CrossRefGoogle ScholarPubMed
Enders, C. K. (2001). The impact of nonnormality on full information maximum-likelihood estimation for structural equation models with missing data. Psychological Methods, 6, 352.CrossRefGoogle ScholarPubMed
Enders, C. K., & Bandalos, D. L. (2001). The relative performance of full information maximum likelihood estimation for missing data in structural equation models. Structural Equation Modeling, 8, 430457.CrossRefGoogle Scholar
Fornell, C. (1983). Issues in the application of covariance structure analysis: A comment. Journal of Consumer Research, 9(4), 443448.CrossRefGoogle Scholar
Freeman, D., Loe, B. S., Chadwick, A., Vaccari, C., Waite, F., Rosebrock, L., … Lambe, S. (2020). COVID-19 vaccine hesitancy in the UK: The Oxford Coronavirus Explanations, Attitudes, and Narratives Survey (OCEANS) II. Psychological Medicine. https://doi.org/10.1017/S0033291720005188.Google Scholar
Freeman, D., Loe, B. S., Yu, L.-M., Freeman, J., Chadwick, A., Vaccari, C., … Lambe, S. (2021). Effects of different types of written vaccination information on COVID-19 vaccine hesitancy in the UK (OCEANS-III): A single-blind, parallel-group, randomised controlled trial. Lancet Public Health, 6, e416e427.CrossRefGoogle ScholarPubMed
Gerbing, D. W., & Anderson, J. C. (1984). On the meaning of within-factor correlated measurement errors. Journal of Consumer Research, 11, 572580.CrossRefGoogle Scholar
Greenland, S., & Drescher, K. (1993). Maximum likelihood estimation of the attributable fraction from logistic models. Biometrics, 49, 865872.CrossRefGoogle ScholarPubMed
Gunzler, D., Chen, T., Wu, P., & Zhang, H. (2013). Introduction to mediation analysis with structural equation modeling. Shanghai Archives of Psychiatry, 25, 390.Google ScholarPubMed
Hiermeier, U. M., & Mofrad, L. (2020). Feasibility of one-session treatment for specific needle phobia in an adult IAPT service – A case series. The Cognitive Behaviour Therapist, 13, e51.CrossRefGoogle Scholar
Holm, S. (1979). A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics, 6, 6570.Google Scholar
Hu, L., & Bentler, P. M. (1999). Cutoff criteria for fit indexes in covariance structure analysis: Conventional criteria versus new alternatives. Structural Equation Modeling, 6, 155.CrossRefGoogle Scholar
Johnson, D. R., Nichol, K., & Lipczynski, K. (2008). Barriers to adult immunization. American Journal of Medicine, 121(S2), s280s35.CrossRefGoogle ScholarPubMed
Jorgensen, T. D., Pornprasertmanit, S., Schoemann, A. M., & Rosseel, Y. (2020). semTools: useful tools for structural equation modeling. R package version 0.5-3. Retrieved from https://CRAN.R-project.org/package=semTools.Google Scholar
Kendler, K., Aggen, S., Werner, M., & Fried, E. (2020). A topography of 21 phobic fears: Network analysis in an epidemiological sample of adult twins. Psychological Medicine. doi: https://doi.org/10.1017/S0033291720004493.Google Scholar
Kline, R. B. (2015). Principles and practice of structural equation modeling. New York: Guilford Publications.Google Scholar
Kline, T. (2005). Psychological testing: A practical approach to design and evaluation. Thousand Oaks, CA: Sage.CrossRefGoogle Scholar
Levin, M. L. (1953). The occurrence of lung cancer in man. Acta-Unio Internationalis Contra Cancrum, 9, 531541.Google ScholarPubMed
Little, R. J., & Rubin, D. B. . (2019). Statistical analysis with missing data (Vol. 793). New Jersey: John Wiley & Sons.Google Scholar
Loken, E., Hettema, J., Aggen, A., & Kendler, K. (2014). The structure of genetic and environmental risk factors for fears and phobias. Psychological Medicine, 44, 23752384.CrossRefGoogle ScholarPubMed
MacKinnon, D. P., Lockwood, C. M., & Williams, J. (2004). Confidence limits for the indirect effect: Distribution of the product and resampling methods. Multivariate Behavioral Research, 39, 99128.CrossRefGoogle ScholarPubMed
McMurtry, C. M., Taddio, A., Noel, M., Antony, M. M., Chambers, C. T., Asmundson, G. J. G., … Scott, J. (2016). Exposure-based interventions for the management of needle fear across the lifespan: A clinical practice guideline and call for further research. Cognitive Behaviour Therapy, 45, 217235.CrossRefGoogle ScholarPubMed
Mithal, P., Simmons, P., Cornelissen, T., Wong, H., Riddell, R. P., McMurtry, C. M., … Taddio, A. (2018). To look or not to look during vaccination: A pilot randomized trial. Canadian Journal of Pain, 2, 18. https://doi.org/10.1080/24740527.2017.1412254.CrossRefGoogle ScholarPubMed
Muris, P., Schmidt, H., & Merckelbach, H. (1999). The structure of specific phobia symptoms among children and adolescents. Behaviour Research and Therapy, 37, 863868.CrossRefGoogle ScholarPubMed
Nikolovski, J., Koldijk, M., Weverling, G. J., Spertus, J., Turakhia, M., Saxon, L., … Navar, A.-M. (2021). Factors indicating intention to vaccinate with a COVID-19 vaccine among older U.S. Adults. PLoS ONE, 16(5), e0251963.CrossRefGoogle ScholarPubMed
Olatunji, B., Ebesutani, C., Sawchuk, C., McKay, D., Lohr, J., & Kleinknecht, R. (2012). Development and initial validation of the Medical Fear Survey–Short Version. Assessment, 19, 318336.CrossRefGoogle ScholarPubMed
Ovanessian, M., Fairbrother, N., Vorstenbosch, V., McCabe, R., Rowa, K., & Antony, M. (2019). Psychometric properties and clinical utility of the specific phobia questionnaire in an anxiety disorders sample. Journal of Psychopathology and Behavioral Assessment, 41, 3652.CrossRefGoogle Scholar
Poole, C. (2015). A history of the population attributable fraction and related measures. Annals of Epidemiology, 25, 147154.CrossRefGoogle ScholarPubMed
Rosseel, Y. (2012). Lavaan: An R package for structural equation modeling and more. Version 0.5-12 (BETA). Journal of Statistical Software, 48, 136.CrossRefGoogle Scholar
Rubin, D. B. (1976). Inference and missing data. Biometrika, 63, 581592.CrossRefGoogle Scholar
Ruiz, J., & Bell, R. (2021). Predictors of intention to vaccinate against COVID-19. Vaccine, 39, 10801086.CrossRefGoogle ScholarPubMed
Saris, W. E., Krosnick, J. A., Revilla, M., & Shae, E. M. (2010). Comparing questions with agree/disagree response options to questions with item-specific response options. Survey Research Methods, 4, 6179.Google Scholar
Sjölander, A., & Vansteelandt, S. (2011). Doubly robust estimation of attributable fractions. Biostatistics (Oxford, England), 12, 112121.CrossRefGoogle ScholarPubMed
Sörbom, D. (1989). Model modification. Psychometrika, 54(3), 371384.CrossRefGoogle Scholar
Stinson, F., Dawson, D., Chou, S. P., Smith, S., Goldstein, R., Ruan, W. J., & Grant, B. (2007). The epidemiology of DSM-IV specific phobia in the USA: Results from the National Epidemiologic Survey on Alcohol and Related Conditions. Psychological Medicine, 37, 10471059.CrossRefGoogle Scholar
Taddio, A., Ipp, M., Thivakaran, S., Jamal, A., Parikh, C., Smart, S., Sovran, J., … Katz, J. (2012). Survey of the prevalence of immunization non-compliance due to needle fears in children and adults. Vaccine, 30, 48074812.CrossRefGoogle ScholarPubMed
Van Houtem, C. M. H. H., Laine, M. L., Boomsma, D. I., Ligthart, L., van Wijk, A. J., & De Jongh, A. (2013). A review and meta-analysis of the heritability of specific phobia subtypes and corresponding fears. Journal of Anxiety Disorders, 27, 379388.CrossRefGoogle ScholarPubMed
Wardenaar, K. J., Lim, C. C. W., Al-Hamzawi, A. O., Alonso, J., Andrade, L. H., Benjet, C., … De Jonge, P. (2017). The cross-national epidemiology of specific phobia in the World Mental Health Surveys. Psychological Medicine, 47, 17441760.CrossRefGoogle ScholarPubMed
Wenzel, A., & Holt, C. (2003). Validation of the multidimensional blood/injury phobia inventory: Evidence for a unitary construct. Journal of Psychopathology and Behavioral Assessment, 25, 203211.CrossRefGoogle Scholar
Yaqub, O., Castle-Clarke, S., Sevdalis, N., & Chataway, J. (2014). Attitudes to vaccination. Social Science and Medicine, 112, 111.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Socio-demographic information (N = 15 014)

Figure 1

Table 2. Socio-demographic information for mediation analysis (N = 14 149)

Figure 2

Table 3. Mediation results using SEM

Figure 3

Fig. 1. Mediational model using SEM for demographic associations with vaccine hesitancy. FOI, fear of injection; VAC.HES, vaccine hesitancy; ETH.B, Black ethnicity; ETH.WO, White other ethnicity; ETH.OM, other mixed background ethnicity; ETH.A, Asian ethnicity; ETH.O, other ethnicity; INC.2040, income group between £20k and £40k; INC.4060, income group between £40k and £60k; INC.G60, income group greater than £60k; INC.PNTS, prefer not to say income group. Reference group for ethnicity is White. Reference group for gender is male.

Figure 4

Fig. 2. Forest plot of associations with vaccine hesitancy (all factors included in the same model). Reference group for gender is male and for ethnicity is White.

Supplementary material: File

Freeman et al. supplementary material

Tables S1-S4

Download Freeman et al. supplementary material(File)
File 22.4 KB