smoothEM: a new approach for the simultaneous assessment of smooth patterns and spikes

06/12/2022
by   Huy Dang, et al.
0

We consider functional data where an underlying smooth curve is composed not just with errors, but also with irregular spikes. We propose an approach that, combining regularized spline smoothing and an Expectation-Maximization algorithm, allows one to both identify spikes and estimate the smooth component. Imposing some assumptions on the error distribution, we prove consistency of EM estimates. Next, we demonstrate the performance of our proposal on finite samples and its robustness to assumptions violations through simulations. Finally, we apply our proposal to data on the annual heatwaves index in the US and on weekly electricity consumption in Ireland. In both datasets, we are able to characterize underlying smooth trends and to pinpoint irregular/extreme behaviors.

READ FULL TEXT

page 18

page 19

page 20

research
11/09/2022

Spline Estimation of Functional Principal Components via Manifold Conjugate Gradient Algorithm

Functional principal component analysis has become the most important di...
research
10/19/2020

A data-driven P-spline smoother and the P-Spline-GARCH-models

Penalized spline smoothing of time series and its asymptotic properties ...
research
02/25/2023

Smooth digital terrain modelling in irregular domain using adaptive finite element thin plate spline smoother

Digital terrain models of geological information occasionally require sm...
research
04/10/2016

Distance for Functional Data Clustering Based on Smoothing Parameter Commutation

We propose a novel method to determine the dissimilarity between subject...
research
10/07/2021

Smooth bootstrapping of copula functionals

The smooth bootstrap for estimating copula functionals in small samples ...
research
04/10/2015

Switching nonparametric regression models for multi-curve data

We develop and apply an approach for analyzing multi-curve data where ea...

Please sign up or login with your details

Forgot password? Click here to reset