Statistical Network Analysis for Functional MRI: Summary Networks and Group Comparisons

Comparing weighted networks in neuroscience is hard, because the topological properties of a given network are necessarily dependent on the number of edges of that network. This problem arises in the analysis of both weighted and unweighted networks. The term density is often used in this context, in order to refer to the mean edge weight of a weighted network, or to the number of edges in an unweighted one. Comparing families of networks is therefore statistically difficult because differences in topology are necessarily associated with differences in density. In this review paper, we consider this problem from two different perspectives, which include (i) the construction of summary networks, such as how to compute and visualize the mean network from a sample of network-valued data points; and (ii) how to test for topological differences, when two families of networks also exhibit significant differences in density. In the first instance, we show that the issue of summarizing a family of networks can be conducted by adopting a mass-univariate approach, which produces a statistical parametric network (SPN). In the second part of this review, we then highlight the inherent problems associated with the comparison of topological functions of families of networks that differ in density. In particular, we show that a wide range of topological summaries, such as global efficiency and network modularity are highly sensitive to differences in density. Moreover, these problems are not restricted to unweighted metrics, as we demonstrate that the same issues remain present when considering the weighted versions of these metrics. We conclude by encouraging caution, when reporting such statistical comparisons, and by emphasizing the importance of constructing summary networks.Comment: 16 pages, 5 figure

On the theta number of powers of cycle graphs

We give a closed formula for Lovasz theta number of the powers of cycle graphs and of their complements, the circular complete graphs. As a consequence, we establish that the circular chromatic number of a circular perfect graph is computable in polynomial time. We also derive an asymptotic estimate for this theta number.Comment: 17 page

Taxonomy of Phanaeus revisited: revised keys to and comments on species of the New World dung beetle genus PhanaeusMacLeay, 1819 (Coleoptera: Scarabaeidae: Scarabaeinae: Phanaeini)

The purpose of this paper is to reassess the taxonomy of Phanaeus MacLeay (Coleoptera: Scarabaeidae) to accommodate new taxa and changes in taxonomic opinion since the publication of Edmonds’ 1994 revision of the genus. The two subgenera and 13 species groups established by Edmonds (1994) remain unchanged. A revised set of keys with accompanying comments and illustrations separates 54 recognized valid species. Seven recently described valid species are incorporated into the revised classification: Phanaeus blackalleri Delgado-Castillo, 1991; P. bordoni Arnaud, 1996; P. changdiazi Kohlmann and Solís, 2001; P. lecourti Arnaud, 2000; P. martinezorum Arnaud, 2000; P. yecoraensis Edmonds, 2004; and P. zapotecus Edmonds, 2006. The new name Phanaeus sororibispinus Edmonds and Zidek replaces Phanaeus alvarengai Arnaud, 1984, a primary junior homonym of P. alvarengai Pereira and d’Andretta, 1955. Three subspecies recognized in 1994 are elevated to species rank, new status: Phanaeus texensis Edmonds, 1994; P. pilatei Harold, 1863; and P. guatemalensis Harold, 1871. Phanaeus obliquans Bates, 1887 is removed from synonymy and given new status as a valid species. Twelve new junior subjective synonyms (bold) are recognized: P. tridens balthasari Arnaud, 2002 (of P. tridens Castelnau, 1840); P. dzidoi Arnaud, 2000 (of P. palaeno Blanchard, 1843); P. genieri Arnaud, 2002 (of P. amethystinus Harold, 1863); P. prasinus jolyi Arnaud, 2001 (of P. prasinus Harold, 1868); P. kirbyi ledezmai Arnaud, 2002 (of P. kirbyi Vigors, 1825); P. achilles lydiae Arnaud, 2000 (of P. achilles Boheman, 1858); P. chalcomelas grossii Arnaud, 2001 (of P. chalcomelas [Perty, 1830]); P pyrois malyi Arnaud, 2002 (of P. pyrois Bates, 1887); P. tridens moroni Arnaud, 2001 (of P. tridens Castenau, 1840); P. lecourti peruanus Arnaud, 2000 (of P. lecourti Arnaud, 2000); P. endymion porioni Arnaud, 2001 (of P. endymion Harold, 1863); P. pseudofurcosus Balthasar, 1939 (of P. tridens Castelnau, 1840); and P. prasinus trinidadensis Arnaud, 2001 (of P. prasinus Harold, 1868). “Phanaeus viridicollis” Olsoufieff, 1924 (sensu Arnaud 2002) is an unavailable name here considered a color variant of P. pyrois Bates, 1887

Group Analysis of Self-organizing Maps based on Functional MRI using Restricted Frechet Means

