Using multivariable ordinary least squares regression to control for hospital-level characteristics , we then estimated the association between the presence of hospitalists and predicted hospital mortality and readmission .