Permutation Inference for Canonical Correlation Analysis

Canonical correlation analysis (CCA) has become a key tool for population neuroimaging, allowing investigation of associations between many imaging and non-imaging measurements. As other variables are often a source of variability not of direct interest, previous work has used CCA on residuals from a model that removes these effects, then proceeded directly to permutation inference. We show that such a simple permutation test leads to inflated error rates. The reason is that residualisation introduces dependencies among the observations that violate the exchangeability assumption. Even in the absence of nuisance variables, however, a simple permutation test for CCA also leads to excess error rates for all canonical correlations other than the first. The reason is that a simple permutation scheme does not ignore the variability already explained by previous canonical variables. Here we propose solutions for both problems: in the case of nuisance variables, we show that transforming the residuals to a lower dimensional basis where exchangeability holds results in a valid permutation test; for more general cases, with or without nuisance variables, we propose estimating the canonical correlations in a stepwise manner, removing at each iteration the variance already explained, while dealing with different number of variables in both sides. We also discuss how to address the multiplicity of tests, proposing an admissible test that is not conservative, and provide a complete algorithm for permutation inference for CCA.Comment: 49 pages, 2 figures, 10 tables, 3 algorithms, 119 reference

Increasing power for voxel-wise genome-wide association studies : the random field theory, least square kernel machines and fast permutation procedures

Imaging traits are thought to have more direct links to genetic variation than diagnostic measures based on cognitive or clinical assessments and provide a powerful substrate to examine the influence of genetics on human brains. Although imaging genetics has attracted growing attention and interest, most brain-wide genome-wide association studies focus on voxel-wise single-locus approaches, without taking advantage of the spatial information in images or combining the effect of multiple genetic variants. In this paper we present a fast implementation of voxel- and cluster-wise inferences based on the random field theory to fully use the spatial information in images. The approach is combined with a multi-locus model based on least square kernel machines to associate the joint effect of several single nucleotide polymorphisms (SNP) with imaging traits. A fast permutation procedure is also proposed which significantly reduces the number of permutations needed relative to the standard empirical method and provides accurate small p-value estimates based on parametric tail approximation. We explored the relation between 448,294 single nucleotide polymorphisms and 18,043 genes in 31,662 voxels of the entire brain across 740 elderly subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Structural MRI scans were analyzed using tensor-based morphometry (TBM) to compute 3D maps of regional brain volume differences compared to an average template image based on healthy elderly subjects. We find method to be more sensitive compared with voxel-wise single-locus approaches. A number of genes were identified as having significant associations with volumetric changes. The most associated gene was GRIN2B, which encodes the N-methyl-d-aspartate (NMDA) glutamate receptor NR2B subunit and affects both the parietal and temporal lobes in human brains. Its role in Alzheimer's disease has been widely acknowledged and studied, suggesting the validity of the approach. The various advantages over existing approaches indicate a great potential offered by this novel framework to detect genetic influences on human brains

In My View

I am pleased that in his review of my book, The Russian Presidency (Naval War College Review, Summer 2001, pp. 165–8), Sergei Khrushchev agreed with me (and we are among a very few in the scholarly community who have agreed on this issue) that presidentialism will serve Russia better than the parliamentary alternative. But elsewhere, Dr. Khrushchev disagrees with me, and I am writing not only to clarify to readers of the Review the nature of that disagreement but also to take issue with places where I believe he has mischaracterized the book

Dynamic filtering of static dipoles in magnetoencephalography

We consider the problem of estimating neural activity from measurements of the magnetic fields recorded by magnetoencephalography. We exploit the temporal structure of the problem and model the neural current as a collection of evolving current dipoles, which appear and disappear, but whose locations are constant throughout their lifetime. This fully reflects the physiological interpretation of the model. In order to conduct inference under this proposed model, it was necessary to develop an algorithm based around state-of-the-art sequential Monte Carlo methods employing carefully designed importance distributions. Previous work employed a bootstrap filter and an artificial dynamic structure where dipoles performed a random walk in space, yielding nonphysical artefacts in the reconstructions; such artefacts are not observed when using the proposed model. The algorithm is validated with simulated data, in which it provided an average localisation error which is approximately half that of the bootstrap filter. An application to complex real data derived from a somatosensory experiment is presented. Assessment of model fit via marginal likelihood showed a clear preference for the proposed model and the associated reconstructions show better localisation

A Unified Approach to Attractor Reconstruction

