Search PubMed⌕ Search

Biomedical subjects

L M C Buydens

Publications and source records attributed to L M C Buydens.

11 recordsLinked to original sources

Interpretation of ANOVA models for microarray data using PCA.

MOTIVATION: ANOVA is a technique, which is frequently used in the analysis of microarray data, e.g. to assess the significance of treatment effects, and to select interesting genes based on P-values. However, it does not give information about what exactly is causing the effect. Our purpose is to improve the interpretation of the results from ANOVA on large microarray datasets, by applying PCA on the individual variance components. Interaction effects can be visualized by biplots, showing genes and variables in one plot, providing insight in the effect of e.g. treatment or time on gene expression. Because ANOVA has removed uninteresting sources of variance, the results are much more interpretable than without ANOVA. Moreover, the combination of ANOVA and PCA provides a simple way to select genes, based on the interactions of interest. RESULTS: It is shown that the components from an ANOVA model can be summarized and visualized with PCA, which improves the interpretability of the models. The method is applied to a real time-course gene expression dataset of mesenchymal stem cells. The dataset was designed to investigate the effect of different treatments on osteogenesis. The biplots generated with the algorithm give specific information about the effects of specific treatments on genes over time. These results are in agreement with the literature. The biological validation with GO annotation from the genes present in the selections shows that biologically relevant groups of genes are selected. AVAILABILITY: R code with the implementation of the method for this dataset is available from http://www.cac.science.ru.nl under the heading "Software".

Algorithms↗

Real-life applications of the MULVADO software package for processing DOSY NMR data.

MULVADO is a newly developed software package for DOSY NMR data processing, based on multivariate curve resolution (MCR), one of the principal multivariate methods for processing DOSY data. This paper will evaluate this software package by using real-life data of materials used in the printing industry: two data sets from the same ink sample but of different quality. Also a sample of an organic photoconductor and a toner sample are analysed. Compared with the routine DOSY output from monoexponential fitting, one of the single channel algorithms in the commercial Bruker software, MULVADO provides several advantages. The key advantage of MCR is that it overcomes the fluctuation problem (non-consistent diffusion coefficient of the same component). The combination of non-linear regression (NLR) and MCR can yield more accurate resolution of a complex mixture. In addition, the data pre-processing techniques in MULVADO minimise the negative effects of experimental artefacts on the results of the data. In this paper, the challenges for analysing polymer samples and other more complex samples will also be discussed.

Magnetic Resonance Spectroscopy↗

Method for the computational comparison of crystal structures.

A new method for assessing the similarity of crystal structures is described. A similarity measure is important in classification and clustering problems in which the crystal structures are the source of information. Classification is particularly important for the understanding of properties of crystals, while clustering can be used as a data reduction step in polymorph prediction. The method described uses a radial distribution function that combines atomic coordinates with partial atomic charges. The descriptor is validated using experimental data from a classification study of clathrate structures of cephalosporins and data from a polymorph prediction run. In both cases, excellent results were obtained.

Journal Article↗

Diagnostic analysis of experimental artefacts in DOSY NMR data by covariance matrix of the residuals.

Multivariate curve resolution (MCR) has been applied to separate pure spectra and pure decay profiles of DOSY NMR data. Given good initial guesses of the pure decay profiles, and combined with the nonlinear least square regression (NLR), MCR can result in good separation of the pure components. Nevertheless, due to the presence of artefacts in experimental data, validation of a MCR model is still necessary. In this paper, the covariance matrix of the residuals (CMR), obtained by postmultiplying the residual matrix with its transpose, is proposed to evaluate the quality of the results of an experimental data set. Plots of the rows of this matrix give a general impression of the covariance in the frequency domain of the residual matrix. Different patterns in the plot indicate possible causes of experimental imperfections. This new criterion can be used as diagnosis in order to improve experimental settings as well as suggest appropriate preprocessing of DOSY NMR data.

Journal Article↗

The use of multivariate MR imaging intensities versus metabolic data from MR spectroscopic imaging for brain tumour classification.

This study investigated the value of information from both magnetic resonance imaging and magnetic resonance spectroscopic imaging (MRSI) to automated discrimination of brain tumours. The influence of imaging intensities and metabolic data was tested by comparing the use of MR spectra from MRSI, MR imaging intensities, peak integration values obtained from the MR spectra and a combination of the latter two. Three classification techniques were objectively compared: linear discriminant analysis, least squares support vector machines (LS-SVM) with a linear kernel as linear techniques and LS-SVM with radial basis function kernel as a nonlinear technique. Classifiers were evaluated over 100 stratified random splittings of the dataset into training and test sets. The area under the receiver operating characteristic (ROC) curve (AUC) was used as a global performance measure on test data. In general, all techniques obtained a high performance when using peak integration values with or without MR imaging intensities. For example for low- versus high-grade tumours, low- versus high-grade gliomas and gliomas versus meningiomas, the mean test AUC was higher than 0.91, 0.94, and 0.99, respectively, when both MR imaging intensities and peak integration values were used. The use of metabolic data from MRSI significantly improved automated classification of brain tumour types compared to the use of MR imaging intensities solely.

Brain Chemistry↗

Improved DOSY NMR data processing by data enhancement and combination of multivariate curve resolution with non-linear least square fitting.

