Accounting for spatial dependency in multivariate spectroscopic data

We examine a hybrid multivariate regression technique to account for the spatial dependency in spectroscopic data due to adjacent measurement locations in the same joint by combining dimension reduction methods and linear mixed effects (LME) modeling. Spatial correlation is a common limitation (assumption of independence) encountered in diagnostic applications involving adjacent measurement locations, such as mapping of tissue properties, and can impede tissue evaluations. Near-infrared spectra were collected from equine joints (n = 5) and corresponding biomechanical (n = 202), compositional (n = 530), and structural (n = 530) properties of cartilage tissue were measured. Subsequently, hybrid regression models for estimating tissue properties from the spectral data were developed in combination with principal component analysis (PCA-LME) scores and least absolute shrinkage and selection operator (LASSO-LME). Performance comparison of PCA-LME and principal component regression, and LASSO-LME and LASSO regression was conducted to evaluate the effects of spatial dependency. A systematic improvement in calibration models’ correlation coefficients and a decrease in cross validation errors were observed when accounting for spatial dependency. Our results indicate that accounting for spatial dependency using a LME-based approach leads to more accurate prediction models

Optimal regression method for near-infrared spectroscopic evaluation of articular cartilage

Abstract Near-infrared (NIR) spectroscopy has been successful in nondestructive assessment of biological tissue properties, such as stiffness of articular cartilage, and is proposed to be used in clinical arthroscopies. Near-infrared spectroscopic data include absorbance values from a broad wavelength region resulting in a large number of contributing factors. This broad spectrum includes information from potentially noisy variables, which may contribute to errors during regression analysis. We hypothesized that partial least squares regression (PLSR) is an optimal multivariate regression technique and requires application of variable selection methods to further improve the performance of NIR spectroscopy-based prediction of cartilage tissue properties, including instantaneous, equilibrium, and dynamic moduli and cartilage thickness. To test this hypothesis, we conducted for the first time a comparative analysis of multivariate regression techniques, which included principal component regression (PCR), PLSR, ridge regression, least absolute shrinkage and selection operator (Lasso), and least squares version of support vector machines (LS-SVM) on NIR spectral data of equine articular cartilage. Additionally, we evaluated the effect of variable selection methods, including Monte Carlo uninformative variable elimination (MC-UVE), competitive adaptive reweighted sampling (CARS), variable combination population analysis (VCPA), backward interval PLS (BiPLS), genetic algorithm (GA), and jackknife, on the performance of the optimal regression technique. The PLSR technique was found as an optimal regression tool ($R^2_{Tissue~thickness} = 75.6\%$, $R^2_{Dynamic~modulus} = 64.9\%$) for cartilage NIR data; variable selection methods simplified the prediction models enabling the use of lesser number of regression components. However, the improvements in model performance with variable selection methods were found to be statistically insignificant. Thus, the PLSR technique is recommended as the regression tool for multivariate analysis for prediction of articular cartilage properties from its NIR spectra

Corrigendum to "Multimodality scoring of chondral injuries in the equine fetlock joint ex vivo" [Osteoarthritis Cartilage 25 (5) (2017 May) 790-798]

The authors have found a systematic error in the Young's modulus values in the manuscript “Multimodality scoring of chondral injuries in the equine fetlock joint ex vivo” published in the May 2017 issue of Osteoarthritis and Cartilage; 25(5):790–98. The formula used to calculate the instantaneous modulus values [Formula presented], and not for Young's modulus [Formula presented] as intended. As a result of this error, the Artscan modulus values [Formula presented] presented in Table 1 are exactly one third of the correct values. This also scales the y-axes of the figures 3C and 3D. This systematic error has no effect on the conclusions, discussion, or any other results/statistics presented in the article. The correct form of the equation with the corrected descriptions is as follows: “Young's modulus was determined based on the Hayes' elastic model of indentation: [Formula presented] where [Formula presented] is the measured indenter force, [Formula presented] is the Poisson's ratio [Formula presented], is the indenter radius of curvature, [Formula presented] is the radius of the indenter, and [Formula presented] and [Formula presented] are theoretical correction factors.” Table I (Corrected). Average ICRS scores (N = 43) and instantaneous moduli [Formula presented] for both surgeons, including SD and 95% CI of each round. Corresponding gold standard values for average histology-based ICRS score and laboratory mechanical testing system based instantaneous modulus [Formula presented] is observed between ICRS grades 0 and 1 based on the average score of multimodal scorings. (C-D) A similar trend is not apparent with Artscan measurements [Formula presented] and ICRS scores (histology) or the average score of multimodal scorings. A single measurement point is not visible (83.1 MPa, at ICRS 1 and ICRS 0) in subfigures C and D, respectively

