Asymptotic robust inferences in the analysis of mean and covariance structures

Structural equation models are widely used in economic, socialand behavioral studies to analyze linear interrelationships amongvariables, some of which may be unobservable or subject to measurementerror. Alternative estimation methods that exploit different distributionalassumptions are now available. The present paper deals with issues ofasymptotic statistical inferences, such as the evaluation of standarderrors of estimates and chi--square goodness--of--fit statistics,in the general context of mean and covariance structures. The emphasisis on drawing correct statistical inferences regardless of thedistribution of the data and the method of estimation employed. A(distribution--free) consistent estimate of $\Gamma$, the matrix ofasymptotic variances of the vector of sample second--order moments,will be used to compute robust standard errors and a robust chi--squaregoodness--of--fit squares. Simple modifications of the usual estimateof $\Gamma$ will also permit correct inferences in the case of multi--stage complex samples. We will also discuss the conditions under which,regardless of the distribution of the data, one can rely on the usual(non--robust) inferential statistics. Finally, a multivariate regressionmodel with errors--in--variables will be used to illustrate, by meansof simulated data, various theoretical aspects of the paper.

Multi-sample analysis of moment-structures: Asymptotic validity of inferences based on second-order moments

In moment structure analysis with nonnormal data, asymptotic valid inferences require the computation of a consistent (under general distributional assumptions) estimate of the matrix $\Gamma$ of asymptotic variances of sample second--order moments. Such a consistent estimate involves the fourth--order sample moments of the data. In practice, the use of fourth--order moments leads to computational burden and lack of robustness against small samples. In this paper we show that, under certain assumptions, correct asymptotic inferences can be attained when $\Gamma$ is replaced by a matrix $\Omega$ that involves only the second--order moments of the data. The present paper extends to the context of multi--sample analysis of second--order moment structures, results derived in the context of (simple--sample) covariance structure analysis (Satorra and Bentler, 1990). The results apply to a variety of estimation methods and general type of statistics. An example involving a test of equality of means under covariance restrictions illustrates theoretical aspects of the paper.

Tectonics of the Simplon massif and Lepontine gneiss dome: deformation structures due to collision between the underthrusting European plate and the Adriatic indenter

this study presents a review of published geological data, combined with original observations on the tectonics of the simplon massif and the Lepontine gneiss dome in the Western Alps. New observations concern the geometry of the Oligocene Vanzone back fold, formed under amphibolite facies conditions, and of its root between Domodossola and Locarno, which is cut at an acute angle by the Miocene, epi- to anchizonal, dextral centovalli strike-slip fault. the structures of the simplon massif result from collision over 50 Ma between two plate boundaries with a different geometry: the underthrusted European plate and the Adriatic indenter. Detailed mapping and analysis of a complex structural interference pattern, combined with observations on the metamorphic grade of the superimposed structures and radiometric data, allow a kinematic model to be developed for this zone of oblique continental collision. the following main Alpine tectonic phases and structures may be distinguished: 1. NW-directed nappe emplacement, starting in the Early Eocene (similar to 50 Ma); 2. W, SW and S- verging transverse folds; 3. transpressional movements on the dextral simplon ductile shear zone since similar to 32 Ma; 4. formation of the Bergell - Vanzone backfolds and of the southern steep belt during the Oligocene, emplacement of the mantle derived 31 - 29 Ma Bergell and Biella granodiorites and porphyritic andesites as well as intrusions of 29-25 Ma crustal aplites and pegmatites; 5. formation of the dextral discrete Rhone-Simplon line and the centovalli line during the Miocene, accompanied by the pull-apart development of the Lepontine gneiss dome - Dent blanche (Valpelline) depression. It is suggested that movements of shortening in fan shaped NW, W and sW directions accompanied the more regular NW- to WNW-directed displacement of the Adriatic indenter during continental collision.

Homologous DNA pairing domain peptides of RecA protein: intrinsic propensity to form beta-structures and filaments.

