Classification Using Generalized Partial Least Squares

The advances in computational biology have made simultaneous monitoring of thousands of features possible. The high throughput technologies not only bring about a much richer information context in which to study various aspects of gene functions but they also present challenge of analyzing data with large number of covariates and few samples. As an integral part of machine learning, classification of samples into two or more categories is almost always of interest to scientists. In this paper, we address the question of classification in this setting by extending partial least squares (PLS), a popular dimension reduction tool in chemometrics, in the context of generalized linear regression based on a previous approach, Iteratively ReWeighted Partial Least Squares, i.e. IRWPLS (Marx, 1996). We compare our results with two-stage PLS (Nguyen and Rocke, 2002A; Nguyen and Rocke, 2002B) and other classifiers. We show that by phrasing the problem in a generalized linear model setting and by applying bias correction to the likelihood to avoid (quasi)separation, we often get lower classification error rates.

Cadaver dogs--a study on detection of contaminated carpet squares

INTRODUCTION: Cadaver dogs are known as valuable forensic tools in crime scene investigations. Scientific research attempting to verify their value is largely lacking, specifically for scents associated with the early postmortem interval. The aim of our investigation was the comparative evaluation of the reliability, accuracy, and specificity of three cadaver dogs belonging to the Hamburg State Police in the detection of scents during the early postmortem interval. MATERIAL AND METHODS: Carpet squares were used as an odor transporting media after they had been contaminated with the scent of two recently deceased bodies (PMI<3h). The contamination occurred for 2 min as well as 10 min without any direct contact between the carpet and the corpse. Comparative searches by the dogs were performed over a time period of 65 days (10 min contamination) and 35 days (2 min contamination). RESULTS: The results of this study indicate that the well-trained cadaver dog is an outstanding tool for crime scene investigation displaying excellent sensitivity (75-100), specificity (91-100), and having a positive predictive value (90-100), negative predictive value (90-100) as well as accuracy (92-100).

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Lipid-poor adenomas on unenhanced CT: does histogram analysis increase sensitivity compared with a mean attenuation threshold?

OBJECTIVE: The purpose of our study was to evaluate the efficacy of CT histogram analysis for further characterization of lipid-poor adenomas on unenhanced CT. MATERIALS AND METHODS: One hundred thirty-two adrenal nodules were identified in 104 patients with lung cancer who underwent PET/CT. Sixty-five nodules were classified as lipid-rich adenomas if they had an unenhanced CT attenuation of less than or equal to 10 H. Thirty-one masses were classified as lipid-poor adenomas if they had an unenhanced CT attenuation greater than 10 H and stability for more than 1 year. Thirty-six masses were classified as lung cancer metastases if they showed rapid growth in 1 year (n = 27) or were biopsy-proven (n = 9). Histogram analysis was performed for all lesions to provide the mean attenuation value and percentage of negative pixels. RESULTS: All lipid-rich adenomas had more than 10% negative pixels; 51.6% of lipid-poor adenomas had more than 10% negative pixels and would have been classified as indeterminate nodules on the basis of mean attenuation alone. None of the metastases had more than 10% negative pixels. Using an unenhanced CT mean attenuation threshold of less than 10 H yielded a sensitivity of 68% and specificity of 100% for the diagnosis of an adenoma. Using an unenhanced CT threshold of more than 10% negative pixels yielded a sensitivity of 84% and specificity of 100% for the diagnosis of an adenoma. CONCLUSION: CT histogram analysis is superior to mean CT attenuation analysis for the evaluation of adrenal nodules and may help decrease referrals for additional imaging or biopsy.

In vivo biochemical 7.0 Tesla magnetic resonance: preliminary results of dGEMRIC, zonal T2, and T2* mapping of articular cartilage