The quality of DOSY NMR data can be improved by careful pre-processing techniques. Baseline drift, peak shift, and phase shift commonly exist in real-world DOSY NMR data. These phenomena seriously hinder the data analysis and should be removed as much as possible. In this paper, a series of preprocessing operations are proposed so that the subsequent multivariate curve resolution can yield optimal results. First, the baseline is corrected according to a method by Golotvin and Williams. Next, frequency and phase shift are removed by a new combination of reference deconvolution (FIDDLE), and a method presented by Witjes et al. that can correct several spectra simultaneously. The corrected data are analysed by the combination of multivariate curve resolution with non-linear least square regression (MCR-NLR). The MCR-NLR method turns out to be more robust and leads to better resolution of the pure components than classic MCR.

Journal Article↗

Exploring field vegetation reflectance as an indicator of soil contamination in river floodplains.

This study investigated the relation between vegetation reflectance and elevated concentrations of the metals Ni, Cd, Cu, Zn and Pb in river floodplain soils. High-resolution vegetation reflectance spectra in the visible to near-infrared (400-1350 nm) were obtained using a field radiometer. The relations were evaluated using simple linear regression in combination with two spectral vegetation indices: the Difference Vegetation Index (DVI) and the Red-Edge Position (REP). In addition, a multivariate regression approach using partial least squares (PLS) regression was adopted. The three methods achieved comparable results. The best R(2) values for the relation between metals concentrations and vegetation reflectance were obtained for grass vegetation and ranged from 0.50 to 0.73. Herbaceous species displayed a larger deviation from the established relationships, resulting in lower R(2) values and larger cross-validation errors. The results corroborate the potential of hyperspectral remote sensing to contribute to the survey of elevated metal concentrations in floodplain soils under grassland using the spectral response of the vegetation as an indicator. Additional constraints will, however, have to be taken into account, as results are resolution- and location-dependent.

Environmental Monitoring↗

Powder pattern indexing using the weighted crosscorrelation and genetic algorithms.

X-ray diffraction is a powerful technique for investigating the structure of crystals and crystalline powders. Unfortunately, for powders, the first step in the structure elucidation process, retrieving the unit cell parameters (indexing), is still very critical. In the present article, an improved approach to powder pattern indexing is presented. The proposed method matches peak positions from experimental X-ray powder patterns with peak positions from trial cells using a recently published method for pattern comparison (weighted crosscorrelation). Trial cells are optimized with Genetic Algorithms. Patterns are not pretreated to remove any existing zero point shift, as this is determined during optimization. Another improvement is the peak assignment procedure. This assignment is needed for determining the similarity between lines from trial cells and experiment. It no longer allows calculated peaks to be assigned twice to different experimental peaks, which is beneficial for the indexing process. The procedure proves to be robust with respect to false peaks and accidental or systematic absensences of reflections, and is successfully applied to powder patterns originating from orthorhombic, monoclinic, and triclinic compounds measured with synchrotron as well as with conventional laboratory X-ray diffractometers.

Journal Article↗

Effect on the partial least-squares prediction of yarn properties combining raman and infrared measurements and applying wavelength selection.

The combination of Raman and infrared spectroscopy on the one hand and wavelength selection on the other hand is used to improve the partial least-squares (PLS) prediction of seven selected yarn properties. These properties are important for on-line quality control during production. From 71 yarn samples, the Raman and infrared spectra are measured and reference methods are used to determine the selected properties. Making separate PLS models for all yarn properties using the Raman and infrared spectra, prior to wavelength selection, reveals that Raman spectroscopy outperforms infrared spectroscopy. If wavelength selection is applied, the PLS prediction error decreases and the correlation coefficient increases for all properties. However, a substantial wavelength selection effect is present for the infrared spectra compared to the Raman spectra. For the infrared spectra, wavelength selection results in PLS prediction errors comparable with the prediction performance of the Raman spectra prior to wavelength selection. Concatenating the Raman and infrared spectra does not enhance the PLS prediction performance, not even after wavelength selection. It is concluded that an infrared spectrometer, combined with a wavelength selection procedure, can be used if no (suitable) Raman instrument is available.

Algorithms↗

Automated correction of unwanted phase jumps in reference signals which corrupt MRSI spectra after eddy current correction.

A commonly applied step in the postprocessing of gradient localized proton MR spectroscopy, is correction for eddy current effects using the water signal as a reference. However, this method can degrade some of the metabolite signals, in particular if applied on proton MR spectroscopic imaging data. This artifact arises from the water reference signal in the presence of a second signal which resonates close to the main water resonance. The interference of both resonances will introduce jumps in the phase of the reference time domain signal. Using this phase for eddy current correction will result in a ringing artifact in the frequency domain of the metabolite signal over the whole frequency range. We propose a moving window correction algorithm, which screens the phase of reference signals and removes phase jumps in time domain caused by interference of signals from multiple spin systems. The phase jumps may be abrupt or gradually distributed over several time data points. Because the correction algorithm only corrects time data points which contain phase jumps, the phase is minimally disrupted. Furthermore, the algorithm is automated for large datasets, correcting only those water reference signals which are corrupted. After correction of the corrupted reference signals, normal eddy current correction may be performed. The algorithm is compared with a method which uses a low-pass filter and tested on simulated data as well as on in vivo proton spectroscopic imaging data from a healthy volunteer and from patients with a brain tumor.

Algorithms↗

On the use of 1H and 13C 1D NMR spectra as QSPR descriptors.

Recently, 1D NMR and IR spectra have been proposed as descriptors containing 3D information. And, as such, said to be suitable for making QSAR and QSPR models where 3D molecular geometries matter, for example, in binding affinities. This paper presents a study on the predictive power of 1D NMR spectra-based QSPR models using simulated proton and carbon 1D NMR spectra. It shows that the spectra-based models are outperformed by models based on theoretical molecular descriptors and that spectra-based models are not easy to interpret. We therefore conclude that the use of such NMR spectra offers no added value.

Journal Article↗