The paper studies estimation of partially linear hazard regression models with varying coefficients for multivariate survival data. A profile pseudopartiallikelihood estimation method is proposed. The estimation of the parameters of the linear part is accomplished via maximization of the profile pseudopartiallikelihood, whereas the varyingcoefficient functions are considered as nuisance parameters that are profiled out of the likelihood. It is shown that the estimators of the parameters are root <i>n</i> consistent and the estimators of the nonparametric coefficient functions achieve optimal convergence rates. Asymptotic normality is obtained for the estimators of the finite parameters and varyingcoefficient functions. Consistent estimators of the asymptotic variances are derived and empirically tested, which facilitate inference for the model. We prove that the varyingcoefficient functions can be estimated as well as if the