Volume 1 Supplement 1

BioSysBio 2007: Systems Biology, Bioinformatics, Synthetic Biology

Open Access

Variance stabilising transformations for NMR metabolomics data

BMC Systems Biology20071(Suppl 1):P22


Published: 8 May 2007

Classifying the fingerprint of an NMR spectrum is a crucial step in many metabolomics experiments. Since many classification techniques such as principal component analysis (PCA) depend upon variance discrepancies, it is important to first maximise any contribution from wanted class variance between biological samples and minimise any contribution from unwanted technical variance arising from the preparation of the samples and measurement of the NMR metabolic fingerprints. The generalised logarithm (glog) transform was developed to stabilise the variance between technical replicates in a two component error model [1] and has also been applied to NMR spectra previously [2]. To increase the effectiveness of the transform on NMR spectra, the glog was extended to include a baseline offset term. This decreases the unwanted noise contribution on the transformed spectra. The extended glog transformation is given as:

z = ln ( ( y y 0 ) + ( y y 0 ) 2 + λ ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWG6bGEcqGH9aqpcyGGSbaBcqGGUbGBdaqadaqaaiabcIcaOiabdMha5jabgkHiTiabdMha5naaBaaaleaacqaIWaamaeqaaOGaeiykaKIaey4kaSYaaOaaaeaacqGGOaakcqWG5bqEcqGHsislcqWG5bqEdaWgaaWcbaGaeGimaadabeaakiabcMcaPmaaCaaaleqabaGaeGOmaidaaOGaey4kaSIaeq4UdWgaleqaaaGccaGLOaGaayzkaaaaaa@45B6@

for z the transformed intensity and y the original intensity of the spectra. y0 and λ are transformation parameters which are to be found.

Here we have applied the extended glog transform to technical replicates of NMR spectra of tissue extracts from marine mussels, to determine the optimised transformation parameters λ and y0. Next we applied the optimised transformation to a data set comprised of two classes of NMR spectra from stressed and unstressed mussels. Following transformation, the results show significantly better separation of the classes on a PCA scores plot than can be achieved with both untransformed data and also data transformed using Pareto scaling, a widely used method in NMR metabolomics [3]. In conclusion, we have demonstrated the value of the extended glog transformation to stabilise the technical variance in an NMR metabolomics dataset and have achieved significantly improved classification of NMR fingerprints from stressed and unstressed animals.

Authors’ Affiliations

School of Biosciences, The University of Birmingham


  1. Rocke D, Lorenzato S: A two-component model of measurement error in analytical chemistry. Technometrics. 1995, 37 (2): 176-184. 10.2307/1269619. 10.2307/1269619View ArticleGoogle Scholar
  2. Purohit P, Rocke D, Viant M, Woodruff D: Discrimination models using variance-stabilizing transformation of metabolomic NMR data. OMICS. 2004, 8 (2): 118-130. 10.1089/1536231041388348PubMedView ArticleGoogle Scholar
  3. Keun H, Ebbels T, Antti H, Bollard M, Beckonert O, Holmes E, Lindon J, Nicholson J: Improved analysis of multivariate data by variable stability scaling: application to NMR-based metabolic profiling. Analytica Chimica Acta. 2003, 490 (1): 265-276. 10.1016/S0003-2670(03)00094-1.View ArticleGoogle Scholar


© Parsons and Viant; licensee BioMed Central Ltd. 2007

This article is published under license to BioMed Central Ltd.