This article considers spline smoothing of variance functions. We focus on selection of smoothing parameters and develop three direct data-driven methods: unbiased risk (UBR), generalized approximate cross validation (GACV) and generalized maximum likelihood (GML). In addition to guaranteed convergence, simulations show that these direct methods perform better than existing indirect UBR, genera...