I would need to use a quite large sparse matrix as panel structure (basically the element s in y ~ x1 + x2 + ... | s ) into the rqpd function https://rdrr.io/r