Asymptotics for Marginal Generalized Linear Models With Sparse Correlations

Marginal generalized linear models can be used for clustered and longitudinal data by fitting a model as if the data were independent and using an empirical estimator of parameter standard errors. We extend this approach to data where the number of observations correlated with a given one grows with sample size and show that parameter estimates are consistent and asymptotically Normal with a slower convergence rate than for independent data, and that an information sandwich variance estimator is consistent. We present two problems that motivated this work, the modelling of patterns of HIV genetic variation and the behavior of clustered data estimators when clusters are large.

Point mutation tRNA(Ser(UCN)) in a child with hearing loss and myoclonus epilepsy

We report on a family with a 12-year-old boy who suffered from a maternally inherited syndrome characterized by a combination of sensorineural hearing loss, myoclonus epilepsy, ataxia, severe psychomotor retardation, short stature, and diabetes mellitus. First, he showed a muscular hypotonia with hearing loss; later, he developed a myoclonus epilepsy, growth failure, and severe psychomotor retardation. At the age of 10 years, he developed diabetes mellitus. After initiation of combined ubiquinone and vitamin C treatment, we observed a progression in psychomotor development. Lactate and pyruvate levels in blood and cerebrospinal fluid were normal. No ragged red fibers or ultrastructural abnormalities were seen in a skeletal muscle biopsy. Biochemical assays of respiratory chain complex activities revealed decreased activity of complexes I and IV. By sequence analysis of mitochondrial DNA encoding transfer ribonucleic acids (RNAs), a homoplasmic T to C substitution at nucleotide position 7512 was found affecting a highly conserved base pair in the tRNA(ser(UCN)) acceptor stem. Asymptomatic family members of the maternal line were heteroplasmic for the mutation in blood samples. Analysis of mitochondrial DNA in patients with hearing loss and myoclonus epilepsy is recommended, even in the absence of laboratory findings. Therapeutically, ubiquinone and antioxidants can be beneficial.

Spastin in the human and mouse central nervous system with special reference to its expression in the hippocampus of mouse pilocarpine model of status epilepticus and temporal lobe epilepsy

In the present in situ hybridization and immunocytochemical studies in the mouse central nervous system (CNS), a strong expression of spastin mRNA and protein was found in Purkinje cells and dentate nucleus in the cerebellum, in hippocampal principal cells and hilar neurons, in amygdala, substantia nigra, striatum, in the motor nuclei of the cranial nerves and in different layers of the cerebral cortex except piriform and entorhinal cortices where only neurons in layer II were strongly stained. Spastin protein and mRNA were weakly expressed in most of the thalamic nuclei. In selected human brain regions such as the cerebral cortex, cerebellum, hippocampus, amygdala, substania nigra and striatum, similar results were obtained. Electron microscopy showed spastin immunopositive staining in the cytoplasma, dendrites, axon terminals and nucleus. In the mouse pilocarpine model of status epilepticus and subsequent temporal lobe epilepsy, spastin expression disappeared in hilar neurons as early as at 2h during pilocarpine induced status epilepticus, and never recovered. At 7 days and 2 months after pilocarpine induced status epilepticus, spastin expression was down-regulated in granule cells in the dentate gyrus, but induced expression was found in reactive astrocytes. The demonstration of widespread distribution of spastin in functionally different brain regions in the present study may provide neuroanatomical basis to explain why different neurological, psychological disorders and cognitive impairment occur in patients with spastin mutation. Down-regulation or loss of spastin expression in hilar neurons may be related to their degeneration and may therefore initiate epileptogenetic events, leading to temporal lobe epilepsy.

Reorganization of CA3 area of the mouse hippocampus after pilocarpine induced temporal lobe epilepsy with special reference to the CA3-septum pathway

We showed that when CA3 pyramidal neurons in the caudal 80% of the dorsal hippocampus had almost disappeared completely, the efferent pathway of CA3 was rarely detectable. We used the mouse pilocarpine model of temporal lobe epilepsy (TLE), and injected iontophoretically the anterograde tracer phaseolus vulgaris leucoagglutinin (PHA-L) into gliotic CA3, medial septum and the nucleus of diagonal band of Broca, median raphe, and lateral supramammillary nuclei, or the retrograde tracer cholera toxin B subunit (CTB) into gliotic CA3 area of hippocampus. In the afferent pathway, the number of neurons projecting to CA3 from medial septum and the nucleus of diagonal band of Broca, median raphe, and lateral supramammillary nuclei increased significantly. In the hippocampus, where CA3 pyramidal neurons were partially lost, calbindin, calretinin, parvalbumin immunopositive back-projection neurons from CA1-CA3 area were observed. Sprouting of Schaffer collaterals with increased number of large boutons in both sides of CA1 area, particularly in the stratum pyramidale, was found. When CA3 pyramidal neurons in caudal 80% of the dorsal hippocampus have almost disappeared completely, surviving CA3 neurons in the rostral 20% of the dorsal hippocampus may play an important role in transmitting hyperactivity of granule cells to surviving CA1 neurons or to dorsal part of the lateral septum. We concluded that reorganization of CA3 area with its downstream or upstream nuclei may be involved in the occurrence of epilepsy.

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.

Differential effect of side of temporal lobe epilepsy on lateralization of hippocampal, temporolateral, and inferior frontal activation patterns during a verbal episodic memory task

The encoding of verbal stimuli elicits left-lateralized activation patterns within the medial temporal lobes in healthy adults. In our study, patients with left- and right-sided temporal lobe epilepsy (LTLE, RTLE) were investigated during the encoding and retrieval of word-pair associates using functional magnetic resonance imaging. Functional asymmetry of activation patterns in hippocampal, inferior frontal, and temporolateral neocortical areas associated with language functions was analyzed. Hippocampal activation patterns in patients with LTLE were more right-lateralized than those in patients with RTLE (P<0.05). There were no group differences with respect to lateralization in frontal or temporolateral regions of interest (ROIs). For both groups, frontal cortical activation patterns were significantly more left-lateralized than hippocampal patterns (P<0.05). For patients with LTLE, there was a strong trend toward a difference in functional asymmetry between the temporolateral and hippocampal ROIs (P=0.059). A graded effect of epileptic activity on laterality of the different regional activation patterns is discussed.

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.