Fine-grained Analysis of Stability and Generalization for Stochastic Bilevel Optimization
Authors: Xuelin Zhang, Hong Chen, Bin Gu, Tieliang Gong, Feng Zheng
Organizations: College of Informatics, Huazhong Agricultural University, Wuhan 430070, China · Engineering Research Center of Intelligent Technology for Agriculture, Ministry of Education, Wuhan 430070, China · School of Artificial Intelligence, Jilin University, Jilin 130012, China · School of Computer Science and Technology, Xi’an Jiaotong University, Xi’an 710049, China · Department of Computer Science and Engineering, Southern University of Science and Technology, Shenzhen 518055, China
Stochastic bilevel optimization (SBO) has been integrated into many machine learning paradigms recently, including hyperparameter optimization, meta learning, and reinforcement learning. Along with the wide range of applications, there have been numerous studies on the computational behavior of SBO. However, the generalization guarantees of SBO methods are far less understood from the lens of statistical learning theory. In this paper, we provide a systematic generalization analysis of the first-order gradient-based bilevel optimization methods. Firstly, we establish the quantitative connections between the on-average argument stability and the generalization gap of SBO methods. Then, we derive the upper bounds of on-average argument stability for single-timescale stochastic gradient descent (SGD) and two-timescale SGD, where three settings (nonconvex-nonconvex (NC-NC), convex-convex (C-C), and strongly-convex-strongly-convex (SC-SC)) are considered respectively. Experimental analysis validates our theoretical findings. Compared with the previous algorithmic stability analysis, our results do not require reinitializing the inner-level parameters at each iteration and are applicable to more general objective functions.
Figures & tables
Algorithms
Stability
SC-SC
C-C
NC-NC
SSGD
l1
O(m1K)
O(m1KC4ln(K))
—
(Theorem 2 )
l2
O(m12(m1+K)K)
O(m12(m1+K)K2C4−1ln2(K))
—
TSGD
l1
O(m1KTC5)
O(m1K2KTC6ln(T))
O(m12KKC2T1+C3T)
(Theorem 3 )
l2
O(m12(m1+K)KT2C5)
O(m12K2(m1+K)2KT2C6ln2(T))
O(m12(m1+K)2KK2C2T1+C3T2)
SSGD
l1
O(m11)
O(m11)
O(m1KC1)
(Proposition 1)
l2
O(m12Km1+K)
O(m12Km1+K)
O(m12(m1+K)K2C1)
Table 1: Summary of the generalization bounds under different settings. For briefly, l1 ( l2 ) represents the l1 ( l2 ) on-average argument stability and C1−C6 are positive constants. m1 is the number of validation samples; K and T are the total numbers of outer and inner iterations. Assume that the output model has a small empirical risk E[RDm1(A(Dm1,Dm2))]=O(m1−1) .
Figure 1: Results of hyperparameter optimization in data reweighting with varying T and K
Figure 2: Results of hyperparameter optimization in data reweighting with varying K and m1
College of Informatics, Huazhong Agricultural University, Wuhan, 430070, China · School of Artificial Intelligence, Jianghan University, Wuhan, 430056, China
Department of Mathematics, The University of Hong Kong, Hong Kong, China · School of Mathematics and Statistics, The University of Sydney, Sydney, Australia