Weighted principal component analysis: a weighted covariance eigendecomposition approach

Abstract

We present a new straightforward principal component analysis (PCA) method based on the diagonalization of the weighted variance-covariance matrix through two spectral decomposition methods: power iteration and Rayleigh quotient iteration. This method allows one to retrieve a given number of orthogonal principal components amongst the most meaningful ones for the case of problems with weighted and/or missing data. Principal coefficients are then retrieved by fitting principal components to the data while providing the final decomposition. Tests performed on real and simulated cases show that our method is optimal in the identification of the most significant patterns within data sets. We illustrate the usefulness of this method by assessing its quality on the extrapolation of Sloan Digital Sky Survey quasar spectra from measured wavelengths to shorter and longer wavelengths. Our new algorithm also benefits from a fast and flexible implementation.

12 pages, 9 figures

References in corpus (6)

• A Tutorial on Principal Component Analysis
• Spectral Classification of Quasars in the Sloan Digital Sky Survey: Eigenspectra; Redshift and Luminosity Effects
• The Sloan Digital Sky Survey quasar catalog: tenth data release
• The Gaia astrophysical parameters inference system (Apsis). Pre-launch description
• Principal Component Analysis with Noisy and/or Missing Data
• A data-driven model for spectra: Finding double redshifts in the Sloan Digital Sky Survey

Cited by in corpus (16)

• Gaia DR3: Apsis III -- Non-stellar content and source classification
• The magnetic field and multiple planets of the young dwarf AU~Mic
• Understanding trained CNNs by indexing neuron selectivity
• The Pale Green Dot: A Method to Characterize Proxima Centauri b using Exo-Aurorae
• $\text{Wapiti}$: a data-driven approach to correct for systematics in RV data -- Application to SPIRou data of the planet-hosting M dwarf GJ 251
• Stellar flare morphology with TESS across the main sequence
• FIESTA II. Disentangling stellar and instrumental variability from exoplanetary Doppler shifts in Fourier domain
• New technique for determining a pulsar period: Waterfall principal component analysis
• Determination of astrophysical parameters of quasars within the Gaia mission
• Gl 725A b: a potential super-Earth detected with SOPHIE and SPIRou in an M dwarf binary system at 3.5 pc
• Statistical uncertainties of the $v_n\{2k\}$ harmonics from Q-cumulants
• Redshift determination through weighted phase correlation: a linearithmic implementation
• The SPIRou Legacy Survey: near-infrared and optical radial velocity analysis of Gl 480 and Gl 382 using SPIRou, HARPS and CARMENES spectrographs
• The SPIRou Legacy Survey: Detection of a nearby world orbiting in the habitable zone of Gl725B achieved by correcting strong telluric contamination in near-infrared radial velocities with WAPITI
• Aspects of Spatially-Correlated Random Fields: Extreme-Value Statistics and Clustering Properties
• Detection and skeletonization of single neurons and tracer injections using topological methods