README.md

In geoFourierFDA: Ordinary Functional Kriging Using Fourier Smoothing and Gaussian Quadrature

geoFourierFDA

In this paper, we have observsed n curves χs1(t), …, χsn(t) in a region, where si = (θi, ηi), i = 1, …, n, and θi is the latitude and ηi is the longitude where the curve χsi was sampled. The goal of this package is to estimate an unsampled curve χs0(t) at s0 ∉ {s1, …, sn}. The Ideia proposed by Giraldo (2011) was simple: the curve χs0(t) is a linear combination of all curves χs1(t), …, χsn(t), i.e., $\widehat{\chi_{\mathbf{s}_0}}(t) = \lambda_1 \chi_{\mathbf{s}_1}(t) + \lambda_2 \chi_{\mathbf{s}_2}(t) + \dots + \lambda_n \chi_{\mathbf{s}_1}(t)$ where λ1, …, λn is solution of the linear system given by

where μ is an constant from the method of Lagrange’s multipliers and the function γ(h) = ∫γ(h; t)dt is called the trace-variogram, where, for each t, γ(h; t) is the semivariogram for the process χ*s1(t), …, χsn(t). More precisely, for each t, a weakly and isotropic spatial process is assumed for χs1(t), …, χsn(t*) and the integration of the semivariogram is carried out. Usually, the integration in the equation (1) is approximated using a modified version of the empirical semivariogram. In this pcackage, we have used the Legendre-Gauss quadrature, which is simple and it explicitly used the definition of the semivariogram.

Installation

This package can be installed using the devtools package.

devtools::install_github("gilberto-sassi/geoFourierFDA")

Examples

In this package, we have used the temperature dataset present in the package fda and in the package geofd. This dataset has temperature measurements from 35 weather stations from Canada. This data can be downloaded at weather.gov.ca. For illustration, we have separated the time series at The Pas station and used all others stations to estimate the curve temperature at The Pas.

How to interpolate a curve at an unmonitored location

# interpolating curve at Halifax using all remaining curves in the functional dataset
data(canada)

# Estimating the temperature at The Pas
geo_fda(canada$m_data, canada$m_coord, canada$ThePas_coord)

Coefficients of smoothing using Fourier series polynomial

# Coefficients of smoothing using Fourier series polynomial
# Coefficients of smoothing at The Pas
data(canada)

coefs <- coef_fourier(canada$ThePas_ts)

Smoothed curve using Fourier series

# Coefficients of smoothing using Fourier series polynomial
# Coefficients of smoothing at The Pas
data(canada)


# coefficients of Fourier series polynomial
coefs <- coef_fourier(canada$ThePas_ts, m)

# points to evaluate curve at interval [-pi, pi]
x <- seq(from = -pi, to = pi, by = 0.01)

# smoothed curve at some points x
y_est <- fourier_b(coefs, x)

References

Giraldo, R, P Delicado, and J Mateu. 2011. “Ordinary Kriging for Function-Valued Spatial Data.” Environmental and Ecological Statistics 18 (3): 411–26.

geoFourierFDA documentation built on Oct. 27, 2021, 5:08 p.m.