Highlights

The dynamic IRT model provides a flexible solution to analyze repeated item responses measuring a latent construct over time

The dynamic IRT model relies on a mixed model to capture the continuoustime nature of the underlying construct

The dynamic IRT model can investigate the item psychometric properties and measurement invariance

Estimation of the dynamic IRT model is made available in the R package lcmm with a companion vignette

The case study describes the depressive symptomatology of patients with endstage renal disease on the transplant waiting list
1 Introduction
Item Response Theory (IRT) models, which exploit the information provided by each individual’s items responses, have received growing interest in health science for capturing latent constructs of interest such as depression, anxiety, fatigue, quality of life, or cognitive functioning Cerou et al. (2019); Gorter et al. (2015); McCall et al. (2021); Abdelhamid et al. (2021). IRT models have interesting properties compared to models coming from classical measurement theory such as Classical Test Theory (CTT) models which aggregate the items into a global score or score per domain. In particular, CTT produces ordinal measurements while IRT generates interval measurements. Hence, with IRT, a unit difference characterizes the same amount when measured from different initial levels on the latent construct scale. IRT also allows for a finer granularity of the level of analysis, done at the item level, which enables a better understanding of the item psychometric properties, and an in depth description of patients’ experience.
In health research, the interest often lies in the longitudinal changes of latent constructs. Examples include the trajectory of anxiety or fatigue in clinical research Rakers et al. (2021); Otto et al. (2021) or the trajectory of functional dependency in epidemiological research on aging Edjolo et al. (2016) based on the observed patients’ repeated responses to questionnaires either selfreported (and named PatientReported Outcomes (PRO)) or reported by the clinician (and named ClinicianReported Outcomes (CRO)). IRT models have been extended to account for repeated item measurements but most of the time, the measurement occasions are necessarily considered as occurring at the exact same time for all patients Cai and Houts (2021)
. This is, however, rarely the case in practice. For instance, in cohort studies, even though visits are planned, the exact timing may differ substantially across individuals. The timescale may also be the time from a specific health event (e.g. diagnosis, registration on a waiting list), independent from the time of study inclusion, so that the timescale becomes per se continuous. This is the case in aging studies where trajectories are assessed according to age, or more generally in all the situations where entry in the study does not correspond to a clearly defined time zero. Models from linear mixed model theory
Laird and Ware (1982) are particularly suited for the analysis of outcomes repeatedly measured over time. They can model the outcome trajectory in continuous time while accounting for the withinsubject correlation. We describe in this paper how IRT modeling can be combined with linear mixed model theory for the analysis of item responses measured repeatedly over time when observation times vary across individuals. We show how to operationalize the latent construct as a latent process defined in continuous time and how to link it with the observed items responses through a measurement model for graded response, the Graded Response Model (GRM) Samejima (1997), at each individual and occasionspecific observation time. Beyond the description of a latent construct trajectory over time from repeated item measurements, this dynamic IRT model may also help to assess the item and scale properties as done with crosssectional IRT methods, and investigate lack of measurement invariance, that is Differential Item Functioning (DIF) between groups, or Response Shift (RS) over more than two measurement times.The proposed approach is contextualized in a clinical example in endstage renal disease, the PREDIALA study Sébille et al. (2016); AuneauEnjalbert et al. (2020). The PREDIALA study aims at studying the experience of patients with endstage renal disease (e.g. quality of life, anxious and depressive symptoms) on the renal transplant waiting list. From a psychological perspective, the waiting list period can be long and anxietyprovoking for patients because of the uncertainty of waiting, the hope of being called for a transplant and the disappointment, and sometimes distress, of not being called. In addition, the onset or worsening of depressive/anxious symptoms can occur as the time on the waiting list increases Tong et al. (2015). Moreover, patients’ experience may also differ according to their renal replacement therapy, that is, whether they are dialyzed or not (i.e. preemptive) AuneauEnjalbert et al. (2020). As patients with chronic diseases have to live and adapt to their illness, its stability or progression over time, they may also understand or interpret the items of the questionnaires differently according to sociodemographic or clinical characteristics (e.g. type of renal replacement therapy) or over time, despite having similar health outcomes. The former situation may induce DIF Holland and Wainer (2009), while the latter may produce RS Sprangers and Schwartz (1999). From a methodological perspective, the relevant timescale was the time elapsed since registration on the waiting list. Since entry in the study occurred after different periods of time spent on the waiting list across individuals, this timescale substantially differed from the classical followup time. It induces large interindividual variations of measurement times between patients and calls for a dynamic IRT model that can handle individualspecific measurement times. Such models can also enable the investigation of measurement invariance between groups (DIF) or over time (RS).
In a sample of subjects, let consider a set of items belonging to the same scale that are measured repeatedly over time, with designating the response to item () for subject () at repeated occasion (). Note that here, we consider the general framework where the number of measurements may differ from one individual to the other (and possibly from one item to the other) so that item measurement is to be associated with its actual time of measurement .
As in classical IRT methodology, we assume that the items measure the same underlying construct called . The major difference is that in a longitudinal IRT setting, the construct is now a latent process defined in continuous time.
Its trajectory over time can be described by a linear mixed model (LMM) to account for the within individual correlation (and betweenindividual variability):
(1) 
where
is a vector of variables including functions of time
, associated with parameters , which describes the shape of trajectory over time of the construct of interest at the population level and its association with covariates; is a vector of variables that almost always includes exclusively functions of time and is associated with the vector of individual randomeffects with ( is usually left unstructured). This second part models the individual departure from the mean trajectory . In some contexts, a Gaussian stochastic process may be added to better reflect the local variations in the individual trajectories. For sake of readability, we do not cover this aspect in the present paper.Identifiability constraint is added to this model to determine the dimension of the latent process. Usually, the intercept is removed from the model (i.e.,
does not include any intercept) so that the mean of the latent process is 0 in the category of reference, and the variance of the first randomeffect is constrained to 1. This first random effect being usually an intercept, this corresponds to assuming that the conditional variance given the covariates is 1 at time
.The latent process is linked to the observations of the items using an itemspecific measurement model. In this work, we focus mainly on binary and ordinal items even though the methodology could also apply to continuous items (see ProustLima et al. ProustLima et al. (2013) for more details). We assume that item is defined with ordinal levels from 0 to .
(2) 
where
is the Gaussian cumulative distribution function,
is a parameter defining the discrimination of item and are the location parameters which correspond to the thresholds defining the change in the successive levels of item . For an ordinal item, we assume .Equation (2) comes from the idea that item takes the level if the underlying construct plus a measurement error of variance lies in the interval as shown below:
(3) 
By considering that the measurement error is Gaussian, it induces that
(4) 
and Equation (2) comes by denoting .
Note that, by considering logistic errors instead of Gaussian errors, one would obtain the logistic ogive model, also very popular in IRT methodology.
Equation (2) defines a model for graded response, also known as Graded Response Model (GRM) Samejima (1997); Baker and Kim (2004). We focus on this type of model in the remaining of the work but acknowledge that any alternative measurement model could be considered instead depending on the distributional assumption and the items type (e.g., binary, continuous). See for instance Saulnier et al. in the same special [REF] issue and Barbieri et al. Barbieri et al. (2017).
The dynamic IRT model defined with equations (1) and (2) assumes that all the items have the same functioning, meaning that the common underlying construct captures all the information of the items and what remains specific to the item is only its location and discrimination/error; this is made clear with equation (3). However, sometimes a different functioning of the items may be suspected according to a covariate at a given time or over time. In IRT methodology, this is called differential item functioning (DIF) or Response Shift (RS) when lack of measurement invariance occurs over time. DIF and RS can also be investigated in the dynamic IRT by completing the measurement model as follows:
(5) 
where is the vector of covariates for which a DIF is suspected and the associated parameters. When is also part of the structural model in Equation (1) (i.e, ), an identifiability constraint is added to the parameters making them correspond to contrasts (i.e., deviations to the mean effect) with . The total effect of becomes the sum of its common effect on the latent process (part of ) and its itemspecific contrast .
With longitudinal data, two sorts of itemspecific functionings may be investigated:

