Generalization Error of Generalized Linear Models in High Dimensions

Total Page:16

File Type:pdf, Size:1020Kb

Generalization Error of Generalized Linear Models in High Dimensions Generalization Error of Generalized Linear Models in High Dimensions Melikasadat Emami 1 Mojtaba Sahraee-Ardakan 1 Parthe Pandit 1 Sundeep Rangan 2 Alyson K. Fletcher 1 3 Abstract a class of generalized linear models (GLMs) of the form 0 At the heart of machine learning lies the question y = φout( x; w ; d); (1) of generalizability of learned rules over previously where x 2 Rp is a vector of input features, y is a scalar out- unseen data. While over-parameterized models 0 p put, w 2 R are weights to be learned, φout(·) is a known based on neural networks are now ubiquitous in 0 machine learning applications, our understanding link function, and d is random noise. The notation x; w of their generalization capabilities is incomplete. denotes an inner product. We use the superscript “0” to This task is made harder by the non-convexity of denote the “true” values in contrast to estimated or postu- the underlying learning problems. We provide lated quantities. The output may be continuous or discrete a general framework to characterize the asymp- to model either regression or classification problems. totic generalization error for single-layer neural We measure the generalization error in a standard man- networks (i.e., generalized linear models) with ner: we are given training data (xi; yi), i = 1;:::;N from arbitrary non-linearities, making it applicable to which we learn some parameter estimate wb via a regularized regression as well as classification problems. This empirical risk minimization of the form framework enables analyzing the effect of (i) over- parameterization and non-linearity during mod- w = argmin Fout(y; Xw) + Fin(w); (2) b w eling; and (ii) choices of loss function, initializa- T tion, and regularizer during learning. Our model where X = [x1 x2 ::: xN ] , is the data matrix, Fout is also captures mismatch between training and test some output loss function, and Fin is some regularizer on distributions. As examples, we analyze a few spe- the weights. We are then given a new test sample, xts, for cial cases, namely linear regression and logistic which the true and predicted values are given by regression. We are also able to rigorously and an- 0 alytically explain the double descent phenomenon yts = φout( xts; w ; dts); ybts = φ(hxts; wbi); (3) in generalized linear models. where dts is the noise in the test sample, and φ(·) is a pos- tulated inverse link function that may be different from the true function φout(·). The generalization error is then de- 1. Introduction fined as the expectation of some expected loss between yts and y of the form A fundamental goal of machine learning is generalization: bts the ability to draw inferences about unseen data from finite E fts(yts; ybts); (4) arXiv:2005.00180v1 [cs.LG] 1 May 2020 training examples. Methods to quantify the generalization error are therefore critical in assessing the performance of for some test loss function fts(·) such as squared error or any machine learning approach. prediction error. This paper seeks to characterize the generalization error for Even for this relatively simple GLM model, the behavior of the generalization error is not fully understood. Recent 1 Department of Electrical and Computer Engineering, Univer- works (Montanari et al., 2019; Deng et al., 2019; Mei & sity of California, Los Angeles, Los Angeles, USA 2Department of Electrical and Computer Engineering, New York University, Montanari, 2019; Salehi et al., 2019) have characterized Brooklyn, New York, USA 3Department of Statistics, University the generalization error of various linear models for clas- of California, Los Angeles, Los Angeles, USA. Correspondence sification and regression in certain large random problem to: Melikasadat Emami <[email protected]>, Mojtaba Sahraee- instances. Specifically, the number of samples N and num- Ardakan <[email protected]>, Parthe Pandit <parthepan- ber of features p both grow without bound with their ratio [email protected]>, Sundeep Rangan <[email protected]>, Alyson K. Fletcher <akfl[email protected]>. satisfying p=N ! β 2 (0; 1), and the samples in the training data xi are drawn randomly. In this limit, the gen- eralization error can be exactly computed. The analysis can Generalization Error of GLMs in High Dimensions explain the so-called double descent phenomena (Belkin using the Convex Gaussian Min-max Theorem in the under- et al., 2019a): in highly under-regularized settings, the test parametrized regime. The results in the current paper can error may initially increase with the number of data samples consider all these models as special cases. Bounds on the N before decreasing. See the prior work section below for generalization error of over-parametrized linear models are more details. also given in (Bartlett et al., 2019; Neyshabur et al., 2018). Although this paper and several other recent works consider Summary of Contributions. Our main result (Theo- only simple linear models and GLMs, much of the motiva- rem1) provides a procedure for exactly computing the tion is to understand generalization in deep neural networks asymptotic value of the generalization error (4) for GLM where classical intuition may not hold (Belkin et al., 2018; models in a certain random high-dimensional regime called Zhang et al., 2016; Neyshabur et al., 2018). In particular, the Large System Limit (LSL). The procedure enables the a number of recent papers have shown the connection be- generalization error to be related to key problem parameters tween neural networks in the over-parametrized regime and including the sampling ratio β = p=N, the regularizer, the kernel methods. The works (Daniely, 2017; Daniely et al., output function, and the distributions of the true weights 2016) showed that gradient descent on over-parametrized and noise. Importantly, our result holds under very general neural networks learns a function in the RKHS correspond- settings including: (i) arbitrary test metrics fts; (ii) arbi- ing to the random feature kernel. Training dynamics of over- trary training loss functions Fout as well as decomposable parametrized neural networks has been studied by (Jacot regularizers Fin; (iii) arbitrary link functions φout; (iv) cor- et al., 2018; Du et al., 2018; Arora et al., 2019; Allen-Zhu related covariates x; (v) underparameterized (β < 1) and et al., 2019), and it is shown that the function learned is in overparameterized regimes (β > 1); and (vi) distributional an RKHS corresponding to the neural tangent kernel. mismatch in training and test data. Section4 discusses in detail the general assumptions on the quantities f , F , ts out Approximate Message Passing. Our key tool to study F , and φ under which Theorem1 holds. in out the generalization error is approximate message passing (AMP), a class of inference algorithms originally developed Prior Work. Many recent works characterize generaliza- in (Donoho et al., 2009; 2010; Bayati & Montanari, 2011) tion error of various machine learning models, including for compressed sensing. We show that the learning problem special cases of the GLM model considered here. For exam- for the GLM can be formulated as an inference problem ple, the precise characterization for asymptotics of predic- on a certain multi-layer network. Multi-layer AMP meth- tion error for least squares regression has been provided in ods (He et al., 2017; Manoel et al., 2018; Fletcher et al., (Belkin et al., 2019b; Hastie et al., 2019; Muthukumar et al., 2018; Pandit et al., 2019) can then be applied to perform 2019). The former confirmed the double descent curve of the inference. The specific algorithm we use in this work (Belkin et al., 2019a) under a Fourier series model and a is the multi-layer vector AMP (ML-VAMP) algorithm of noisy Gaussian model for data in the over-parameterized (Fletcher et al., 2018; Pandit et al., 2019) which itself builds regime. The latter also obtained this scenario under both on several works (Opper & Winther, 2005; Fletcher et al., linear and non-linear feature models for ridge regression 2016; Rangan et al., 2019; Cakmak et al., 2014; Ma & Ping, and min-norm least squares using random matrix theory. 2017). The ML-VAMP algorithm is not necessarily the most Also, (Advani & Saxe, 2017) studied the same setting for computationally efficient procedure for the minimization deep linear and shallow non-linear networks. (2). For our purposes, the key property is that ML-VAMP enables exact predictions of its performance in the large sys- The analysis of the the generalization for max-margin linear tem limit. Specifically, the error of the algorithm estimates classifiers in the high dimensional regime has been done in in each iteration can be predicted by a set of deterministic re- (Montanari et al., 2019). The exact expression for asymp- cursive equations called the state evolution or SE. The fixed totic prediction error is derived and in a specific case for points of these equations provide a way of computing the two-layer neural network with random first-layer weights, asymptotic performance of the algorithm. In certain cases, the double descent curve was obtained. A similar dou- the algorithm can be proven to be Bayes optimal (Reeves, ble descent curve for logistic regression as well as linear 2017; Gabrie´ et al., 2018; Barbier et al., 2019). discriminant analysis has been reported by (Deng et al., 2019). Random feature learning in the same setting has This approach of using AMP methods to characterize the also been studied for ridge regression in (Mei & Monta- generalization error of GLMs was also explored in (Barbier nari, 2019). The authors have, in particular, shown that et al., 2019) for i.i.d. distributions on the data. The explicit highly over-parametrized estimators with zero training error formulae for the asymptotic mean squared error for the are statistically optimal at high signal-to-noise ratio (SNR).
Recommended publications
  • A Simple Algorithm for Semi-Supervised Learning with Improved Generalization Error Bound
    A Simple Algorithm for Semi-supervised Learning with Improved Generalization Error Bound Ming Ji∗‡ [email protected] Tianbao Yang∗† [email protected] Binbin Lin\ [email protected] Rong Jiny [email protected] Jiawei Hanz [email protected] zDepartment of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA yDepartment of Computer Science and Engineering, Michigan State University, East Lansing, MI, 48824, USA \State Key Lab of CAD&CG, College of Computer Science, Zhejiang University, Hangzhou, 310058, China ∗Equal contribution Abstract which states that the prediction function lives in a low In this work, we develop a simple algorithm dimensional manifold of the marginal distribution PX . for semi-supervised regression. The key idea It has been pointed out by several stud- is to use the top eigenfunctions of integral ies (Lafferty & Wasserman, 2007; Nadler et al., operator derived from both labeled and un- 2009) that the manifold assumption by itself is labeled examples as the basis functions and insufficient to reduce the generalization error bound learn the prediction function by a simple lin- of supervised learning. However, on the other hand, it ear regression. We show that under appropri- was found in (Niyogi, 2008) that for certain learning ate assumptions about the integral operator, problems, no supervised learner can learn effectively, this approach is able to achieve an improved while a manifold based learner (that knows the man- regression error bound better than existing ifold or learns it from unlabeled examples) can learn bounds of supervised learning. We also veri- well with relatively few labeled examples.
    [Show full text]
  • Statistical Mechanics Methods for Discovering Knowledge from Production-Scale Neural Networks
    Statistical Mechanics Methods for Discovering Knowledge from Production-Scale Neural Networks Charles H. Martin∗ and Michael W. Mahoney† Tutorial at ACM-KDD, August 2019 ∗Calculation Consulting, [email protected] †ICSI and Dept of Statistics, UC Berkeley, https://www.stat.berkeley.edu/~mmahoney/ Martin and Mahoney (CC & ICSI/UCB) Statistical Mechanics Methods August 2019 1 / 98 Outline 1 Prehistory and History Older Background A Very Simple Deep Learning Model More Immediate Background 2 Preliminary Results Regularization and the Energy Landscape Preliminary Empirical Results Gaussian and Heavy-tailed Random Matrix Theory 3 Developing a Theory for Deep Learning More Detailed Empirical Results An RMT-based Theory for Deep Learning Tikhonov Regularization versus Heavy-tailed Regularization 4 Validating and Using the Theory Varying the Batch Size: Explaining the Generalization Gap Using the Theory: pip install weightwatcher Diagnostics at Scale: Predicting Test Accuracies 5 More General Implications and Conclusions Outline 1 Prehistory and History Older Background A Very Simple Deep Learning Model More Immediate Background 2 Preliminary Results Regularization and the Energy Landscape Preliminary Empirical Results Gaussian and Heavy-tailed Random Matrix Theory 3 Developing a Theory for Deep Learning More Detailed Empirical Results An RMT-based Theory for Deep Learning Tikhonov Regularization versus Heavy-tailed Regularization 4 Validating and Using the Theory Varying the Batch Size: Explaining the Generalization Gap Using the Theory: pip install weightwatcher Diagnostics at Scale: Predicting Test Accuracies 5 More General Implications and Conclusions Statistical Physics & Neural Networks: A Long History 60s: I J. D. Cowan, Statistical Mechanics of Neural Networks, 1967. 70s: I W. A. Little, “The existence of persistent states in the brain,” Math.
    [Show full text]
  • Multitask Learning with Local Attention for Tibetan Speech Recognition
    Hindawi Complexity Volume 2020, Article ID 8894566, 10 pages https://doi.org/10.1155/2020/8894566 Research Article Multitask Learning with Local Attention for Tibetan Speech Recognition Hui Wang , Fei Gao , Yue Zhao , Li Yang , Jianjian Yue , and Huilin Ma School of Information Engineering, Minzu University of China, Beijing 100081, China Correspondence should be addressed to Fei Gao; [email protected] Received 22 September 2020; Revised 12 November 2020; Accepted 26 November 2020; Published 18 December 2020 Academic Editor: Ning Cai Copyright © 2020 Hui Wang et al. *is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. In this paper, we propose to incorporate the local attention in WaveNet-CTC to improve the performance of Tibetan speech recognition in multitask learning. With an increase in task number, such as simultaneous Tibetan speech content recognition, dialect identification, and speaker recognition, the accuracy rate of a single WaveNet-CTC decreases on speech recognition. Inspired by the attention mechanism, we introduce the local attention to automatically tune the weights of feature frames in a window and pay different attention on context information for multitask learning. *e experimental results show that our method improves the accuracies of speech recognition for all Tibetan dialects in three-task learning, compared with the baseline model. Furthermore, our method significantly improves the accuracy for low-resource dialect by 5.11% against the specific-dialect model. 1. Introduction recognition, dialect identification, and speaker recognition in a single model.
    [Show full text]
  • Lecture 9: Generalization
    Lecture 9: Generalization Roger Grosse 1 Introduction When we train a machine learning model, we don't just want it to learn to model the training data. We want it to generalize to data it hasn't seen before. Fortunately, there's a very convenient way to measure an algorithm's generalization performance: we measure its performance on a held-out test set, consisting of examples it hasn't seen before. If an algorithm works well on the training set but fails to generalize, we say it is overfitting. Improving generalization (or preventing overfitting) in neural nets is still somewhat of a dark art, but this lecture will cover a few simple strategies that can often help a lot. 1.1 Learning Goals • Know the difference between a training set, validation set, and test set. • Be able to reason qualitatively about how training and test error de- pend on the size of the model, the number of training examples, and the number of training iterations. • Understand the motivation behind, and be able to use, several strate- gies to improve generalization: { reducing the capacity { early stopping { weight decay { ensembles { input transformations { stochastic regularization 2 Measuring generalization So far in this course, we've focused on training, or optimizing, neural net- works. We defined a cost function, the average loss over the training set: N 1 X L(y(x(i)); t(i)): (1) N i=1 But we don't just want the network to get the training examples right; we also want it to generalize to novel instances it hasn't seen before.
    [Show full text]
  • Issues in Using Function Approximation for Reinforcement Learning
    To appear in: Proceedings of the Fourth Connectionist Models Summer School Lawrence Erlbaum Publisher, Hillsdale, NJ, Dec. 1993 Issues in Using Function Approximation for Reinforcement Learning Sebastian Thrun Anton Schwartz Institut fÈur Informatik III Dept. of Computer Science UniversitÈat Bonn Stanford University RÈomerstr. 164, D±53225 Bonn, Germany Stanford, CA 94305 [email protected] [email protected] Reinforcement learning techniques address the problem of learning to select actions in unknown, dynamic environments. It is widely acknowledged that to be of use in complex domains, reinforcement learning techniques must be combined with generalizing function approximation methods such as arti®cial neural networks. Little, however, is understood about the theoretical properties of such combinations, and many researchers have encountered failures in practice. In this paper we identify a prime source of such failuresÐnamely, a systematic overestimation of utility values. Using Watkins' Q-Learning [18] as an example, we give a theoretical account of the phenomenon, deriving conditions under which one may expected it to cause learning to fail. Employing some of the most popular function approximators, we present experimental results which support the theoretical ®ndings. 1 Introduction Reinforcement learning methods [1, 16, 18] address the problem of learning, through experimentation, to choose actions so as to maximize one's productivity in unknown, dynamic environments. Unlike most learning algorithms that have been studied in the ®eld of Machine Learning, reinforcement learning techniques allow for ®nding optimal action sequences in temporal decision tasks where the external evaluation is sparse, and neither the effects of actions, nor the temporal delay between actions and its effects on the learner's performance is known to the learner beforehand.
    [Show full text]
  • Effective Dimensionality Revisited
    Under review as a conference paper at ICLR 2021 RETHINKING PARAMETER COUNTING: EFFECTIVE DIMENSIONALITY REVISITED Anonymous authors Paper under double-blind review ABSTRACT Neural networks appear to have mysterious generalization properties when using parameter counting as a proxy for complexity. Indeed, neural networks often have many more parameters than there are data points, yet still provide good general- ization performance. Moreover, when we measure generalization as a function of parameters, we see double descent behaviour, where the test error decreases, increases, and then again decreases. We show that many of these properties become understandable when viewed through the lens of effective dimensionality, which measures the dimensionality of the parameter space determined by the data. We relate effective dimensionality to posterior contraction in Bayesian deep learning, model selection, width-depth tradeoffs, double descent, and functional diversity in loss surfaces, leading to a richer understanding of the interplay between parameters and functions in deep models. We also show that effective dimensionality compares favourably to alternative norm- and flatness- based generalization measures. 1 INTRODUCTION Parameter counting pervades the narrative in modern deep learning. “One of the defining properties of deep learning is that models are chosen to have many more parameters than available training data. In light of this capacity for overfitting, it is remarkable that simple algorithms like SGD reliably return solutions with low test error” (Dziugaite and Roy, 2017). “Despite their massive size, successful deep artificial neural networks can exhibit a remarkably small difference between training and test performance” (Zhang et al., 2017). “Increasing the number of parameters of neural networks can give much better prediction accuracy” (Shazeer et al., 2017).
    [Show full text]
  • Generalizing to Unseen Domains: a Survey on Domain Generalization
    PREPRINT 1 Generalizing to Unseen Domains: A Survey on Domain Generalization Jindong Wang, Cuiling Lan, Member, IEEE, Chang Liu, Yidong Ouyang, Wenjun Zeng, Fellow, IEEE, Tao Qin, Senior Member, IEEE Abstract—Machine learning systems generally assume that the Sketch Cartoon Art painting Photo training and testing distributions are the same. To this end, a key requirement is to develop models that can generalize to unseen distributions. Domain generalization (DG), i.e., out-of- distribution generalization, has attracted increasing interests in recent years. Domain generalization deals with a challenging setting where one or several different but related domain(s) are given, and the goal is to learn a model that can generalize to Training set Test set an unseen test domain. Great progress has been made in the area of domain generalization for years. This paper presents Fig. 1. Examples from the dataset PACS [1] for domain generalization. The the first review of recent advances in this area. First, we training set is composed of images belonging to domains of sketch, cartoon, provide a formal definition of domain generalization and discuss and art paintings. DG aims to learn a generalized model that performs well several related fields. We then thoroughly review the theories on the unseen target domain of photos. related to domain generalization and carefully analyze the theory behind generalization. We categorize recent algorithms into three classes: data manipulation, representation learning, and learning datasets) that will generalize well on unseen testing domains. strategy, and present several popular algorithms in detail for For instance, given a training set consisting of images coming each category.
    [Show full text]
  • Batch Policy Learning Under Constraints
    Batch Policy Learning under Constraints Hoang M. Le 1 Cameron Voloshin 1 Yisong Yue 1 Abstract deed, many such real-world applications require the primary objective function be augmented with an appropriate set of When learning policies for real-world domains, constraints (Altman, 1999). two important questions arise: (i) how to effi- ciently use pre-collected off-policy, non-optimal Contemporary policy learning research has largely focused behavior data; and (ii) how to mediate among dif- on either online reinforcement learning (RL) with a focus on ferent competing objectives and constraints. We exploration, or imitation learning (IL) with a focus on learn- thus study the problem of batch policy learning un- ing from expert demonstrations. However, many real-world der multiple constraints, and offer a systematic so- settings already contain large amounts of pre-collected data lution. We first propose a flexible meta-algorithm generated by existing policies (e.g., existing driving behav- that admits any batch reinforcement learning and ior, power grid control policies, etc.). We thus study the online learning procedure as subroutines. We then complementary question: can we leverage this abundant present a specific algorithmic instantiation and source of (non-optimal) behavior data in order to learn se- provide performance guarantees for the main ob- quential decision making policies with provable guarantees jective and all constraints. As part of off-policy on both primary objective and constraint satisfaction? learning, we propose a simple method for off- We thus propose and study the problem of batch policy policy policy evaluation (OPE) and derive PAC- learning under multiple constraints.
    [Show full text]
  • Learning Better Structured Representations Using Low-Rank Adaptive Label Smoothing
    Published as a conference paper at ICLR 2021 LEARNING BETTER STRUCTURED REPRESENTATIONS USING LOW-RANK ADAPTIVE LABEL SMOOTHING Asish Ghoshal, Xilun Chen, Sonal Gupta, Luke Zettlemoyer & Yashar Mehdad {aghoshal,xilun,sonalgupta,lsz,mehdad}@fb.com Facebook AI ABSTRACT Training with soft targets instead of hard targets has been shown to improve per- formance and calibration of deep neural networks. Label smoothing is a popular way of computing soft targets, where one-hot encoding of a class is smoothed with a uniform distribution. Owing to its simplicity, label smoothing has found wide- spread use for training deep neural networks on a wide variety of tasks, ranging from image and text classification to machine translation and semantic parsing. Complementing recent empirical justification for label smoothing, we obtain PAC- Bayesian generalization bounds for label smoothing and show that the generaliza- tion error depends on the choice of the noise (smoothing) distribution. Then we propose low-rank adaptive label smoothing (LORAS): a simple yet novel method for training with learned soft targets that generalizes label smoothing and adapts to the latent structure of the label space in structured prediction tasks. Specifi- cally, we evaluate our method on semantic parsing tasks and show that training with appropriately smoothed soft targets can significantly improve accuracy and model calibration, especially in low-resource settings. Used in conjunction with pre-trained sequence-to-sequence models, our method achieves state of the art performance on four semantic parsing data sets. LORAS can be used with any model, improves performance and implicit model calibration without increasing the number of model parameters, and can be scaled to problems with large label spaces containing tens of thousands of labels.
    [Show full text]
  • Learning Universal Graph Neural Network Embeddings with Aid of Transfer Learning
    Learning Universal Graph Neural Network Embeddings With Aid Of Transfer Learning Saurabh Verma Zhi-Li Zhang Department of Computer Science Department of Computer Science University of Minnesota Twin Cities University of Minnesota Twin Cities [email protected] [email protected] Abstract Learning powerful data embeddings has become a center piece in machine learn- ing, especially in natural language processing and computer vision domains. The crux of these embeddings is that they are pretrained on huge corpus of data in a unsupervised fashion, sometimes aided with transfer learning. However currently in the graph learning domain, embeddings learned through existing graph neural networks (GNNs) are task dependent and thus cannot be shared across different datasets. In this paper, we present a first powerful and theoretically guaranteed graph neural network that is designed to learn task-independent graph embed- dings, thereafter referred to as deep universal graph embedding (DUGNN). Our DUGNN model incorporates a novel graph neural network (as a universal graph encoder) and leverages rich Graph Kernels (as a multi-task graph decoder) for both unsupervised learning and (task-specific) adaptive supervised learning. By learning task-independent graph embeddings across diverse datasets, DUGNN also reaps the benefits of transfer learning. Through extensive experiments and ablation studies, we show that the proposed DUGNN model consistently outperforms both the existing state-of-art GNN models and Graph Kernels by an increased accuracy of 3% − 8% on graph classification benchmark datasets. 1 Introduction Learning powerful data embeddings has become a center piece in machine learning for producing superior results. This new trend of learning embeddings from data can be attributed to the huge success of word2vec [31, 35] with unprecedented real-world performance in natural language processing (NLP).
    [Show full text]
  • 1 Introduction 2 the Learning Problem
    9.520: Statistical Learning Theory and Applications February 8th, 2010 The Learning Problem and Regularization Lecturer: Tomaso Poggio Scribe: Sara Abbasabadi 1 Introduction Today we will introduce a few basic mathematical concepts underlying conditions of the learning theory and generalization property of learning algorithms. We will dene empirical risk, expected risk and empirical risk minimization. Then regularization algorithms will be introduced which forms the bases for the next few classes. 2 The Learning Problem Consider the problem of supervised learning where for a number of input data X, the output labels Y are also known. Data samples are taken independently from an underlying distribution µ(z) on Z = X × Y and form the training set S: (x1; y1);:::; (xn; yn) that is z1; : : : ; zn. We assume independent and identically distributed (i.i.d.) random variables which means the order does not matter (this is the exchangeability property). Our goal is to nd the conditional probability of y given a data x. µ(z) = p(x; y) = p(yjx) · p(x) p(x; y) is xed but unknown. Finding this probabilistic relation between X and Y solves the learning problem. XY P(y|x) P(x) For example suppose we have the data xi;Fi measured from a spring. Hooke's law F = ax gives 2 2 p(F=x) = δ(F −ax). For additive Gaussian noise Hooke's law F = ax gives p(F=x) = e−(F −ax) =2σ . 2-1 2.1 Hypothesis Space Hypothesis space H should be dened in every learning algorithm as the space of functions that can be explored.
    [Show full text]
  • Principles and Algorithms for Forecasting Groups of Time Series: Locality and Globality
    Principles and Algorithms for Forecasting Groups of Time Series: Locality and Globality Pablo Montero-Manso Rob J Hyndman Discipline of Business Analytics Department of Econometrics and Business Statistics University of Sydney Monash University Darlington, NSW 2006, Australia Clayton, VIC 3800, Australia [email protected] [email protected] March 30, 2021 Abstract Forecasting of groups of time series (e.g. demand for multiple products offered by a retailer, server loads within a data center or the number of completed ride shares in zones within a city) can be approached locally, by considering each time series as a separate regression task and fitting a function to each, or globally, by fitting a single function to all time series in the set. While global methods can outperform local for groups composed of similar time series, recent empirical evidence shows surprisingly good performance on heterogeneous groups. This suggests a more general applicability of global methods, potentially leading to more accurate tools and new scenarios to study. However, the evidence has been of empirical nature and a more fundamental study is required. Formalizing the setting of forecasting a set of time series with local and global methods, we provide the following contributions: • We show that global methods are not more restrictive than local methods for time series forecasting, a result which does not apply to sets of regression problems in general. Global and local methods can produce the same forecasts without any assumptions about similarity of the series in the set, therefore global models can succeed in a wider range of problems than previously thought.
    [Show full text]