Open access peer-reviewed chapter

# On Statistical Assessments of Racial/Ethnic Inequalities in Cigarette Purchase Price among Daily Smokers in the United States: Non-Hispanic Whites Pay Least

Written By

Julia N. Soulakova and Trung Ha

Submitted: 07 July 2020 Reviewed: 14 July 2020 Published: 21 August 2020

DOI: 10.5772/intechopen.93380

From the Edited Volume

## Recent Advances in Numerical Simulations

Edited by Francisco Bulnes and Jan Peter Hessling

Chapter metrics overview

View Full Metrics

## Abstract

We discuss statistical methods suitable for comparing multiple populations versus one reference population and consider two common problems: (1) detecting all significant mean differences and (2) demonstrating that all mean differences are significant. Discussed methods include the Bonferroni approach (both problems), Min test (problem 2), and Strassburger-Bretz-Hochberg (SBH) confidence interval for estimating the smallest mean difference (problem 2). They illustrate the methods using the pooled 2010–2015 Tobacco Use Supplement to the Current Population Survey (TUS-CPS) data on the cigarette purchase price (per pack) reported by adult daily smokers (n = 34,728). The goal was to show that among seven considered racial/ethnic groups of daily smokers, non-Hispanic (NH) Whites paid least for cigarettes (on average). We used the design-based multiple linear regression to derive the estimates and raw p-values. The Min test supported the study goal. Likewise, SBH lower 95% confidence interval bound was $0.08, indicating that the other racial/ethnic groups of daily smokers paid at least eight cents more for a pack of cigarettes (on average) than did non-Hispanic Whites. However, Bonferroni method (that was originally proposed for problem 1) failed to support the study goal. The study highlights the importance of choosing the right statistical method for a given problem. ### Keywords • balanced repeated replications • complex survey • multiple comparisons • statistical multiple-testing problems ## 1. Introduction In this chapter, we discuss statistical methods for comparing multiple populations relative to one population (termed “reference”). These types of multiple comparisons commonly arise in behavioral science, for example, when multiple racial/ethnic groups are compared to non-Hispanic (NH) White smokers in terms of tobacco-use-related behaviors [1, 2, 3, 4]. When the statistical parameter of interest is the mean difference, the most common study goal is one of the following two goals. Goal 1 is to detect all significant mean differences among the considered populations (versus the reference population), that is, to draw an individual conclusion regarding significance of each mean difference. Goal 2 is to demonstrate that all mean differences among the considered ones are significant. Note that if one assessed Goal 1 and concluded that each mean difference was significant then s/he has (indirectly) assessed Goal 2 as well. Other more intricate study goals, such as the ones arising in pharmaceutical statistics which involve a hierarchical structure among the primary and secondary end points, were addressed elsewhere and are outside of the scope of this chapter [5, 6, 7, 8, 9, 10, 11]. We discuss how Goals 1 and 2 can be assessed in a study of racial and ethnic disparities, where Hispanic (H) population and five non-Hispanic populations such as American Indian/Alaska Native (AIAN), Asian (ASIAN), Black/African American (BAA), Hawaiian/Pacific Islander (HPI), and Multiracial (MULT), are compared to non-Hispanic White (W) population in terms of the mean differences. μAIANμW,μASIANμW,μBAAμW,μHμW,μHPIμW,andμMULTμW,E1 where μAIAN, μASIAN,μBAA, μH, μHPI, μMULT, and μW denote, respectively, the mean responses for AIAN, ASIAN, BAA, H, HPI, MULT, and W populations. Furthermore, suppose that each positive mean difference in Eq. (1) corresponds to a significant result, for example, the first difference being positive implies that the mean response among AIANs is greater than the mean response among Ws. Then the null and alternative hypotheses corresponding to theithdifference, where i denotes AIAN, ASIAN, BAA, H, HPI, and MULT, can be stated as H0i:μiμW0andHai:μiμW>0.E2 Finally, let pi(i=AIAN,ASIAN,BAA,H,HPI,MULT) denote a p-value corresponding to testing H0i versus Hai. As a result, we have six pairs of hypotheses and six p-values. Suppose the overall error rate for assessing Goal 1 (Goal 2) is fixed at α-level. Then to assess Goal 1 (to detect all significant mean differences), we should first rescale each p-value pi, for example, via Bonferroni, Holm, or Hochberg approaches [6, 12, 13, 14]. This rescaling is essential to control the overall error rate at the nominal α-level. For example, in our case with six null hypotheses, the p-values rescaled via Bonferroni method are given as 6pi (i.e., we multiply each original p-value by six). Second, we compare each rescaled p-value with α. If piα, then we reject H0i and conclude that the ith difference is significant (positive); if pi>α, then we accept H0i and conclude that the ith difference is not significant. As a result, we draw an individual conclusion regarding significance of each mean difference. Alternatively to the above hypothesis testing, one could construct the lower 1001α/6% confidence intervals for the mean differences in Eq. (1) and use the lower bounds to differentiate between the significant and insignificant mean differences; this approach was discussed elsewhere [15]. To assess Goal 2 (to demonstrate that all differences are significant), one can use the Min test that is an intersection-union test [16, 17, 18, 19, 20, 21]. The p-value for the Min test, denoted by p, is simply the largest p-value among the six individual p-values: p=maxpAIANpASIANpBAApHpHPIpMULT.E3 If pα, then we reject all H0is and conclude that all mean differences are significant (positive); if p>α, then we fail to reject all H0is and conclude that there is at least one insignificant mean difference. Note that we cannot comment on the significance of an individual mean difference, because we tested whether all mean differences are significant (i.e., whether the smallest mean difference is significant). Nonetheless, the Min test is more suitable for assessing Goal 2 than Bonferroni approach or another approach proposed for assessing Goal 1. A statistical method is usually proposed for a specific problem and thus, the methods should be used accordingly: the union-intersection hypothesis (Goal 1) should be tested via Bonferroni or another union-intersection test, while the intersection-union hypothesis (Goal 2) should be tested via the Min or another intersection-union test [12, 22]. Alternatively to the Min test, we can use the Strassburger-Bretz-Hochberg (SBH) confidence interval approach as follows [23, 24]. First, we compute the lower 1001α% confidence intervals for the mean differences in Eq. (1). Let these bounds be denoted as LAIAN, LASIAN, LBAA, LH, LHPI, and LMULT. Second, let L denote the smallest bound among these bounds, that is, L=minLAIANLASIANLBAALHLHPILMULT.E4 Then the SBH lower 1001α% confidence interval for the smallest mean difference is given by L,+.If L>0, we conclude that all mean differences are significant, and if L0, then we conclude that there is at least one insignificant mean difference. We note that one needs to identify the appropriate statistical method to compute the individual p-values and confidence bounds. The choice depends on the study design, probability distributions, and other statistical considerations. The Min test and the SBH interval were discussed for parallel and factorial designs, where sample mean responses followed normal distributions with known variances or unknown (common) variance, as well as Binomial and several other distributions [20, 21, 23, 24, 25, 26, 27]. In addition, one needs to decide whether the analyses should adjust for explanatory factors, for example, sociodemographic characteristics [28, 29, 30]. Such adjustments may help reduce the effect of confounding factors and therefore, improve estimation [31, 32]. For example, Golden et al. examined how much smokers pay for a pack of cigarettes, on average, in the United States using data from the 2010–2011 Tobacco Use Supplement to the Current Population Survey (TUS-CPS) [1]. Among several design-based multiple linear regression models for the mean purchase price per pack (PPP) used in the study, one model adjusted for smokers’ sociodemographic and smoking-related characteristics, cigarette purchase attributes, and the survey wave [1]. Despite availability and benefits of the Min test and SBH interval, these methods have not received much attention in behavioral sciences. We illustrate benefits of using these methods over Bonferroni method and simplicity of applications of these methods. We consider a study of racial and ethnic disparities in cigarette purchase prices conducted to demonstrate that W daily smokers, on average, purchase cigarettes at lower prices than do AIAN, ASIAN, BAA, H, HPI, and MULT daily smokers in the United States. This goal was motivated by results of a prior study revealing that BAA, H, and ASIAN/HPI (ASIAN and HPI combined) smokers paid higher PPP, on average, relative to W smokers, in the United States in the period from 2010 to 2011 [1]. Advertisement ## 2. Methods and results ### 2.1 Using data to derive the p-values and lower confidence interval bounds We used the pooled 2010–2011 and 2014–2015 TUS-CPS data for adult daily smokers (n = 34,728) who reported the price of the last self-purchased pack or carton of cigarettes. The reported prices were used to compute the (average) PPP. The overall cohort was representative of about 23,370,261 adult daily smokers, where 12% were 18–24 years old, 38% were 25–44 years old, and 50% were 45+ years old, and 54% were men and 47% were women. The racial/ethnic representation was as follows: 76% were W, 11% were BAA, 8% were H, 2% were MULT, 2% were ASIAN, 1% were AIAN, and less than 1% were HPI. All racial/ethnic groups were well represented in the sample: the smallest number of respondents (96) corresponded to HPI daily smokers. Additional sample characteristics have been described in a prior study of purchasing cigarettes on Indian reservations [33]. We fixed the overall error rate at α = 5% and fitted a design-based multiple linear regression (R2 ≈ 30%, F(25, 160) ≈ 257, p < 0.0001) to model the mean PPP as a function of daily smokers’ characteristics, location of the purchase (on/off Indian reservation), survey mode (phone, in-person), and survey period (2010–2011, 2014–2015). The daily smokers’ characteristics included race/ethnicity, age, sex, marital status, education, employment record, region of residency (West, South, Midwest, and Northeast), metropolitan area of residency (metro, nonmetro), and heavy smoking indicator. The analysis incorporated statistical methods recommended in the methodological guidelines for analysis of the CPS and CPS supplements [34, 35]. Specifically, because the CPS incorporates complex sampling, we estimated variance using balanced repeated replications [36]. The main and 160 replicate weights for this approach have been made available for public use by the U.S. Census Bureau [34, 35]. The analysis was performed using SAS®9.4 software [37]; the SAS®9.4 Survey Package procedures suitable for analysis of TUS-CPS have been discussed elsewhere [38]. Table 1 depicts the estimated model coefficients and their standard errors for all covariates. As is shown in Table 1, smokers’ sex and survey mode (phone, in-person) were not significant. FactorEstimated coefficientStandard errorp-Value* Intercept3.640.12* Race/ethnicity (reference group is W) AIAN versus W0.610.13* ASIAN versus W0.620.09* BAA versus W0.510.04* H versus W0.610.06* HPI versus W0.830.220.0002 MULT versus W0.220.080.0087 Age (reference group is 45+ years old) 18–24 years old versus 45+ years old0.190.04* 25–44 years old versus 45+ years old0.200.02* Sex Female versus male0.000.020.9052 Marital status (reference group is widowed/divorced/separated) Married (living with a spouse)0.020.020.4878 Never married0.220.03* Highest level of education (reference group is some college/Bachelor’s degree) Graduate degree0.110.080.1632 High school/equivalent−0.140.02* Less than high school−0.150.03* Employment status (reference group is unemployed) Employed (at work or absent) versus unemployed0.160.03* Not in labor force versus unemployed−0.150.04* Place where cigarettes were purchased (reference group is “on Indian reservation”) (reference group is “yes) No versus yes1.570.12* U.S. region of residency Midwest versus West0.090.030.0091 Northeast versus West1.750.05* South versus West−0.620.03* Metropolitan area of residency Metropolitan area versus nonmetropolitan area0.320.03* Heavy smoking indicator Heavy (20+ cigarettes per day) versus non-heavy smoker−0.200.02* Survey mode Personal interview versus phone interview−0.010.020.5205 Survey period 2010–2011 versus 2014–2015−0.380.02* ### Table 1. Design-based multiple linear regression for the mean cigarette price per pack. p-value < 0.0001. Table 1 presents the individual p-values for comparisons of racial/ethnic populations of daily smokers versus W daily smokers (based on the model): pAIAN<0.0001, pASIAN<0.0001, pBAA<0.0001, pH<0.0001, pHPI=0.0002, and pMULT=0.0087.The individual lower 95% confidence interval bounds for the mean PPP difference for each racial/ethnic population of daily smokers relative to W daily smokers were computed using the formula: Li=d^it0.95,df=160SEd^i,i=AIAN,ASIAN,BAA,H,HPI,MULT,E5 where d^i denotes the estimated mean PPP difference relative to W daily smokers (the point estimate for the ith mean difference), SEd^i is the standard error of the estimate (computed using the balanced repeated replications), and t0.05,df=160=1.6544 is the 95th percentile of the central t-distribution with 160 degrees of freedom (the number of degrees of freedom matches the number of the replicate weights) [34, 35, 36]. We note that there are alternative methods to construct the lower bounds, for example, using the standard normal distribution instead of the central t-distribution [34, 36]. Figure 1 depicts the lower bounds Lis and the estimated mean differences d^is for all racial/ethnic populations (relative to the W population). These bounds were computed using proc surveyreg procedure with lsmestimate statements (with “cl,” “e,” “upper,” and “alpha = 0.05” options) when fitting the model using SAS software. Alternatively, we could use the lsmeans statement (with “adj = bon,” “cl,” and “alpha = 0.1” options), and select the comparisons of interest out of all 21 pair-wise comparisons reported and note the lower bound of the two-sided 90% confidence interval reported in the output. ### 2.2 Demonstrating the study goal via the min test and SBH confidence interval The p-value for the Min test is p=0.0087, indicating that at 5% significance level we reject the null hypothesis in favor of the alternative. The corresponding SBH lower 95% confidence interval bound for the mean PPP difference is$0.08 (see Figure 1). Therefore, all six racial/ethnic groups of daily smokers paid, on average, higher PPP relative to W daily smokers in the United States in the periods from 2010–2011 to 2014–2015.

