Prev Chronic DisPreventing Chronic Disease1545-1151Centers for Disease Control and Prevention220056303222908PCDv86_10_0296Original ResearchPeer ReviewedA Latent Class Modeling Approach to Evaluate Behavioral Risk Factors and Health-Related Quality of LifeJiangYongwenPhDCenter for Health Data and Analysis, Rhode Island Department of Health. Dr Jiang is also affiliated with the Brown University School of Medicine, Providence, Rhode Island
3 Capitol Hill, Providence, RI 02908401-222-5797Yongwen.Jiang@health.ri.gov
ZackMatthew M.MD, MPHCenters for Disease Control and Prevention, Atlanta, Georgia
1120111510201186A1372011Introduction

The Behavioral Risk Factor Surveillance System (BRFSS) monitors multiple health indicators related to 4 domains: risky behaviors, health conditions, health care access, and use of preventive services. When evaluating the effect of these indicators on health-related quality of life (HRQOL), conventional analytical methods focus only on individual risks and thus are not ideally suited for analyzing complex relationships among many health indicators. The objectives of this study were to 1) summarize and group multiple related health indicators within a health domain by using latent class modeling and 2) analyze how 24 health indicators in 4 health domains were associated with 2 HRQOL outcomes to identify Rhode Island adult populations at highest risk for poor HRQOL.

Methods

The 2008 Rhode Island BRFSS, a population-based, random-digit–dialed telephone survey, collected responses from 4,786 adults aged 18 years or older. We used latent class modeling to assign 24 health indicators to high-, intermediate-, and low-risk groups within 4 domains. The effects of all risks on HRQOL were then assessed with logistic regression modeling.

Results

The latent class model with 3 classes fitted the 4 domains best. Respondents with more health conditions and limited health care access were more likely to have frequent physical distress. Those with more health conditions, risky behaviors, and limited health care access were more likely to have frequent mental distress. Use of preventive health services did not affect risk for frequent physical or mental distress.

Conclusion

The latent class modeling approach can be applied to identifying high-risk subpopulations in Rhode Island for which interventions may have the most substantial effect on HRQOL.

Introduction

The 2 overarching goals of Healthy People 2010 are "to increase quality and years of healthy life" and "to eliminate health disparities" (1). Tracking population health-related quality of life (HRQOL) can monitor progress toward national health objectives (2) and helps identify health disparities. Besides asking HRQOL questions, population-based health surveys such as the Behavioral Risk Factor Surveillance System (BRFSS) monitor multiple health domains, including risky behaviors, health conditions, health care access, and use of preventive services, which are all associated with HRQOL. Each of these domains has many indicators, most of which are highly correlated. In epidemiologic studies, adjusting for multiple, often correlated, confounders is necessary to draw valid inferences about associations. Because traditional analytical methods such as linear regression (3) or logistic regression (4-12) focus only on individual risks, they are not ideally suited to analyzing complex relationships involving multiple health domains. Indicators measure only special aspects of each domain; therefore, traditional methods must include several indicators to assess these domains comprehensively. However, including highly correlated items can lead to collinearity and problems with estimating and evaluating the association of these domains with outcomes such as HRQOL. This study shows how latent class modeling can summarize indicators within a health domain into latent classes and how logistic regression can then associate these latent classes with HRQOL. Specifically, we applied these methods to Rhode Island 2008 BRFSS data to analyze associations among 24 indicators in 4 health domains with 2 HRQOL outcomes to identify Rhode Island populations at highest risk for poor HRQOL.

MethodsData source

BRFSS is a telephone survey administered in all 50 states and 4 US territories with funding and specifications from the Centers for Disease Control and Prevention (CDC). BRFSS monitors behavioral health risks that contribute to the leading causes of disease and death among adults aged 18 years or older. It also monitors access to health care and certain health conditions. Rhode Island has participated in BRFSS since 1984. From the 2008 Rhode Island BRFSS, a population-based, random-digit–dialed telephone survey of 4,786 respondents, we selected 24 dichotomous indicators of 4 health domains — risky behaviors, health conditions, health care access, and use of preventive services (13) — and 2 HRQOL outcomes. Although CDC's institutional review board (IRB) approved the survey administration, IRB approval was unnecessary for these secondary analyses.

Indicators

