A Dynamic Factor Model for Multivariate Counting Process Data

We propose a dynamic multiplicative factor model for process data, which arise from complex problem-solving items, an emerging testing mode in large-scale educational assessment. The proposed model can be viewed as an extension of the classical frailty models developed in survival analysis for multivariate recurrent event times, but with two important distinctions: (i) the factor (frailty) is of primary interest; (ii) covariates are internal and embedded in the factor. It allows us to explore low dimensional structure with meaningful interpretation. We show that the proposed model is identifiable and that the maximum likelihood estimators are consistent and asymptotically normal. Furthermore, to obtain a parsimonious model and to improve interpretation of parameters therein, variable selection and estimation for both fixed and random effects are developed through suitable penalisation. The computation is carried out by a stochastic EM combined with the Metropolis algorithm and the coordinate descent algorithm. Simulation studies demonstrate that the proposed approach provides an effective recovery of the true structure. The proposed method is applied to analysing the log-file of an item from the Programme for the International Assessment of Adult Competencies (PIAAC), where meaningful relationships are discovered.