If instead of the Min test we used the Bonferroni approach, then the adjusted p-values would be less than 0.0006 for four comparisons (AIAN versus W, ASIAN versus W, BAA versus W, and H versus W), 0.0012 for one comparison (HPI versus W), and 0.0522 for one comparison (MULT versus W). Therefore, we would conclude that only AIAN, ASIAN, BAA, H, and HPI daily smokers pay higher PPP, on average, than do W daily smokers; and would fail to demonstrate that all six considered racial/ethnic groups of daily smokers pay higher PPP, on average, relative to W daily smokers.

## 3. Discussion

The choice of the reference group as “W daily smokers” was based on the study goal and prior studies of cigarette purchasing behaviors of smokers [1, 33]. The choice of the reference group as well as the statistical methods should always align with the study goal and should be made prior to the data analysis. Specifically, when examining racial/ethnic disparities, using “W” as the reference group could be logical in some studies but not logical in the other studies. For example, if the study goal is to show that purchasing cigarettes on Indian reservations is most prevalent among AIAN smokers, then “AIAN smokers” should be chosen as the reference group. In addition, while both Bonferroni method and the Min test are simple to use, in practice, only Bonferroni method results in individual conclusions regarding each comparison. However, Bonferroni method is less powerful than the Min test when applied to an intersection-union problem (to assess Goal 2) [6, 12].

