CaltechAUTHORS
  A Caltech Library Service

Analysis of Incomplete Climate Data: Estimation of Mean Values and Covariance Matrices and Imputation of Missing Values

Schneider, Tapio (2001) Analysis of Incomplete Climate Data: Estimation of Mean Values and Covariance Matrices and Imputation of Missing Values. Journal of Climate, 14 (5). pp. 853-871. ISSN 0894-8755. https://resolver.caltech.edu/CaltechAUTHORS:SCHNjc01b

[img]
Preview
PDF
See Usage Policy.

215Kb

Use this Persistent URL to link to this item: https://resolver.caltech.edu/CaltechAUTHORS:SCHNjc01b

Abstract

Estimating the mean and the covariance matrix of an incomplete dataset and filling in missing values with imputed values is generally a nonlinear problem, which must be solved iteratively. The expectation maximization (EM) algorithm for Gaussian data, an iterative method both for the estimation of mean values and covariance matrices from incomplete datasets and for the imputation of missing values, is taken as the point of departure for the development of a regularized EM algorithm. In contrast to the conventional EM algorithm, the regularized EM algorithm is applicable to sets of climate data, in which the number of variables typically exceeds the sample size. The regularized EM algorithm is based on iterated analyses of linear regressions of variables with missing values on variables with available values, with regression coefficients estimated by ridge regression, a regularized regression method in which a continuous regularization parameter controls the filtering of the noise in the data. The regularization parameter is determined by generalized cross-validation, such as to minimize, approximately, the expected mean-squared error of the imputed values. The regularized EM algorithm can estimate, and exploit for the imputation of missing values, both synchronic and diachronic covariance matrices, which may contain information on spatial covariability, stationary temporal covariability, or cyclostationary temporal covariability. A test of the regularized EM algorithm with simulated surface temperature data demonstrates that the algorithm is applicable to typical sets of climate data and that it leads to more accurate estimates of the missing values than a conventional noniterative imputation technique.


Item Type:Article
ORCID:
AuthorORCID
Schneider, Tapio0000-0001-5687-2287
Additional Information:© Copyright by American Meteorological Society 2001 (Manuscript received December 3, 1999, in final form March 27, 2000) This research project was supported by a NASA Earth System Science Fellowship. I thank Keith Dixon for providing the ensemble of climate simulations; Isaac Held, John Lanzante, Michael Mann, and Arnold Neumaier for comments on drafts of this paper; and Heidi Swanson for editing the manuscript.
Issue or Number:5
Record Number:CaltechAUTHORS:SCHNjc01b
Persistent URL:https://resolver.caltech.edu/CaltechAUTHORS:SCHNjc01b
Alternative URL:http://dx.doi.org/10.1175/1520-0442(2001)014<0853:AOICDE>2.0.CO;2
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:3973
Collection:CaltechAUTHORS
Deposited By: Archive Administrator
Deposited On:20 Jul 2006
Last Modified:02 Oct 2019 23:08

Repository Staff Only: item control page