On estimation of covariate-specific residual time quantiles under the proportional hazards model