The study indicated that W daily smokers paid significantly less for cigarettes, on average, than the other six racial/ethnic groups of daily smokers in the United States in the period from 2010–2011 to 2014–2015. The earlier reported finding (see model 6 in [1]) was that non-Hispanic White smokers, on average, paid significantly less for cigarettes than did BAA, AIAN, ASIAN/HPI (combined), and H smokers, and paid similar prices to the prices paid by “other non-Hispanic” smokers [1]. While the results might seem to disagree, the direct comparisons between these two findings are problematic, because the studies concerned different populations of smokers (daily smokers in our study, and daily and occasional smokers in the prior study) and time periods (overall 2010–2011 and 2014–2015 in our study, and 2010–2011 in the prior study). Moreover (though, the authors did not mention the method they used to adjust for multiple comparisons, if any), the authors considered the union-intersection problem that is conceptually different from the intersection-union problem addressed in our study [1].

Our study has several potential limitations. First, we considered the population of daily smokers, and thus, results should not be generalized to other populations of smokers such as occasional smokers. Indeed, daily and occasional smokers have very different cigarette purchasing behaviors, for example, daily smokers are more likely to purchase cigarettes in cartons rather than packs and travel to another state or Indian reservations to purchase cigarettes at lower prices [1, 39, 40]. Second, the analysis was based on a certain regression model where the mean PPP was modeled as a function of smokers’ characteristics, location of the purchase, survey mode, and survey period. Another model could potentially lead to a different conclusion, for example, only two out of six models indicated significantly higher mean PPP for AIAN smokers relative to W smokers [1]. Another potential limitation is a lack of a theoretical proof that the SBH interval for the smallest mean PPP difference has indeed confidence level of 1001α%. The probability coverage of the SBH confidence interval depends on the probability coverage of the individual confidence intervals for the mean differences [23]. Because we used the statistical methods outlined in the CPS methodological guidelines for constructing the individual intervals, we believe that the resulting SBH interval has the probability coverage close to 1001α% level.