In the analysis of complex, nonlinear time series, scientists in a variety of disciplines have relied on a time delayed embedding of their data, i.e. attractor reconstruction. The process has focused primarily on heuristic and empirical arguments for selection of the key embedding parameters, delay and embedding dimension. This approach has left several long-standing, but common problems unresolved in which the standard approaches produce inferior results or give no guidance at all. We view the current reconstruction process as unnecessarily broken into separate problems. We propose an alternative approach that views the problem of choosing all embedding parameters as being one and the same problem addressable using a single statistical test formulated directly from the reconstruction theorems. This allows for varying time delays appropriate to the data and simultaneously helps decide on embedding dimension. A second new statistic, undersampling, acts as a check against overly long time delays and overly large embedding dimension. Our approach is more flexible than those currently used, but is more directly connected with the mathematical requirements of embedding. In addition, the statistics developed guide the user by allowing optimization and warning when embedding parameters are chosen beyond what the data can support. We demonstrate our approach on uni- and multivariate data, data possessing multiple time scales, and chaotic data. This unified approach resolves all the main issues in attractor reconstruction.Comment: 22 pages, revised version as submitted to CHAOS. Manuscript is currently under review. 4 Figures, 31 reference

False Discovery Rate and Localizing Power

False discovery rate (FDR) is commonly used for correction for multiple testing in neuroimaging studies. However, when using two-tailed tests, making directional inferences about the results can lead to vastly inflated error rate, even approaching 100% in some cases. This happens because FDR only provides weak control over the error rate, meaning that the proportion of error is guaranteed only globally over all tests, not within subsets, such as among those in only one or another direction. Here we consider and evaluate different strategies for FDR control with two-tailed tests, using both synthetic and real imaging data. Approaches that separate the tests by direction of the hypothesis test, or by the direction of the resulting test statistic, more properly control the directional error rate and preserve FDR benefits, albeit with a doubled risk of errors under complete absence of signal. Strategies that combine tests in both directions, or that use simple two-tailed p-values, can lead to invalid directional conclusions, even if these tests remain globally valid. To enable valid thresholding for directional inference, we suggest that imaging software should allow the possibility that the user sets asymmetrical thresholds for the two sides of the statistical map. While FDR continues to be a valid, powerful procedure for multiple testing correction, care is needed when making directional inferences for two-tailed tests, or more broadly, when making any localized inference

Invasive Haemophilus influenzae Disease in Adults ≥65 Years, United States, 2011.

BackgroundSince the introduction of the Haemophilus influenzae serotype b vaccine, H influenzae epidemiology has shifted. In the United States, the largest burden of disease is now in adults aged ≥65 years. However, few data exist on risk factors for disease severity and outcome in this age group.MethodsA retrospective case-series review of invasive H influenzae infections in patients aged ≥65 years was conducted for hospitalized cases reported to Active Bacterial Core surveillance in 2011.ResultsThere were 299 hospitalized cases included in the analysis. The majority of cases were caused by nontypeable H influenzae, and the overall case fatality ratio (CFR) was 19.5%. Three or more underlying conditions were present in 63% of cases; 94% of cases had at least 1. Patients with chronic heart conditions (congestive heart failure, coronary artery disease, and/or atrial fibrillation) (odds ratio [OR], 3.27; 95% confidence interval [CI], 1.65-6.46), patients from private residences (OR, 8.75; 95% CI, 2.13-35.95), and patients who were not resuscitate status (OR, 2.72; 95% CI, 1.31-5.66) were more likely to be admitted to the intensive care unit (ICU). Intensive care unit admission (OR, 3.75; 95% CI, 1.71-8.22) and do not resuscitate status (OR, 12.94; 95% CI, 4.84-34.55) were significantly associated with death.ConclusionsWithin this age group, burden of disease and CFR both increased significantly as age increased. Using ICU admission as a proxy for disease severity, our findings suggest several conditions increased risk of disease severity and patients with severe disease were more likely to die. Further research is needed to determine the most effective approach to prevent H influenzae disease and mortality in older adults

Occurrence and Reproduction of the Alabama Shad, Alosa alabamae Jordan and Evermann, in the Ouachita River System of Arkansas

The anadromous Alabama shad, Alosa alabamae, has drastically declined in abundance in recent decades throughout its historic range and has previously been reported in Arkansas from only five localities. Three of those locality records are pre- 1900. Sampling by seine in the Ouachita River drainage system of southern Arkansas in July and August of 1997 and 1998 produced more than 300 juvenile A. alabamae from two localities on the Little Missouri River and four localities on the Ouachita River. One record of an adult Alabama shad, taken on 4 April1997 by an angler below Remmel Dam on the Ouachita River, was also documented. Adults apparently ascended the Ouachita River and spawned successfully in 1997 and 1998 despite the construction of four locks and dams on that river in Louisiana and Arkansas in the 1980s. The Ouachita River drainage and a few streams in east-central Missouri are currently the only known spawning areas for A. alabamae in noncoastal regions of the entire Mississippi River basin. Continued survival of the Alabama shad in Arkansas depends on protecting critical spawning and nursery habitats in the Ouachita River system from deleterious alteration and on preserving that migratory species\u27 access to those habitats