classical DIF with including timeindependent covariates only; this explores how differently level parameters of a specific item differ according to individual characteristics;

item response shift with including functions of time ; this explores how level parameters of a specific item change over time.
We consider here a maximum likelihood framework for the estimation of the dynamic IRT model.
Let denote the total vector of parameters defined in the structural part of the model described in (1) and in the K measurement equations described in (2). This vector includes:

the fixed effects (except the intercept for identifiability).

the parameters specifying the VarianceCovariance matrix of the randomeffects. To ensure that is positive definite, we consider the parameters of the Cholesky upper triangular transformation C (i.e., ) with first element fixed at 1 for identifiability.

the discrimination parameters. We consider parameters so that the item discriminations .

the vector of item locations of each item . To account for the constraint that , we consider the vector so that and for .

In case of differential effect of covariates on items (Equation (5)), the vector of parameters with .
Let denote all the repeated item information of subject . The contribution of subject to the likelihood is
(6) 
where the integration over the distribution of the pvector of random effects is obtained by Quasi MonteCarlo approximation following proposals of Philipson et al. Philipson et al. (2020). We systematically considered 1000 points in this work.
The maximum likelihood estimators of are obtained by maximizing the loglikelihood . This is achieved with the MarquardtLevenberg algorithm, a robust Newtonlike algorithm, with stringent convergence criteria on the parameters, the loglikelihood and the first and second derivatives of the loglikelihood (see Philipps et al. Philipps et al. (2021) for details).
The maximum likelihood estimates are denoted and their variance, obtained by the inverse of the Hessian, is denoted .
The dynamic IRT model can be estimated with the multlcmm function of R package lcmm ProustLima et al. (2017). A package vignette provides a tutorial that fully describes the present dynamic IRT model estimation and posterior computations on a simulated dataset that mimics the PREDIALA data.
We can compute several posterior quantities from the estimates
, and confidence intervals around these quantities can be obtained by approximating the posterior distribution by MonteCarlo simulations using the asymptotic distribution of the parameters
. We consider for this 2000 random draws. In the following, we omit or more generally in the equations for better readability.Predicted trajectories of the construct can be computed either at the population level (i.e., marginally to the randomeffects) or at the individual level (i.e., conditionally to the individual randomeffects). The predicted trajectory at the population level is computed for a profile of covariates :
The predicted trajectory at the individual level is computed given the individual covariates , and all the information on the items :
where the expected randomeffect is approximated by the mode of the posterior distribution .
With binary items, the Item Characteristic Curve (ICC) describes the probability of the highest item level according to the underlying construct level. With ordinal items, ICC translates into two curves:

