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.

Posterior Simulation in the Generalized Linear Model with Semiparmetric Random Effects

Generalized linear mixed models with semiparametric random effects are useful in a wide variety of Bayesian applications. When the random effects arise from a mixture of Dirichlet process (MDP) model, normal base measures and Gibbs sampling procedures based on the Pólya urn scheme are often used to simulate posterior draws. These algorithms are applicable in the conjugate case when (for a normal base measure) the likelihood is normal. In the non-conjugate case, the algorithms proposed by MacEachern and Müller (1998) and Neal (2000) are often applied to generate posterior samples. Some common problems associated with simulation algorithms for non-conjugate MDP models include convergence and mixing difficulties. This paper proposes an algorithm based on the Pólya urn scheme that extends the Gibbs sampling algorithms to non-conjugate models with normal base measures and exponential family likelihoods. The algorithm proceeds by making Laplace approximations to the likelihood function, thereby reducing the procedure to that of conjugate normal MDP models. To ensure the validity of the stationary distribution in the non-conjugate case, the proposals are accepted or rejected by a Metropolis-Hastings step. In the special case where the data are normally distributed, the algorithm is identical to the Gibbs sampler.

Gauss-Seidel Estimation of Generalized Linear Mixed Models with Application to Poisson Modeling of Spatially Varying Disease Rates

Generalized linear mixed models (GLMMs) provide an elegant framework for the analysis of correlated data. Due to the non-closed form of the likelihood, GLMMs are often fit by computational procedures like penalized quasi-likelihood (PQL). Special cases of these models are generalized linear models (GLMs), which are often fit using algorithms like iterative weighted least squares (IWLS). High computational costs and memory space constraints often make it difficult to apply these iterative procedures to data sets with very large number of cases. This paper proposes a computationally efficient strategy based on the Gauss-Seidel algorithm that iteratively fits sub-models of the GLMM to subsetted versions of the data. Additional gains in efficiency are achieved for Poisson models, commonly used in disease mapping problems, because of their special collapsibility property which allows data reduction through summaries. Convergence of the proposed iterative procedure is guaranteed for canonical link functions. The strategy is applied to investigate the relationship between ischemic heart disease, socioeconomic status and age/gender category in New South Wales, Australia, based on outcome data consisting of approximately 33 million records. A simulation study demonstrates the algorithm's reliability in analyzing a data set with 12 million records for a (non-collapsible) logistic regression model.

Numerical solutions of elliptic inverse problems via the equation error method

To estimate a parameter in an elliptic boundary value problem, the method of equation error chooses the value that minimizes the error in the PDE and boundary condition (the solution of the BVP having been replaced by a measurement). The estimated parameter converges to the exact value as the measured data converge to the exact value, provided Tikhonov regularization is used to control the instability inherent in the problem. The error in the estimated solution can be bounded in an appropriate quotient norm; estimates can be derived for both the underlying (infinite-dimensional) problem and a finite-element discretization that can be implemented in a practical algorithm. Numerical experiments demonstrate the efficacy and limitations of the method.

Inverse thermodilution with conventional pulmonary artery catheters for the assessment of cerebral, hepatic, renal, and femoral blood flow

Assessment of regional blood flow changes is difficult in the clinical setting. We tested whether conventional pulmonary artery catheters (PACs) can be used to measure regional venous blood flows by inverse thermodilution (ITD). Inverse thermodilution was tested in vitro and in vivo using perivascular ultrasound Doppler (USD) flow probes as a reference. In anesthetized pigs, PACs were inserted in jugular, hepatic, renal, and femoral veins, and their measurements were compared with simultaneous USD flow measurements from carotid, hepatic, renal, and femoral arteries and from portal vein. Fluid boluses were injected through the PAC's distal port, and temperature changes were recorded from the proximally located thermistor. Injectates of 2 and 5 mL at 22 degrees C and 4 degrees C were used. Flows were altered by using a roller pump (in vitro), and infusion of dobutamine and induction of cardiac tamponade, respectively. In vitro: At blood flows between 400 mL . min-1 and 700 mL . min-1 (n = 50), ITD and USD correlated well (r = 0.86, P < 0.0001), with bias and limits of agreement of 3 +/- 101 mL . min-1. In vivo: 514 pairs of measurements had to be excluded from analysis for technical reasons, and 976 were analyzed. Best correlations were r = 0.87 (P < 0.0001) for renal flow and r = 0.46 (P < 0.0001) for hepatic flow. No significant correlation was found for cerebral and femoral flows. Inverse thermodilution using conventional PAC compared moderately well with USD for renal but not for other flows despite good in vitro correlation in various conditions. In addition, this method has significant technical limitations.

SPATIAL UNDERSTANDING OF MATRIX INVERSION FOR INVERSE FORCE ESTIMATION

