“…To alleviate the computation burden for smoothing splines, various basis selection methods have been developed. These methods are of similar nature to the subsampling methods, which are widely used in large-scale data analysis [Mahoney, 2011, Drineas et al, 2012, Ma et al, 2015b, Ma and Sun, 2015, Meng et al, 2017, Zhang et al, 2018b, Ai et al, 2021, Xie et al, 2019, Ma et al, 2020, Yu et al, 2020, Ai et al, 2020, Meng et al, 2020a, Li and Meng, 2021.…”