- Open Access
Item response theory and factor analysis as a mean to characterize occurrence of response shift in a longitudinal quality of life study in breast cancer patients
Health and Quality of Life Outcomes volume 12, Article number: 32 (2014)
The occurrence of response shift (RS) in longitudinal health-related quality of life (HRQoL) studies, reflecting patient adaptation to disease, has already been demonstrated. Several methods have been developed to detect the three different types of response shift (RS), i.e. recalibration RS, 2) reprioritization RS, and 3) reconceptualization RS. We investigated two complementary methods that characterize the occurrence of RS: factor analysis, comprising Principal Component Analysis (PCA) and Multiple Correspondence Analysis (MCA), and a method of Item Response Theory (IRT).
Breast cancer patients (n = 381) completed the EORTC QLQ-C30 and EORTC QLQ-BR23 questionnaires at baseline, immediately following surgery, and three and six months after surgery, according to the “then-test/post-test” design. Recalibration was explored using MCA and a model of IRT, called the Linear Logistic Model with Relaxed Assumptions (LLRA) using the then-test method. Principal Component Analysis (PCA) was used to explore reconceptualization and reprioritization.
MCA highlighted the main profiles of recalibration: patients with high HRQoL level report a slightly worse HRQoL level retrospectively and vice versa. The LLRA model indicated a downward or upward recalibration for each dimension. At six months, the recalibration effect was statistically significant for 11/22 dimensions of the QLQ-C30 and BR23 according to the LLRA model (p ≤ 0.001). Regarding the QLQ-C30, PCA indicated a reprioritization of symptom scales and reconceptualization via an increased correlation between functional scales.
Our findings demonstrate the usefulness of these analyses in characterizing the occurrence of RS. MCA and IRT model had convergent results with then-test method to characterize recalibration component of RS. PCA is an indirect method in investigating the reprioritization and reconceptualization components of RS.
Health-related quality of life (HRQoL) is a subjective clinical endpoint that has been increasingly important in health outcomes research and particularly in cancer clinical trials over the past two decades  as well as in breast cancer . Although overall survival is still considered as the primary objective and the primary endpoint in many studies, most clinical trials now integrate HRQoL as an endpoint in order to investigate the clinical benefit for the patient.
One major objective of measuring HRQoL over time is determining the extent to which treatment toxicities or disease progression can affect patients’ HRQoL level. However, self-assessment of HRQoL is subjective, i.e. it is dependent on the patient's internal standards and definition of HRQoL [3–5]. As patients can adapt to disease and treatment toxicities, their health and HRQoL expectations can also change over time. These changes could result in a response shift (RS) effect [6–8].
RS can be defined as “a change in the meaning of one’s self-evaluation of a target construct as a result of: (a) a change in the respondent’s internal standards of measurement (i.e. scale recalibration); (b) a change in the respondent’s values (i.e. the importance of component domains constituting the target construct, [reprioritization]) or (c) a redefinition of the target construct (i.e. reconceptualization)” .
Different methods have been proposed to assess RS [10–13]. The most widely used is the “Then-test” method, which assesses patients’ pre-test HRQoL levels retrospectively. The test involves asking patients post-treatment to provide their current levels (post-test) but also their pre-test levels in retrospect (then-test). This method is based on the assumption that patients rate their HRQoL post-test and pre-test levels with the same criteria, since the assessments occur at the same time point. The recalibration component of RS should thus be taken into account when comparing post-test and then-test scores. Comparing the mean of the pre-test and then-test scores explores recalibration component of RS .
Statistical methods have also been investigated to detect RS. First, factor analyses have been explored to detect RS [14, 15]. An alternative to investigate RS with factor analysis is the use of structural equation modeling (SEM) [11, 16, 17]. These models can evaluate all types of RS if they are experienced by a substantial part of individual in the population analyzed . These models are based on means and covariance structures and rely on observed scores. At this time and to our knowledge, these models have never been applied to an European Organization for Research and Treatment of Cancer (EORTC) HRQoL questionnaire in order to highlight RS effect. Principal Component Analysis (PCA) is a special case of SEM. Item Response Theory (IRT) could also be considered to explore RS effect but up to now these models remain less applied and mainly through differential item functioning [18, 19]. Contrary to SEM, IRT models are not based on the observed scores but directly on items answers. In fact, in SEM, the raw score is assumed to be a good representation of the latent trait (i.e. HRQoL), while in IRT the items responses play a key role and the relation between the items responses and the latent trait are not linear in IRT.
While occurrence of RS in HRQoL studies has been demonstrated , approaches that could reinforce the proof that each component of RS occurs should be investigate to complement the results from other methods. All methods that highlight RS have their strengths and weaknesses then similar trend obtained from different tools should increase accuracy of the results characterizing RS occurrence and the confidence of the results. Moreover, the studies to detect the occurrence of RS are generally performed with two measurement time points while in oncology clinical trials more than two assessments is usually planned. Then we need also tools for the longitudinal analysis of the potential occurrence of RS.
The intent of this study was thus to investigate statistical methods to characterize the occurrence of RS for HRQoL in breast cancer (BC) patients.
The primary objective was to assess if Multiple Correspondence Analysis (MCA), which is a factor analysis and a model of IRT named the Linear Logistic model with Relaxed Assumptions (LLRA), had convergent results with then-test method to characterize recalibration component of RS.
The secondary objective was to assess if Principal Component Analysis (PCA), which is another factor analysis model, could be a valuable tool to longitudinally identify the reconceptualization and reprioritization components of RS independently of the occurrence of the recalibration component of RS.
Patients and eligibility criteria
A prospective, multicenter, randomized cohort study was performed in the cancer care centers at Dijon, Nancy, and the university hospitals of Strasbourg and Reims (cities of France). It is a collaboration between different teams with complementary skills and an interest in the topic and all these teams are involved in quality of life research. All women initially hospitalized between February 2006 and February 2008 for diagnosis or treatment of primary or suspected BC were eligible for inclusion. We anticipated that patients with no confirmed BC could constitute a control group. Nevertheless, due to the low effective (less than 10% of patients included) they could not constitute a larger control group then they were excluded from the analyses. Women with cancer other than BC, already undergoing BC treatment, or with a previous history of cancer were excluded. Written informed consent was obtained from every participant and the protocol was approved by Dijon University Hospital Ethics committee .
Health-related quality of life assessment
HRQoL was evaluated using the EORTC QLQ-C30 and EORTC QLQ-BR23 BC specific tool at four time points: at baseline (initial examination or initial hospitalization), at discharge following initial hospitalization, at three months (M3) and six months (M6) [21, 22]. The QLQ-C30 and its BC module BR23 are validated tools in assessing HRQoL in cancer, specifically in BC [21, 22]. The QLQ-C30 comprises 30 items and measures five functional scales (physical, role, emotional, cognitive and social functioning), global health status (GHS), financial difficulties and eight symptom scales (fatigue, nausea and vomiting, pain, dyspnea, insomnia, appetite loss, constipation, diarrhea) . The BR23 module comprises 23 items that generate four functional scales (body image, sexual functioning, sexual enjoyment, future perspective) and four symptom scales (systemic therapy side effects (STSE), breast symptoms, arm symptoms, upset by hair loss) .
Response categories vary from 1 to 4 on a Likert scale for the QLQ-C30 and BR23 questionnaires, with 1 corresponding to the best state for functional scales or no symptoms, and 4 corresponding to the worst state for functional scales or the highest symptomatic level. For sexual dimensions, the response categories are reversed. Scores are generated according to the EORTC Scoring Manual . These scores vary from 0 (worst) to 100 (best) for the functional dimensions and GHS, and from 0 (best) to 100 (worst) for the symptom dimensions.
A five-point difference in EORTC HRQoL scores is considered as the minimal clinically important difference (MCID) .
Then test assessment
In this study, the retrospective pre-test/post-test design was used to detect recalibration . At each follow-up time point, one prospective and one retrospective measurement were performed. The retrospective assessments at the end of initial hospitalization and at M3 refer to baseline HRQoL. At M6, the retrospective measurement refers to HRQoL at M3. The order of the then-test and post-test of HRQoL questionnaires was randomized with a 1:1 allocation and stratification by center to assess the impact of the order on RS occurrence and estimate. In arm A, the order of the questionnaires was post-test/then-test. In arm B, the order was then-test/post-test. Authorization was obtained from the EORTC HRQoL unit to adapt the HRQoL questionnaires (EORTC QLQ-C30 and module BR23) to the then-test assessment. The impact of the retrospective or prospective administration of the questionnaire on RS occurrence has already been analyzed in a previous paper showing no order effect and is not treated in the present paper .
Treatments as well as clinical and sociodemographic variables were recorded at inclusion.
Studied population and missing data
Variables collected at baseline were described with median and range for continuous variables and percentage for qualitative variables, with percentage of missing data. No imputation was performed on missing items. Scores were calculated if at least half of the items were answered according to the recommendations of the EORTC scoring manual . No imputation was performed on missing scores.
MCA and LLRA were both performed on patients with all items of the studied dimension (each dimension of the QLQ-C30 and the QLQ-BR23) filled out at the then-test and the pre-test measurement time points and with a MCID between then-test and pre-test of at least 5 points for the given dimension. This selection was done in order to retain a clinically meaningful difference of the recalibration occurrence.
PCA were performed on patients with all scores available at the four prospective measurement times for one questionnaire (QLQ-C30 or BR23).
For each analysis, patients retained were compared to those excluded according to baseline characteristics in order to check the random missing data profile and then a possible selection bias.
For each score, the mean difference (MD) between each then-test and the corresponding pre-test was calculated and described as mean (SD). The existence of a significant recalibration was tested with a Wilcoxon matched pairs test. The effect size was calculated in order to assess the magnitude of RS effect and was defined as the mean change score between the then-test and the corresponding pre-test dividing by the standard deviation of patients at the prospective measurement time.
The primary objective was to assess if MCA and the LLRA model of IRT had convergent results with the then-test method to characterize the recalibration component of RS.
Firstly, recalibration was thus explored by MCA . MCA is a factor analysis dedicated to qualitative variables and can identify links between categories of polytomous variables. This method is thus well adapted to the items constructed on a Likert scale. This analysis was applied to items of each dimension according to the Then-test method, i.e. with pre-test and then-test measures of the same HRQoL. Only recalibration was explored with this method since only one dimension was included. Therefore, recalibration was confirmed by a correlation between two different response categories of the same item measured at pre-test and at then-test measurement time . The study was limited to the first two axes.
IRT and Classical Test Theory differ in terms of score calculation. Classical Test Theory is mainly based on observed scores while in IRT, item responses play the key role: IRT models the item responses to the latent trait by a probabilistic model. The raw score is thus not considered as a good representation of the latent trait but the response to each item is considered directly. The relationship between the observed score and the latent trait is no longer linear. They are generally linked and modeled by a logistic function. The IRT models introduce the concepts of item easiness parameters and person parameters.
The person parameter corresponds to the level of the patient on the latent trait (e.g. the level of HRQoL). The item parameter is the location of the item on the latent trait and corresponds to a level of difficulty or easiness in this model.
The LLRA requires neither items’ unidimensionality nor distributional assumptions about the population of subjects . In addition, the LLRA can fit with polytomous responses and was developed in order to measure the change occurring between several measurement time points . To give up the unidimensionality of the items, items have to be measured at two measurement time points or more [32, 33].
The main idea of the LLRA model is not to consider longitudinal change as a change in person parameters, but rather as a change in item parameters. In this way, person parameters are fixed over time and only item parameters vary. Since person parameters are nuisance parameters, we can estimate the item parameter trend instead of the person parameter trend by conditional maximum likelihood . Indeed, fewer parameters have to be estimated and they did not depend on the sample considered.
One item I with parameter βi evaluated twice on an individual can be seen as a pair of virtual items I*1 and I*2 with two item parameters β*i1 and β*i2 respectively. For the pre-test, β*i1 = βi while for the then-test β*i2 = βi + τ where τ is the upward or downward trend effect of item easiness parameter. This parameter is targeted by LLRA . In cases of polytomous items, for each item with (m + 1) response categories there are m category parameters. The trend parameter τ is the same for each category parameter. The design matrix was constructed such that there is one trend parameter for each item. If possible, the trend was generalized for all items of a dimension. The general form of LLRA, a longitudinal IRT model adapted to polytomous items, is based on the partial credit approach .
A positive (or negative) trend τ for one item implies that the item easiness parameter increases (or decreases) at the time of the then-test measurement compared to the pre-test measurement. Patients choose higher (or lower) response categories in the retrospective then-test measure than in the prospective one for this item. In this way, recalibration would be indicated by a significant positive or negative trend for one dimension.
Convergent results between MCA and IRT would correspond to:
- a significant positive trend parameter for IRT and some upward recalibration profiles highlighted by MCA (i.e. patients choose upper response categories at the then-test assessment as compared to the prospective measurement time) more than some downward recalibration profiles (patients choose lower response categories at the then-test assessment as compared to the prospective measurement time).
- a significant negative trend parameter for IRT and some downward recalibration profiles highlighted by MCA (patients choose lower response categories at the then-test assessment) rather than some upward recalibration profiles (patients choose upper response categories at the then-test assessment).
- an insignificant trend parameter for IRT and well-balanced recalibration profiles highlighted by MCA (as many patients choose higher than lower response categories at the retrospective measurement time as compared to the prospective measure).
GHS was excluded from MCA and LLRA because of the high number of response categories. There are seven response categories for both items measuring GHS. To apply a longitudinal model of IRT, all seven categories have to be represented at each measurement time point, which was not the case in the present study. GHS was excluded from MCA in order to be consistent with IRT.
Reprioritization and reconceptualization
The secondary objective was to assess if PCA could be a valuable tool to longitudinally identify the reconceptualization and reprioritization components of RS independently of the occurrence of recalibration component of RS.
PCA was performed on patients with all scores available at all prospective measurement times and for one questionnaire (QLQ-C30 or BR23) on the scores generated for all dimensions of each prospective questionnaire [12, 14, 15, 36]. PCA was performed only for one questionnaire in order to have clear and understandable graphs. Reprioritization was indicated by a change in scores generating the first two principal components: scales generating the first principal component are a priority to patients while those generating the second principal component are secondary. Changes occurring at the first principal component are considered as major and those occurring at the second principal component as minor. In this way, changes were qualified in the first axis of “major reprioritization” and in the second axis of “secondary reprioritization”. The study was limited to the first two principal components, according to the Scree test . Reconceptualization was reflected by a change in the structure of the graph of correlations between scores and principal components, as well as in the connection or opposition of some scores. Concerning the module BR23, sexual enjoyment and hair loss were excluded from the analysis given the number of missing values.
The statistical significance level was reduced to p = 0.002 for all analyses in order to prevent false positive results due to the number of multiple comparisons performed (alpha risk 0.05 divided by the number of dimensions analyzed).
Between February 2006 and February 2008, 381 patients were included in the four participating centers. Mean age was 58.4 (standard deviation = 11) years. Three hundred and forty (89%) patients had confirmed BC. Complete clinical and pathologic features of the population are given in Table 1.
HRQoL questionnaires completion and missing data
Table 2 describes the number of completed QLQ-C30 and BR23 questionnaires at each measurement time.
(93%) patients had at least one HRQoL score at baseline, 311 (91%) on discharge following initial hospitalization (i.e. after surgery), 304 (89%) at M3 and 290 (85%) at M6.
Median time for HRQoL assessments between baseline and the discharge following initial hospitalization was 6 days, range [1.5; 81.5].
Patients retained for MCA and LLRA with a 5-point MCID were similar to those excluded according to baseline characteristics for each analysis (data not shown). Patients retained for PCA with all the four prospective measurement times were similar to those excluded except that they seem to be older (data not shown).
After surgery (Table 3), the recalibration effect was statistically and clinically significant for emotional functioning (MD = 5.36) and future perspectives (MD = 7.41) dimensions with a moderate effect size (0.21 and 0.24 respectively). At M3, the recalibration effect was statistically and clinically significant for role (MD = -6.50), emotional (MD = 6.97) and social functioning (MD = -5.01), insomnia (MD = -6.93), body image (MD = -8.16) and future perspectives (MD = 6.95) dimensions.
At M6, the recalibration effect was statistically and clinically significant for physical (MD = 5.10), role (MD = 8.55) and social functioning (MD = 6.02) and for fatigue (MD = -11.03), pain (MD = -6.02), insomnia (MD = -5.64), body image (MD = 7.78) and breast symptoms (MD = -7.28).
Recalibration and MCA
All results obtained on the QLQ-C30 and QLQ-BR23 are summarized in Table 4 and in Table 5, respectively. Qi_k (resp. Ri_k) refers to the k-th response category of the i-th item of a prospective (resp. retrospective) questionnaire on the graph.
Figure 1 presents the graph obtained for baseline and the then-test performed after surgery for role functioning. 272 patients answered the items 6 (Were you limited in doing either your work or other daily activities?) and 7 (Were you limited in pursuing your hobbies or other leisure time activities?) measuring the role functioning scale at baseline and at the retrospective measurement after surgery referring to the baseline HRQoL. Response categories are coded “1/2/3/4” respectively for “Not at all/A little/Quite a bit/Very much”. Among these patients, 84 (31%) had a MCID of at least 5 points between the two measures. Figure 1 highlights two main patterns of recalibration: patients who had reported an excellent role functioning at baseline (i.e. had chosen response category 1 for both items 6 and 7 at baseline) and who had declared a slightly worse role functioning level when they reevaluated this dimension retrospectively after surgery (i.e. chose response category 2 for both items 6 and 7 at the retrospective measurement time), and vice versa (i.e. had chosen response category 2 for both items 6 and 7 at baseline and had chosen response category 2 for both items 6 and 7 at the retrospective measurement time). The first profile is suggested by an association between Q6_1, Q7_1, R6_2 and R7_2. The reverse profile corresponds to the association between Q6_2, Q7_2, R6_1 and R7_1. Recalibration profiles are less explicit for patients who had reported a low role functioning at baseline (i.e. had chosen response category 3 or 4 for both items 6 and 7 at baseline). Indeed, these patients are fewer and they did not follow a unique recalibration profile. Patients who had reported a relatively low role functioning at baseline (i.e. had chosen response category 3 for both items 6 and 7 at baseline) either tended to revise their opinion upwards or downwards by choosing either response category 4 or response category 2 for both items 6 and 7 at the retrospective assessment after surgery.
Recalibration and IRT using LLRA
Positive trend parameters (τ = +) indicated that at the pre-test measurement, patients had overestimated their functional level or had underestimated their symptomatic or sexual level.
Based on the first retrospective reassessment of their baseline HRQoL (Table 6), patients had significantly underestimated their emotional (τ = -0.62, p < 0.001) and cognitive functioning (τ = -1.15, p < 0.001) and their level of arm symptoms (τ = 0.64, p < 0.001). Patients also had overestimated the presence of insomnia (τ = -0.49, p = 0.001) and diarrhea (τ = -1.34, p < 0.001).
Based on the second retrospective reassessment of their baseline HRQoL, patients had significantly overestimated their role (τ = 0.71, p < 0.001) and social functioning (τ = 0.66, p < 0.001), their body image (τ = 0.83, p < 0.001) and insomnia level (τ = -0.49, p < 0.001) and had underestimated their level of pain (τ = 0.37, p = 0.001), and arm symptoms (τ = 0.86, p < 0.001).
Regarding HRQoL at M3, patients had significantly underestimated their physical (τ = -0.84, p < 0.001), role (τ = -0.60, p < 0.001), cognitive (τ = -0.53, p < 0.001) and social functioning (τ = -0.51, p < 0.001) as well as their body image (τ = -0.66, p < 0.001). Patients also had overestimated their emotional functioning (τ = 0.22) and their levels of fatigue (τ = -0.54), pain (τ = -0.54), arm (τ = -0.37) and breast (τ = -0.76) symptoms (p ≤ 0.001).
To summarize, after surgery, the recalibration effect was statistically significant for 6/22 dimensions of the QLQ-C30 and BR23 according to the IRT model while for the then-test method it was only clinically significant for 2 of these dimensions (emotional functioning and future perspective). At M3, the recalibration effect was statistically significant for 7/22 dimensions of the QLQ-C30 and BR23 according to the IRT model and to the then-test method except for the arm symptoms (MD = 2.43, p = 0.011). At M6, the recalibration effect was statistically significant for 11/22 dimensions of the QLQ-C30 and BR23 according to the IRT model. The same results were observed for the then-test method except for the emotional (p = 0.054) and cognitive functioning (p = 0.004) and the arm symptoms (p = 0.010). A significant and clinically recalibration was also observed according to the classical then-test method for insomnia (MD = -5.64, p = 0.002) and not according to the IRT model (p = 005).
Reprioritization and reconceptualization with PCA
PCA was performed on each prospective measure performed at baseline, post-surgery, and at M3 and M6, on patients with all scores available at the prospective measurement time for one questionnaire (QLQ-C30 or QLQ-BR23).
For the QLQ-C30 (respectively QLQ-BR23), 192 (50.4%) patients (respectively 154 (40.4%)) were retained with all scores available at the four prospective measurement times.
Concerning the QLQ-C30 (Figure 2), functional scales became more interrelated and related to the first principal component, reflecting a strong positive correlation between these scales (Table 7). This is observed at each measurement time point. Fatigue and pain remained strongly correlated at each measurement time point, a little less at M3. Diarrhea and financial difficulties were correlated just after surgery (Figure 2B). Nausea and vomiting were correlated to appetite loss at M6 (Figure 2D).
Reprioritization was mainly secondary, as it mostly affected the second principal component: fatigue and pain were still priority symptoms at each prospective measure, since they were still highly correlated with the first principal component at each prospective measure. These symptoms mainly affected physical, social and role functioning as well as GHS. At M6, all functional scales were affected by these symptoms. The second axis highlighted secondary symptoms, namely insomnia at baseline, diarrhea after surgery, nausea and vomiting at M3 and M6 (Table 7).
Concerning the QLQ-BR23, STSE affected the patients’ body image since these dimensions remained strongly negatively correlated (Figure 3). These scales were highly correlated with the first principal component (Table 7). Post-surgery, arm symptoms as well as body image and STSE were significantly correlated with the first principal component. Thus, arm symptoms were equally important to body image and STSE regarding patient HRQoL level. At M3 and M6, future perspectives became significantly correlated with the first principal component and thus gained importance, highlighting a major reprioritization. Breast symptoms and sexual functioning were not significantly associated with the two first axes at baseline: they were minor dimensions. After surgery, only sexual functioning was a relevant factor, while at M3 and M6, only breast symptoms were relevant.
Reconceptualization is illustrated by changes in correlations (positive or negative) between variables at each measurement time point. At each measurement time point, body image score was opposed to STSE score. Post-surgery, STSE was associated with arm symptoms. At M3 and M6, a high body image score was associated with high future perspective.
The present study demonstrates that response shift effect occurred in patients with primary breast cancer, just after surgery, as well as at 3 and 6 months. The intent of this study was to investigate statistical methods to characterize the occurrence of response shift in breast cancer patients.
Our primary objective was to assess if MCA and IRT model had convergent results with the then-test method to characterize recalibration component of RS.
Both methods explored are convergent to the then-test method. When the then-test method highlighted a clinically significant recalibration, MCA highlighted a general trend to overestimate or underestimate their HRQoL level choosing higher or lower response categories according to the direction of the recalibration effect. IRT model showed a statistically significant general trend (positive or negative) of item easiness parameter with the exception of insomnia at 6 months for which a recalibration is not detected by IRT at the alpha level p = 0.002 but borderline (p = 0.005). When the mean difference between the then-test and the pre-test is not significant, i.e. no clinically significant recalibration occurs, MCA highlighted as many patients recalibrate upward than downward their HRQoL level and then there is not a general trend to overestimate or underestimate their HRQoL level. The IRT model also highlighted that the trend of item easiness parameter is not significant. However, some discrepancies are observed: for 4/6 dimensions for which a recalibration was detected by IRT and not significant according to the then-test method at the first retrospective assessment, 1/7 at the second one, and 3/11 at the last retrospective assessment. Thus, the IRT model detects more recalibration effect than the classical then-test method.
MCA and IRT models highlight convergent results. Based on the retrospective assessment of baseline HRQoL after surgery and according to the LLRA, the trend of item easiness parameter is insignificant for social functioning. The corresponding MCA shows readjustment between response categories 1 and 2 and between response categories 3 and 4. In this way, as many patients had chosen higher than lower response categories at the retrospective measurement time as compared to the baseline measure, which is consistent with the LLRA. Based on retrospective assessment of baseline HRQoL at 3 months and according to the LLRA, patients had overestimated their body image with a positive trend of item easiness parameter. Regarding the corresponding MCA, it highlights a readjustment from response categories 2 to response categories 3 and from response categories 3 to response categories 4. In this way, patients had chosen higher response categories at the retrospective measurement time compared to the prospective measure indicating an overestimation of body image at baseline. Based on the retrospective assessment of the three months HRQoL at 6 months, patients had overestimated their pain level with a negative trend of item easiness parameter according to the LLRA. Regarding the MCA performed on pain at the same measurement times, it shows a recalibration between response categories 1 and 2, and only from response categories 4 to 3, not from 3 to 4.
The secondary objective was to assess if PCA could be a valuable tool to longitudinally identify the reconceptualization and reprioritization components of RS independently of the occurrence of recalibration component of RS. PCA indicated a reprioritization of the HRQoL domains as evaluated by the QLQ-C30. Patients’ anxiety probably related to the diagnosis of cancer and surgery seemed to be a major concern at baseline before the start of treatment, along with insomnia, which generated the second principal component, after fatigue and pain, which generated the first principal component. After surgery, diarrhea symptoms increased in importance, reflecting the impact of treatment. These results underline how patients adapt to their disease. At 3 months and 6 months, nausea and vomiting were more important as compared to diarrhea, also reflecting the toxicities of cancer treatment, especially chemotherapy. Regarding the QLQ-BR23, patients with a high level of systemic therapy side effects after surgery also tended to report a high level of arm symptoms, which can be due to the recent surgery. From 3 months, arm symptoms become less important, while future perspectives gained importance for primary breast cancer patients. Our results indicate that there is no correlation between breast symptoms and sexual functioning at 3 and 6 months.
No reprioritization was observed for the QLQ-C30 and QLQ-BR23 between the measures at M3 and M6. Patients seemed to assess their HRQoL with the same relative importance at 6 months as they did at 3 months suggesting that after treatment initiation they have a more “stabilized” appreciation of HRQoL dimensions.
The reprioritization of symptomatic scales enables interpretation of HRQoL levels and changes and impact of treatments and disease on HRQoL. Then based on these results we suggest that the occurrence of the reprioritization component of RS should be taken into account in the interpretation of the results of the longitudinal analysis. Deterioration of a scale, which becomes more important over time for the patient, could have a strong impact on patient’s overall HRQoL level and could indicate priority for care. Conversely, deterioration of a scale, which for the patient loses importance over time, could have a minor impact on patient’s general HRQoL level.
Reconceptualization is reflected by changes in connections and contrasts between variables, and more generally by changes in graph structure of PCA. The functional scales of the QLQ-C30 became increasingly interrelated. When one functional scale is affected by cancer treatment or disease progression, then it is likely that all the other functional scales are affected. Moreover, patients had associated nausea and vomiting to appetite loss at 6 months.
These results suggest that PCA is an indirect method in investigating the reprioritization and reconceptualization components of RS.
The main limitation of this work is the use of the Then-test as the standard method to explore recalibration. The Then-test method is increasingly called into question [41–43], mainly because it can induce a recall bias . Indeed, the second reassessment of baseline HRQoL was three months after baseline and the reassessment of HRQoL at M3 was three months after the prospective measure so it may induce a recall bias.
Schwartz et al. have proposed some guidelines to improve the stringency of the Then-test method . In their paper, Schwartz et al. recommended to include a control group, which would not susceptible to RS. As RS is a treatment-dependent phenomenon, we tried to constitute a control group including patients with only a suspicion of BC. However, the number of patients with no confirmed BC was not sufficient to constitute a control group. Others explanations of detection of RS effect could then also be formulated as social desirability responding. This hypothesis cannot be verified since a control group could not be constituted. As it was recommended by Schwartz et al., we reported effect sizes for recalibration in order to assess the magnitude of RS effect. A Bonferroni correction of type I error rate was performed in order to minimize false positive conclusions. The guidelines also recommended to use internal or external validation approaches as performance-based, perception-based, and evaluation-based items/subscales for internal validation of then-test results and clinical measures indicating health state at baseline and follow-up for external validation. However, we failed to include such approaches in this present study. The instructions of the retrospective questionnaires clearly indicate the patients to think back to the referent time as advisable by Schwartz et al. Moreover, the nomenclature used in this paper to characterize the recalibration component of RS is those recommended . The Then-test method is based on the assumption that patients rate their HRQoL post-test and pre-test levels with the same criteria, since the assessments occur at the same time point. A test of the measurement invariance of the Then-test method would be necessary in this study in order to validate the then-test and to assess the possible recall bias due to its retrospective nature. This would be planned in another analysis using the Oort’s procedure [16, 44].
Based on this study, substituting the then-test with the LLRA and MCA to explore the recalibration component of RS cannot be recommended at this time. Nevertheless, IRT using LLRA could reinforce the Then-test method because of the improved interpretation of recalibration. This model is effective, and the results are clearer, more explicit and easy to summarize and to interpret. These methods should be used in other studies to validate their ability to reinforce the then-test method.
SEM is often used nowadays to demonstrate RS [11, 16, 17, 45–48]. These models are not dependent on the Then-test method. However, they are based on the raw score and not on the items. In this way, IRT as compared to SEM could be more informative. Moreover, at this time, SEM has never been applied to the EORTC HRQoL questionnaires in order to highlight occurrence of the response shift effect.
It would be interesting to compare the statistical method described in the present study (factor analysis and IRT) to SEM applied on prospective measure in another paper in order to check their ability to capture all the three components of RS. There is a need to investigate all these methods using simulated data in order to establish differences using these three methods.
Factor analysis presents the advantage of graphically exploring all the components of RS. This visual representation is interesting in order to explore reconceptualization, which is the most conceptual component of the RS effect. Moreover, at this time, few methods have been proposed to identify this component  and in our point of view no gold standard has emerged. In addition, no additional questionnaires are required for exploring reconceptualization and reprioritization. Thus, the use of PCA on the scores of the main questionnaires seems to be adequate in exploring these components. SEM is also often used to assess these components. However, our objective was to investigate the PCA method already used in the past [14, 15, 36] and not to apply SEM. Indeed, PCA is a special case of SEM.
IRT models and factor analysis are mostly used in the development and validation of HRQoL questionnaires [49–55]. However, several studies have begun to use IRT in longitudinal studies of HRQoL [29, 56–61], underscoring the potential of these models in longitudinal analyses. Moreover, longitudinal IRT model was used in order to characterize recalibration component of RS. Few studies have investigated RS using IRT while differential item functioning based on IRT was also proposed as alternative approach [18, 62, 63].
Finally, PCA were performed on patients with all scores available at all the prospective measurement times. Only 40% to 50% of patients were thus retained in the analysis but these patients were comparable to those excluded according to baseline characteristics except there was an age effect which may reflect a selection bias.
The data presented in this article confirm the potential of IRT models in longitudinal HRQoL studies, especially their ability to characterize more precisely the recalibration component of RS. Our data also underline the interest of PCA to characterize reprioritization and reconceptualization components of RS. These results confirm the need to take recalibration into account when comparing longitudinal HRQoL data between patient groups and the need to explore the other components in order to better interpret results [64, 65]. The items of these questionnaire are prone to response shift effect since they are evaluation-based items. Then an objective assessment by the patient cannot be made. Despite the fact that items of these questionnaires are prone to RS effect. Some work is still needed to provide both a longitudinal analysis method easy to understand for the clinician and to extract the potential measurement bias due to the occurrence of a response shift effect. Another solution would be to develop or use other questionnaires not prone to response shift effect with more performance-based items . Future studies should investigate the ability of these statistical methods to capture all components of RS without the then-test method.
Global Health Status
European Organization for Research and Treatment of Cancer
Health-related quality of life
Item response theory
Linear logistic model with relaxed assumptions
Multiple Correspondence Analyses
Minimal clinically important difference
Principal component analyses
Structural equation modeling
Systemic therapy side effects.
Osoba D: Health-related quality of life and cancer clinical trials. Ther Adv Med Oncol 2011, 3: 57–71. 10.1177/1758834010395342
Montazeri A: Health-related quality of life in breast cancer patients: a bibliographic review of the literature from 1974 to 2007. J Exp Clin Cancer Res 2008, 27: 32. 10.1186/1756-9966-27-32
Ubel PA, Peeters Y, Smith D: Abandoning the language of “response shift”: a plea for conceptual clarity in distinguishing scale recalibration from true changes in quality of life. Qual Life Res 2010, 19: 465–471. 10.1007/s11136-010-9592-x
Wiklund I: Assessment of patient-reported outcomes in clinical trials: the example of health-related quality of life. Fundam Clin Pharmacol 2004, 18: 351–363. 10.1111/j.1472-8206.2004.00234.x
Bullinger M: Assessing health related quality of life in medicine. An overview over concepts, methods and applications in international research. Restor Neurol Neurosci 2002, 20: 93–101.
Gibbons FX: Social comparison as a mediator of response shift. Soc Sci Med 1999, 48: 1517–1530. 10.1016/S0277-9536(99)00046-5
Howard GS, Dailey PR, Gulianick NA: The feasibility of informed pretests in attenuating response shift bias. Appl Psychol Meas 1979, 3: 481–494. 10.1177/014662167900300406
Howard GS, Ralph KM, Gulianick NA, Maxwell SE, Nance SW, Gerber SK: Internal invalidity in pretest-posttest self-report evaluations and a reevaluation of retrospective pretests. Appl Psychol Meas 1979, 3: 1–23. 10.1177/014662167900300101
Sprangers MA, Schwartz CE: Integrating response shift into health related quality of life research: a theoretical model. Soc Sci Med 1999, 48: 1507–1515. 10.1016/S0277-9536(99)00045-3
Korfage IJ, de Koning HJ, Essink-Bot ML: Response shift due to diagnosis and primary treatment of localized prostate cancer: a then-test and a vignette study. Qual Life Res 2007, 16: 1627–1634. 10.1007/s11136-007-9265-6
Oort FJ, Visser MR, Sprangers MA: An application of structural equation modeling to detect response shifts and true change in quality of life data from cancer patients undergoing invasive surgery. Qual Life Res 2005, 14: 599–609. 10.1007/s11136-004-0831-x
Schwartz CE, Sprangers MA: Methodological approaches for assessing response shift in longitudinal health-related quality-of-life research. Soc Sci Med 1999, 48: 1531–1548. 10.1016/S0277-9536(99)00047-7
Sprangers MA, Van Dam FS, Broersen J, Lodder L, Wever L, Visser MR, Oosterveld P, Smets EM: Revealing response shift in longitudinal research on fatigue–the use of the thentest approach. Acta Oncol 1999, 38: 709–718. 10.1080/028418699432860
Schmitt N: The use of analysis of covariance structures to assess beta and gamma change. Multivar Behav Res 1982, 17: 343–358. 10.1207/s15327906mbr1703_3
Ahmed S, Mayo NE, Corbiere M, Wood-Dauphinee S, Hanley J, Cohen R: Change in quality of life of people with stroke over time: true change or response shift? Qual Life Res 2005, 14: 611–627. 10.1007/s11136-004-3708-0
Oort FJ: Using structural equation modeling to detect response shifts and true change. Qual Life Res 2005, 14: 587–598. 10.1007/s11136-004-0830-y
King-Kallimanis BL, Oort FJ, Nolte S, Schwartz CE, Sprangers MA: Using structural equation modeling to detect response shift in performance and health-related quality of life scores of multiple sclerosis patients. Qual Life Res 2011, 20: 1527–1540. 10.1007/s11136-010-9844-9
Craig S: Measuring Change Retrospectively: An Examination Based on Item Response Theory. In “Measuring Behavioral Change: Methodological Considerations.” Symposium Presented at the Annual Conference of the Society for Industrial and Organizational Psychology Edited by: Martineau J. 2000.
Meade AW, Ellington JK, Craig SB: Exploratory Measurement Invariance: A New Method Based on Item Response Theory. Symposium Presented at the 19th Annual Conference of the Society for Industrial and Organizational Psychology, Chicago, IL 2004.
Dabakuyo TS, Guillemin F, Conroy T, Velten M, Jolly D, Mercier M, Causeret S, Cuisenier J, Graesslin O, Gauthier M, Bonnetain F: Response shift effects on measuring post-operative quality of life among breast cancer patients: a multicenter cohort study. Qual Life Res 2013, 22: 1–11.
Aaronson NK, Ahmedzai S, Bergman B, Bullinger M, Cull A, Duez NJ, Filiberti A, Flechtner H, Fleishman SB, de Haes JC, Kaasa S, Klee M, Osoba D, Ravasi D, Robe PB, Schraub S, Sneeuw K, Sullivan M, Takeda F: The European Organization for Research and Treatment of Cancer QLQ-C30: a quality-of-life instrument for use in international clinical trials in oncology. J Natl Cancer Inst 1993, 85: 365–376. 10.1093/jnci/85.5.365
Sprangers MA, Groenvold M, Arraras JI, Franklin J, te Velde A, Muller M, Franzini L, Williams A, de Haes HC, Hopwood P, Cull A, Aaronson NK: The European Organization for Research and Treatment of Cancer breast cancer-specific quality-of-life questionnaire module: first results from a three-country field study. J Clin Oncol 1996, 14: 2756–2768.
Fayers PM, Aaronson NK, Bjordal K, Groenvold M, Curran D: Bottomley AobotEQoLG. In EORTC QLQ-C30 Scoring Manual. 3rd edition. EORTC; 2001. edn; 2001
Osoba D, Rodrigues G, Myles J, Zee B, Pater J: Interpreting the significance of changes in health-related quality-of-life scores. J Clin Oncol 1998, 16: 139–144.
Hoffman DL, De Leeuw J: Interpreting multiple correspondence analysis as an multidimensional scaling method. Mark Lett 1992, 3: 259–272. 10.1007/BF00994134
Greenacre MJ: Interpreting multiple correspondence analysis. Appl Stochastic Model Data Anal 2006, 7: 195–210.
Fischer GH: Some Probabilistic Models for Measuring Change. In Advances in Psychological and Educational Measurement. Edited by: de Gruijter DNM, van der Kamp LJT. New York: Wiley; 1976.
Fischer GH: Some Latent Trait Models for Measuring Change in Qualitative Observations. In New Horizons in Testing. Latent Trait Theory and Computerized Adaptive Testing. Edited by: Weiss DJ. New York: Academic Press; 1983.
Fischer GH, Ponocny I: An extension of the patial credit model with an application to the measurement of change. Psychometrika 1994, 59: 177–192. 10.1007/BF02295182
Fischer G, Parzer P: An extension of the rating scale model with an application to the measurement of change. Psychometrika 1991, 56: 637–651. 10.1007/BF02294496
Rasch G: Probabilistic Models for Some Intelligence and Attainment Tests. The Danish Institute of Educational Research, Copenhagen: MESA Press; 1960.
Fischer GH: Linear Logistic Models for Change. Rasch Models. Foundations, Recent Developments and Applications. New York: Springer; 1995.
Van der Linden WJ, Hambleton RK: Handbook of Modern Item Response Theory. New York: Springer Verlag; 1997.
Mair P, Hatzinger R: CML based estimation of extended Rasch models with the eRm package in R. Psychol Sci Q 2007, 49: 26–43.
Mair P, Hatzinger R: Extended rasch modeling: the eRm package for the application of IRT models in R. J Stat Softw 2007, 20: 1–20.
Barclay-Goddard R, Epstein JD, Mayo NE: Response shift: a brief overview and proposed research priorities. Qual Life Res 2009, 18: 335–346. 10.1007/s11136-009-9450-x
Costello AB, Osborne JW: Best practices in exploratory factor analysis: four recommendations for getting the most from your analysis. Pract Asses Res Eval 2005, 10: 173–178.
Development Core Team R: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3–900051–07–0. URL http://www.R-project.org/
Lê S, Josse J, Husson F: FactoMineR: an R package for multivariate analysis. J Stat Softw 2008, 25: 1–18.
Hatzinger R, Rusch T: IRT models with relaxed assumptions in eRm: a manual-like instruction. Psychol Sci Q 2009, 51: 87–120.
Schwartz CE, Sprangers MA: Guidelines for improving the stringency of response shift research using the thentest. Qual Life Res 2010, 19: 455–464. 10.1007/s11136-010-9585-9
Schwartz CE, Rapkin BA: Understanding appraisal processes underlying the thentest: a mixed methods investigation. Qual Life Res 2012, 21: 381–388. 10.1007/s11136-011-0023-4
Visser MR, Oort FJ, Sprangers MA: Methods to detect response shift in quality of life data: a convergent validity study. Qual Life Res 2005, 14: 629–639. 10.1007/s11136-004-2577-x
Nolte S, Elsworth GR, Sinclair AJ, Osborne RH: Tests of measurement invariance failed to support the application of the “then-test”. J Clin Epidemiol 2009, 62: 1173–1180. 10.1016/j.jclinepi.2009.01.021
Ahmed S, Bourbeau J, Maltais F, Mansour A: The Oort structural equation modeling approach detected a response shift after a COPD self-management program not detected by the Schmitt technique. J Clin Epidemiol 2009, 62: 1165–1172. 10.1016/j.jclinepi.2009.03.015
Gandhi PK, Ried LD, Huang IC, Kimberlin CL, Kauf TL: Assessment of response shift using two structural equation modeling techniques. Qual Life Res 2013, 22: 461–471. 10.1007/s11136-012-0171-1
Donaldson GW: Structural equation models for quality of life response shifts: promises and pitfalls. Qual Life Res 2005, 14: 2345–2351. 10.1007/s11136-005-3977-2
King-Kallimanis BL, Oort FJ, Visser MR, Sprangers MA: Structural equation modeling of health-related quality-of-life data illustrates the measurement and conceptual perspectives on response shift. J Clin Epidemiol 2009, 62: 1157–1164. 10.1016/j.jclinepi.2009.04.004
Edelen MO, Reeve BB: Applying item response theory (IRT) modeling to questionnaire development, evaluation, and refinement. Qual Life Res 2007, 16(Suppl 1):5–18.
Floyd FJ, Widaman KF: Factor analysis in the development and refinement of clinical assessment instruments. Psychol Assess 1995, 7: 286–299.
Hambleton RK: Emergence of item response modeling in instrument development and data analysis. Med Care 2000, 38: II60-II65.
Lai JS, Crane PK, Cella D: Factor analysis techniques for assessing sufficient unidimensionality of cancer related fatigue. Qual Life Res 2006, 15: 1179–1190. 10.1007/s11136-006-0060-6
McLachlan SA, Devins GM, Goodwin PJ: Factor analysis of the psychosocial items of the EORTC QLQ-C30 in metastatic breast cancer patients participating in a psychosocial intervention study. Qual Life Res 1999, 8: 311–317. 10.1023/A:1008939818096
Reise SP, Widaman KF, Pugh RH: Confirmatory factor analysis and item response theory: two approaches for exploring measurement invariance. Psychol Bull 1993, 114: 552–566.
Smith AB, Wright P, Selby PJ, Velikova G: A rasch and factor analysis of the functional assessment of cancer therapy-general (FACT-G). Health Qual Life Outcomes 2007, 5: 19. 10.1186/1477-7525-5-19
Liu L, Hedeker D, Mermelstein R: Modeling nicotine dependence: an application of a longitudinal IRT model for the analysis of adolescent nicotine dependence syndrome scale. Nicotine Tob Res 2013, 15: 326–333. 10.1093/ntr/nts125
Blanchin M, Hardouin JB, Le Neel T, Kubis G, Blanchard C, Mirallié E, Sébille V: Comparaison of CTT and rasch-based approaches for the analysis of longitudinal patient reported outcomes. Stat Med 2011, 30: 825–838.
Swartz RJ, Schwartz C, Basch E, Cai L, Fairclough DL, McLeod L, Mendoza TR, Rapkin B: The king’s foot of patient-reported outcomes: current practices and new developments for the measurement of change. Qual Life Res 2011, 20: 1159–1167. 10.1007/s11136-011-9863-1
van Nispen RM, Knol DL, Neve HJ, van Rens GH: A multilevel item response theory model was investigated for longitudinal vision-related quality-of-life data. J Clin Epidemiol 2010, 63: 321–330. 10.1016/j.jclinepi.2009.06.012
Douglas JA: Item response models for longitudinal quality of life data in clinical trials. Stat Med 1999, 18: 2917–2931. 10.1002/(SICI)1097-0258(19991115)18:21<2917::AID-SIM204>3.0.CO;2-N
Glas CA, Geerlings H, van de Laar MA, Taal E: Analysis of longitudinal randomized clinical trials using item response models. Contemp Clin Trials 2009, 30: 158–170. 10.1016/j.cct.2008.12.003
Craig S: Implicit theories and beta change in longitudinal evaluations of training effectiveness: an investigation using item response theory. 2002.
Schwartz C, Martha J, Swaim B, Bode R, Kim D: Detecting response shift in using rasch analysis of then test data. Qual Life Res 2007., A-15(suppl):
Bonnetain F, Dahan L, Maillard E, Ychou M, Mitry E, Hammel P, Legoux JL, Rougier P, Bedenne L, Seitz JF: Time until definitive quality of life score deterioration as a means of longitudinal analysis for treatment trials in patients with metastatic pancreatic adenocarcinoma. Eur J Cancer 2010, 46: 2753–2762. 10.1016/j.ejca.2010.07.023
Hamidou Z, Dabakuyo TS, Mercier M, Fraisse J, Causeret S, Tixier H, Padeano MM, Loustalot C, Cuisenier J, Sauzedde JM, Smail M, Combier JP, Chevillote P, Rosbuger C, Arveux P, Bonnetain F: Time to deterioration in quality of life score as a modality of longitudinal analysis in patients with breast cancer. Oncologist 2011, 16: 1458–1468. 10.1634/theoncologist.2011-0085
We thank Holly Sandu for correcting the manuscript. This work was supported by a grant from the “Institut National du Cancer”. The study sponsor had no role in the conception, the design of the study, the data acquisition and analysis or in the manuscript preparation.
The authors declare that they have no competing interests.
AA performed the statistical analyses and interpretation and written the manuscript, CBM interpreted the data and drafted the manuscript, TC, FG, MV, DJ, MM, SC, JC, OG designed the study, SC, JC, OG included the patients, ZH interpreted the data, FB designed the study, written protocol, managed the statistical analyses, interpreted the data and review the draft. All authors read and approved the final manuscript.
About this article
Cite this article
Anota, A., Bascoul-Mollevi, C., Conroy, T. et al. Item response theory and factor analysis as a mean to characterize occurrence of response shift in a longitudinal quality of life study in breast cancer patients. Health Qual Life Outcomes 12, 32 (2014). https://doi.org/10.1186/1477-7525-12-32
- Health-related quality of life
- Longitudinal analysis
- Response shift
- Factor analysis
- Item Response Theory