We assessed 5 risky behaviors: current smoking (smokes cigarettes daily or some days), binge drinking (≥5 drinks for men, ≥4 drinks for women on 1 occasion during the past 30 d), heavy drinking (≥3 drinks/d for men, ≥2 drinks/d for women), drinking and driving in the past 30 days, and not always wearing a seatbelt (Figure 1).

Indicators, health domains, outcomes, and covariates for latent class modeling approach to analyze 2008 data from the Rhode Island Behavioral Risk Factor Surveillance System. Abbreviations: BMI, body mass index; PSA, prostate-specific antigen; Pap, Papanicolaou.

FigureThis diagram shows how the indicators, domains, outcomes, and sociodemographic characteristics are related. Twenty-four manifest dichotomous health indicators (represented by rectangles) are grouped into 4 domains (circles) that are related to 2 outcomes after adjustment for sociodemographic characteristics. The first 5 indicators measure the domain risky behaviors; the next 8 indicators, health conditions; the next 5 indicators, health care access; and the remaining 6 indicators, use of preventive services. The outcomes, frequent physical distress and frequent mental distress, are then regressed on these 4 domains and on age group, sex, race/ethnicity, education level, and employment status.

Respondents reported ever having had a diagnosis of the following 8 health conditions: asthma (ever being told by a health professional that they had asthma and they have asthma now), diabetes, obesity (body mass index ≥30 kg/m2, calculated from self-reported height and weight), disability (limited in any activities or use special equipment), loss of 6 or more teeth because of decay or gum disease, angina/coronary heart disease, heart attack/myocardial infarction, or stroke.

Respondents answered 5 questions about access to health care: being uninsured (no health care coverage), not seeing a physician because of cost (despite needing to, during the past year), having no regular health care provider, no routine checkup in the past year, or no dental visit in the past year.

We assessed use of preventive services with 6 questions: never had a pneumonia vaccination, had no influenza vaccination in the past 12 months, no prostate-specific antigen (PSA) test in the past 2 years (men ≥40 y), never had a sigmoidoscopy or colonoscopy (men and women ≥50 y), had no mammogram in past 2 years (women ≥40 y), or had no Papanicolaou (Pap) test in the past 3 years (women with an intact cervix).

Outcomes: 2 HRQOL items

The 2 HRQOL items asked respondents about their recent physically unhealthy days ("Now thinking about your physical health, which includes physical illness and injury, for how many days during the past 30 days was your physical health not good?") and mentally unhealthy days ("Now thinking about your mental health, which includes stress, depression, and problems with emotions, for how many days during the past 30 days was your mental health not good?"). These variables were dichotomized at 14 or more days (frequent physical or mental distress) and 0 to 13 days (no physical or mental distress). Others have used this cutpoint to identify frequent distress (5,10).

We included as potential confounders in the multivariable analyses age group (18-44 y, 45-64 y, ≥65 y), sex, race/ethnicity (non-Hispanic white, Hispanic, non-Hispanic other [which includes African Americans because they are rare in the adult Rhode Island population]), education level, and employment status (2,12,14). We excluded income from the models because 14% of respondents lacked income information; we instead used education level and employment status as indicators of socioeconomic status because both highly correlate with income.

Analyses

The latent class model collapses many discrete indicators into a few meaningful latent classes, categories that refer to levels of a domain, and estimates the probability that a particular person belongs to a specific latent class. The unit of analysis is the response pattern (15,16). In our study, the 5 indicators of risk behaviors have a total of 32 (25) possible response patterns (eg, current smoker vs not current smoker); the 8 indicators of health conditions, 256 (28) possible response patterns; the 5 indicators of access to health care, 32 (25) possible response patterns; and the 6 indicators of use of preventive services, 16 (24) (men) and 32 (25) (women) possible response patterns (17-19). For example, for the 5 indicators of risky behavior, latent class modeling can summarize the 32 possible response patterns into an interpretable number of latent classes (eg, 2-4).

The latent class model provides the prevalence of each latent class (marginal probabilities), and the class-specific response probabilities of each indicator (conditional probabilities) (14,20-23). Formulae and parameters for latent classes appear in the Appendix.

To determine the most parsimonious model, we sequentially fitted models from 1 to more latent classes (16,18,24), compared successive models using the Bayesian information criterion (BIC), and chose the model with the smallest BIC values (25,26) (Table 1).

