We first start with a base case scenario which represents the status quo. In the base case, the rate of switching to SFP use is representative of a population of AS with their existing perceptions of nicotine harm. We compare the base case to four nicotine perception scenarios in which the rate of switching to SFP use is adjusted based on stratified transition probabilities associated with the responses to the question “Do you believe nicotine is the chemical that causes most of the cancer caused by smoking cigarettes?” from Wave 4 data of the Population Assessment of Tobacco and Health (PATH) Study (R04_AC9120, Four-item scale from “Definitely not” to “Definitely yes”). Each nicotine perception scenario modeled the population impact expected to occur if all AS held the same level of perception of nicotine harm as those in each response group from “Definitely not” to “Definitely yes”. We modeled nicotine perception scenarios corresponding to each response rather than dichotomizing them into correct versus incorrect (9, 11, 12) to highlight the trend in population health impact associated with each response on this perception scale. The public health impact of nicotine perceptions is estimated as the difference in all-cause mortality and cigarette prevalence among AS, between the base case and each nicotine perception scenario.
In this research, we use a validated Agent-Based Population Model (ABM) as previously described (15–17) which projects future cigarette prevalence and all-cause mortality beginning in the year 2000. The ABM begins by initializing a hypothetical population of 2.81 million agents (1/100th of the year 2000 U.S. population) that is representative of the U.S. population for age, sex and tobacco use status. The model is iterated 100 times to simulate the entire U.S. population. The initial population mirrors the age and sex distribution data from the year 2000 U.S. Census. Each agent in the initial population is assigned one of three tobacco use statuses representative of people who have never smoked (NT), people who currently smoke (CS) and people who formerly smoked (FS). Tobacco use status is assigned by sex and ages ≥ 18 using information from the National Health Interview Survey (NHIS) Sample Adult Questionnaire data for the year 2000. In our analysis of NHIS data, CS are defined as those who reported having smoked at least 100 cigarettes in their lifetime and are currently smoking every day or some days now. FS were defined as those who reported having smoked at least 100 cigarettes in their lifetime and currently not smoking at the time of interview. NT are defined as those who have never smoked or have not reached 100 cigarettes in their lifetime. These definitions are commonly used by the Centers for Disease Control and Prevention (CDC) to estimate tobacco use prevalence for the U.S. adult population (18). Since NHIS does not provide tobacco use information for ages < 18 years, tobacco use status assigned to the younger U.S. population, ages 10–17, by sex were estimated from the 2000 National Youth Tobacco Survey (NYTS). The NYTS is a nationally representative survey of middle and high school students focused exclusively on patterns of tobacco use. In analysis of the NYTS, we define youth who have never smoked or have not smoked 100 cigarettes in their lifetime as NT; those who reported having smoked 100 cigarettes in their lifetime and having smoked in past 30 days as CS; and those who reported having smoked 100 cigarettes in their lifetime but not having smoked in past 30 days as FS.
Each agent in the initial population is assigned tobacco use history which is updated over the 100-year simulation timeframe from year 2001 to 2100. Agents who initialize as CS or FS are assigned with their associated years of smoking and/or years stopped smoking and the age(s) at which the agent initiated and/or stopped smoking. Age and sex-specific probabilities from U.S. birth cohort smoking history data developed by Jeon et al. were used to assign when a CS or FS in the model’s starting population initiated or stopped smoking (19). The age and sex specific cigarette smoking initiation and cessation probabilities were generated by researchers who used NHIS surveys administered from 1964 to 2015 to estimate birth cohort smoking histories (20) and details of the methodological approach and the resulting data are available on the Cancer Intervention and Surveillance Modelling Network website (CISNET, https://resources.cisnet.cancer.gov/projects/ - Publication Support and Modeling Resources).
Once the initial population is generated, the following algorithms are executed in 1-year time intervals throughout the simulation time frame:
Mortality Sub-model:
A mortality sub-model is used to estimate the survival probability of each agent based on their age, sex and current or former tobacco use history. The mortality sub-model was developed using data from a Kaiser Permanente (KP) Medical Care Program Cohort study (21), which included number of deaths, person- years, smoking status, age, sex, years smoked, and years since quitting smoking. The KP Study data were adjusted using the Human Mortality Database (HMD) to be representative of the U.S. population in the Year 2000. Mortality rates throughout the simulation time frame are further adjusted to account for expected age-specific changes in mortality over time using the methodology described by Carter et al., (22).
Transition sub-model:
At each time interval within a simulation, agents are provided with an option to change or maintain their current tobacco use status. These decisions are governed by the agent’s defined current status, age and sex specific transition probabilities.
Population Update:
The age of agents who survive at each time interval is increased by 1-year increments. New agents are added to the population each year to account for birth and net immigration based on U.S. Census projections (23). Agents who enter the population via immigration were assigned tobacco use status and history similar to that used for the initial population.
ABM scenarios are projected through the year 2100. We quantify the public health impact of varying perceptions of nicotine harm as the difference in adult prevalence and cumulative all-cause deaths for ages 35–85 between the base case and nicotine perception scenarios. The 35–85 age range is used since most smoking mortalities are expected to occur within this range (1).
Base Case Scenario Inputs
The base case scenario models AS prevalence and all-cause mortality where cigarette smoking continues to be the predominant tobacco use behavior. Transition probabilities for initiation from NT to CS and cessation of CS to FS by sex and age are taken from CISNET, where the smoking history of U.S. birth cohorts is reported using NHIS data (19). CISNET transition probabilities are available by age (0–99), sex and year through 2015. CISNET initiation probabilities by sex and age are updated each year in the model from 2001 through 2015, at which point they are held constant for the remaining simulation timeframe. Cessation probabilities in CISNET are based on at least two years of successful smoking cessation. This is an important consideration since minimal relapse has been reported (24) after a period of two years which minimizes the relapse transition between FS and CS in our model. The model updates the CISNET smoking cessation probabilities from CISNET each year between 2001–2013 with the corresponding CISNET estimates. Beginning in the year 2014, the FS state is split to differentiate individuals that successfully quit smoking and do not use SFP from those who successfully quit smoking but are using SFP now. The SFP use state represents the proportion of FS who completely switched to smoke-free product use. Figure 1 provides a diagram of the base case tobacco use states and transitions.
To differentiate FS from SFP use, beginning in 2014 transitions from CS to FS and CS to SFP use are estimated using data from Wave 1 (W1) 2013/2014 to Wave 5 (W5) 2018/2019 of the PATH study which is funded by the FDA Center for Tobacco Products (CTP) and administered by the National Institute on Drug Abuse (NIDA). The PATH study was designed to generate longitudinal epidemiologic data on tobacco-use behavior and health in the U.S. population (25). As of 2022, five waves of PATH data are publicly available for analyses (26). In our analysis of PATH, we identified AS (n = 6,349) in W1 as those who reported smoking at least 100 cigarettes in their lifetime and currently smoke every day or some days. FS were identified as W1 AS who no longer smoked in subsequent waves and did not use ENDS (defined in PATH to include products such as e-cigarettes, e-hookahs, e-cigars, e-pipes, personal vaporizers, vape pens, and hookah pens) or smokeless tobacco (ST, defined in PATH to include moist loose snus, snuff, dip and spit or chewing tobacco or snus every day or some days). SFP use was identified as W1 AS who no longer smoked but also indicated every day or someday use of ENDS, ST or snus in subsequent waves. We calculated successful quitting and switching as people who smoked in W1 who transitioned to either FS or SFP use in W2 (2015/2016) and maintained their FS or SFP use status through W5 (an approximate 3-year of follow up). W1 AS who indicated they had not used ENDS, ST or snus every day or some days in W2 but had used in any subsequent wave, W3 through W5, were excluded from our analysis. This was done to ensure the transition rates reflected long-term sustained behavior. Estimates from five waves of PATH were used to obtain successful quitting/switching probabilities, which follows a similar methodology as CISNET successful smoking cessation calculation. We observed higher successful cigarette cessation probabilities based on the PATH study data, which aligns with observations from existing research that indicates a recent increasing trend in cigarette cessation (27). We applied PATH cessation probabilities to account for recent changes in smoking prevalence observed between 2014 and 2018 since they are not captured in the CISNET cessation rates.
Based on our analysis of PATH W1 to W5, 3.92% male and 4.35% female AS transitioned to FS in W2 and remained FS through W5. Adults who smoked in W1 transitioned to SFP use in W2 and remained SFP users through W5 at a rate of 1.32% and 0.44%, male and female respectively.
The all-cause mortality probability assigned to CS that switch to SFP use is based on the excess relative risk (ERR) of SFP use compared to cigarette smoking. We applied an ERR of 10%, which reflects a reasonable aggregate estimate for the various SFPs. An ERR of 10% is slightly higher than the estimate of 9% for current ST product use relative to smoking (28). Applying a multi-criteria decision analysis (MCDA) model developed by an international expert panel convened by the Independent Scientific Committee on Drugs (29), the authors assigned the relative importance of different types of harm related to the use of nicotine-containing products compared to cigarette smoking. The panel estimated a relative risk of 5% for snus and 4% for ENDS products relative to smoking. An ERR value of 10% represents a conservative choice for SFP use in the model given that this tobacco use state is representative of exclusive or other combinations of ST, snus and ENDS use. Additionally, we do not allow for transition to former SFP use, therefore agents entering the SFP use state carry their former CS risk in addition to the ERR for current SFP use throughout the remaining simulation timeframe.
Nicotine Perception Scenarios
Nicotine perception scenarios were simulated in which all AS transition to SFP use at rates associated with different level of perception of nicotine harm. PATH W4 (2016–2018) data were used to assess AS perceptions of nicotine harm based on their responses to the question “Do you believe nicotine is the chemical that causes most of the cancer caused by smoking cigarettes?” (R04_AC9120, Four-item scale from “Definitely not” to “Definitely yes”). We choose W4 as the baseline wave to analyze the levels of nicotine perception because, in addition to also being a true cross-sectional wave due to sample replenishment, it is more recent compared to W1 (2013–2014) and may reflect the current state of nicotine perception among AS. Table 1 shows the demographics characteristics of W4 AS overall and by perception of nicotine harm. We note that striking differences were observed among the different age and race/ethnicity subgroups in response to the nicotine perceptions question. While the nicotine misperceptions of harm were high across all age subgroups, they trended to be higher among the 45 + age subgroup. More than 70% of AS in the 45 + age subgroup responded “Definitely yes” or “Probably yes” regarding harm from nicotine, compared to ~ 50–55% for those in the lower age groups. Additionally, a higher proportion (~ 76%) of AS in the Non-Hispanic Black and Hispanic subgroups indicated higher perceptions of nicotine harm compared to the Non-Hispanic Whites (~ 48%). The data in Table 1 are important as the observed differences may result in different rates of switching to SFP use between demographic subgroups. Modeling can theoretically be used to investigate population health impacts by demographic characteristics however, the sample sizes needed to provide robust input estimates have limited most models to include only age and sex variables. We focused our simulations on the overall group of W4 established AS (18+) due to sample size and model limitations (i.e., current model is not capable of utilizing population inputs by race or economic status) however it is important to consider the data in Table 1 when interpreting the results.
Table 1
AS demographics by nicotine perception held
Characteristic
|
Total
|
Perception of nicotine harm: “Do you believe nicotine is the chemical that causes most of the cancer caused by smoking cigarettes?”
|
Definitely yes
|
Probably yes
|
Probably not
|
Definitely not
|
|
7007
|
1338
|
2943
|
2090
|
515
|
Sex
|
Female
|
46.38%
|
19.01%
|
45.99%
|
28.95%
|
6.04%
|
Male
|
53.62%
|
18.06%
|
43.19%
|
31.64%
|
7.11%
|
Age
|
18–24
|
8.15%
|
13.79%
|
36.36%
|
36.85%
|
12.99%
|
25–44
|
45.47%
|
16.84%
|
39.97%
|
35.52%
|
7.67%
|
45+
|
46.38%
|
21.02%
|
50.36%
|
24.18%
|
4.44%
|
Race/Ethnicity
|
Non-Hispanic White
|
68.78%
|
13.47%
|
44.42%
|
35.02%
|
7.09%
|
Non-Hispanic Black
|
13.08%
|
33.94%
|
43.08%
|
17.57%
|
5.41%
|
Non-Hispanic Others
|
5.42%
|
19.06%
|
43.02%
|
31.21%
|
6.7%
|
Hispanic
|
12.71%
|
29.63%
|
46.89%
|
18.26%
|
5.21%
|
Education
|
LT College
|
57.39%
|
22.45%
|
46.92%
|
25.19%
|
5.44%
|
Some College
|
31.70%
|
14.08%
|
41.74%
|
36.58%
|
7.59%
|
College Grad
|
10.91%
|
9.64%
|
40.32%
|
40.33%
|
9.7%
|
Smoking Behavior
|
Every Day
|
76.32%
|
18.51%
|
45.52%
|
29.84%
|
6.13%
|
Some Day
|
23.68%
|
18.58%
|
41.12%
|
32.12%
|
8.18%
|
Nicotine perception scenario inputs were developed by first estimating the overall rate of switching from CS to SFP use based on longitudinal analysis of PATH W4 (2016–2018) and W5 (2018–2019). The analysis included W4 established AS (18+) that responded to the nicotine perception question (i.e., “Definitely not”, “Probably not”, “Probably yes” or “Definitely yes”) (n = 6886) following the same definitions used to define CS and SFP users as used in the base case. The overall transition rate for W4 AS who completely switched to SFP use by W5 was estimated to be 3.94%. The overall switching rate was then stratified by W4 AS perceptions of nicotine harm. We did not estimate switching rates by age or sex to preserve sample size for each of the nicotine perception response groups. Figure 2 shows the proportions of W4 AS that comprise each response group and their associated transition rates to SFP use in W5. As shown in Fig. 2, 6.61% of W4 AS (n = 515 shown in Table 1) responded “Definitely not”. Of those AS who responded “Definitely not”, 8.39% quit smoking and switched to SFP use by W5, more than twice the overall rate. This is in contrast with 18.52% (n = 1338 in Table 1) of AS who responded “Definitely yes”, of which 2.59% switched to exclusive SFP use at W5. The error bars shown in Fig. 2 represent 95% confidence intervals for each response group switching rate.
A total of four nicotine perception scenarios were simulated corresponding to each of the four responses to the question “Do you believe nicotine is the chemical that causes most of the cancer caused by smoking cigarettes?”. The model inputs corresponding to each nicotine perception scenario were estimated by adjusting the base case sustained switching rates with the relative percent change between the W4 to W5 overall and stratified transition rates corresponding to each different level of nicotine perception. Relative percent changes are calculated as the stratified rate associated with one level of the nicotine harm perception minus the overall rate divided by the overall2 rate of CS to SFP use between W4 and W5. Relative percent changes were used to adjust the base case rates by sex to avoid over estimation of switching behavior. Switching rates between W4 to W5 would not be expected to be comparable to the sustained 1-year switching rate used in the base case since W5 is a follow up study conducted approximately two years after W4 data collection and may not reflect sustained switching. Therefore, we only applied relative percentage changes between the overall and stratified switching rates to the base case to evaluate the population impact associated with various levels of nicotine harm perception. Table 2 provides the base case transition rate inputs by sex, relative percent change adjustments and the final rates used in each of the four nicotine perception scenarios. As shown in Table 2, AS who responded “Definitely not” switched to SFP use at a rate 113% greater (more than twice) than the overall observed rate of switching while those that responded “Definitely yes” switched at a rate 34% lower than the overall. Also included in Table 2 are relative percent change adjustments corresponding to 95% confidence limits of the transition rates calculated as the ratio of the upper and lower 95% confidence limit of the transition rate associated with one level of the nicotine harm perception over the overall transition rate of CS to SFP use between W4 and W5. Sensitivity scenarios were conducted based on these 95% confidence limits of the stratified transition rates. Each scenario assumes that all AS will exhibit the switching behavior associated with the specific response groups.
Table 2
Nicotine perception scenario CS to SFP switch rates derived from adjustment of the base case rate using the relative percent change between W4 – W5 overall and each nicotine perception response group.
Nicotine Perception Scenario
|
Base case CS to SFP Switch Rates
|
Relative % change Adjustment*
|
CS to SFP Switch Rates corresponding to each Nicotine Perception Scenario**
|
Male
|
Female
|
Male
|
Female
|
Definitely not
|
1.32
|
0.44
|
113% (47%,204%)
|
2.81
|
0.94
|
Probably not
|
28% (1%, 60%)
|
1.69
|
0.56
|
Probably yes
|
-21% (-1%, -37%)
|
1.04
|
0.35
|
Definitely yes
|
-34% (-5%, -55%)
|
0.87
|
0.29
|
*Relative percent changes corresponding to 95% confidence limits of the stratified switch rates are shown. |
**CS to SFP switch rate calculated as Base Rate + Base Rate * Relative % change adjustment. |
The difference in overall prevalence and all-cause mortality between the base case and nicotine perception scenarios is used to quantify the population health impact corresponding to AS perceptions of nicotine harm.