Use of posterior predictive assessments to evaluate model fit in multilevel logistic regression

Assessing the fit of a model is an important final step in any statistical analysis, but this is not straightforward when complex discrete response models are used. Cross validation and posterior predictions have been suggested as methods to aid model criticism. In this paper a comparison is made between four methods of model predictive assessment in the context of a three level logistic regression model for clinical mastitis in dairy cattle; cross validation, a prediction using the full posterior predictive distribution and two “mixed” predictive methods that incorporate higher level random effects simulated from the underlying model distribution. Cross validation is considered a gold standard method but is computationally intensive and thus a comparison is made between posterior predictive assessments and cross validation. The analyses revealed that mixed prediction methods produced results close to cross validation whilst the full posterior predictive assessment gave predictions that were over-optimistic (closer to the observed disease rates) compared with cross validation. A mixed prediction method that simulated random effects from both higher levels was best at identifying the outlying level two (farm-year) units of interest. It is concluded that this mixed prediction method, simulating random effects from both higher levels, is straightforward and may be of value in model criticism of multilevel logistic regression, a technique commonly used for animal health data with a hierarchical structure.

Organisational justice and smoking: the Finnish Public Sector Study

Study Objective: To examine the extent to which justice of decision-making procedures and interpersonal relations is associated with smoking. Setting: Ten municipalities and 21 hospitals in Finland. Design and Participants: Cross-sectional data derived from the Finnish Public Sector Study were analysed with logistic regression analysis models with generalized estimating equations. Analyses of smoking status were based on 34 021 employees. Separate models for heavy smoking (>20 cigarettes per day) were calculated for 6295 current smokers. Main results: After adjustment for age, education, socio-economic position, marital status, job contract, and negative affectivity, smokers who reported low procedural justice were about 1.4 times more likely to smoke >20 cigarettes per day compared with their counterparts with high justice. In a similar way, after adjustments, low justice in interpersonal treatment was significantly associated with an elevated prevalence of heavy smoking (odds ratio (OR) = 1.35, 95% CI = 1.03 to 1.77 for men and OR = 1.41, 95% CI = 1.09 to 1.83 for women). Further adjustment for job strain and effort-reward imbalance had little effect on these results. There were no associations between justice components and smoking status or ex-smoking. Conclusions: The extent to which employees are treated with justice in the workplace seems to be associated with smoking intensity independently of established stressors at work.

Work stress, smoking status and smoking intensity: an observational study of 46 190 employees

Study objective: To examine the relationship between work stress, as indicated by the job strain model and the effort-reward imbalance model, and smoking. Setting: Ten municipalities and 21 hospitals in Finland. Design and Participants: Binary logistic regression models for the prevalence of smoking were related to survey responses of 37 309 female and 8881 male Finnish public sector employees aged 17-65. Separate multinomial logistic regression models were calculated for smoking intensity for 8130 smokers. In addition, binary logistic regression models for ex-smoking were fitted among 16 277 former and current smokers. In all analyses, adjustments were made for age, basic education, occupational status, type of employment and marital status. Main results: Respondents with high effort-reward imbalance or lower rewards were more likely to be smokers. Among smokers, an increased likelihood of higher intensity of smoking was associated with higher job strain and higher effort-reward imbalance and their components such as low job control and low rewards. Smoking intensity was also higher in active jobs in women, in passive jobs and among employees with low effort expenditure. Among former and current smokers, high job strain, high effort-reward imbalance and high job demands were associated with a higher likelihood of being a current smoker. Lower effort was associated with a higher likelihood of ex-smoking. Conclusions: This evidence suggests an association between work stress and smoking and implies that smoking cessation programs may benefit from the taking into account the modification of stressful features of work environment. Key words: effort-reward imbalance; job strain; smoking. Abbreviations: OR, odds ratio; CI, confidence interval; SES, socioeconomic status

Effort-reward imbalance and sedentary lifestyle: an observational study in a large occupational cohort

