Estimation of trace-variogram using Legendre-Gauss quadrature

Functional Data Analysis is known for its application in several fields of science. In some cases, functional datasets are constituted by spatially indexed curves. The primary goal of this paper is to supply a straightforward and precise approach to interpolate these curves, that is, the aim is to estimate a curve at an unmonitored location. It is proven that the best linear unbiased estimator for this unsampled curve is the solution of a linear system, where the coefficients and the constant terms of the system are formed using a function called trace-variogram. In this paper, we propose using Legendre-Gauss quadrature to estimate the trace-variogram. This estimator's suitable numerical properties are shown in simulation studies for normal and non-normal datasets. Simulation results indicated that the proposed methodology outperforms the established estimation procedure. An R package was built and is available at the CRAN repository. The novel estimation methodology is illustrated with a real dataset on temperature curves from 35 weather stations in Canada.