the item category probability curve also known as category characteristic curve which describes the probability of a response in a given item category according to the underlying construct level:
(7) 
The item score expectation curve which can also be computed as a function of the underlying construct level:
(8) 
These two curves allow representing the items and their properties. The items locations describe where the items function along the construct level while the steepness of the item score expectation characterizes the items discriminations. For example, the steeper the curve for an item, the better it can discriminate between two different construct levels.
The predicted trajectory over time of each item can be computed for a profile of covariates as follows:
(9) 
where is computed as in equation (8) with , and the integral over the randomeffects distribution is obtained by Quasi MonteCarlo approximation.
The Fisher information provides a quantification of the level of information brought by each item, and each item level. It is computed using the second derivatives of the item level probability denoted for item and level . The item information function for category is defined as follows (calculations are detailed in Section 1 of the supplementary material):
(10) 
The information curve provides the summary at the item level as follows:
(11) 
We applied the dynamic IRT model to analyze the repeated measures of depressive symptomatology in the PREDIALA study. The HADS (Hospital Anxiety and Depression scale) was used to measure anxiety and depression disorders Zigmond and Snaith (1983). The HADS consists of 14 items on a 4point Likert scale, seven of which are related to anxiety symptoms and seven to depression symptoms. Only the depression symptom domain is presented here. The 7 items rated from 0 (total agreement) to 3 (total disagreement) are as follows:

Item 2 “I still enjoy the things I used to enjoy” (Enjoy)

Item 4 “I can laugh and see the funny side of things” (Laugh)

Item 6 “I feel cheerful” (Cheerful)

Item 8 “I feel I am slowed down” (Slow)

Item 10 “I have lost interest in my appearance” (Appearance)

Item 12 “I look forward with enjoyment to things” (Looking forward)