To examine further the influence of high-risk classes on health, we performed 2 logistic regressions to assess the association of all these latent classes with the 2 outcomes, physical health and mental health. These models estimate odds ratios with 95% confidence intervals for frequent physical and mental distress of the higher risk classes relative to the lower risk classes, adjusting for age, sex, race/ethnicity, education, and employment status.

We used Latent Gold version 4.5 (Statistical Innovations, Inc, Belmont, Massachusetts) for latent class modeling because this software can accommodate the complex sampling design of BRFSS and can also handle missing data. We performed all other analyses with SAS version 9.1 (SAS Institute, Inc, Cary, North Carolina) and accounted for the complex sampling design of BRFSS. Significance was set at P < .05.

Results

The latent class model with 3 latent classes provided the best fit for all 4 health domains (Table 1). Nine percent of Rhode Island's adult population belonged to the highest risk class for the 5 risky behaviors (Figure 2): 29% or more in this class reported currently smoking cigarettes, binge drinking, heavy drinking, drinking and driving, or not wearing seatbelts. Approximately 31% of the population made up an intermediate class who do not always wear seatbelts, are current smokers, and drink alcohol in binges but do not drink alcohol heavily or drink and drive. In the remaining 60% of the population, less than 8% reported any of these risky behaviors.

Probabilities of risky behaviors by latent class, Rhode Island Behavioral Risk Factor Surveillance System, 2008. Behaviors refer to the past 30 days. Binge drinking defined as ≥5 drinks for men or ≥4 drinks for women on 1 occasion. Heavy drinking defined as ≥3 drinks/d for men or ≥2 drinks/d for women.

Figure
BehaviorClass 3: High Risk (9.1% of Population), %Class 2: Intermediate Risk (31.0% of Population), %Class 1: Low Risk (59.9% of Population), %
Does not always wear seatbelt28.640.17.2
Current smoking33.931.07.8
Binge drinking99.517.45.1
Heavy drinking53.82.52.1
Drinking and driving29.53.40.4

Six percent of the adult population belonged to the highest risk class for the 8 health conditions (Figure 3): at least 15% in this class reported having asthma, diabetes, obesity, a disability, 6 or more teeth lost, angina or coronary heart disease, a myocardial infarction, or a stroke. In 18% of the population, 20% or more reported having asthma, diabetes, obesity, a disability, or 6 or more teeth lost, but less than 6% reported the cardiovascular conditions — angina, heart attack, or stroke. In the remaining three-fourths of the population, 10% or less reported any of these conditions except obesity (16%).

Probabilities of health conditions by latent class, Rhode Island Behavioral Risk Factor Surveillance System, 2008. Health conditions defined as having ever been told by a health professional that they had the condition. Obesity defined as body mass index ≥30 kg/m2.

Figure
ConditionClass 3: Multiple Conditions (5.7% of Population), %Class 2: Intermediate Conditions (18.1% of Population), %Class 1: Few Conditions (76.2% of Population), %
Asthma15.620.08.0
Diabetes33.020.42.4
Obesity41.039.616.5
Disability58.755.59.7
≥6 teeth removed54.437.54.9
Angina63.32.10.8
Heart attack64.30.00.4
Stroke17.25.20.5

Ten percent of Rhode Island's adult population belonged to the highest risk class for difficulty accessing the 5 kinds of health care services (Figure 4): at least 45% of these reported having no health care coverage, not being able to see a physician because of cost, having no regular health care provider, having no health checkup in the past year, and having no dental visits in the past year. In 7% of the adult population, at least 78% said they had visited the dentist, had health coverage, and could afford to see a doctor, but at least half had not had a routine checkup in the past year or a regular provider. Less than 16% of the remaining 83% of the adult population reported difficulty accessing any of these kinds of health care.

Probabilities of health care access by latent class, Rhode Island Behavioral Risk Factor Surveillance System, 2008.

Figure
Access IndicatorClass 3: Limited Access (10.2% of Population), %Class 2: Intermediate Access (7.3% of Population), %Class 1: Has Access (82.5% of Population), %
No regular provider44.297.20.0
No routine checkup in past year61.355.812.9
No dental visit in past year58.322.116.2
No health coverage68.319.22.9
Does not see doctor due to cost73.70.14.8