Multimodality scoring of chondral injuries in the equine fetlock joint ex vivo

Objective: We investigate the potential of a prototype multimodality arthroscope, combining ultrasound, optical coherence tomography (OCT) and arthroscopic indentation device, for assessing cartilage lesions, and compare the reliability of this approach with conventional arthroscopic scoring ex vivo. Design: Areas of interest (AIs, N = 43) were selected from equine fetlock joints (N = 5). Blind-coded AIs were independently scored by two equine surgeons employing International Cartilage Repair Society (ICRS) scoring system via conventional arthroscope and multimodality arthroscope, in which high-frequency ultrasound and OCT catheters were attached to an arthroscopic indentation device. In addition, cartilage stiffness was measured with the indentation device, and lesions in OCT images scored using custom-made automated software. Measurements and scorings were performed twice in two separate rounds. Finally, the scores were compared to histological ICRS scores. Results: OCT and arthroscopic examinations showed the highest average agreements (55.2%) between the scoring by surgeons and histology scores, whereas ultrasound had the lowest (50.6%). Average intraobserver agreements of surgeons and interobserver agreements between rounds were, respectively, for conventional arthroscope (68.6%, 69.8%), ultrasound (68.6%, 68.6%), OCT (65.1%, 61.7%) and automated software (65.1%, 59.3%). Conclusions: OCT imaging supplemented with the automated software provided the most reliable lesion scoring. However, limited penetration depth of light limits the clinical potential of OCT in assessing human cartilage thickness; thus, the combination of OCT and ultrasound could be optimal for reliable diagnostics. Present findings suggest imaging and quantitatively analyzing the entire articular surface to eliminate surgeon-related variation in the selection of the most severe lesion to be scored

Combination of optical coherence tomography and near infrared spectroscopy enhances determination of articular cartilage composition and structure

Abstract Conventional arthroscopic evaluation of articular cartilage is subjective and poorly reproducible. Therefore, implementation of quantitative diagnostic techniques, such as near infrared spectroscopy (NIRS) and optical coherence tomography (OCT), is essential. Locations (n = 44) with various cartilage conditions were selected from mature equine fetlock joints (n = 5). These locations and their surroundings were measured with NIRS and OCT (n = 530). As a reference, cartilage proteoglycan (PG) and collagen contents, and collagen network organization were determined using quantitative microscopy. Additionally, lesion severity visualized in OCT images was graded with an automatic algorithm according to International Cartilage Research Society (ICRS) scoring system. Artificial neural network with variable selection was then employed to predict cartilage composition in the superficial and deep zones from NIRS data, and the performance of two models, generalized (including all samples) and condition-specific models (based on ICRS-grades), was compared. Spectral data correlated significantly (p < 0.002) with PG and collagen contents, and collagen orientation in the superficial and deep zones. The combination of NIRS and OCT provided the most reliable outcome, with condition-specific models having lower prediction errors (9.2%) compared to generalized models (10.4%). Therefore, the results highlight the potential of combining both modalities for comprehensive evaluation of cartilage during arthroscopy

Dataset on equine cartilage near infrared spectra, composition, and functional properties

Abstract Near infrared (NIR) spectroscopy is a well-established technique that is widely employed in agriculture, chemometrics, and pharmaceutical engineering. Recently, the technique has shown potential in clinical orthopaedic applications, for example, assisting in the diagnosis of various knee-related diseases (e.g., osteoarthritis) and their pathologies. NIR spectroscopy (NIRS) could be especially useful for determining the integrity and condition of articular cartilage, as the current arthroscopic diagnostics is subjective and unreliable. In this work, we present an extensive dataset of NIRS measurements for evaluating the condition, mechanical properties, structure, and composition of equine articular cartilage. The dataset contains NIRS measurements from 869 different locations across the articular surfaces of five equine fetlock joints. A comprehensive library of reference values for each measurement location is also provided, including results from a mechanical indentation testing, digital densitometry imaging, polarized light microscopy, and Fourier transform infrared spectroscopy. The published data can either be used as a model of human cartilage or to advance equine veterinary research