Item 14 “I can enjoy a good book or radio or TV programme” (Leisure)
Responses to items 8 and 10 are reversed so that higher levels systematically indicate more intense symptoms.
Our objective was to describe the trajectory of depressive symptomatology over time from registration on the waiting list for a renal transplant, and to describe the possible differences according to patients’ renal replacement therapy at inclusion in the study, i.e. patients either dialyzed or not (preemptive). Indeed, being on the transplant waiting list may be experienced differently between dialyzed and preemptive patients. It can be hypothesized that, for example, the depressive symptoms experienced by dialyzed patients are more pronounced as compared to preemptive patients due to their experience with dialysis and expectations of associated complications. It is therefore possible that the need for clinical and psychological support is not the same for all patients.
A secondary objective was to assess whether the functioning of some items differed according to the group (preemptive or dialyzed) or shifted with time on the waiting list. Indeed, patients may perceive the items differently according to their renal replacement therapy and over time, despite having similar depression levels.
We included in the analysis all the patients from PREDIALA who entered the study within 48 months following their registration on the waiting list. They were either in the dialyzed or preemptive group at entry, and had at least one measure for each of the 7 items of the HADS before the end of the study. The end of the study was defined by either a switch in group (from preemptive to dialyzed status), a clinical event (mainly transplantation) or the administrative censoring. From the initial 577 patients included in PREDIALA study, this selection lead to a final sample of 561 patients and 1136 repeated visits. Among them 356 (63.5% ) were men and 288 (51.3%) were under dialysis. The median age at entry was 59 years (range 1967 years), and the patients had been on the waiting list for a very variable time ranging from 0.1 to 43.1 months (median 5.1 months) at entry in the cohort. This leads to substantial variability in measurement timings across patients at entry and during followup as shown in Figure 1. This continuous distribution of the measurement times which would have been ignored using standard IRT methods is naturally handled in the dynamic IRT model thanks to the definition of the underlying construct as a latent process in continuous time.
The dynamic model was defined following equation (1) for the trajectory of underlying depressive symptomatology and equation (2) for the 7 itemspecific measurement model. As we did not have any assumption regarding the shape of trajectory over time of depressive symptomatology, we used a basis of natural cubic splines with 2 internal knots placed at tertiles of the measurement time distribution, that is 7 and 15 months, and boundary knots placed at 0 and 60 months. Each of the four functions of time (intercept and 3 splines functions) was associated with fixed effects specific to the group (Preemptive/Dialyzed) to assess the mean trajectories, and individual correlated random effects to account for the correlation within repeated measures of each individual. For the measurement models, we assumed in the main analysis that all items functioned similarly, i.e., no DIF and no response shift occurred. We then explored in secondary analyses whether some items functioned differently by group (adding an itemspecific contrast on group), and whether some items were affected by response shift over time (adding an itemspecific contrast of the 3 time functions). The parameter estimates of these three models are provided in Table S1 of the supplementary material.
In the following description, 1 unit of depressive symptomatology, called 1 SD, corresponds to the interindividual variability at registration in the dialyzed group.
The mean predicted trajectory over time of depressive symptomatology, displayed in Figure 2, varied over time and according to the group. In the preemptive group, the level of depressive symptomatology increased during the first year on the waiting list by 0.243 (0.012,0.498) SD and then remained stable. The level of depressive symptomatology was higher in the dialyzed group compared to the preemptive group at the time of registration (difference of 0.482 (0.814,0.149) SD). It then slightly decreased during the first year to reach a similar level as in the preemptive group, and then increased again after approximately 2 years on the waiting list by a mean annual rate of 0.245 (0.053,0.438) SD (computed from 2 to 6 years).
We exploited the dynamic IRT model to assess the HADS depressive symptomatology items characteristics. To help appreciating the item characteristics, we assessed the range of the distribution of the underlying depressive symptomatology based on the estimates of the dynamic IRT model and an hypothetical population of 100000 preemptive patients and 100000 dialyzed patients with measures every month from registration up to 72 months. The resulting 95% prediction interval of the underlying depressive symptomatology was [6.10,5.90] with the 10% and 90% percentiles of the distribution at 3.00 and 2.74, respectively.
Table 1 provides the estimated locations and discrimination while Figure 3 shows the curves of item expectations (top) and curves of item information (bottom) according to the underlying depressive symptomatology. Figure S1 and S2 of the supplementary material further display for each item category the probability curve and the information function, respectively.
The easiest items, in terms of their difficulty parameters, were Item 8 (Slow) and Item 2 (Enjoy) while the most difficult one was Item 14 (Leisure). This means that the level of depression required to respond to the most unfavorable response categories (i.e. indicative of higher depressive symptoms) of items 2 and 8 (e.g., “Sometimes” for item 8 “I feel as if I am slowed down”) was lower than the level of depression required to respond to the most unfavorable response categories of item 14 (e.g., “Not often” for item 14 “I can enjoy a good book or radio or TV program”). The most discriminant items with the steepest curves, representing their ability to discriminate patients with different levels of depression, were items 2, 4 and 12 concerning the ability to enjoy, laugh and look forward, respectively. The estimated curve of the Fisher information plotted in Figure 3 (bottom) also underlines the major role of these 3 items compared to the others. In contrast, item 14 about leisure does not bring much information in this population as it seems to measure much higher levels of depression than the other items.
Item  category 0  1  category 1  2  category 2  3  discrimination  