Twenty-one percent of Rhode Island's adult population belonged to the highest risk class for not using the 6 kinds of preventive health care services (Figure 5): at least 75% of these had never had a pneumonia vaccination and had not had an influenza vaccination in the past year; at least 77% of those aged 50 or older had never had a sigmoidoscopy or a colonoscopy; 99% of the men aged 40 or older had not had a PSA test in the past 2 years; and at least 48% of women had not had a Pap test in the past 3 years or a mammogram in the past 2 years for women aged 40 or older. In 52% of the population, more than 70% had not had a pneumonia vaccination or an influenza vaccination in the past year but at least 65% had used 1 or more of the other preventive services. In the remaining 27% of the population, at least 70% used all of the different kinds of preventive health care services.

Probabilities of use of preventive services by latent class, Rhode Island Behavioral Risk Factor Surveillance System, 2008. The prostate-specific antigen (PSA) test question was asked only of men. The mammogram and Papanicolaou (Pap) test questions were asked only of women.

Figure
ServiceClass 3: Lower Use (21.1% of Population), %Class 2: Intermediate Use (51.6% of Population), %Class 1: Higher Use (27.3% of Population), %
No pneumonia vaccine ever85.394.028.5
No flu vaccine in past year74.970.920.5
No PSA test in past 2 y, ages ≥40 y98.833.311.5
No sigmoidoscopy or colonoscopy ever, ages ≥50 y76.633.312.4
No mammogram in past 2 y, ages ≥40 y81.43.26.0
No Pap test in past 3 y48.40.113.6

Approximately 11% of Rhode Island adults reported frequent physical or mental distress. Frequent physical distress was more common in older adults, those with a high school education or less, the retired, the unemployed, and those unable to work (Table 2). Frequent mental distress was more common in younger adults, women, those with some college or less education, homemakers, students, the unemployed, and those unable to work.

Respondents with more health conditions and limited health care access were more likely to have frequent physical distress (Table 3). Those with more risky behaviors, health conditions, and limited health care access were more likely to have frequent mental distress. Use of preventive health services did not affect risk for frequent physical or mental distress.

Discussion

Our results are consistent with previous findings (8-11,27) that indicators of individual health conditions, such as obesity (8,9) and asthma (27), or individual health care access, such as no health care coverage (11), are associated with frequent physical and mental distress. Individual risky behaviors such as binge drinking (10,11), are only related to frequent mental distress. Our methods allow us to go beyond these results by including many correlated items as indicators for 4 health domains in assessing 2 outcomes. Other authors have either enumerated patterns of adherence to specific health behaviors (28) or summed the number of risk behaviors or health conditions to provide a summary score for analysis (29). The latter approach has been criticized because it weighs risk behaviors or health conditions equally regardless of their severity (30). Latent class modeling is qualitatively more informative than a summary score. For instance, the highest-risk latent class of risky behaviors for Rhode Island adults differs most from lower-risk classes by their larger percentages who binge drink, drink heavily, or drink and drive rather than who do not always use seatbelts or currently smoke.

Combining latent class modeling with logistic regression provides 3 advantages in analyzing data from surveys like BRFSS. First, by grouping multiple correlated health indicators into a few health domains, latent class modeling allows these indicators to become associated with the outcomes through these domains and minimizes multicollinearity. To avoid this problem, many studies (8-11,27) evaluate only some indicators that partially reflect only a specific aspect of each health domain, making their analyses subject to confounding from other excluded indicators. Second, combining latent class models and logistic regression models allows several indicators to contribute to each health domain, thus improving the reliability of these domains by reducing the variation in the health domain due to any single indicator. Finally, once the number of latent classes is selected, latent class models do not require specifying arbitrary cutoffs to distinguish these classes.

These methods analyze men and women together rather than separating them. When analyzed by sex, the patterns of the 4 health domains for men and women were similar. Although some of the indicators in the "preventive services" domain are sex-specific, the Latent Gold software identifies these indicators as structurally missing and treats them as missing at random and totally dependent on the covariates (31). Although men generally use preventive services less often than women, the latent class patterns for the 3 preventive service indicators that are not sex-specific are similar among men and women. Finally, only about 1,690 men and 3,096 women responded to the 2008 Rhode Island BRFSS, making it difficult to carry out sex-specific analyses with so many indicators while avoiding artifacts.

