Juliette Mukangango, Amanda Muyskens, Benjamin W. Priest
{"title":"A Robust Approach to Gaussian Processes Implementation","authors":"Juliette Mukangango, Amanda Muyskens, Benjamin W. Priest","doi":"arxiv-2409.11577","DOIUrl":null,"url":null,"abstract":"Gaussian Process (GP) regression is a flexible modeling technique used to\npredict outputs and to capture uncertainty in the predictions. However, the GP\nregression process becomes computationally intensive when the training spatial\ndataset has a large number of observations. To address this challenge, we\nintroduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest\nneighbor and leave-one-out cross-validation during training. This approach\nenables the evaluation of large spatial datasets with state-of-the-art accuracy\nand speed in certain spatial problems. Despite these advantages, conventional\nquadratic loss functions used in the MuyGPs optimization such as Root Mean\nSquared Error(RMSE), are highly influenced by outliers. We explore the behavior\nof MuyGPs in cases involving outlying observations, and subsequently, develop a\nrobust approach to handle and mitigate their impact. Specifically, we introduce\na novel leave-one-out loss function based on the pseudo-Huber function (LOOPH)\nthat effectively accounts for outliers in large spatial datasets within the\nMuyGPs framework. Our simulation study shows that the \"LOOPH\" loss method\nmaintains accuracy despite outlying observations, establishing MuyGPs as a\npowerful tool for mitigating unusual observation impacts in the large data\nregime. In the analysis of U.S. ozone data, MuyGPs provides accurate\npredictions and uncertainty quantification, demonstrating its utility in\nmanaging data anomalies. Through these efforts, we advance the understanding of\nGP regression in spatial contexts.","PeriodicalId":501215,"journal":{"name":"arXiv - STAT - Computation","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Computation","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.11577","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Gaussian Process (GP) regression is a flexible modeling technique used to
predict outputs and to capture uncertainty in the predictions. However, the GP
regression process becomes computationally intensive when the training spatial
dataset has a large number of observations. To address this challenge, we
introduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest
neighbor and leave-one-out cross-validation during training. This approach
enables the evaluation of large spatial datasets with state-of-the-art accuracy
and speed in certain spatial problems. Despite these advantages, conventional
quadratic loss functions used in the MuyGPs optimization such as Root Mean
Squared Error(RMSE), are highly influenced by outliers. We explore the behavior
of MuyGPs in cases involving outlying observations, and subsequently, develop a
robust approach to handle and mitigate their impact. Specifically, we introduce
a novel leave-one-out loss function based on the pseudo-Huber function (LOOPH)
that effectively accounts for outliers in large spatial datasets within the
MuyGPs framework. Our simulation study shows that the "LOOPH" loss method
maintains accuracy despite outlying observations, establishing MuyGPs as a
powerful tool for mitigating unusual observation impacts in the large data
regime. In the analysis of U.S. ozone data, MuyGPs provides accurate
predictions and uncertainty quantification, demonstrating its utility in
managing data anomalies. Through these efforts, we advance the understanding of
GP regression in spatial contexts.