A basic approach to study a NVH problem is to break down the system in three basic elements – source, path and receiver. While the receiver (response) and the transfer path can be measured, it is difficult to measure the source (forces) acting on the system. It becomes necessary to predict these forces to know how they influence the responses. This requires inverting the transfer path. Singular Value Decomposition (SVD) method is used to decompose the transfer path matrix into its principle components which is required for the inversion. The usual approach to force prediction requires rejecting the small singular values obtained during SVD by setting a threshold, as these small values dominate the inverse matrix. This assumption of the threshold may be subjected to rejecting important singular values severely affecting force prediction. The new approach discussed in this report looks at the column space of the transfer path matrix which is the basis for the predicted response. The response participation is an indication of how the small singular values influence the force participation. The ability to accurately reconstruct the response vector is important to establish a confidence in force vector prediction. The goal of this report is to suggest a solution that is mathematically feasible, physically meaningful, and numerically more efficient through examples. This understanding adds new insight to the effects of current code and how to apply algorithms and understanding to new codes.

Generalized dilating diathesis in patients with popliteal arterial aneurysm

BACKGROUND: Aim of this study was to analyse the relationship between popliteal artery aneurysm (PAA) and generalized arteriomegaly. PATIENTS AND METHODS: In this consecutive serie, thirty-three patients (1 woman, mean age 69.7 +/- 9.6 years) undergoing PAA repair between 1996 and 2000 agreed to participate in a duplex screening program to assess the diameters of the infrarenal abdominal aorta, common and external iliac, common and superficial femoral and contralateral popliteal arteries as well as common carotid and brachial arteries. RESULTS: The prevalence of arteriomegaly and aneurysmal disease, respectively, was as follows: abdominal aorta 15/33 (45.5%) and 8/33 (24.2%), common iliac artery 34/66 (51.5%) and 23/66 (34.8%), common femoral artery 55/66 (83.3%) and 7/66 (10.6%) as well as contralateral popliteal artery 7/33 (21.2%) 15/33 (45.5%). Significantly larger carotid artery diameters were found comparing PAA patients with age- and body surface adjusted healthy controls (p < 0.001). Furthermore, patients with multiple peripheral arterial aneurysms had significantly larger diameters of the brachial (p < 0.02) and external iliac (p < 0.005). CONCLUSIONS: Our findings support the hypothesis of a diathesis for a generalized arteriomegaly with a predilection for further aneurysms of the abdominal aorta, iliac arteries, femoral and contralateral popliteal arteries in patients with PAA.

Chronic stimulation of the globus pallidus internus for treatment of non-dYT1 generalized dystonia and choreoathetosis: 2-year follow up

OBJECT: The authors studied the long-term efficacy of deep brain stimulation (DBS) of the posteroventral lateral globus pallidus internus up to 2 years postoperatively in patients with primary non-DYT1 generalized dystonia or choreoathetosis. The results are briefly compared with those reported for DBS in DYT1 dystonia (Oppenheim dystonia), which is caused by the DYT1 gene. METHODS: Enrollment in this prospective expanded pilot study was limited to adult patients with severely disabling, medically refractory non-DYT1 generalized dystonia or choreoathetosis. Six consecutive patients underwent follow-up examinations at defined intervals of 3 months, 1 year, and 2 years postsurgery. There were five women and one man, and their mean age at surgery was 45.5 years. Formal assessments included both the Burke-Fahn-Marsden dystonia scale and the recently developed Unified Dystonia Rating Scale. Two patients had primary generalized non-DYT1 dystonia, and four suffered from choreoathetosis secondary to infantile cerebral palsy. Bilateral quadripolar DBS electrodes were implanted in all instances, except in one patient with markedly asymmetrical symptoms. There were no adverse events related to surgery. The Burke-Fahn-Marsden scores in the two patients with generalized dystonia improved by 78 and 71% at 3 months, by 82 and 69% at 1 year, and by 78 and 70% at 2 years postoperatively. This was paralleled by marked amelioration of disability scores. The mean improvement in Burke-Fahn-Marsden scores in patients with choreoathetosis was 12% at 3 months, 29% at 1 year, and 23% at 2 years postoperatively, which was not significant. Two of these patients thought that they had achieved marked improvement at 2 years postoperatively, although results of objective evaluations were less impressive. In these two patients there was a minor but stable improvement in disability scores. All patients had an improvement in pain scores at the 2-year follow-up review. Medication was tapered off in both patients with generalized dystonia and reduced in two of the patients with choreoathetosis. All stimulation-induced side effects were reversible on adjustment of the DBS settings. Energy consumption of the batteries was considerably higher than in patients with Parkinson disease. CONCLUSIONS: Chronic pallidal DBS is a safe and effective procedure in generalized non-DYT1 dystonia, and it may become the procedure of choice in patients with medically refractory dystonia. Postoperative improvement of choreoathetosis is more modest and varied, and subjective ratings of outcome may exceed objective evaluations.