Chiropractic & Manual Therapies 2023 (Aug 25); 31: 32 ~ FULL TEXT
Maja Stupar • Pierre Côté • Linda J. Carroll •
Robert J. Brison
Eleanor Boyle • Heather M. Shearer • J. David Cassidy
Institute for Disability and Rehabilitation Research,
Ontario Tech University,
Objective: Few clinical prediction models are available to clinicians to predict the recovery of patients with post-collision neck pain and associated disorders. We aimed to develop evidence-based clinical prediction models to predict (1) self-reported recovery and (2) insurance claim closure from neck pain and associated disorders (NAD) caused or aggravated by a traffic collision.
Methods: The selection of potential predictors was informed by a systematic review of the literature. We used Cox regression to build models in an incident cohort of Saskatchewan adults (n = 4,923). The models were internally validated using bootstrapping and replicated in participants from a randomized controlled trial conducted in Ontario (n = 340). We used C-statistics to describe predictive ability.
Results: Participants from both cohorts (Saskatchewan and Ontario) were similar at baseline. Our prediction model for self-reported recovery included prior traffic-related neck injury claim, expectation of recovery, age, percentage of body in pain, disability, neck pain intensity and headache intensity (C = 0.643; 95% CI 0.634–0.653). The prediction model for claim closure included prior traffic-related neck injury claim, expectation of recovery, age, percentage of body in pain, disability, neck pain intensity, headache intensity and depressive symptoms (C = 0.637; 95% CI 0.629–0.648).
Conclusions: We developed prediction models for the recovery and claim closure of NAD caused or aggravated by a traffic collision. Future research needs to focus on improving the predictive ability of the models.
Keywords: Clinical prediction model; Disability; Health recovery; Neck pain; Rehabilitation; Whiplash.
From the FULL TEXT Article:
Predicting prognosis of neck pain and associated disorders (NAD) related to a traffic collision is challenging for clinicians, insurers and policy makers. This is due, in part, to the lack of high-quality evidence about the prognosis of NAD.  To date, the best evidence on the prognosis of NAD suggests that half of those with NAD secondary to traffic collisions will improve within three months and recover within six months.  Moreover, prognostic factors associated with recovery include post-crash psychological factors such as expectation of poor recovery, pain-related depression, anxiety, fear, frustration or anger and poor coping.  Improving the ability of clinicians to predict the outcomes of patients with post-collision neck pain is important because it is estimated that 50% of persons with this condition continue to experience symptoms one-year post-collision. [1, 2] Therefore, our ability to effectively predict the course of NAD is limited.
Clinical prediction models are tools developed to assist clinicians predict clinical outcomes using a combination of patient, clinical and other variables.  Using valid clinical prediction models improves prognosis prediction compared to clinical judgement alone.  Several models aimed at predicting recovery from NAD secondary to traffic collisions have been developed. [4–15] However, few of these models have been internally  or externally validated [5–7, 15, 16] and only one study developed a clinical prediction rule.  Therefore, there is a need to develop validated clinical prediction models for NAD secondary to traffic collisions.
Our objective was to develop evidence-based clinical prediction models to predict recovery from NAD secondary to traffic collision.
We used the following steps to meet our objective:
We developed prediction models in a population-based inception-cohort of Saskatchewan adults injured in traffic collisions.
We internally validated the models using bootstrapping methods.
We repeated the models in a sample of adults from Ontario enrolled in a randomized controlled trial.
We used the TRIPOD Checklist for Prediction Model Development to report our study (Appendix A).
Selection of potential prognostic Factors/Independent variables
Best-evidence synthesis methods to identify prognostic factors
We conducted a systematic review of the literature to identify potential prognostic factors to be included in our development model. 
In the systematic review, the prognostic factors were classified according to three levels of evidence:
(a) Phase I are exploratory studies reporting unadjusted associations;
(b) Phase II are exploratory studies reporting associations measured in multivariable models; and
(c) Phase III are confirmatory studies that test the independence of associations. [1, 17]
Therefore, evidence of association from Phase III studies can be used with more confidence than evidence from Phase I and Phase II studies.
Summary of the systematic review on prognostic factors
Prognostic factors identified from Phase III studies
Based on confirmatory evidence, the following factors were deemed to be independent predictors of delayed recovery or health improvement: prior NAD secondary to a traffic collision, poor expectation of recovery, passive coping and greater levels of depression, anxiety, fear, frustration and/or anger related to the post-collision pain. 
Prognostic factors identified from phase I and II studies
Based on preliminary evidence, older age may be associated with poorer disability recovery. In addition, poor or delayed recovery may be associated with: prior generalized pain, kinesiophobia, cervical radiculopathy, post-crash cold pain threshold, pain intensity and disability, number and severity of symptoms, symptoms of acute stress disorder and post-traumatic stress disorder, and early post-crash onset of depressive symptoms. 
Other prognostic factors
We also considered factors that were not identified as predictors in the systematic review but that clinicians may consider as associated with recovery based on clinical experience (i.e. sociodemographic and clinical factors) [List of variables in Appendix B].
Development and Internal Validation: the Saskatchewan dataset
We developed the model using data from a population-based, inception cohort of traffic injury from Saskatchewan.  Participants with acute NAD were identified by answering “yes” to the question, “Did the accident cause neck or shoulder pain?”  We restricted our analysis to participants who provided baseline information within 21 days of the traffic collision. Eligible for the cohort were all Saskatchewan residents who received treatment for a road traffic injury from a regulated health professional (physician, physical therapist, chiropractor or massage therapist) or made an insurance claim for treatment of their injuries to Saskatchewan Government Insurance (SGI) between December 1, 1997 and November 30, 1999. Participants were interviewed approximately 6, 12, 26, 36 and 52 weeks post-injury.
Replicating: repeating the Models in the Ontario dataset
We replicated our models using data from a randomized clinical trial that investigated the effectiveness of physiotherapy and physician-based care for the management of NAD.  The trial was conducted across five primary care clinics in Ontario between February 2008 and April 2011.  Similar to the development cohort, participants were recruited within 21 days of their collision. These participants were also interviewed approximately 6, 12, 26, 36 and 52 weeks post-injury. Blind assessment, sample size and the flow of participants are described elsewhere. 
Primary outcome: self-reported global recovery
In both studies, recovery was assessed by asking “How well do you feel you are recovering from your injuries?” The response options included:
(1) “all better (cured),”
(2) “feeling quite a bit improved,”
(3) “feeling some improvement,”
(4) “feeling no improvement,”
(5) “getting a little worse,” or
(6) “getting much worse.” [18, 20]
The self-reported global recovery question is reliable and valid in individuals with NAD secondary to a traffic collision. [18, 20–22] We defined recovery as feeling “all better (cured)” or “feeling quite a bit improved”. [18, 20] A participant was considered recovered on the first interview when recovery was reported.
Secondary outcomes: claim Closure
Claim closure has been used as an indicator of recovery in previous research. [1, 23–25] We defined claim closure as the number of days from the date of the injury to the date on which the claim was closed (i.e. payments ceased and a final agreement was reached between the insurer and the claimant).  Closure usually coincides with the end of treatment, the attainment of maximal medical improvement, or with the end of income-replacement payments. Individuals with NAD who close their claims have lower neck pain intensity, better physical functioning and fewer depressive symptoms than those who keep their claims open. [26, 27]
We used Cox regression to derive and replicate the prediction models.  We performed a complete case analysis. We described the association between predictors and time-to-recovery/time-to-claim-closure using hazard rate ratios (HRR) and 95% confidence intervals (CI). All analyses were conducted using SAS Version 9.3 and Stata MP 12. In the Saskatchewan cohort, we defined recovery as the time to first report of recovery with no subsequent report of having failed to recover. In the Ontario trial, we defined recovery as the time to first report of recovery. Therefore, recovery could occur at any follow-up period (6 weeks, 3, 6, 9 or 12 months) for self-reported recovery and at any time within the study period for claim closure. Observations on participants who did not recover were censored at the last report of non-recovery or at the 12–month follow-up. We handled ties using Efron approximation. 
We examined the crude associations between each prognostic factor and the outcome using Cox regression. The statistical significance (p < 0.05) was assessed with the Wald Chi-square test.  Factors significantly associated with the outcome were included in the multivariable model. The proportionality assumptions were tested with Kaplan-Meier (categorical variables) curves and using time-interactions. Continuous measures (age, percent of body in pain, WDQ) were also categorized and tested as categorical variables.
Development of the Prediction Model
Prior to inclusion of potential prognostic factors in the multivariable models, we tested all variables to determine which were associated with outcomes of self-reported recovery and claim closure and we tested for multicollinearity. Pearson correlations were used to assess normally distributed variables and Spearman rank correlations to assess nominal and ordinal variables or variables with skewed distributions. Pairs of variables with correlations of r = 0.7 or more were considered as highly associated for inclusion in the same predictive model.  Second, we assessed the variance inflation factors (VIF) of each variable in each model. We considered VIF greater than 2.5 to designate multicollinearity. 
We used three stages to develop our models. In Stage 1, we built a model that included Phase III prognostic factors available in the Saskatchewan dataset. Our datasets did not specifically ask about pain-related anxiety or worry; however, we made the assumption that participants answered about anxiety or worry related to their NAD pain-related symptoms. In Stage 2, we added prognostic factors supported by Phase I and II evidence. Finally, in Stage 3, we tested the contribution of additional clinically relevant variables. A lack of significant association was identified by HRR 95% confidence intervals that cross one.  We removed factors in a backward fashion, starting with the least statistically significant, until all associations remaining in the model were significant. We did not test for interactions.
We determined the predictive ability of each model using Harrell’s concordance C-statistic.  The C-statistic estimates the predictive ability of the model by computing the area under a receiver-operating curve (AUC).  The C-statistic ranges from 0 to 1; a value of 0.5 or less denotes that the prediction is due to chance while values approaching 1 demonstrate better predictive ability. 
Internal validation of the Prediction Model
We tested the internal validity of our predictive models using bootstrapping.  We performed 3000 repetitions per model. Our model was deemed internally valid if the standard error of the point estimate following bootstrapping did not change significantly compared to the error obtained at the development stage.
Repeating the Prediction Models in the Ontario Population
We tested our models by repeating them using data from a randomized controlled trial conducted in Ontario. [3, 19]We replicated the models developed in the Saskatchewan population by testing the Stage 1 and 2 models with the same or similar variables available in the Ontario population [Table 1]. In the Saskatchewan sample, we used a combination variable for disability that included “yes” responses for single questions of functional limitations related to work, study or other activities. In the Ontario dataset, disability was tested using the continuous baseline scores of the Whiplash Disability Questionnaire.  In the Saskatchewan sample, we used a percent body in pain variable developed from the body pain diagram and ranging from 0 to 100. In the Ontario dataset, body pain was represented by a 3–category variable stating whether the person had muscle, bone or joint pain. To test the impact of the body pain variable on the prediction of each model, we repeated the prediction models excluding this variable since its construct differed between the two populations.
The Saskatchewan and Ontario samples included 4923 and 340 participants, respectively. Participants in both samples had similar baseline characteristics [Table 1]. Two thirds of participants were female (66. 3% in Saskatchewan; 66.7% in Ontario), most were employed (84.1% for both) and most presented shortly post-collision (9 days [range 0–21]; 6 days [range 0–25]). Participants had a mean age of 38.3 (s.d. 15.1) years in Saskatchewan and 40.5 (s.d. 13.2) years in Ontario. The mean baseline neck pain was 6.5 (s.d. 2.1) in Saskatchewan and 5.7 (s.d. 13.2) in Ontario and the 12–month follow-up rate (84.4%; 78.8%) was also similar. Finally, median time-to-recovery based on the global recovery question was similar (95 days in Saskatchewan and 98 days in Ontario). In the Saskatchewan (development) sample, 4,810/4,923 participants were analyzed in the final model. In the Ontario sample, 311/340 observations were analyzed in the final model.
The proportionality assumption was met for all models. All potential prognostic factors identified in our systematic review were significantly associated with the outcomes.
No multicollinearity was present in the development dataset. VIF values in all models ranged between 1 and 1.7.
Stage 1: prediction model with phase III prognostic factors
All prognostic factors with Phase III evidence were significantly associated with time-to-recovery and the predictive ability of the model based on the C-statistic was higher than chance [Table 2]. Predictors of recovery included in the model: prior neck injury claim with SGI; anxiety/worry (as proxies of pain-related or NAD-related anxiety and fear); depression (Centre for Epidemiological Studies Depression Scale: CES-D > 27) and expectation of recovery at baseline [Table 3].
Stage 2: prediction models with phase I and II prognostic factors
The second stage model included age (in years), prior neck injury claim, percentage of body in pain, baseline neck pain and headache intensity, disability and expectation of recovery as statistically significant predictors of time-to-recovery [Table 3]. The C-statistic for this model was slightly superior to the previous model [Table 2].
Stage 3: prediction models with clinically relevant prognostic factors
The addition of clinically relevant variables did not improve the model built in Stage II (C-statistic ranged from 0.64 to 0.65).
The final model
Based on the Saskatchewan sample, the predictive ability of the Stage 2 model was higher than the Stage 1 model. Therefore, the final model demonstrated that older age, prior neck injury claim, greater percentage of body in pain, greater baseline neck pain, greater baseline headache intensity, disability and poor expectation of recovery are predictive of greater time-to-recovery [Tables 2 and 3].
Similarly in the Saskatchewan sample, we found that all prognostic factors with Phase III evidence remained significantly associated with time-to-claim-closure [Table 4], the predictive ability of the model based on the C-statistic was higher than chance [Table 2] and clinically relevant variables did not improve the model built in Stage II (C-statistic ranged from 0.64 to 0.65). Therefore, the final model predicting longer time-to-claim closure included: older age, prior neck injury claim, greater percentage of body in pain, greater intensity of baseline neck pain and headache intensity, disability, depression and poor expectation of recovery [Table 4].
Internal validation of the Prediction Models
Bootstrapped models demonstrated similar predictive ability as the development models for time-to-recovery and for time-to-claim closure [Table 2].
Repeating prediction models in the Ontario Population
The final models were repeated in the Ontario dataset [Table 2]. Baseline expectation category of ‘never get better’ was excluded from the repeated model due to the small sample size (n = 8); therefore, the validity results for that model are reported without values for that predictor response. Excluding two variables that differed in construct between study samples did not cause a significant change in the predictive ability of the models; with disability removed (Saskatchewan recovery model: C = 0.643, 95% CI 0.633–0.652; Ontario recovery model: C = 0.637, 95% CI 0.599–0.675) or percent of body in pain removed (Saskatchewan recovery model: C = 0.641, 95% CI 0.631–0.651; Ontario recovery model: C = 0.649, 95% CI 0.613–0.686).
Summary of evidence
We developed and validated a predictive model of recovery using strong methodology and results from a best-evidence synthesis.  Our final model demonstrated that expectation of recovery, age, prior neck injury claim, percentage of body in pain, baseline neck pain and headache intensity, and disability are predictive of longer time-to-recovery. In addition, depressive symptoms were also predictive of longer time-to-claim closure. Models for both self-reported recovery and insurance claim closure included the same predictors with the exception of depression being an extra variable in the prediction of claim closure. However, these models had limited predictive ability.
Comparison with previous literature
We identified twelve previous predictive models for the recovery of NAD secondary to traffic collisions. The previous studies used similar outcomes and follow-up periods. The most common outcomes reported in these studies were disability [4–7, 9, 10, 13], persistent neck pain [12, 14] and self-reported recovery. [8, 11, 15] Most studies assessed their outcomes at one-year post-injury. [4–6, 9, 10, 12, 13] Participants were recruited from various sources (e.g., emergency departments, general practitioner practices, rehabilitation clinics, insurance companies and advertisements). In addition, prognostic factors tested in the models and methods used to develop and validate the models varied greatly among studies.
Our results differ from previous predictive models due to methodological differences. Previous models were developed in small samples that limited the number of modelled prognostic factors. Small sample sizes can lead to model uncertainty and unreliable predictions.  We developed our models in a large, population-based inception cohort of patients with incident traffic-related NAD (N = 4,923) and replicated them using data from a smaller randomized controlled trial with a similar patient population (N = 340).
Important differences between our models and the ones proposed by Atherton et al.:
(1) method of analysis (logistic regression);
(2) ill-defined method of selection of prognostic factors; and
(3) failure to report predictive ability statistics;
(4) use of all incident injury cases. 
Only one previous study internally validated their model  and five externally validated their model. [5–7, 15, 16] Three external validation studies included prognostic factors not available to us (i.e. measure of sympathetic vasoconstriction, range of motion assessments, hyperarousal subscale). [6, 7, 16] One study did not use an appropriate sample to validate their model (floor effect).  Finally, one study’s population may differ from ours since they recruited through emergency departments only. 
We derived and validated a clinical prediction model for patients with post-collision neck pain. Our models predict self-reported global recovery and time to claim closure beyond chance; however, further investigation is needed before the models are recommended for clinical practice. Specifically, our models need to be validated in different populations and clinical settings, and an impact analysis needs to be conducted.  Impact analysis are necessary to determine the effect of a clinical prediction model on clinician behaviour and patient outcomes.
Study strengths and limitations
Our study has strengths. We used a large, population-based inception-cohort to develop our models which gives robust model estimates and statistical certainty.  We replicated our models in Ontario, a distinct population with similar demographic and clinical characteristics.  Both cohorts included incident cases only. Finally, we used evidence-based methods to identify prognostic factors and applied data-based methods to identify additional clinically relevant prognostic factors.
Our study also has limitations. Although the Saskatchewan and Ontario samples had similar baseline characteristics, the methods of recruitment for the two studies were different. Our dataset provided a predetermined set of potential prognostic factors. Therefore, we were not able to test additional potential factors including passive coping, cervical radiculopathy, post-crash cold pain threshold, symptoms of Acute Stress Disorder/Post-Traumatic Stress Disorder and kinesiophobia. It is possible that adding these variables to our models could improve their predictive ability. Our anxiety or worry questions were not specifically pain-related and may have included generalized anxiety or worry. However, we assumed that participants answered about their NAD pain-related anxiety or worry since these questions were asked as part of a questionnaire about their NAD-related experience. Our models demonstrated a limited predictive ability; however, this limited predictive ability is predictive beyond chance. Considering that predictive ability of clinical judgement alone is unknown, a model that provides predictive capacity beyond chance adds some certainty to clinical practice. Therefore, we recommend the use of these predictive models until models with higher predictive ability are identified likely including predictive factors not included in our models. Expectation of recovery is the most predictive factor identified by our models and this is a factor that clinicians can focus on with their patients immediately to improve recovery.
Future prediction models should consider potential prognostic factors identified by the best-evidence synthesis that were not available in our datasets. Furthermore, future research should consider factors that have not been studied-to-date that could add to the limited predictive ability of current predictive models.
We developed clinical prediction models that predict recovery and claim closure in individuals with NAD following traffic collisions. Prognostic factors included expectation of recovery, age, having a prior neck injury claim, percentage of body in pain, baseline neck pain and headache intensity, and disability. In addition, depressive symptoms remained predictive in the model predicting claim closure. Our models have limited predictive ability and require an impact analysis before being used in clinical settings.
Supplementary Material 1: TRIPOD Checklist
Supplementary Material 2: Additional variables tested in Stage 3 of the clinical prediction model development
area under a receiver-operating curve
Centre for Epidemiological Studies Depression Scale
95% confidence intervals
Hazard rate ratios
neck pain and associated disorders
Saskatchewan Government Insurance
Variance inflation factors
Whiplash Disability Questionnaire
MS contributed to the design of the work; conducted the analysis; contributed to the interpretation of the results and drafted the work. PC and EB contributed to the design of the work; acquired the data; oversaw the analysis; contributed to the interpretation of the results and substantively revised the paper. LC, HS, and JDC contributed to the design of the work; acquired the data; contributed to the interpretation of the results and substantively revised the paper. RJB contributed to the design of the work; contributed to the interpretation of the results and substantively revised the paper. All authors have approved the submitted version and have agreed both to be personally accountable for the author’s own contributions.
Ethical approval and Consent to participate
The study was approved by the Ontario Tech University Research Ethics Board (#13–107).
PC reports research grants unrelated to this work from Canadian Institutes of Health Research, Canada Research Chair Program, Canadian Chiropractic Association, Canadian Chiropractic Research Foundation, College of Chiropractors of British Columbia; travel expenditures unrelated to this work from Eurospine; fees medical-legal expertise from the Canadian Chiropractic Protective Association and NCMIC. No other conflicts to declare.
Carroll, LJ, Holm, LW, Hogg-Johnson, S et al.
Course and Prognostic Factors for Neck Pain in Whiplash-associated
Disorders (WAD): Results of the Bone and Joint Decade 2000-2010
Task Force on Neck Pain and Its Associated Disorders
Spine (Phila Pa 1976). 2008 (Feb 15); 33 (4 Suppl): S83–92
Shearer HM, Carroll LJ, Côté P, Randhawa K, Southerst D, Varatharajan S et al.
The Course and Factors Associated with Recovery of Whiplash-associated
Disorders: An Updated Systematic Review by the Ontario Protocol for
Traffic Injury Management (OPTIMa) Collaboration
European J Physiotherapy 2020; 23 (5): 1–16
Clinical prediction models a practical approach to
development, validation, and updating.
New York, NY: Springer; 2009.
Ritchie C, Hendrikz J, Kenardy J, Sterling M.
Derivation of a clinical prediction rule to identify both chronic
moderate/severe disability and full recovery following whiplash injury.
Asenlof P, Bring A, Soderlund A.
The clinical course over the first year of Whiplash Associated
Disorders (WAD): Pain-related disability predicts
outcome in a mildly affected sample.
BMC Musculoskelet Disord. 2013;14.
Sterling M, Hendrikz J, Kenardy J, Kristjansson E, Dumas JP, Niere K, et al.
Assessment and validation of prognostic models for poor
functional recovery 12 months after whiplash injury:
a multicentre inception cohort study.
Sterling M, Jull G, Kenardy J.
Physical and psychological factors maintain
long-term predictive capacity post-whiplash injury.
Bohman T, Cote P, Boyle E, Cassidy JD, Carroll LJ, Skillgate E.
Prognosis of patients with whiplash-associated disorders
consulting physiotherapy: development of
a predictive model for recovery.
BMC Musculoskelet Disord. 2012;13.
Kasch H, Qerama E, Kongsted A, Bach FW, Bendix T, Jensen TS.
The risk assessment score in acute whiplash injury
predicts outcome and reflects biopsychosocial factors.
Kasch H, Bach FW, Jensen TS.
Handicap after acute whiplash injury:
a 1-year prospective study of risk factors.
Dufton JA, Kopec JA, Wong H, Cassidy JD, Quon J, McIntosh G, et al.
Prognostic factors associated with minimal improvement
following acute whiplash-associated disorders.
Atherton K, Wiles NJ, Lecky FE, Hawes SJ, Silman AJ, Macfarlane GJ, et al.
Predictors of persistent neck pain after whiplash injury.
Emerg Med J. 2006;23:195–201.
Hendriks EJM, Scholten-Peeters GGM, Van Der Windt DAWM.
Prognostic factors for poor recovery in
acute whiplash patients.
Hartling L, Pickett W, Brison RJ.
Derivation of a clinical decision rule for whiplash associated
disorders among individuals involved in rear-end collisions.
Accid Anal Prev. 2002;34:531–9.
Rydman E, Ponzer S, Ottosson C, Jarnbert-Pettersson H.
Predicting nonrecovery among whiplash patients in the
emergency room and in an insurance company setting.
Eur Spine J 2016:1–8.
Ritchie C, Hendrikz J, Jull G, Elliott J, Sterling M.
External validation of a clinical prediction rule to predict full
recovery and ongoing moderate/severe disability following
acute whiplash injury.
J Orthop Sports Phys Ther. 2015;45:242–50.
Côté P, Cassidy JD, Carroll LJ, Frank JW, Bombardier C.
A systematic review of the prognosis of acute whiplash and
a new conceptual framework to synthesize the literature.
Cassidy JD, Carroll LJ, Côté P, Frank J.
Does multidisciplinary rehabilitation benefit whiplash recovery?:
results of a population-based incidence cohort study.
Côté P, Boyle E, Shearer HM, Stupar M, Jacobs C, Cassidy JD, et al.
Is a Government-regulated Rehabilitation Guideline More Effective
than General Practitioner Education or Preferred-provider
Rehabilitation in Promoting Recovery from Acute
Whiplash-associated Disorders? A Pragmatic
Eandomised Controlled Trial
BMJ Open. 2019 (Jan 24); 9 (1): e021283
Ngo T, Stupar M, Côté P, Boyle E, Shearer H.
A study of the test-retest reliability of the self-perceived general
recovery and self-perceived change in neck pain questions in
patients with recent whiplash-associated disorders.
Eur Spine J. 2010;19:957–62.
Ottosson C, Pettersson H, Johansson SE, Nyren O, Ponzer S.
Recovery after minor traffic injuries: a randomized controlled trial.
PLoS Clin Trials. 2007;2:e14.
Carroll LJ, Jones DC, Ozegovic D, Cassidy JD.
How well are you recovering? The association between a simple question
about recovery and patient reports of pain intensity and pain
disability in whiplash-associated disorders.
Disabil Rehabilitation. 2012;34:45–52.
Côté P, Hogg-Johnson S, Cassidy JD, Carroll L, Frank JW, Bombardier C.
Early aggressive care and delayed recovery from whiplash:
isolated finding or reproducible result?
Arthr Rhuem. 2007;57:861–8.
Côté P, Hogg-Johnson S, Cassidy JD, Carroll LJ, Frank JW, Bombardier C.
Initial patterns of clinical care and recovery from whiplash injuries:
a population-based cohort study.see comment.
Arch Intern Med. 2005;165:2257–63.
Cassidy JD, Carroll LJ, Côté P, Lemstra M, Berglund A, Nygren A.
Effect of eliminating compensation for pain and suffering
on the outcome of insurance claims for whiplash injury.
N Engl J Med. 2000;342:1179–86.
Côté P, Hogg-Johnson S, Cassidy JD, Carroll L, Frank JW.
The association between neck pain intensity, physical functioning,
depressive symptomatology and time-to-claim-closure after whiplash.
J Clin Epidemiol. 2001;54:275–86.
Boyle E, Cassidy JD, Cote P, Carroll LJ.
The relationship between insurance claim closure and recovery after
traffic injuries for individuals with whiplash associated disorders.
Disabil Rehabil. 2017;39:889–96.
Survival analysis using the SAS system: a practical guide.
Cary, N.C.: SAS Institute; 1995.
Allison PD, Chatterjee S, Hadi AS.
Logistic regression using the SAS System:
theory and application + regression. Analysis by Example:
Pinfold M, Niere KR, O’Leary EF, Hoving JL, Green S, Buchbinder R.
Validity and internal consistency of a whiplash-specific disability measure.
Cowley LE, Farewell DM, Maguire S, et al.
Methodological standards for the development and evaluation
of clinical prediction rules: a review of the literature.
Diagn Progn Res. 2019;3:16.
Return to WHIPLASH
Return to CHRONIC NECK PAIN