The linearization coefficient $\mathcal{L}(L_{n_1}(x)\dots L_{n_k}(x))$ of classical Laguerre polynomials $L_n(x)$ is known to be equal to the number of $(n_1,\dots,n_k)$-derangements, which are permutations with a certain condition. Kasraoui, Stanton and Zeng found a $q$-analog of this result using $q$-Laguerre polynomials with two parameters $q$ and $y$. Their formula expresses the linearization coefficient of $q$-Laguerre polynomials as the generating function for $(n_1,\dots,n_k)$-derangements with two statistics counting weak excedances and crossings. In this paper their result is proved by constructing a sign-reversing involution on marked perfect matchings.