The 20 amino acid residue peptides derived from RecA loop L2 have been shown to be the pairing domain of RecA. The peptides bind to ss- and dsDNA, unstack ssDNA, and pair the ssDNA to its homologous target in a duplex DNA. As shown by circular dichroism, upon binding to DNA the disordered peptides adopt a beta-structure conformation. Here we show that the conformational change of the peptide from random coil to beta-structure is important in binding ss- and dsDNA. The beta-structure in the DNA pairing peptides can be induced by many environmental conditions such as high pH, high concentration, and non-micellar sodium dodecyl sulfate (6 mM). This behavior indicates an intrinsic property of these peptides to form a beta-structure. A beta-structure model for the loop L2 of RecA protein when bound to DNA is thus proposed. The fact that aromatic residues at the central position 203 strongly modulate the peptide binding to DNA and subsequent biochemical activities can be accounted for by the direct effect of the aromatic amino acids on the peptide conformational change. The DNA-pairing domain of RecA visualized by electron microscopy self-assembles into a filamentous structure like RecA. The relevance of such a peptide filamentous structure to the structure of RecA when bound to DNA is discussed.

Tertiary Himalayan structures and metamorphism in the Kulu Valley (Mandi-Khoksar transect of the Western Himalaya) - Shikar-Beh-nappe and crystalline nappe