Future research may target development and implementation of procedures for the Min test and SBH interval. Specifically, the software packages developed for analysis of complex survey data currently offer just a few multiple comparison methods. For example, the SAS Survey Package offers a built-in procedure for Bonferroni adjustments but lacks procedures for the multiple testing (interval estimation) such as the Min test (SBH interval). Availability of the “Min test” and “SBH interval” procedures would enable researchers to incorporate these methods directly in their analyses of complex survey data.

## 4. Conclusion

In our study, results of the Min test (and SBH interval) were different from the results of the Bonferroni method. Specifically, using the Min test (and SBH interval), we demonstrated that all six racial/ethnic groups of daily smokers paid, on average, higher PPP relative to W daily smokers in the United States in the periods from 2010–2011 to 2014–2015. However, using the Bonferroni method, we failed to demonstrate this claim. This discrepancy highlights the importance of choosing the appropriate statistical method for assessing the minimum among multiple mean differences (relative to one reference population). Availability of the “Min test” and “SBH interval” procedures in survey packages would help facilitate application of these methods in behavioral research.

## Acknowledgments

The authors are thankful to Richard Pack, B.S., for providing editing comments.

## Conflict of interest

There is no conflict of interest.

## Funding

Research reported in this publication was supported by the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Number R01MD009718. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