Studies of functional MRI data are increasingly concerned with the estimation of differences in spatio-temporal networks across groups of subjects or experimental conditions. Unsupervised clustering and independent component analysis (ICA) have been used to identify such spatio-temporal networks. While these approaches have been useful for estimating these networks at the subject-level, comparisons over groups or experimental conditions require further methodological development. In this paper, we tackle this problem by showing how self-organizing maps (SOMs) can be compared within a Frechean inferential framework. Here, we summarize the mean SOM in each group as a Frechet mean with respect to a metric on the space of SOMs. We consider the use of different metrics, and introduce two extensions of the classical sum of minimum distance (SMD) between two SOMs, which take into account the spatio-temporal pattern of the fMRI data. The validity of these methods is illustrated on synthetic data. Through these simulations, we show that the three metrics of interest behave as expected, in the sense that the ones capturing temporal, spatial and spatio-temporal aspects of the SOMs are more likely to reach significance under simulated scenarios characterized by temporal, spatial and spatio-temporal differences, respectively. In addition, a re-analysis of a classical experiment on visually-triggered emotions demonstrates the usefulness of this methodology. In this study, the multivariate functional patterns typical of the subjects exposed to pleasant and unpleasant stimuli are found to be more similar than the ones of the subjects exposed to emotionally neutral stimuli. Taken together, these results indicate that our proposed methods can cast new light on existing data by adopting a global analytical perspective on functional MRI paradigms.Comment: 23 pages, 5 figures, 4 tables. Submitted to Neuroimag

Amélioration des performances d'un modèle stochastique de génération de hyétogrammes horaires: application au pourtour méditerranéen français