We summarized 24 indicators in 4 health domains to predict 2 outcomes. Some may object that we did not include some important indicators, that some health behaviors were misclassified as health conditions, or that the health domains provide little explanatory power beyond a summary score of risky behaviors and health conditions. BRFSS does not ask the same questions every year. By choosing 2008, we included questions on mammography screening, Pap tests, PSA tests, colorectal screening, oral health, and drinking and driving in the past 30 days. We also excluded questions asked only in odd years about physical activity, eating fruits and vegetables, the burden of arthritis, and gastrointestinal diseases. We excluded questions about HIV and AIDS because very few respondents reported either condition. We also excluded 2 questions on emotional support and life satisfaction that appeared to form a fifth latent health domain because this domain complicated interpretation while not adding much to the 4 existing domains. We included obesity as an indicator for a health condition rather than for a health behavior because we view obesity as an outcome of other behaviors (ie, inadequate diet and physical inactivity) (28).

Besides not being able to include some health indicators, this study has several limitations. BRFSS is cross-sectional, making it difficult to determine whether the risky health behaviors, health conditions, and access to health care increased risk for frequent physical and mental distress or whether this distress increased risk for these behaviors, conditions, and health care access. Because BRFSS respondents reported information, not being able to validate this information and possible recall bias may have affected the observed associations. Latent class modeling has its own limitations that require interpretation. Each class should have a substantial number of observations because small classes may be artifactual outliers. The latent classes themselves should be interpretable in light of current theory and other empirical evidence about the domain they represent.

Despite these limitations, other surveys with multiple health indicators such as the Youth Risk Behavior Survey and the Pregnancy Risk Assessment Monitoring System may benefit from the combined use of latent class modeling and logistic regression in analyzing outcomes of interest. The methods developed here can be applied to other outcomes (eg, diabetes, asthma) and subpopulations (eg, racial/ethnic minorities, people with disabilities).

We conducted separate latent class models for the 4 domains because incorporating all the risk indicators into 1 latent class model was too complicated to allow model convergence. Exploring larger models would be worthwhile. Moreover, when the probability of a pattern of indicators does not differ substantially for different latent classes within a domain, these patterns do not distinguish these latent classes; thus, their labels may not accurately describe all their patterns. For example, those who currently smoke but do not indulge in of any of the other 4 risky behaviors (4% of the sample) are classified into the low-risk class because more in this class have this pattern (42%) than in the intermediate-risk class (35%) and high-risk class (23%).

This study has identified high-risk subpopulations in Rhode Island where health service and health policy interventions may improve HRQOL — among residents with high probabilities of risky behaviors, multiple health conditions, or limited health care access.

We thank our colleagues in the Rhode Island Department of Health for their comments and suggestions as the manuscript evolved. We also thank William W. Thompson, PhD, and Kurt Greenlund, PhD, Centers for Disease Control and Prevention, for reviewing and commenting on the final draft. This work and the 2008 Rhode Island Behavioral BRFSS were supported in part by the Chronic Disease Prevention and Health Promotion Programs Cooperative Agreement no. 5U58DP122791-05.

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.

Suggested citation for this article: Jiang Y, Zack MM. A latent class modeling approach to evaluate behavioral risk factors and health-related quality of life. Prev Chronic Dis 2011;8(6) http://www.cdc.gov/pcd/issues/2011/nov/10_0296.htm. Accessed [date].