INTRODUCTION: Ultra-high-field whole-body systems (7.0 T) have a high potential for future human in vivo magnetic resonance imaging (MRI). In musculoskeletal MRI, biochemical imaging of articular cartilage may benefit, in particular. Delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) and T2 mapping have shown potential at 3.0 T. Although dGEMRIC, allows the determination of the glycosaminoglycan content of articular cartilage, T2 mapping is a promising tool for the evaluation of water and collagen content. In addition, the evaluation of zonal variation, based on tissue anisotropy, provides an indicator of the nature of cartilage ie, hyaline or hyaline-like articular cartilage.Thus, the aim of our study was to show the feasibility of in vivo dGEMRIC, and T2 and T2* relaxation measurements, at 7.0 T MRI; and to evaluate the potential of T2 and T2* measurements in an initial patient study after matrix-associated autologous chondrocyte transplantation (MACT) in the knee. MATERIALS AND METHODS: MRI was performed on a whole-body 7.0 T MR scanner using a dedicated circular polarization knee coil. The protocol consisted of an inversion recovery sequence for dGEMRIC, a multiecho spin-echo sequence for standard T2 mapping, a gradient-echo sequence for T2* mapping and a morphologic PD SPACE sequence. Twelve healthy volunteers (mean age, 26.7 +/- 3.4 years) and 4 patients (mean age, 38.0 +/- 14.0 years) were enrolled 29.5 +/- 15.1 months after MACT. For dGEMRIC, 5 healthy volunteers (mean age, 32.4 +/- 11.2 years) were included. T1 maps were calculated using a nonlinear, 2-parameter, least squares fit analysis. Using a region-of-interest analysis, mean cartilage relaxation rate was determined as T1 (0) for precontrast measurements and T1 (Gd) for postcontrast gadopentate dimeglumine [Gd-DTPA(2-)] measurements. T2 and T2* maps were obtained using a pixelwise, monoexponential, non-negative least squares fit analysis; region-of-interest analysis was carried out for deep and superficial cartilage aspects. Statistical evaluation was performed by analyses of variance. RESULTS: Mean T1 (dGEMRIC) values for healthy volunteers showed slightly different results for femoral [T1 (0): 1259 +/- 277 ms; T1 (Gd): 683 +/- 141 ms] compared with tibial cartilage [T1 (0): 1093 +/- 281 ms; T1 (Gd): 769 +/- 150 ms]. Global mean T2 relaxation for healthy volunteers showed comparable results for femoral (T2: 56.3 +/- 15.2 ms; T2*: 19.7 +/- 6.4 ms) and patellar (T2: 54.6 +/- 13.0 ms; T2*: 19.6 +/- 5.2 ms) cartilage, but lower values for tibial cartilage (T2: 43.6 +/- 8.5 ms; T2*: 16.6 +/- 5.6 ms). All healthy cartilage sites showed a significant increase from deep to superficial cartilage (P < 0.001). Within healthy cartilage sites in MACT patients, adequate values could be found for T2 (56.6 +/- 13.2 ms) and T2* (18.6 +/- 5.3 ms), which also showed a significant stratification. Within cartilage repair tissue, global mean values showed no difference, with 55.9 +/- 4.9 ms for T2 and 16.2 +/- 6.3 ms for T2*. However, zonal assessment showed only a slight and not significant increase from deep to superficial cartilage (T2: P = 0.174; T2*: P = 0.150). CONCLUSION: In vivo T1 dGEMRIC assessment in healthy cartilage, and T2 and T2* mapping in healthy and reparative articular cartilage, seems to be possible at 7.0 T MRI. For T2 and T2*, zonal variation of articular cartilage could also be evaluated at 7.0 T. This zonal assessment of deep and superficial cartilage aspects shows promising results for the differentiation of healthy and affected articular cartilage. In future studies, optimized protocol selection, and sophisticated coil technology, together with increased signal at ultra-high-field MRI, may lead to advanced biochemical cartilage imaging.

T2 mapping in the knee after microfracture at 3.0 T: correlation of global T2 values and clinical outcome - preliminary results

OBJECTIVE: The aim of our study was to correlate global T2 values of microfracture repair tissue (RT) with clinical outcome in the knee joint. METHODS: We assessed 24 patients treated with microfracture in the knee joint. Magnetic resonance (MR) examinations were performed on a 3T MR unit, T2 relaxation times were obtained with a multi-echo spin-echo technique. T2 maps were obtained using a pixel wise, mono-exponential non-negative least squares fit analysis. Slices covering the cartilage RT were selected and region of interest analysis was done. An individual T2 index was calculated with global mean T2 of the RT and global mean T2 of normal, hyaline cartilage. The Lysholm score and the International Knee Documentation Committee (IKDC) knee evaluation forms were used for the assessment of clinical outcome. Bivariate correlation analysis and a paired, two tailed t test were used for statistics. RESULTS: Global T2 values of the RT [mean 49.8ms, standards deviation (SD) 7.5] differed significantly (P<0.001) from global T2 values of normal, hyaline cartilage (mean 58.5ms, SD 7.0). The T2 index ranged from 61.3 to 101.5. We found the T2 index to correlate with outcome of the Lysholm score (r(s)=0.641, P<0.001) and the IKDC subjective knee evaluation form (r(s)=0.549, P=0.005), whereas there was no correlation with the IKDC knee form (r(s)=-0.284, P=0.179). CONCLUSION: These findings indicate that T2 mapping is sensitive to assess RT function and provides additional information to morphologic MRI in the monitoring of microfracture.

