Variance function estimation in high-dimensions