2000 Healthy People 2010: understanding and improving health 2nd edition Washington (DC) US Department of Health and Human Services ZackMMMoriartyDGStroupDFFordESMokdadAH 119 5 2004 493 505 Worsening trends in adult health-related quality of life and self-rated health — United States, 1993-2001 Public Health Rep 15313113 GoinsRTSpencerSMKrummelDA 96 6 2003 552 557 Effect of obesity on health-related quality of life among Appalachian elderly South Med J 12938781 BarrettDHBoehmerTKBootheVLFlandersWDBarrettDH 168 11 2003 941 947 Health-related quality of life of US military personnel: a population-based study Mil Med 14680052 LiCFordESMokdadAHBalluzLSBrownDWGilesWH 11 4 2008 689 699 Clustering of cardiovascular disease risk factors and health-related quality of life among US adults Value Health 18194400 DominickKLAhernFMGoldCHHellerDA 2004 2 5 Health-related quality of life among older adults with arthritis Health Qual Life Outcomes 14720300 BrownDWBalluzLSHeathGWMoriartyDGFordESGilesWHMokdadAH 37 5 2003 520 528 Associations between recommended levels of physical activity and health-related quality of life. Findings from the 2001 Behavioral Risk Factor Surveillance System (BRFSS) survey Prev Med 14572437 FordESMoriartyDGZackMMMokdadAHChapmanDP 9 1 2001 21 31 Self-reported body mass index and health-related quality of life: findings from the Behavioral Risk Factor Surveillance System Obes Res 11346664 HassanMKJoshiAVMadhavanSSAmonkarMM 27 10 2003 1227 1232 Obesity and health-related quality of life: a cross-sectional analysis of the US population Int J Obes Relat Metab Disord 14513071 OkoroCABrewerRDNaimiTSMoriartyDGGilesWHMokdadAH 26 3 2004 230 233 Binge drinking and health-related quality of life: do popular perceptions match reality? Am J Prev Med 15026103 StrineTWBalluzLChapmanDPMoriartyDGOwensMMokdadAH 26 3 2004 213 216 Risk behaviors and healthcare coverage among adults by frequent mental distress status, 2001 Am J Prev Med 15026100 JiangYHesserJE 2006 4 14 Associations between health-related quality of life and demographics and health risks. Results from Rhode Island's 2002 behavioral risk factor survey Health Qual Life Outcomes 16515690 DuBardCAGizliceZ 98 11 2008 2021 2028 Language spoken and differences in health status, access to care, and receipt of preventive services among US Hispanics Am J Public Health 18799780 JiangYHesserJE 2008 6 49 Patterns of health-related quality of life and patterns associated with health risks among Rhode Island adults Health Qual Life Outcomes 18620582 FlahertyBP 2002 68 Suppl 1 S7 20 Assessing reliability of categorical substance use measures with latent class analysis Drug Alcohol Depend 12324172 StraussSMRindskopfDMAstone-TwerellJMDes JarlaisDCHaganH 83 1 2006 15 24 Using latent class analysis to identify patterns of hepatitis C service provision in drug-free treatment programs in the US Drug Alcohol Depend 16289523 McCutcheonAL 1987 Latent class analysis Newbury Park (CA) Sage Publications ReboussinBASongEYShresthaALohmanKKWolfsonM 83 3 2006 199 209 A latent class analysis of underage problem drinking: evidence from a community sample of 16-20 year olds Drug Alcohol Depend 16359829 van LangNDFerdinandRFOrmelJVerhulstFC 44 6 2006 849 860 Latent class analysis of anxiety and depressive symptoms of the Youth Self-Report in a general population sample of young adolescents Behav Res Ther 16122697 DaytonCMMacreadyGB 83 401 1988 173 178 Concomitant-variable latent-class models J Am Stat Assoc XueQLBandeen-RocheK 58 1 2002 110 120 Combining complete multivariate outcomes with incomplete covariate information: a latent class approach Biometrics 11890305 Bandeen-RocheKMigliorettiDLZegerSLRathouzPJ 92 440 1997 1375 1386 Latent variable regression for multiple discrete outcomes J Am Stat Assoc LarsenK 60 1 2004 85 92 Joint analysis of time-to-event and multiple binary indicators of latent classes Biometrics 15032777 MongaNRehmJFischerBBrissetteSBruneauJEl-GuebalyN 88 1 2007 1 8 Using latent class analysis (LCA) to analyze patterns of drug use in a population of illegal opioid users Drug Alcohol Depend 17049753 CarlsonRGWangJFalckRSSiegalHA 79 2 2005 167 179 Drug use practices among MDMA/ecstasy users in Ohio: a latent class analysis Drug Alcohol Depend 16002026 GuoJWallMAmemiyaY 7 1 2006 145 163 Latent class regression on latent factors Biostatistics 16079163 FordESManninoDMHomaDMGwynnCReddSCMoriartyDGMokdadAH 123 1 2003 119 127 Self-reported asthma and health-related quality of life: findings from the Behavioral Risk Factor Surveillance System Chest 12527612 BerriganDDoddKTroianoRPKrebs-SmithSMBarbashRB 36 5 2003 615 623 Patterns of health behavior in US adults Prev Med 12689807 LiCFordESZhaoGMokdadAH 49 2-3 2009 253 259 Associations of health risk factors and chronic illnesses with life dissatisfaction among US adults: the Behavioral Risk Factor Surveillance System, 2006 Prev Med 19501613 SackerAWigginsRDClarkePBartleyM 25 3 2003 215 222 Making sense of symptom checklists: a latent class approach to the first 9 years of the British Household Panel Survey J Public Health Med 14575196 LittleRJARubinDB 1987 New York (NY) John Wiley Statistical analysis with missing data Appendix: Parameters for Latent Classes