Mean 20-year followup of Bernese periacetabular osteotomy

The goal of the Bernese periacetabular osteotomy is to correct the deficient acetabular coverage in hips with developmental dysplasia to prevent secondary osteoarthrosis. We determined the 20-year survivorship of symptomatic patients treated with this procedure, determined the clinical and radiographic outcomes of the surviving hips, and identified factors predicting poor outcome. We retrospectively evaluated the first 63 patients (75 hips) who underwent periacetabular osteotomy at the institution where this technique was developed. The mean age of the patients at surgery was 29 years (range, 13-56 years), and preoperatively 24% presented with advanced grades of osteoarthritis. Four patients (five hips) were lost to followup and one patient (two hips) died. The remaining 58 patients (68 hips) were followed for a minimum of 19 years (mean, 20.4 years; range, 19-23 years) and 41 hips (60%) were preserved at last followup. The overall mean Merle d'Aubigné and Postel score decreased in comparison to the 10-year value and was similar to the preoperative score. We observed no major changes in any of the radiographic parameters during the 20-year postoperative period except the osteoarthritis score. We identified six factors predicting poor outcome: age at surgery, preoperative Merle d'Aubigné and Postel score, positive anterior impingement test, limp, osteoarthrosis grade, and the postoperative extrusion index. Periacetabular osteotomy is an effective technique for treating symptomatic developmental dysplasia of the hip and can maintain the natural hip at least 19 years in selected patients. LEVEL OF EVIDENCE: Level III, prognostic study.

Parametric combustion modeling for ethanol-gasoline fuelled spark ignition engines

Ethanol-gasoline fuel blends are increasingly being used in spark ignition (SI) engines due to continued growth in renewable fuels as part of a growing renewable portfolio standard (RPS). This leads to the need for a simple and accurate ethanol-gasoline blends combustion model that is applicable to one-dimensional engine simulation. A parametric combustion model has been developed, integrated into an engine simulation tool, and validated using SI engine experimental data. The parametric combustion model was built inside a user compound in GT-Power. In this model, selected burn durations were computed using correlations as functions of physically based non-dimensional groups that have been developed using the experimental engine database over a wide range of ethanol-gasoline blends, engine geometries, and operating conditions. A coefficient of variance (COV) of gross indicated mean effective pressure (IMEP) correlation was also added to the parametric combustion model. This correlation enables the cycle combustion variation modeling as a function of engine geometry and operating conditions. The computed burn durations were then used to fit single and double Wiebe functions. The single-Wiebe parametric combustion compound used the least squares method to compute the single-Wiebe parameters, while the double-Wiebe parametric combustion compound used an analytical solution to compute the double-Wiebe parameters. These compounds were then integrated into the engine model in GT-Power through the multi-Wiebe combustion template in which the values of Wiebe parameters (single-Wiebe or double-Wiebe) were sensed via RLT-dependence. The parametric combustion models were validated by overlaying the simulated pressure trace from GT-Power on to experimentally measured pressure traces. A thermodynamic engine model was also developed to study the effect of fuel blends, engine geometries and operating conditions on both the burn durations and COV of gross IMEP simulation results.

The mean anatomical shape of the tibial plateau at the knee arthroplasty resection level: an investigation using MRI

Clinical assessments after Total Knee Arthroplasty (TKA) show persisting pain after implantation in over 20% of patients. Impingement of soft tissue around the knee, due to imprecise geometry of the tibial implant, can be one reason for persisting ailment. Two hundred and thirty seven MRI scans were evaluated using an active contour detection algorithm (snake) to obtain a high-resolution mean anatomical shape of the tibial plateau. Differences between female and male, older and younger (40) and left and right averaged shapes were determined. The shapes obtained were asymmetric throughout. Absolute differences between the subgroups fell short of inter-individual variations represented by calculated one-sigma confidence intervals. Our results indicate that a differentiation in TKA tibial plateau design by gender, age, or side is of minor relevance.