est  SE  est  SE  est  SE  est  SE  
2  Enjoy  0.46  0.13  0.77  0.14  1.52  0.18  1.29  0.13 
4  Laugh  0.26  0.12  0.74  0.13  1.91  0.21  1.56  0.16 
6  Cheerful  0.48  0.14  1.58  0.20  3.34  0.38  0.85  0.09 
8  Slow  1.51  0.19  0.40  0.13  1.69  0.20  0.95  0.10 
10  Appearance  0.05  0.12  1.00  0.16  2.27  0.26  0.88  0.10 
12  Looking Forward  0.32  0.13  0.72  0.13  1.82  0.20  1.46  0.15 
14  Leisure  0.83  0.17  3.18  0.42  4.11  0.54  0.56  0.07 
Estimate (and associated standard error (SE) obtained by
method) of locations and discrimination of the 7 items of HADS measuring Depressive symptomatologyWe first explored any differential item functioning on the group (dialyzed versus preemptive). Estimates are provided in Table 2 along with those of the main model that ignores DIF. Overall, the Chisquare test assessing simultaneously the 6 contrasts on group did not reject the null “no DIF on group” assumption (p=0.266). However, taken individually, the difference between groups for item 2 (Enjoy) was significantly larger than for the other items (itemspecific effect of preemptive group on the underlying level estimated at 0.150 (0.269,0.0314) SD). This suggests that this item is more difficult for preemptive patients (location parameter shifted to +0.150) at the same underlying level of depressive symptomatology; preemptive patients tend to respond more readily to more favorable response categories than patients under dialysis, despite having similar depression levels. In addition, accounting or not for DIF impacted the conclusions: the group effect on the underlying depressive symptomatology was not significant anymore when accounting for DIF suggesting that the difference between groups in the model without DIF was mainly carried by item 2.
coef  SE  p  coef  SE  p  
intercept  0.000      0.000     
ns1  0.305  0.172  0.075  0.304  0.229  0.184 
ns2  0.039  0.222  0.862  0.038  0.590  0.949 
ns3  0.538  0.250  0.031  0.538  0.310  0.083 
group  0.482  0.170  0.004  0.458  0.413  0.267 
ns1:preemptive  0.477  0.223  0.032  0.473  0.265  0.075 
ns2:preemptive  0.428  0.335  0.201  0.429  0.491  0.382 
ns3:preemptive  0.390  0.292  0.182  0.391  0.338  0.247 
Contrasts on preemptive:  (global p=0.266)  
Item 2  0.150  0.061  0.013  
Item 4  0.031  0.059  0.600  
Item 6  0.020  0.077  0.800  
Item 8  0.053  0.100  0.593  
Item 10  0.040  0.158  0.802  
Item 12  0.037  0.061  0.543  
** Item 14  0.138  0.124  0.266 
** coefficient not estimated but obtained as minus the sum of the others.
We secondly explored any item responses shift over time by adding itemspecific contrasts on the 3 natural cubic splines functions of time (ns1, ns2, ns3). The resulting predicted trajectories of each item in this model compared to the model assuming no response shift are given in Figure 4. No clear response shift was identified here although the overall test assessing simultaneously ns1, ns2 and ns3 contrasts in the model with RS suggested some potential for a different behavior of item 2 over time compared to the others (p=0.073). Each item behavior over time was close to the one of the underlying construct showed in Figure 2 with some slight differences in the intensity of change in the first months after entry on the waiting list. Note that these slight differences do not reflect directly the statistical test which only focused on the model assuming RS.
We have described how to combine the item response theory with the linear mixed model theory for itemlevel analysis of longitudinal PRO or CRO data when measurement times may vary across individuals. Using a real case example with the PREDIALA study, we have shown that this dynamic IRT model can describe the latent construct trajectory over time and its determinants, while simultaneously assessing the item and scale properties, and exploring lack of measurement invariance between groups (DIF) and over time (response shift).
This analysis of PREDIALA data helped to better understand the experience of patients with endstage renal disease on the renal transplant waiting list in terms of depressive symptoms. DIF was highlighted on item 2 (Enjoy) indicating that patients under dialysis had more difficulty in reporting having enjoyment than preemptive patients, despite having similar levels of depression. This may reveal that the need for clinical and psychological support may not be the same for all patients, according to their renal replacement therapy. Response shift was not significantly evidenced despite a trend on this same item 2. Adjusting for DIF and response shift, the level of depressive symptomatology of preemptive patients tended to slightly increase during the first year and to remain stable afterwards. The level of depressive symptomatology of patients under dialysis tended to be close to the one of preemptive patients after 2 years on the waiting list and to increase afterwards. Although it has been reported that the time on waiting list should be reduced to limit depressive symptoms Corruble et al. (2010) and improve healthrelated quality of life Ong et al. (2013), the shortage of grafts unfortunately often makes this difficult to achieve.
The dynamic IRT model we described here unites the strengths of IRT and LMM theories. On the one hand, the use of a structural mixed model makes it possible to operationalize the latent construct as a latent process defined in continuous time and thus takes into account that most health phenomena intrinsically evolve in continuous time. On the other hand, the use of IRT methodology to define the measurement scale at each individual and occasionspecific observation time enables a precise modeling of the items constituting the measurement scale, and their properties.
However, as for all methodologies, this method is not without some limitations and some further path of research may be put forward. First, we focused here on a specific measurement model, the GRM, which translates the discretization of the underlying latent process into ordinal categories as shown in Equation (3) Commenges et al. (2015). However, coming from IRT models, this measurement model does not possess the specific objectivity property as Rasch Measurement Theory (RMT) models do Andrich (2011). It would thus be interesting to adapt the proposed methodology to RMT measurement models (see Blanchin et al. [REF] in this special issue and Barbieri et al. Barbieri et al. (2017)). Of note, changing the measurement model does not impact the estimation procedure nor the structural part of the model, and the current version of the program already handles different measurement models for continuous items in addition to GRM for ordinal items.
Second, by relying on the mixed model theory and the maximum likelihood estimation, the dynamic IRT model relies on the missing at random assumption both for monotonic and for intermittent missingness Little (1995). In the presence of informative dropout, a joint model for the repeated item responses and the time to dropout, as described by Saulnier et al. in this special issue (REF), should be favored. This joint model combines the dynamic IRT model with a survival model that captures the association between the underlying latent construct and the dropout (or any other event of interest). Third, the methodology is fully parametric, so that analytic choices are systematically made. For instance, to simplify the application setting, we only included timeindependent covariates even though timedependent covariates could also be considered, such as a changing of group during the followup. We also globally tested the lack of measurement invariance over the items’ parameters across the overall followup although a more precise assessment could be done regarding which function of time is affected or at which time the lack of measurement invariance occurs.
To conclude, by extending the IRT methodology to longitudinal data, and considering the time as continuous, our methodology provides a versatile and flexible approach for modeling item responses measured repeatedly over time as encountered in numerous longitudinal health studies.
Acknowledgements The authors gratefully thank the coinvestigators of the study: Magali Giral, Aurélie Meurette, Emmanuel Morelon, and Laetitia Albano. The authors express sincere thanks to Elodie FaurelPaul and Astrid Fleury for their assistance for the study, as well as all the participants for their contribution to the study. The authors wish to thank members of the clinical research assistant team and DIVAT Consortium Collaborators (Medical Doctors, Surgeons, HLA Biologists). Nantes: Gilles Blancho, Julien Branchereau, Diego Cantarovich, Agnès Chapelet, Jacques Dantal, Clément Deltombe, Lucile Figueres, Claire Garandeau, Magali Giral, Caroline GourraudVercel, Maryvonne Hourmant, Georges Karam, Clarisse Kerleau, Aurélie Meurette, Simon Ville, Christine Kandell, Anne Moreau, Karine Renaudin, Anne Cesbron, Florent Delbos, Alexandre Walencik, Anne Devis; Lyon E. Hériot : Lionel Badet, Maria Brunet, Fanny Buron, Rémi Cahen, Sameh Daoud, Coralie Fournie, Arnaud Grégoire, Alice Koenig, Charlène Lévi, Emmanuel Morelon, Claire PouteilNoble, Thomas Rimmelé, Olivier Thaunat.
Funding This work was funded by the French National Research Agency (Project DyMES  ANR18C36000401) and the French Ministry of Health (PHRC130224, 2013).
Ethical approval The PREDIALA study is part of the PreKitQoL study which is registered on the ClinicalTrials.gov Registry (RC14_0078, NCT02154815). It has obtained approval from the ethical Committee for Persons’ Protection (CPP, Tours, 2014S8), and from the advisory committee on research data and information in health (CCTIRS, Paris, 14.314).
Conflict of interest All authors have no competing interest.
Supplementary Material Supplementary material can be found on the journal website.
Data availability The software is openly available in the R package lcmm (on https://github.com/CecileProustLima/lcmm and on cran). The R script along with a dataset mimicking the PREDIALA study are also provided as a package vignette. The raw PREDIALA data can not be shared.
2 Methods
3 Application to PREDIALA study
4 Conclusions
References
 Assessing cognitive abilities using the waisiv: an item response theory approach. International Journal of Environmental Research and Public Health 18 (13), pp. 6835. Note: PMID: 34202249 PMCID: PMC8297006 External Links: Document Cited by: §1.
 Rating scales and rasch measurement. Expert Review of Pharmacoeconomics & Outcomes Research 11 (5), pp. 571–585. Note: PMID: 21958102 External Links: Document Cited by: §1.
 Comparison of longitudinal quality of life outcomes in preemptive and dialyzed patients on waiting list for kidney transplantation. Quality of Life Research: An International Journal of Quality of Life Aspects of Treatment, Care and Rehabilitation 29 (4), pp. 959–970 (eng). External Links: ISSN 15732649, Document Cited by: §1.
 Item Response Theory. Parameter Estimation Techniques. 2nd edition, Statistics: Textbooks & Monographs, Marcel Dekker, New York. Note: 20/03/2009ID  3773LB  DIM2585 / C. ProustLima (Biostat) External Links: ISBN 9780824758257 Cited by: §1.
 Item response models for the longitudinal analysis of healthrelated quality of life in cancer clinical trials. BMC Medical Research Methodology 17 (1), pp. 148. External Links: ISSN 14712288, Link, Document Cited by: §1, §1.
 Longitudinal Analysis of PatientReported Outcomes in Clinical Trials: Applications of Multilevel and Multidimensional Item Response Theory. Psychometrika 86 (3), pp. 754–777 (eng). External Links: ISSN 18600980, Document Cited by: §1.
 Application of item response theory to model disease progression and agomelatine effect in patients with major depressive disorder. The AAPS journal 22 (1), pp. 4. Note: PMID: 31720897 External Links: Document Cited by: §1.
 Chapter 5. Extensions of mixed models. In Dynamical Biostastical Models, (en). Cited by: §1.
 Progressive Increase of Anxiety and Depression in Patients Waiting for a Kidney Transplantation. Behavioral Medicine 36 (1), pp. 32–36. Note: Publisher: Taylor & Francis _eprint: https://doi.org/10.1080/08964280903521339 External Links: ISSN 08964289, Link, Document Cited by: §1.
 Natural History of Dependency in the Elderly: A 24Year PopulationBased Study Using a Longitudinal Item Response Theory Model. American Journal of Epidemiology 183 (4), pp. 277–285 (eng). Note: Number: 4 External Links: ISSN 14766256, Document Cited by: §1.
 Why item response theory should be used for longitudinal questionnaire data analysis in medical research. BMC Medical Research Methodology 15 (1), pp. 55. External Links: ISSN 14712288, Link, Document Cited by: §1.
 Differential Item Functioning. Routledge, NewYork, NY. Cited by: §1.
 Randomeffects models for longitudinal data. Biometrics 38 (4), pp. 963–74 (eng). Note: ET  1982/12/01ID  31Laird, N M Ware, J H United states Biometrics Biometrics. 1982 Dec;38(4):96374. External Links: ISSN 0006341X (Print) 0006341X (Linking) Cited by: §1.
 Modeling the DropOut Mechanism in RepeatedMeasures Studies. Journal of the American Statistical Association 90 (431), pp. 1112–1121. External Links: ISSN 01621459, Link, Document Cited by: §1.
 Examining suicide assessment measures for research use: using item response theory to optimize psychometric assessment for research on suicidal ideation in major depressive disorder. Suicide & LifeThreatening Behavior. Note: PMID: 34237156 External Links: Document Cited by: §1.
 What factors really matter? Healthrelated quality of life for patients on kidney transplant waiting list. Annals of the Academy of Medicine, Singapore 42 (12), pp. 657–666 (eng). External Links: ISSN 03044602 Cited by: §1.
 Illness representations, coping and anxiety among men with localized prostate cancer over an 18months period: A parallel vs. levelcontrast mediation approach. PsychoOncology (eng). External Links: ISSN 10991611, Document Cited by: §1.
 Robust and efficient optimization using a MarquardtLevenberg algorithm with R package marqLevAlg. R journal in press. Cited by: §1.
 Faster Monte Carlo estimation of joint models for timetoevent and multivariate longitudinal data. Computational Statistics & Data Analysis 151, pp. 107010 (en). External Links: ISSN 01679473, Link, Document Cited by: §1.
 Analysis of multivariate mixed longitudinal data: a flexible latent process approach. The British journal of mathematical and statistical psychology 66 (3), pp. 470–487 (eng). External Links: ISSN 20448317, Document Cited by: §1.
 Estimation of Extended Mixed Models Using Latent Classes and Latent Processes: The R Package lcmm. Journal of Statistical Software, Articles 78 (2), pp. 1–56. External Links: ISSN 15487660, Link, Document Cited by: §1.
 Trajectories of Fatigue, Psychological Distress, and Coping Styles After Mild Traumatic Brain Injury: A 6Month Prospective Cohort Study. Archives of Physical Medicine and Rehabilitation, pp. S0003–9993(21)00462–7 (eng). External Links: ISSN 1532821X, Document Cited by: §1.
 Graded Response Model. In Handbook of Modern Item Response Theory, W. J. van der Linden and R. K. Hambleton (Eds.), pp. 85–100 (en). External Links: ISBN 9781475726916, Link, Document Cited by: §1, §1.
 Prospective, multicenter, controlled study of quality of life, psychological adjustment process and medical outcomes of patients receiving a preemptive kidney transplant compared to a similar population of recipients after a dialysis period of less than three years–The PreKitQoL study protocol. BMC nephrology 17, pp. 11 (eng). External Links: ISSN 14712369, Document Cited by: §1.
 Integrating response shift into healthrelated quality of life research: a theoretical model. Social Science & Medicine (1982) 48 (11), pp. 1507–1515 (eng). External Links: ISSN 02779536, Document Cited by: §1.
 ’Suspended in a paradox’patient attitudes to waitlisting for kidney transplantation: systematic review and thematic synthesis of qualitative studies. Transplant International: Official Journal of the European Society for Organ Transplantation 28 (7), pp. 771–787 (eng). External Links: ISSN 14322277, Document Cited by: §1.
 The hospital anxiety and depression scale. Acta Psychiatrica Scandinavica 67 (6), pp. 361–370 (eng). External Links: ISSN 0001690X, Document Cited by: §1.
Comments
There are no comments yet.