The Crystalline Nappe of the High Himalayan Crystalline has been examined along the Kulu Valley and its vicinity (Mandi-Khoksar transect). This nappe was believed to have undergone deformation related only to its transport towards the SW essentially during the `'Main Central Thrust event''. New data has led to the conclusion that during the Himalayan orogeny, two distinctive phases, related to two opposite transport directions, characterize the evolution of this part of the chain, before the creation of the late NE-vergent backfolding. The first phase corresponds to an early NE-vergent folding and thrusting, creating the Tandi Syncline and the NE-oriented Shikar Beh Nappe stack, with a displacement amplitude of about 50 km. Two schistosities, together with a strong stretching lineation are developed at a deep tectonic level under amphibolite facies conditions (kyanite-staurolite-garnet-two mica schists). At a higher tectonic level and in the southern part of the section (Tandy Syncline and southern Kulu Valley between Kulu and Mandi) one or two schistosities are developed in the greenschist facies grade rocks (garnet-biotite and biotite schists). These structures and the associated Barrovian type metamorphism are all related to the NE-verging Shikar Beh Nappe. The creation of the NE-verging Shikar Beh Nappe may be explained by the reactivation of a SW dipping listric normal fault of the N Indian flexural passive margin, during the early stages of the Himalayan orogeny. In the second phase, the still hot metamorphic rocks of the Shikar Beh Nappe were folded and thrust towards the SW (mainly along the MBT and the MCT with a displacement in excess of 100 km) onto the cold, low-grade metamorphic rocks of the Larji-Kulu-Rampur Window or, near Mandi, on the non-metamorphic sandstones of the Ganges Molasse (Siwaliks). Sense of shear criteria and a strong NE-SW stretching-lineation indicate that the Crystalline Nappe has been overthrusted towards the SW. Thermometry on synkinematically crystallised garnet-biotite and garnet-hornblende pairs reveals the lower amphibolite facies temperature conditions related to the Crystalline Nappe formation. From the muscovite and biotite Rb-Sr cooling ages, the Shikar Beh Nappe emplacement occurred before 32 Ma and the southwestward thrusting of the Crystalline Nappe began before 21 Ma. Our model involving two opposite directions of thrusting goes against the conventional idea of only one main SW-oriented transport direction in the High Himalayan Crystalline Nappes.

A methodology for the classification of convective structures using meteorological radar: Application to heavy rainfall events on the Mediterranean coast of the Iberian Peninsula

During the period 1996-2000, forty-three heavy rainfall events have been detected in the Internal Basins of Catalonia (Northeastern of Spain). Most of these events caused floods and serious damage. This high number leads to the need for a methodology to classify them, on the basis of their surface rainfall distribution, their internal organization and their physical features. The aim of this paper is to show a methodology to analyze systematically the convective structures responsible of those heavy rainfall events on the basis of the information supplied by the meteorological radar. The proposed methodology is as follows. Firstly, the rainfall intensity and the surface rainfall pattern are analyzed on the basis of the raingauge data. Secondly, the convective structures at the lowest level are identified and characterized by using a 2-D algorithm, and the convective cells are identified by using a 3-D procedure that looks for the reflectivity cores in every radar volume. Thirdly, the convective cells (3-D) are associated with the 2-D structures (convective rainfall areas). This methodology has been applied to the 43 heavy rainfall events using the meteorological radar located near Barcelona and the SAIH automatic raingauge network.

Analysis of warm season thunderstorms using an object-oriented tracking method based on radar and total lightning data

Monitoring thunderstorms activity is an essential part of operational weather surveillance given their potential hazards, including lightning, hail, heavy rainfall, strong winds or even tornadoes. This study has two main objectives: firstly, the description of a methodology, based on radar and total lightning data to characterise thunderstorms in real-time; secondly, the application of this methodology to 66 thunderstorms that affected Catalonia (NE Spain) in the summer of 2006. An object-oriented tracking procedure is employed, where different observation data types generate four different types of objects (radar 1-km CAPPI reflectivity composites, radar reflectivity volumetric data, cloud-to-ground lightning data and intra-cloud lightning data). In the framework proposed, these objects are the building blocks of a higher level object, the thunderstorm. The methodology is demonstrated with a dataset of thunderstorms whose main characteristics, along the complete life cycle of the convective structures (development, maturity and dissipation), are described statistically. The development and dissipation stages present similar durations in most cases examined. On the contrary, the duration of the maturity phase is much more variable and related to the thunderstorm intensity, defined here in terms of lightning flash rate. Most of the activity of IC and CG flashes is registered in the maturity stage. In the development stage little CG flashes are observed (2% to 5%), while for the dissipation phase is possible to observe a few more CG flashes (10% to 15%). Additionally, a selection of thunderstorms is used to examine general life cycle patterns, obtained from the analysis of normalized (with respect to thunderstorm total duration and maximum value of variables considered) thunderstorm parameters. Among other findings, the study indicates that the normalized duration of the three stages of thunderstorm life cycle is similar in most thunderstorms, with the longest duration corresponding to the maturity stage (approximately 80% of the total time).

Climate oscillations and species interactions: large-scale congruence but regional differences in the phylogeographic structures of an alpine plant and its monophagous insect

Aim. To predict the fate of alpine interactions involving specialized species, using a monophagous beetle and its host-plant as a case study. Location. The Alps. Methods. We investigated genetic structuring of the herbivorous beetle Oreina gloriosa and its specific host-plant Peucedanum ostruthium. We used genome fingerprinting (in the insect and the plant) and sequence data (in the insect) to compare the distribution of the main gene pools in the two associated species and to estimate divergence time in the insect, a proxy for the temporal origin of the interaction. We quantified the similarity in spatial genetic structures by performing a Procrustes analysis, a tool from the shape theory. Finally, we simulated recolonization of an empty space analogous to the deglaciated Alps just after ice retreat by two lineages from two species showing unbalanced dependence, to examine how timing of the recolonization process, as well as dispersal capacities of associated species, could explain the observed pattern. Results. Contrasting with expectations based on their asymmetrical dependence, patterns in the beetle and plant were congruent at a large scale. Exceptions occurred at a regional scale in areas of admixture, matching known suture zones in Alpine plants. Simulations using a lattice-based model suggested these empirical patterns arose during or soon after recolonization, long after the estimated origin of the interaction c. 0.5 million years ago. Main conclusions. Species-specific interactions are scarce in alpine habitats because glacial cycles have limited opportunities for coevolution. Their fate, however, remains uncertain under climate change. Here we show that whereas most dispersal routes are paralleled at large scale, regional incongruence implies that the destinies of the species might differ under changing climate. This may be a consequence of the host-dependence of the beetle that locally limits the establishment of dispersing insects.

Frequency resolved admittance spectroscopy measurements on In0.52Al0.48As/InxGa1¿xAs/In0.52Al0.48As single quantum well structures

In this work a new admittance spectroscopy technique is proposed to determine the conduction band offset in single quantum well structures (SQW). The proposed technique is based on the study of the capacitance derivative versus the frequency logarithm. This method is found to be less sensitive to parasitic effects, such as leakage current and series resistance, than the classical conductance analysis. Using this technique, we have determined the conduction band offset in In0.52Al0.48As/InxGa1¿xAs/In0.52Al0.48As SQW structures. Two different well compositions, x=0.53, which corresponds to the lattice¿matched case and x=0.60, which corresponds to a strained case, and two well widths (5 and 25 nm) have been considered. The average results are ¿Ec=0.49±0.04 eV for x=0.53 and ¿Ec =0.51±0.04 eV for x=0.6, which are in good agreement with previous reported data.

Raman scattering and photoluminescence analysis of silicon on insulator structures obtained by single and multiple oxygen implants

An analysis of silicon on insulator structures obtained by single and multiple implants by means of Raman scattering and photoluminescence spectroscopy is reported. The Raman spectra obtained with different excitation powers and wavelengths indicate the presence of a tensile strain in the top silicon layer of the structures. The comparison between the spectra measured in both kinds of samples points out the existence in the multiple implant material of a lower strain for a penetration depth about 300 nm and a higher strain for higher penetration depths. These results have been correlated with transmission electron microscopy observations, which have allowed to associate the higher strain to the presence of SiO2 precipitates in the top silicon layer, close to the buried oxide. The found lower strain is in agreement with the better quality expected for this material, which is corroborated by the photoluminescence data.

Raman microstructural analysis of silicon-on-insulator formed by high dose oxygen ion implantation: As-implanted structures

A microstructural analysis of silicon-on-insulator samples obtained by high dose oxygen ion implantation was performed by Raman scattering. The samples analyzed were obtained under different conditions thus leading to different concentrations of defects in the top Si layer. The samples were implanted with the surface covered with SiO2 capping layers of different thicknesses. The spectra measured from the as-implanted samples were fitted to a correlation length model taking into account the possible presence of stress effects in the spectra. This allowed quantification of both disorder effects, which are determined by structural defects, and residual stress in the top Si layer before annealing. These data were correlated to the density of dislocations remaining in the layer after annealing. The analysis performed corroborates the existence of two mechanisms that generate defects in the top Si layer that are related to surface conditions during implantation and the proximity of the top Si/buried oxide layer interface to the surface before annealing.

BIO-INSPIRED COMPUTATIONAL TECHNIQUES APPLIED TO THE CLUSTERING AND VISUALIZATION OF SPATIO-TEMPORAL GEOSPATIAL DATA

The coverage and volume of geo-referenced datasets are extensive and incessantly¦growing. The systematic capture of geo-referenced information generates large volumes¦of spatio-temporal data to be analyzed. Clustering and visualization play a key¦role in the exploratory data analysis and the extraction of knowledge embedded in¦these data. However, new challenges in visualization and clustering are posed when¦dealing with the special characteristics of this data. For instance, its complex structures,¦large quantity of samples, variables involved in a temporal context, high dimensionality¦and large variability in cluster shapes.¦The central aim of my thesis is to propose new algorithms and methodologies for¦clustering and visualization, in order to assist the knowledge extraction from spatiotemporal¦geo-referenced data, thus improving making decision processes.¦I present two original algorithms, one for clustering: the Fuzzy Growing Hierarchical¦Self-Organizing Networks (FGHSON), and the second for exploratory visual data analysis:¦the Tree-structured Self-organizing Maps Component Planes. In addition, I present¦methodologies that combined with FGHSON and the Tree-structured SOM Component¦Planes allow the integration of space and time seamlessly and simultaneously in¦order to extract knowledge embedded in a temporal context.¦The originality of the FGHSON lies in its capability to reflect the underlying structure¦of a dataset in a hierarchical fuzzy way. A hierarchical fuzzy representation of¦clusters is crucial when data include complex structures with large variability of cluster¦shapes, variances, densities and number of clusters. The most important characteristics¦of the FGHSON include: (1) It does not require an a-priori setup of the number¦of clusters. (2) The algorithm executes several self-organizing processes in parallel.¦Hence, when dealing with large datasets the processes can be distributed reducing the¦computational cost. (3) Only three parameters are necessary to set up the algorithm.¦In the case of the Tree-structured SOM Component Planes, the novelty of this algorithm¦lies in its ability to create a structure that allows the visual exploratory data analysis¦of large high-dimensional datasets. This algorithm creates a hierarchical structure¦of Self-Organizing Map Component Planes, arranging similar variables' projections in¦the same branches of the tree. Hence, similarities on variables' behavior can be easily¦detected (e.g. local correlations, maximal and minimal values and outliers).¦Both FGHSON and the Tree-structured SOM Component Planes were applied in¦several agroecological problems proving to be very efficient in the exploratory analysis¦and clustering of spatio-temporal datasets.¦In this thesis I also tested three soft competitive learning algorithms. Two of them¦well-known non supervised soft competitive algorithms, namely the Self-Organizing¦Maps (SOMs) and the Growing Hierarchical Self-Organizing Maps (GHSOMs); and the¦third was our original contribution, the FGHSON. Although the algorithms presented¦here have been used in several areas, to my knowledge there is not any work applying¦and comparing the performance of those techniques when dealing with spatiotemporal¦geospatial data, as it is presented in this thesis.¦I propose original methodologies to explore spatio-temporal geo-referenced datasets¦through time. Our approach uses time windows to capture temporal similarities and¦variations by using the FGHSON clustering algorithm. The developed methodologies¦are used in two case studies. In the first, the objective was to find similar agroecozones¦through time and in the second one it was to find similar environmental patterns¦shifted in time.¦Several results presented in this thesis have led to new contributions to agroecological¦knowledge, for instance, in sugar cane, and blackberry production.¦Finally, in the framework of this thesis we developed several software tools: (1)¦a Matlab toolbox that implements the FGHSON algorithm, and (2) a program called¦BIS (Bio-inspired Identification of Similar agroecozones) an interactive graphical user¦interface tool which integrates the FGHSON algorithm with Google Earth in order to¦show zones with similar agroecological characteristics.

Improved segmentation of deep brain grey matter structures using magnetization transfer (MT) parameter maps.

Basal ganglia and brain stem nuclei are involved in the pathophysiology of various neurological and neuropsychiatric disorders. Currently available structural T1-weighted (T1w) magnetic resonance images do not provide sufficient contrast for reliable automated segmentation of various subcortical grey matter structures. We use a novel, semi-quantitative magnetization transfer (MT) imaging protocol that overcomes limitations in T1w images, which are mainly due to their sensitivity to the high iron content in subcortical grey matter. We demonstrate improved automated segmentation of putamen, pallidum, pulvinar and substantia nigra using MT images. A comparison with segmentation of high-quality T1w images was performed in 49 healthy subjects. Our results show that MT maps are highly suitable for automated segmentation, and so for multi-subject morphometric studies with a focus on subcortical structures.

From protein sequences to 3D-structures and beyond: the example of the UniProt knowledgebase.

With the dramatic increase in the volume of experimental results in every domain of life sciences, assembling pertinent data and combining information from different fields has become a challenge. Information is dispersed over numerous specialized databases and is presented in many different formats. Rapid access to experiment-based information about well-characterized proteins helps predict the function of uncharacterized proteins identified by large-scale sequencing. In this context, universal knowledgebases play essential roles in providing access to data from complementary types of experiments and serving as hubs with cross-references to many specialized databases. This review outlines how the value of experimental data is optimized by combining high-quality protein sequences with complementary experimental results, including information derived from protein 3D-structures, using as an example the UniProt knowledgebase (UniProtKB) and the tools and links provided on its website ( http://www.uniprot.org/ ). It also evokes precautions that are necessary for successful predictions and extrapolations.