Depuis quelques années, un modèle stochastique de génération de hyétogrammes horaires est développé au groupement d'Aix-en-Provence du Cemagref, pour être couplé à une modélisation de la pluie en débit, fournissant ainsi une multitude de scénarios de crues analysés statistiquement et utilisés en prédétermination des débits de crues. L'extension de la zone d'application du modèle de pluies horaires au-delà de sa zone de conception, a fait apparaître une hétérogénéité dans les résultats. Ce constat a entraîné certaines modifications du modèle comme : la recherche d'une loi de probabilité théorique peu sensible aux problèmes d'échantillonnage pour une variable du modèle (intensité d'une averse), la prise en compte originale de la dépendance observée entre deux variables du modèle (durée et intensité d'une averse), et la modélisation de la persistance des averses au sein d'une même période pluvieuse. Ces différentes modifications apportées au modèle initial ont entraîné une très nette amélioration de ses performances sur la cinquantaine de postes pluviographiques du pourtour méditerranéen français. On obtient ainsi un outil beaucoup plus robuste et validé sur une zone étendue, capable de fournir de multiples formes de hyétogrammes, couvrant toute la gamme des fréquences, permettant ainsi de s'affranchir des pluies de projet uniques. On aborde aussi une nouvelle approche du comportement à l'infini des distributions de fréquences des pluies qui semble parfois supérieur à une tendance strictement exponentielle. De plus, l'étude de plusieurs événements par an dont chacun présente plusieurs réalisations des différentes variables du modèle augmente la taille des échantillons analysés, semblant rendre la méthode plus rapidement fiable qu'une approche statistique classique basée par exemple sur l'ajustement de valeurs maximales annuelles.A stochastic model for generating hourly hyetographs has been recently developed, in the Cemagref of Aix-en-Provence, to be coupled with a rainfall runoff conversion modelling. Thus, by simulation of very long periods (1000 years for example), we obtain a large number of hourly hyetographs and flood scenarios that are statistically studied and used in flood predetermination problems. The rainfall model studied is based on the theory that rainfall can be linked to a random and intermittent process whose evolution is described by stochastic laws. It is also based on the hypothesis of independence between variables describing hyetographs and on the hypothesis of the stationary nature of the phenomenon studied. Generating a rainfall time series involves two steps : descriptive study of the phenomenon (nine independent variables are chosen to describe the phenomenon and these variables are defined by a theoretical law of probability fitted to the observations) and creation of a rainfall time series using descriptive variables generated randomly from their law of probability. Initially developed on the Réal Collobrier watershed data, the model has been applied to fifty raingauges located on the Mediterranean French seaboard. The extension of the model applying area has shown heterogeneousness in the results. Therefore, modifications have been made to the model to improve its performances. Among these modifications, three of them have presented notable improvements. A study of the sensitivity of the parameters has been made. Parameters of shape variables and of some other variables had only a slight influence on depth of generated rainfalls. But, the law of mean rainfall intensities clearly differentiates the stations. Then, a theoretical probability distribution for the storm intensity variable, less sensitive to the sampling problems, has been searched. An exponential distribution is fitted to the value smaller than four times the mean of the variable. A slope breakage was then introduced to generate all the values beyond this limit. The breakage at the value four times the mean of the variable and modelling this breakage were based on a study of so-called "regional" distributions of the storm intensity variable. These distributions were designed by clustering the variable's homogenized values for all 50 studied stations. A second modification has been made to develop new model for the observed dependence between two variables (duration and intensity of the storm). The study of this dependence has been considered directly based on the cumulative frequency of the two variables. Then, an additional parameter was defined to model the dependence between the probabilities of the two variables. This parameter characterises the cumulative frequency curve of the sum of the probabilities of the two variables. This point, neglected during a long time, has been very important in the improvement of the model. Finally, the modelling of storm persistence in a same rainfall episode has been studied to generate some high 24 hours maximum rainfalls. Persistence modelling is entirely justified by the fact that "ordinary storms" cluster together around the "main storm" (the "main storm" is the greatest storm of an episode and the "ordinary storms" are the other storms of the episode). When the study of this phenomenon is extended, it can be observed that there is a certain positive dependency between occurrence probability of the "main storm" and occurrence probability of storms which come before or after it. Two combined effects occur : within one rainy episode, the strongest "ordinary storms" are preferentially clustered together around the "main storm", and considering the number of "ordinary storms" throughout all the episodes, the strongest storms close to the "main storm" are preferentially associated with the strongest "main storms" and vice versa. This modification improves the performances of the altitude raingauges, which are characterised by high daily rainfall accumulations. The different modifications added to the initial model, give very important improvements on the calibration of the fifty raingauges studied on the French Mediterranean seaboard. Its aptitude to generate rains observed in Mediterranean climate, strongly variables, consolidates us in the idea of its application on a zone much larger. The generation of hyetographs makes it possible to use the maximum the temporal information of the rain. Thus, we obtain a reliable tool, validated on a large area, for simulating hyetographs and hourly flood scenarios at all frequencies, and used instead of a unique design storm and design flood. The approach allows a new cumulative probability curve extrapolation, which seems sometimes greater than an exponential behaviour. Moreover, the study of many events per year, with many occurrences of the different variables of the model, increase the analysed sample size and seems to make the method more reliable than a statistical approach simply based, for example, on the fitting of annual maximum values

La modélisation stochastique des pluies horaires et leur transformation en débits pour la prédétermination des crues

Pour étudier les distributions de fréquences des variables hydrologiques (pluies et débits) au pas de temps horaire, une méthodologie associant un générateur de chroniques de pluies horaires et un modèle conceptuel global de transformation de la pluie en débit a été développée. Sur une période de simulation donnée, la méthode génère une collection de scénarios de crues vraisemblables utilisée en prédétermination des risques hydrologiques. Les distributions de fréquences des variables hydrologiques sont construites empiriquement à partir des événements de pluies et de crues générés. L'extrapolation des distributions de fréquences des variables hydrologiques vers les fréquences rares se fait de façon empirique en augmentant la période de simulation, et non plus sur l'ajustement direct des distributions observées. Le principe de cette méthode (appelée SHYPRE : Simulation d'HYdrogrammes pour la PREdétermination) est donc d'utiliser les observations pour décrire le phénomène, afin de le reproduire statistiquement et de s'affranchir ainsi du manque d'observation. Son utilisation permet une estimation originale des quantiles de crues de fréquences courantes à rares et présente l'intérêt d'obtenir une information temporelle complète sur ces crues. De plus, on montre que l'approche fournit une estimation de quantiles de crues bien plus robuste que les ajustements statistiques des distributions observées, même pour les événements de fréquences courantes. Cette robustesse provient d'une meilleure prise en compte de l'information pluviométrique et de la stabilité de la paramétrisation du modèle pluie-débit.A statistical approach encompassing a stochastic model to generate hourly rainfall and rainfall runoff was used to study frequency distributions of hydrologic variables. The method generates numerous different flood events over a given simulation period to evaluate hydrologic risks. Entitled Simulated HYdrographs for flood PRobability Estimation (SHYPRE), it makes use of observed values to describe hydrological phenomena and successfully reproduces observed-value statistics. Frequency distributions of hydrologic variables are built empirically from model-generated rainfall and flood events. Extrapolation of these frequency distributions to rare frequencies is performed by simulation over longer periods, rather than by direct fit of theoretical probability distributions over observed values. This approach yields different estimations of flood quantiles for common to rare frequencies as well as complete temporal flood data. Moreover, SHYPRE estimates of flood quantiles are more stable than statistical distributions fitted onto observed values, even for frequent events. The improvement stems from better use of rainfall data and from the parametric stability of the rainfall model and rainfall-runoff model

A taxonomic review of the neotropical genus Coprophanaeus Olsoufieff, 1924 (Coleoptera: Scarabaeidae, Scarabaeinae)

The Neotropical genus Coprophanaeus Olsoufieff (1924), as classified here, comprises 38 species distributed among three subgenera (Megaphanaeus Olsoufieff, Metallophanaeus Olsoufieff, and Coprophanaeus s. str.) and eight species groups. Keys presented help to identify supraspecific and species taxa, all of which are illustrated and diagnosed. Lectotypes are designated for Phanaeus ignecinctus Felsche and Phanaeus ohausi Felsche. Coprophanaeus corythus (Harold), formerly regarded as a subspecies of C. telamon (Erichson), assumes species status. Coprophanaeus magnoi Arnaud, described as a subspecies of C. milon (Blanchard), is raised to species status. New taxonomic interpretations result in 10 new subjective synonymies (junior synonym listed first): Phanaeus machadoi Pereira and d’Andretta = Coprophanaeus saphirinus (Perty); Phanaeus costatus Olsoufieff = Coprophanaeus cyanescens (Olsoufieff); Phanaeus worontzowi Pessôa and Lane = Coprophanaeus cyanescens (Olsoufieff); Coprophanaeus kohlmanni Arnaud = Coprophanaeus morenoi Arnaud; Coprophanaeus pluto nogueirai Arnaud = Coprophanaeus pluto (Harold); Coprophanaeus edmondsi Arnaud = Coprophanaeus conocephalus (Olsoufieff); Coprophanaeus uhleri Malý and Pokorný = Coprophanaeus chiriquensis (Olsoufieff); Coprophanaeus henryi Malý and Pokorný = Coprophanaeus gilli Arnaud; Phanaeus perseus Harold = Coprophanaeus corythus (Harold); Coprophanaeus telamon nevinsoni Arnaud and Gámez = Coprophanaeus corythus; and Coprophanaeus florenti Arnaud = Coprophanaeus ohausi (Felsche). The status of the following names remains unresolved: Phanaeus strandi Balthasar; Coprophanaeus rigoutorum Arnaud; C. terrali Arnaud; C. lichyi Arnaud; C. lecromi Arnaud; C. larseni Arnaud; and C. vazdemeloi Arnaud

Pertinence d'une méthode de prédétermination des crues basée sur la simulation stochastique de pluies horaires.

International audienceExtreme events are rarely observed, so their analysis is generally based on observations of more frequent values. The relevance of flood frequency analysis (FFA) method depends on its capability to estimate the frequency of extreme values with reasonable accuracy using extrapolation. A FFA method based on stochastic simulation of flood event is assessed based on its reliability and stability. For such an assessment, different training/testing decompositions are performed for a set of data from more than 1,000 gauging stations. We showed that the method enables relevant "predictive" estimates, e.g. by assigning correct return periods to the record values that are systematically absent in calibration data sets. The model is also highly stable vis-a-vis the sampling. This characteristic is linked to the use of regional statistical rainfall data and a simple rainfall-runoff model that requires calibrating only one parameterLes événements extrêmes sont par nature rarement observés, c'est pourquoi leur estimation est généralement basée sur l'observation de valeurs plus courantes. La pertinence d'une méthode de prédétermination des événements extrêmes dépend donc de sa capacité à raisonnablement extrapoler les distributions de fréquences vers les valeurs extrêmes. Dans cette étude, une méthode de prédétermination de crues basée sur la simulation de scénarios de pluies horaires, est évaluée sur sa capacité à produire des estimations justes et stables. Cette évaluation s'appuie sur différents tests d'échantillonnage sur les périodes de calage et de validation, appliqués sur un jeu de données conséquent (plus de 1000 stations). Nous montrons que la méthode est capable de fournir une estimation pertinente sur les événements extrêmes bien que ceux-ci soit systématiquement ôtés de la période de calage. La méthode montre aussi une grande stabilité face à l'échantillonnage. Cette caractéristique est liée à l'utilisation d'une information statistique régionale sur la pluie et à la simplicité de la modélisation hydrologique paramétrée par un seul paramètre

High-temperature-grown buffer layer boosts electron mobility in epitaxial La-doped BaSnO3_3/SrZrO3_3 heterostructures

By inserting a SrZrO$_3$ buffer layer between the film and the substrate, we demonstrate a significant reduction of the threading dislocation density with an associated improvement of the electron mobility in La:BaSnO$_3$ films. A room temperature mobility of 140 cm$^2$ V$^{-1}\text{s}^{-1}$ is achieved for 25-nm-thick films without any post-growth treatment. The density of threading dislocations is only $4.9\times 10^{9}$ cm$^{-2}$ for buffered films prepared on (110) TbScO$_3$ substrates by pulsed laser deposition.Comment: 5 pages, 4 figure