For a latent class model with 3 categorical indicators (c1, c2 and c3) and 5 indicators (Y1, Y2, Y3, Y4, Y5), the posterior probability of class c  under pattern of y1, y2, y3, y4 and y5 is

where

= P(y1|c)P(y2|c)P(y3|c)P(y4|c)P(y5|c)

Therefore we need to calculate P(c) and P(yk|c) first.

P(c) can be calculated from parameters of Model for Clusters as

, c=1,2,3

When Yk is a binary indicator, with 0 for not present and 1 present, P(yk | c) can be calculated from parameters for Models for Indicators and Overall Intercepts, ie,

, k=1,2,3,4,5; c=1,2,3.

1. Parameters for risk behaviors

Model for Clusters

c1c2c3
γ1γ2γ3
0−0.6589−1.8837

Models for Indicators, βkc, k = 1 to 5; c = 1 to 3

Risk behaviorCluster
c1c2c3
Does not always wear seatbelt02.15801.6457
Current smoking01.67391.8032
Binge drinking01.35638.1758
Heavy drinking00.14683.9802
Drinking and driving04.37926.8685

Overall Intercepts

Does not always wear seatbeltα1−2.5588
Current smokingα2−2.4719
Binge drinkingα3−2.9152
Heavy drinkingα4−3.8289
Drinking and drivingα5−7.7384

2. Parameters for health conditions

Model for Clusters

c1c2c3
γ1γ2γ3
0−1.4362−2.5869

Models for Indicators, βkc, k = 1 to 8; c = 1 to 3

Health indicatorsCluster
c1c2c3
Asthma01.06160.7619
Diabetes02.33192.9862
Obesity01.20441.2616
Disability02.45852.5869
≥6 Teeth removed02.44633.1331
Angina00.95345.3237
Heart attack0−2.51086.0319
Stroke02.38553.7092

Overall Intercepts

Asthmaα1−2.4489
Diabetesα2−3.6929
Obesityα3−1.6249
Disabilityα4−2.2358
≥6 Teeth removedα5−2.9565
Anginaα6−4.7777
Heart attackα7−5.4432
Strokeα8−5.2806

3. Parameters for health care access

Model for Clusters

c1c2c3
γ1γ2γ3
0−2.0890−2.4236

Models for Indicators, βkc, k = 1 to 8; c = 1 to 3

Health care acccessCluster
c1c2c3
No regular provider07.555311.3359
No routine checkup in past year02.37222.1454
No dental visit in past year01.97870.3852
No health coverage04.28712.0795
Unable to see a doctor because of cost04.0284−4.3380

Overall Intercepts

No regular providerα1−7.7890
No routine checkup in past yearα2−1.9135
No dental visit in past yearα3−1.6427
No health coverageα4−3.5175
Unable to see a doctor because of costα5−2.9970

4. Parameters for preventive service

Model for Clusters

c1c2c3
γ1γ2γ3
0−0.6369−0.8973

Models for Indicators, βkc, k = 1 to 8; c = 1 to 3

Preventive serviceCluster
c1c2c3
No pneumonia vaccine ever0−3.6641−0.9908
No flu vaccine in past year0−2.24580.1996
No PSA test in past 2 years, ages ≥400−1.34155.1220
No sigmoidoscopy or colonoscopy ever, ages ≥500−1.26121.8784
No mammogram in past 2 years, ages ≥4000.66064.8812
No Pap test in past 3 years05.45417.2361

Overall Intercepts

No pneumonia vaccine everα12.7465
No flu vaccine in past yearα20.8924
No PSA test in past 2 years, ages ≥40α3−0.6957
No sigmoidoscopy or colonoscopy ever, ages ≥50α4−0.6932
No mammogram in past 2 years, ages ≥40α5−3.4063
No Pap test in past 3 yearsα6−7.3014