Meningococcal Disease in Patients With Human Immunodeficiency Virus Infection: A Review of Cases Reported Through Active Surveillance in the United States, 2000-2008.

BackgroundAlthough human immunodeficiency virus (HIV) infection is an established risk factor for several bacterial infections, the association between HIV infection and meningococcal disease remains unclear.MethodsExpanded chart reviews were completed on persons with meningococcal disease and HIV infection reported from 2000 through 2008 from 9 US sites participating in an active population-based surveillance system for meningococcal disease. The incidence of meningococcal disease among patients meeting Centers for Disease Control and Prevention acquired immune deficiency syndrome (AIDS) surveillance criteria was estimated using data from the National HIV Surveillance System for the participating sites.ResultsThirty-three cases of meningococcal disease in individuals with HIV infection were reported from participating sites, representing 2.0% of all reported meningococcal disease cases. Most (75.8%) persons with HIV infection were adult males aged 25 to 64 years old. Among all meningococcal disease cases aged 25 to 64 years old, case fatality ratios were similar among HIV-infected and HIV-uninfected persons (13.3% vs 10.6%; P = .6). The cumulative, mean incidence of meningococcal disease among patients aged 25 to 64 years old with HIV infection ever classified as AIDS was 3.5 cases per 100000 person years (95% confidence interval [CI], 2.1-5.6), compared with 0.3 cases per 100000 person years (95% CI, 0.3-0.3) for persons of the same age group not reported to have AIDS (relative risk = 12.9; 95% CI, 7.9-20.9).ConclusionsIndividuals with HIV infection meeting the AIDS surveillance case definition have a higher incidence of meningococcal disease compared with the general adult population

Metagenetic analysis of patterns of distribution and diversity of marine meiobenthic eukaryotes

AimMeiofaunal communities that inhabit the marine benthos offer unique opportunities to simultaneously study the macroecology of numerous phyla that exhibit different life-history strategies. Here, we ask: (1) if the macroecology of meiobenthic communities is explained mainly by dispersal constraints or by environmental conditions; and (2) if levels of meiofaunal diversity surpass existing estimates based on morphological taxonomy. LocationUK and mainland European coast. MethodsNext-generation sequencing techniques (NGS; Roche 454 FLX platform) using 18S nuclear small subunit ribosomal DNA (rDNA) gene. Pyrosequences were analysed using AmpliconNoise followed by chimera removal using Perseus. ResultsRarefaction curves revealed that sampling saturation was only reached at 15% of sites, highlighting that the bulk of meiofaunal diversity is yet to be discovered. Overall, 1353 OTUs were recovered and assigned to 23 different phyla. The majority of sampled sites had c. 60-70 unique operational taxonomic units (OTUs) per site, indicating high levels of beta diversity. The environmental parameters that best explained community structure were seawater temperature, geographical distance and sediment size, but most of the variability (R-2=70%-80%) remains unexplained. Main conclusionsHigh percentages of endemic OTUs suggest that meiobenthic community composition is partly niche-driven, as observed in larger organisms, but also shares macroecological features of microorganisms by showing high levels of cosmopolitanism (albeit on a much smaller scale). Meiobenthic communities exhibited patterns of isolation by distance as well as associations between niche, latitude and temperature, indicating that meiobenthic communities result from a combination of niche assembly and dispersal processes. Conversely, isolation-by-distance patterns were not identified in the featured protists, suggesting that animals and protists adhere to radically different macroecological processes, linked to life-history strategies.Natural Environment Research Council (NERC) [NE/E001505/1, NE/F001266/1, MGF-167]; Portuguese Foundation for Science and Technology (FCT) [SFRH/BD/27413/2006, SFRH/BPD/80447/2014]; EPSRC [EP/H003851/1]; BBSRC CASE studentship; Unilever; Biotechnology and Biological Sciences Research Council [987347]; Engineering and Physical Sciences Research Council [EP/H003851/1]; Natural Environment Research Council [NE/F001290/1, NE/F001266/1, NE/E001505/1, NBAF010002]info:eu-repo/semantics/publishedVersio