Objectives: To investigate the association between effort-reward imbalance (ERI) at work and sedentary lifestyle. Methods: Cross-sectional data from the ongoing Finnish Public Sector Study related to 30 433 women and 7718 men aged 17-64 were used (n = 35 918 after exclusion of participants with missing values in covariates). From the responses to a questionnaire, an aggregated mean score for ERI in a work unit was assigned to each participant. The outcome was sedentary lifestyle defined as <2.00 metabolic equivalent task (MET) hours/day. Logistic regression with generalized estimating equations was used as an analysis method to include both individual and work unit level predictors in the models. Adjustments were made for age, marital status, occupational status, job contract, smoking, and heavy drinking. Results: Twenty five percent of women and 27% of men had a sedentary lifestyle. High individual level ERI was associated with a higher likelihood of sedentary lifestyle both among women (odds ratio (OR) = 1.08, 95% CI 1.01 to 1.16) and men (OR = 1.17, 95% CI 1.02 to 1.33). These associations were not explained by relevant confounders and they were also independent of work unit level job strain measured as a ratio of job demands and control. Conclusions: A mismatch between high occupational effort spent and low reward received in turn seems to be associated with an elevated risk of sedentary lifestyle, although this association is relatively weak.

Effort-reward imbalance at work and the co-occurrence of lifestyle risk factors: cross-sectional survey in a sample of 36,127 public sector employees

Background In occupational life, a mismatch between high expenditure of effort and receiving few rewards may promote the co-occurrence of lifestyle risk factors, however, there is insufficient evidence to support or refute this hypothesis. The aim of this study is to examine the extent to which the dimensions of the Effort-Reward Imbalance (ERI) model – effort, rewards and ERI – are associated with the co-occurrence of lifestyle risk factors. Methods Based on data from the Finnish Public Sector Study, cross-sectional analyses were performed for 28,894 women and 7233 men. ERI was conceptualized as a ratio of effort and rewards. To control for individual differences in response styles, such as a personal disposition to answer negatively to questionnaires, occupational and organizational -level ecological ERI scores were constructed in addition to individual-level ERI scores. Risk factors included current smoking, heavy drinking, body mass index ≥25 kg/m2, and physical inactivity. Multinomial logistic regression models were used to estimate the likelihood of having one risk factor, two risk factors, and three or four risk factors. The associations between ERI and single risk factors were explored using binary logistic regression models. Results After adjustment for age, socioeconomic position, marital status, and type of job contract, women and men with high ecological ERI were 40% more likely to have simultaneously ≥3 lifestyle risk factors (vs. 0 risk factors) compared with their counterparts with low ERI. When examined separately, both low ecological effort and low ecological rewards were also associated with an elevated prevalence of risk factor co-occurrence. The results obtained with the individual-level scores were in the same direction. The associations of ecological ERI with single risk factors were generally less marked than the associations with the co-occurrence of risk factors. Conclusion This study suggests that a high ratio of occupational efforts relative to rewards may be associated with an elevated risk of having multiple lifestyle risk factors. However, an unexpected association between low effort and a higher likelihood of risk factor co-occurrence as well as the absence of data on overcommitment (and thereby a lack of full test of the ERI model) warrant caution in regard to the extent to which the entire ERI model is supported by our evidence.

Relationship between hyperacute blood pressure and outcome after ischemic stroke: data from the VISTA Collaboration

Background and Purpose—High blood pressure (BP) is associated independently with poor outcome after acute ischemic stroke, although in most analyses “baseline” BP was measured 24 hours or more postictus, and not during the hyperacute period. Methods—Analyses included 1722 patients in hyperacute trials (recruitment 8 hours) from the Virtual Stroke International Stroke Trial Archive (VISTA) Collaboration. Data on BP at enrolment and after 1, 2, 16, 24, 48, and 72 hours, neurological impairment at 7 days (NIHSS), and functional outcome at 90 days (modified Rankin scale) were assessed using logistic regression models, adjusted for confounding variables; results are for 10-mm Hg change in BP. Results—Mean time to enrolment was 3.7 hours (range 1.0 to 7.9). High systolic BP (SBP) was significantly associated with increased neurological impairment (odds ratio, OR 1.06, 95% confidence interval, 95% CI 1.01 to 1.12), and poor functional outcome; odds ratios for both increased with later BP measurements made at up to 24 hours poststroke. Smaller (versus larger) declines in SBP over the first 24 hours were significantly associated with poor NIHSS scores (OR 1.16, 95% CI 1.05 to 1.27) and functional outcome (OR 1.23, 95% CI 1.13 to 1.34). A large variability in SBP was also associated with poor functional outcome. Conclusions—High SBP and large variability in SBP in the hyperacute stages of ischemic stroke are associated with increased neurological impairment and poor functional outcome, as are small falls in SBP over the first 24 hours.