P-splines quantile regression estimation in varying coefficient models