Exploring functional data analysis and wavelet principal component analysis on ecstasy (MDMA) wastewater data

Abstract

BackgroundWastewater-based epidemiology (WBE) is a novel approach in drug use epidemiology which aims to monitor the extent of use of various drugs in a community. In this study, we investigate functional principal component analysis (FPCA) as a tool for analysing WBE data and compare it to traditional principal component analysis (PCA) and to wavelet principal component analysis (WPCA) which is more flexible temporally.MethodsWe analysed temporal wastewater data from 42 European cities collected daily over one week in March 2013. The main temporal features of ecstasy (MDMA) were extracted using FPCA using both Fourier and B-spline basis functions with three different smoothing parameters, along with PCA and WPCA with different mother wavelets and shrinkage rules. The stability of FPCA was explored through bootstrapping and analysis of sensitivity to missing data.ResultsThe first three principal components (PCs), functional principal components (FPCs) and wavelet principal components (WPCs) explained 87.5-99.6 % of the temporal variation between cities, depending on the choice of basis and smoothing. The extracted temporal features from PCA, FPCA and WPCA were consistent. FPCA using Fourier basis and common-optimal smoothing was the most stable and least sensitive to missing data.ConclusionFPCA is a flexible and analytically tractable method for analysing temporal changes in wastewater data, and is robust to missing data. WPCA did not reveal any rapid temporal changes in the data not captured by FPCA. Overall the results suggest FPCA with Fourier basis functions and common-optimal smoothing parameter as the most accurate approach when analysing WBE data.Electronic supplementary materialThe online version of this article (doi:10.1186/s12874-016-0179-2) contains supplementary material, which is available to authorized users.

Genetic determinants of glucose levels in pregnancy: genetic risk scores analysis and GWAS in the Norwegian STORK cohort

“…1B). A common smoothing parameter was used in both smoothing procedures (21). These 497 and 506 individually fitted curves formed the basis for two separate functional principal component analyses (FPCAs), in which curve shape information was extracted.…”

Functional principal component analysis for identifying the child growth pattern using longitudinal birth cohort data

“…Independent functional principal component curves describe the important modes of temporal variability in growth across the individual fitted curves. FPCA also reduces the dimensions of the problem by representing functions in terms of a finite set of functions and further functional linear model was used to assess the association between factors and trajectories [ 20 , 22 – 24 , 32 – 36 ]. The conditional kernel density estimators plot was used to identify the subgroup of the growth functions and, the proportion of children contributing to each subgroup were estimated.…”

Erratum to: Exploring functional data analysis and wavelet principal component analysis on ecstasy (MDMA) wastewater data

“…After publication of the original article [ 1 ], it came to the authors’ attention that there were errors in Fig. 3 , Fig.…”