## References

1. 1. Golden SD, Kong AY, Ribisl KM. Racial and ethnic differences in what smokers report paying for their cigarettes. Nicotine & Tobacco Research. 2016;18(7):1649-1655
2. 2. Soulakova JN, Li J, Crockett LJ. Race/ethnicity and intention to quit cigarette smoking. Preventive Medicine Reports. 2017;5:160-165. Available from: http://linkinghub.elsevier.com/retrieve/pii/S2211335516301632
3. 3. Cokkinides VE, Halpern MT, Barbeau EM, Ward E, Thun MJ. Racial and ethnic disparities in smoking-cessation interventions: Analysis of the 2005 National Health Interview Survey. American Journal of Preventive Medicine. 2008;34(5):404-412
4. 4. Tran S-TT, Rosenberg KD, Carlson NE. Racial/ethnic disparities in the receipt of smoking cessation interventions during prenatal care. Maternal and Child Health Journal. 2010;14(6):901-909. Available from: http://link.springer.com/10.1007/s10995-009-0522-x [Accessed: 10 December 2019]
5. 5. Dmitrienko A, Soulakova JN, Millen BA. Three methods for constructing parallel gatekeeping procedures in clinical trials. Journal of Biopharmaceutical Statistics. 2011;21(4):768-786
6. 6. Hochberg Y, Tamhane A. Multiple Comparison Procedures. New York: Wiley; 1987
7. 7. Dmitrienko A, Offen W, Wang O, Xiao D. Gatekeeping procedures in dose–response clinical trials based on the Dunnett test. Pharmaceutical Statistics. 2006;5(1):19-28. Available from: http://doi.wiley.com/10.1002/pst.190 [Accessed: 22 February 2018]
8. 8. Dmitrienko A, Offen WW, Westfall PH. Gatekeeping strategies for clinical trials that do not require all primary effects to be significant. Statistics in Medicine. 2003;22(15):2387-2400. Available from: http://doi.wiley.com/10.1002/sim.1526 [Accessed: 22 February 2018]
9. 9. Soulakova JN. Comparison of gatekeeping and other testing methods for identifying superior drug combinations in bifactorial designs with isotonic parameters. Journal of Biopharmaceutical Statistics. 2011;21(4):635-649
10. 10. Chen X, Luo X, Capizzi T. The application of enhanced parallel gatekeeping strategies. Statistics in Medicine. 2005;24(9):1385-1397. Available from: http://doi.wiley.com/10.1002/sim.2005 [Accessed: 22 February 2018]
11. 11. Sugitani T, Morikawa T. Gatekeeping strategies and graphical approaches in clinical trials with hierarchically structured study objectives: A review. Japanese Journal of Biometrics. 2017;38(1):41-78
12. 12. Hsu JC. Multiple Comparisons: Theory and Methods. New York: CRC Press; 1996
13. 13. Soulakova JN. On identifying effective and superior drug combinations via Holm’s procedure based on the Min tests. Journal of Biopharmaceutical Statistics. 2009;19(2):280-291
14. 14. Holm S. A simple sequentially rejective multiple test procedure on JSTOR. Scandinavian Journal of Statistics. 1979;6(2):65-70. Available from: https://www.jstor.org/stable/4615733?seq=1#page_scan_tab_contents [Accessed: 25 July 2018]
15. 15. Johnson RA, Wichern DW. Applied Multivariate Statistical Analysis. Upper Saddle River, New Jersey: Prentice Hall; 2002. pp. 232-234
16. 16. Lehmann EL. Testing multiparameter hypotheses. Annals of Mathematical Statistics. 1952;23:541-552
17. 17. Berger RL. Multiparameter hypothesis testing and acceptance sampling. Technometrics. 1982;24(4):295-300. Available from: http://www.jstor.org/stable/1267823?origin=crossref [Accessed: 20 February 2018]
18. 18. Berger RL, Hsu JC. Bioequivalence trials, intersection-union tests and equivalence confidence sets. Statistical Science. 1996;11(4):283-319. Available from: https://projecteuclid.org/download/pdf_1/euclid.ss/1032280304 [Accessed: 20 February 2018]
19. 19. Berger RL. Likelihood Ratio Tests and Intersection-Union Tests. Raleigh, North Carolina: Institute of Statistics Mimeo Series Number 2288; 1996. Available from: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.17.4041&rep=rep1&type=pdf [Accessed: 20 February 2018]
20. 20. Laska EM, Meisner MJ. Testing whether an identified treatment is best. Biometrics. 1989;45(4):1139. Available from: http://www.jstor.org/stable/2531766?origin=crossref [Accessed: 20 February 2018]
21. 21. Laska EM, Tang D-I, Meisner MJ. Testing hypotheses about an identified treatment when there are multiple endpoints. Journal of the American Statistical Association. 1992;87(419):825. Available from: http://www.jstor.org/stable/2290221?origin=crossref [Accessed: 20 February 2018]
22. 22. Casella G, Berger RL. Statistical Inference. 2nd ed. Pacific Grove, CA: Duxbury Press; 2001
23. 23. Strassburger K, Bretz F, Hochberg Y. Compatible Confidence Intervals for Intersection Union Tests Involving Two Hypotheses. Institute of Mathematical Statistics; 2004. Vol. 47. pp. 129-142. Available from: http://projecteuclid.org/euclid.lnms/1196285631 [Accessed: 20 February 2018]
24. 24. Soulakova JN. Generalized confidence intervals compatible with the Min test for simultaneous comparisons of one subpopulation to several other subpopulations. Communications in Statistics - Theory and Methods. 2017;46(19):9441-9449
25. 25. Hung HMJ. Global tests for combination drug studies in factorial trials. Statistics in Medicine. 1996;15(3):233-247 Available from: http://doi.wiley.com/10.1002/%28SICI%291097-0258%2819960215%2915%3A3%3C233%3A%3AAID-SIM167%3E3.0.CO%3B2-H [Accessed: 22 February 2018]
26. 26. Hung HMJ, Chi GYH, Lipicky RJ. Testing for the existence of a desirable dose combination. Biometrics. 1993;49(1):85. Available from: http://www.jstor.org/stable/2532604?origin=crossref [Accessed: 22 February 2018]
27. 27. Westfall P, Ho S-Y, Prillaman B. Properties of multiple intersection-union tests for multiple endpoints in combination therapy trials. Journal of Biopharmaceutical Statistics. 2001;11(3):125-138
28. 28. Cartmell KB, Miner C, Carpenter MJ, Vitoc CS, Biggers S, Onicescu G, et al. Secondhand smoke exposure in young people and parental rules against smoking at home and in the car. Public Health Reports. 2011;126(4):575-582
29. 29. Johnson T, Mott J. The reliability of self-reported age of onset of tobacco, alcohol and illicit drug use. Addiction. 2001;96(8):1187-1198.Available from: http://onlinelibrary.wiley.com/doi/10.1046/j.1360-0443.2001.968118711.x/full [Accessed: 15 December 2016]
30. 30. Soulakova JN, Crockett LJ. Unassisted quitting and smoking cessation methods used in the United States: Analyses of 2010-2011 tobacco use supplement to the current population survey data. Nicotine & Tobacco Research. 2017;20(1):30-39
31. 31. Jiang J. Linear and Generalized Linear Mixed Models and their Applications. New York, USA: Springer Science & Business Media, LLC New York; 2007
32. 32. Charles EM, Shayle RS, John MN. Generalized, Linear, and Mixed Models. 2nd ed. Hoboken, New Jersey, USA: John Wiley and Sons Inc.; 2008
33. 33. Soulakova JN, Pack R, Ha T. Patterns and correlates of purchasing cigarettes on Indian reservations among daily smokers in the United States. Drug and Alcohol Dependence. 2018;192:88-93. Available from: http://www.ncbi.nlm.nih.gov/pubmed/30243144 [Accessed: 12 October 2018]
34. 34. U.S. Bureau of Labor Statistics, U.S. Census Bureau. Design and Methodology: Current Population Survey, Technical Paper 66 [Internet]. 2006. Available from: https://www.census.gov/prod/2006pubs/tp-66.pdf
35. 35. U.S. Department of Commerce, U.S. Census Bureau. National Cancer Institute and Food and Drug Administration co-sponsored Tobacco Use Supplement to the Current Population Survey. 2014-15. Technical Documentation [Internet]. 2016. Available from: https://www.census.gov/programs-surveys/cps/technical-documentation/complete.html
36. 36. Wolter KM. Introduction to Variance Estimation. New York, USA: Springer; 2007. p. 447
37. 37. SAS Institute Inc. SAS® 9.4 Product Documentation. Cary, NC, USA: SAS Institute Inc; 2013
38. 38. Ha T, Soulakova JN. Statistical Analyses of Public Health Surveys Using SAS® Survey Package. In: SESUG Paper 189 [Internet]. 2017. Available from: http://analytics.ncsu.edu/sesug/2017/SESUG2017_Paper-189_Final_PDF.pdf
39. 39. Cornelius ME, Driezen P, Hyland A, Fong GT, Chaloupka FJ, Cummings KM. Trends in cigarette pricing and purchasing patterns in a sample of US smokers: Findings from the ITC US surveys (2002-2011). Tobacco Control. 2015;24(Supplement 3):iii4-iii10. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24917617 [Accessed: 22 February 2018]
40. 40. Pesko MF, Kruger J, Hyland A. Cigarette price-minimization strategies by U.S. smokers. American Journal of Public Health. 2012;102(9):19-21

Written By

Julia N. Soulakova and Trung Ha

Submitted: 07 July 2020 Reviewed: 14 July 2020 Published: 21 August 2020