In this paper we develop a novel procedure of variable screening for multivariate additive random-effects model, based on B-spline function approximations. With these approximations, the so-called signal-to-noise ratio (SNR) can be defined to inform importance each covariate in model. Then, SNR-based forward filtering is conducted covariates by using iterative projections multiple response data...