Fit of Latent Class Models Using Bayesian Information Criterion (BIC)a for Health Domains of the Rhode Island Behavioral Risk Factor Surveillance System, 2008

No. of Classes in ModelHealth Domain, BIC Value

Risky BehaviorHealth ConditionHealth Care AccessPreventive Service Use
117019.1723700.5019911.7718451.38
216142.5122412.9818306.6417836.85
316119.4022318.9018171.8717661.99
416130.1122354.4118180.7317698.54

Models were sequentially fitted from 1 to more latent classes (16,18,24) and compared to successive models using BIC. The smallest BIC values (25,26) indicate the most parsimonious model.

Prevalence of Frequent Physical and Mental Distress, by Demographic Characteristics, Rhode Island Adults, Behavioral Risk Factor Surveillance System, 2008

CharacteristicNo. of Respondents, % (n = 4,786)aFrequent Physical Distress, % (95% CI) (n = 635)Frequent Mental Distress, % (95% CI) (n = 500)
Age, y
18-441,331 (48)7.2 (5.4-9.0)12.2 (9.7-14.7)
45-641,966 (33)12.0 (10.3-13.6)10.4 (8.9-11.9)
≥651,446 (18)18.1 (15.9-20.3)7.0 (5.5-8.4)
Sex
Male1,690 (48)9.4 (7.7-11.1)8.8 (6.6-11.0)
Female3,096 (52)12.0 (10.6-13.5)12.2 (10.6-13.8)
Race/ethnicity
Non-Hispanic white4,148 (86)10.5 (9.3-11.8)10.4 (9.0-11.9)
Hispanic320 (9)13.1 (9.1-17.1)11.6 (7.1-16.0)
Non-Hispanic other251 (5)11.0 (7.0-15.1)11.4 (6.5-16.3)
Education
High school or less1,826 (37)15.8 (13.6-18.1)13.9 (11.3-16.4)
Some college1,115 (25)9.9 (7.7-12.0)12.0 (9.2-14.7)
College graduate1,829 (38)6.6 (5.3-7.9)6.5 (4.8-8.3)
Employment status
Employed2,554 (61)5.9 (4.6-7.1)8.3 (6.9-9.8)
Homemaker/student376 (11)7.2 (3.5-10.8)14.6 (8.1-21.1)
Retired1,253 (16)17.3 (14.9-19.6)6.1 (4.6-7.7)
Unemployed271 (7)17.3 (11.6-23.0)19.2 (11.9-26.5)
Unable to work318 (4)56.2 (49.4-63.0)35.2 (28.8-41.7)

Abbreviation: CI, confidence interval.

Some categories do not add to the total because of missing responses. Percentages are weighted.

Odds Ratios of Frequent Physical and Mental Distress, by Health Domain, Rhode Island Adults, Behavioral Risk Factor Surveillance System, 2008a

Health DomainFrequent Physical DistressFrequent Mental Distress

OR (95% CI)P ValueOR (95% CI)P Value
Risk behavior
Low risk1 [Reference]NA1 [Reference]NA
Intermediate risk1.21 (0.88-1.65).242.32 (1.68-3.21)<.001
High risk1.16 (0.47-2.89).753.85 (2.13-6.98)<.001
Health condition
Few conditions1 [Reference]NA1 [Reference]NA
Intermediate conditions2.40 (1.69-3.41)<.0012.52 (1.74-3.64)<.001
Multiple conditions3.14 (2.13-4.62)<.0012.34 (1.41-3.87)<.001
Health care access
Has access1 [Reference]NA1 [Reference]NA
Intermediate access0.99 (0.54-1.83).980.86 (0.44-1.68).66
Limited access2.10 (1.14-3.86).022.40 (1.38-4.17).002
Use of preventive service
High use1 [Reference]NA1 [Reference]NA
Intermediate use1.05 (0.74-1.49).771.24 (0.84-1.83).27
Low use1.01 (0.67-1.53).950.98 (0.62-1.54).92

Abbreviations: OR, odds ratio; CI, confidence interval; NA, not applicable.

Analyses were adjusted for age, sex, race/ethnicity, education, employment status, and all other latent variables in the table. Frequent physical or mental distress was defined as reporting 14 or more physically or mentally unhealthy days in the last 30 days.