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.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Ultra-High Field (7 Tesla) MRI for Gamma Knife Surgery Targeting of the Ventro-Intermediate Nucleus: A Pilot in Vivo Study

Aim: Gamma Knife surgery (GKS) is a non-invasive neurosurgical stereotactic procedure, increasingly used as an alternative to open functional procedures. This includes the targeting of the ventro-intermediate (Vim) nucleus of the thalamus for tremor. We currently perform an indirect targeting, using the "quadrilatere of Guyot," as the Vim nucleus is not visible on current 3 Tesla (T) MRI acquisitions. The primary objective of the current study was to enhance anatomic imaging for Vim GKS using high-field (7 T) MRI, with the aim of refining the visualization and precision of anatomical targeting. Method: Five young healthy subjects (mean age 23 years) were scanned both on 3 and 7 T MRI in Lausanne University Hospital (CHUV) and Center for Biomedical Imaging (CIBM). Classical T1-weighted MPRAGE, T2 CISS sequences (replacing former ventriculography) and diffusion tensor imaging were acquired at 3T. We obtained high-resolution susceptibility weighted images (SWI) at 7T for the visualization of thalamic subparts. SWI was further integrated for the first time into Leksell Gamma Plan® (LGP) software and co-registered with the 3T images. A simulation of targeting of the Vim was done using the "quadrilatere of Guyot" methodology on the 3T images. Furthermore, a correlation with the position of the found target on SWI was performed. The atlas of Morel et al. was used to confirm the findings on a detailed computer analysis outside LGP. Also, 3T and 7T MRI of one patient undergoing GKS Vim thalamotomy, were obtained before and 2 years after the procedure, and studied similarly. Results: The use of SWI provided a superior resolution and improved image contrast within the central gray matter. This allowed visualization and direct delineation of groups of thalamic nuclei in vivo, including the Vim. The position of the target, as assessed with the "quadrilatere of Guyot" method on 3 T, perfectly matched with the supposed one of the Vim on the SWI. Furthermore, a 3-dimensional model of the Vim target area was created on the basis of 3T and 7T images. Conclusion: This is the first report of the integration of SWI high-field MRI into the LGP in healthy subjects and in one patient treated GKS Vim thalamotomy. This approach aims at the improvement of targeting validation and further direct targeting of the Vim in tremor. The anatomical correlation between the direct visualization on 7T and the current targeting methods on 3T seems to show a very good anatomical matching.

Revisiting Field Capacity (FC): variation of definition of FC and its estimation from pedotransfer functions

Taking into account the nature of the hydrological processes involved in in situ measurement of Field Capacity (FC), this study proposes a variation of the definition of FC aiming not only at minimizing the inadequacies of its determination, but also at maintaining its original, practical meaning. Analysis of FC data for 22 Brazilian soils and additional FC data from the literature, all measured according to the proposed definition, which is based on a 48-h drainage time after infiltration by shallow ponding, indicates a weak dependency on the amount of infiltrated water, antecedent moisture level, soil morphology, and the level of the groundwater table, but a strong dependency on basic soil properties. The dependence on basic soil properties allowed determination of FC of the 22 soil profiles by pedotransfer functions (PTFs) using the input variables usually adopted in prediction of soil water retention. Among the input variables, soil moisture content θ (6 kPa) had the greatest impact. Indeed, a linear PTF based only on it resulted in an FC with a root mean squared residue less than 0.04 m³ m-3 for most soils individually. Such a PTF proved to be a better FC predictor than the traditional method of using moisture content at an arbitrary suction. Our FC data were compatible with an equivalent and broader USA database found in the literature, mainly for medium-texture soil samples. One reason for differences between FCs of the two data sets of fine-textured soils is due to their different drainage times. Thus, a standardized procedure for in situ determination of FC is recommended.

Development of LRFD Procedures for Bridge Pile Foundations in Iowa: Volume II: Field Testing of Steel Piles in Clay, Sand, and Mixed Soils and Data Analysis, September 2011

In response to the mandate on Load and Resistance Factor Design (LRFD) implementations by the Federal Highway Administration (FHWA) on all new bridge projects initiated after October 1, 2007, the Iowa Highway Research Board (IHRB) sponsored these research projects to develop regional LRFD recommendations. The LRFD development was performed using the Iowa Department of Transportation (DOT) Pile Load Test database (PILOT). To increase the data points for LRFD development, develop LRFD recommendations for dynamic methods, and validate the results ofLRFD calibration, 10 full-scale field tests on the most commonly used steel H-piles (e.g., HP 10 x 42) were conducted throughout Iowa. Detailed in situ soil investigations were carried out, push-in pressure cells were installed, and laboratory soil tests were performed. Pile responses during driving, at the end of driving (EOD), and at re-strikes were monitored using the Pile Driving Analyzer (PDA), following with the CAse Pile Wave Analysis Program (CAPWAP) analysis. The hammer blow counts were recorded for Wave Equation Analysis Program (WEAP) and dynamic formulas. Static load tests (SLTs) were performed and the pile capacities were determined based on the Davisson’s criteria. The extensive experimental research studies generated important data for analytical and computational investigations. The SLT measured loaddisplacements were compared with the simulated results obtained using a model of the TZPILE program and using the modified borehole shear test method. Two analytical pile setup quantification methods, in terms of soil properties, were developed and validated. A new calibration procedure was developed to incorporate pile setup into LRFD.

Bilateral Sequential NAION: A Comparison of Visual Outcomes in Fellow Eyes Using Quantitative Analysis of Goldmann Visual Fields

Purpose: Previous studies of the visual outcome in bilateral non-arteritic anterior ischemic optic neuropathy (NAION) have yielded conflicting results, specifically regarding congruity between fellow eyes. Prior studies have used measures of acuity and computerized perimetry but none has compared Goldmann visual field outcomes between fellow eyes. In order to better define the concordance of visual loss in this condition, we reviewed our cases of bilateral sequential NAION, including measures of visual acuity, pupillary function and both pattern and severity of visual field loss.Methods: We performed a retrospective chart review of 102 patients with a diagnosis of bilateral sequential NAION. Of the 102 patients, 86 were included in the study for analysis of final visual outcome between the affected eyes. Visual function was assessed using visual acuity, Goldmann visual fields, color vision and RAPD. A quantitative total visual field score and score per quadrant was analyzed for each eye using the numerical Goldmann visual field scoring method previously described by Esterman and colleagues. Based upon these scores, we calculated the total deviation and pattern deviation between fellow eyes and between eyes of different patients. Statistical significance was determined using nonparametric tests.Results: A statistically significant correlation was found between fellow eyes for multiple parameters, including logMAR visual acuity (P = 0.0101), global visual field (P = 0.0001), superior visual field (P = 0.0001), and inferior visual field (P = 0.0001). In addition, the mean deviation of both total (P = 0.0000000007) and pattern (P = 0.000000004) deviation analyses was significantly less between fellow eyes ("intra"-eyes) than between eyes of different patients ("inter"-eyes).Conclusions: Visual function between fellow eyes showed a fair to moderate correlation that was statistically significant. The pattern of vision loss was also more similar in fellow eyes than between eyes of different patients. These results may help allow better prediction of visual outcome for the second eye in patients with NAION. These findings may also be useful for evaluating efficacy of therapeutic interventions.

Cost-Benefit Analysis of Sequential Warning Lights in Nighttime Work Zone Tapers, TPF-5(081), 2011

Improving safety at nighttime work zones is important because of the extra visibility concerns. The deployment of sequential lights is an innovative method for improving driver recognition of lane closures and work zone tapers. Sequential lights are wireless warning lights that flash in a sequence to clearly delineate the taper at work zones. The effectiveness of sequential lights was investigated using controlled field studies. Traffic parameters were collected at the same field site with and without the deployment of sequential lights. Three surrogate performance measures were used to determine the impact of sequential lights on safety. These measures were the speeds of approaching vehicles, the number of late taper merges and the locations where vehicles merged into open lane from the closed lane. In addition, an economic analysis was conducted to monetize the benefits and costs of deploying sequential lights at nighttime work zones. The results of this study indicates that sequential warning lights had a net positive effect in reducing the speeds of approaching vehicles, enhancing driver compliance, and preventing passenger cars, trucks and vehicles at rural work zones from late taper merges. Statistically significant decreases of 2.21 mph mean speed and 1 mph 85% speed resulted with sequential lights. The shift in the cumulative speed distributions to the left (i.e. speed decrease) was also found to be statistically significant using the Mann-Whitney and Kolmogorov-Smirnov tests. But a statistically significant increase of 0.91 mph in the speed standard deviation also resulted with sequential lights. With sequential lights, the percentage of vehicles that merged earlier increased from 53.49% to 65.36%. A benefit-cost ratio of around 5 or 10 resulted from this analysis of Missouri nighttime work zones and historical crash data. The two different benefitcost ratios reflect two different ways of computing labor costs.

Generation mean analysis in relation to polygenic systems with epistasis and fixed genes

Epistatic effects involving genic combinations of fixed and non fixed genes are shown to contribute to the genotypic mean of any population. These effects define specific additive x additive and additive x dominant epistatic components. As such components are not estimable, their relative importance cannot be assessed. These epistatic effects can cause bias in the estimates of the additive and dominance components to which they are confounded. The magnitude of the bias depends on the relative values of the epistatic effects, comparatively to deviations d and h, type of prevailing epistasis and direction of dominance.

International Equity on Greenhouse Gas Emissions and World Levels: an integrated analysis through distributive welfare indices

Up until now, analyses of the international distribution of pollutant emissions have not paid sufficient attention to the implications that, in terms of social welfare, the combined evolution of the global world average entails. In this context, this paper proposes the use of environmental welfare indices, taken and adapted from the literature on social welfare and inequality, in order to make a comprehensive examination of the international equity factor and the mean factor in this field. The proposed methodology is implemented empirically in order to explore the evolution in distributive-based environmental welfare on a global level for the three main pollutants with greenhouse gas effects: CO2, CH4 and NO, both globally and for selected years during the period of 1990- 2005. The main results found are as follows: firstly, typically, the environmental welfare associated with the overall greenhouse gases decreased significantly over the period, due primarily to the role of CO2; secondly, in contrast, the global welfare associated with CH4 and NO improved; and thirdly, typically, the evolutions can be attributed to a greater extent to the mean component than to the distributive component, although there are exceptions. These results would seem to be relevant in policy terms. JEL codes: D39; Q43; Q56. Keywords: environmental welfare: greenhouse gases; environmental equity.

Frontshear and backshear instabilities of the mean longshore current

An analytical model based on Bowen and Holman [1989] is used to prove the existence of instabilities due to the presence of a second extremum of the background vorticity at the front side of the longshore current. The growth rate of the so-called frontshear waves depends primarily upon the frontshear but also upon the backshear and the maximum and the width of the current. Depending on the values of these parameters, either the frontshear or the backshear instabilities may dominate. Both types of waves have a cross-shore extension of the order of the width of the current, but the frontshear modes are localized closer to the coast than are the backshear modes. Moreover, under certain conditions both unstable waves have similar growth rates with close wave numbers and angular frequencies, leading to the possibility of having modulated shear waves in the alongshore direction. Numerical analysis performed on realistic current profiles confirm the behavior anticipated by the analytical model. The theory has been applied to a current profile fitted to data measured during the 1980 Nearshore Sediment Transport Studies experiment at Leadbetter Beach that has an extremum of background vorticity at the front side of the current. In this case and in agreement with field observations, the model predicts instability, whereas the theory based only on backshear instability fai led to do so.

Development of LRFD Design Procedures for Bridge Piles in Iowa – Field Testing of Steel H-Piles in Clay, Sand, and Mixed Soils and Data Analysis (Volume II), TR-583, 2011

In response to the mandate on Load and Resistance Factor Design (LRFD) implementations by the Federal Highway Administration (FHWA) on all new bridge projects initiated after October 1, 2007, the Iowa Highway Research Board (IHRB) sponsored these research projects to develop regional LRFD recommendations. The LRFD development was performed using the Iowa Department of Transportation (DOT) Pile Load Test database (PILOT). To increase the data points for LRFD development, develop LRFD recommendations for dynamic methods, and validate the results of LRFD calibration, 10 full-scale field tests on the most commonly used steel H-piles (e.g., HP 10 x 42) were conducted throughout Iowa. Detailed in situ soil investigations were carried out, push-in pressure cells were installed, and laboratory soil tests were performed. Pile responses during driving, at the end of driving (EOD), and at re-strikes were monitored using the Pile Driving Analyzer (PDA), following with the CAse Pile Wave Analysis Program (CAPWAP) analysis. The hammer blow counts were recorded for Wave Equation Analysis Program (WEAP) and dynamic formulas. Static load tests (SLTs) were performed and the pile capacities were determined based on the Davisson’s criteria. The extensive experimental research studies generated important data for analytical and computational investigations. The SLT measured load-displacements were compared with the simulated results obtained using a model of the TZPILE program and using the modified borehole shear test method. Two analytical pile setup quantification methods, in terms of soil properties, were developed and validated. A new calibration procedure was developed to incorporate pile setup into LRFD.

Total hip replacement with a collarless polished cemented anatomic stem: clinical and gait analysis results at ten years follow-up.

PURPOSE: The aim of this study was to determine outcomes of total hip replacement (THR) with the Lemania cemented femoral stem. METHODS: A total of 78 THR patients were followed and compared to 17 "fit", healthy, elderly and 72 "frail" elderly subjects without THR, using clinical outcome measures and a portable, in-field gait analysis device at five and ten years follow-up. RESULTS: Forty-one patients (53%), mean age 83.4 years, available at ten years follow-up, reported very good to excellent satisfaction. Mean Harris Hip and Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) scores were 81.2 and 10.5 points, respectively, with excellent radiological preservation of proximal femur bone stock. Spatial and temporal gait parameters were close to the fit group and better than the frail group. CONCLUSIONS: Lemania THR demonstrated very good, stable clinical and radiological results at ten years in an older patient group, comparable to other cemented systems for primary THR. Gait analysis confirmed good walking performance in a real-life environment.

Structural analysis, clay mineralogy and K-Ar dating of fault gouges from Centovalli Line (Central Alps) for reconstruction of their recent activity

Between the cities of Domodossola and Locarno, the complex ``Centovalli Line'' tectonic zone of the Central Alps outlines deformation phases over a long period of time (probably starting similar to 30 Ma ago) and under variable P-T conditions. The last deformation phases developed gouge-bearing faults with a general E-W trend that crosscuts the roots of the Alpine Canavese zone and the Finero ultramafic body. Kinematic indicators show that the general motion was mainly dextral associated with back thrusting towards the S. The <2 mu m clay fractions of fault gouges from Centovalli Line consist mainly of illite, smectite and chlorite with varied illite-smectite, chlorite-smectite and chlorite-serpentine mixed-layers. Constrained with the illite crystallinity index, the thermal conditions induced by the tectonic activity show a gradual trend from anchizonal to diagenetic conditions. The <2 and <0.2 mu M clay fractions, and hydrothermal K-feldspar separates all provide K-Ar ages between 14.2 +/- 2.9 Ma and roughly 0 Ma, with major episodes at about 12,8, 6 and close to 0 Ma These ages set the recurrent tectonic activity and the associated fluid circulations between Upper Miocene and Recent. On the basis of the K-Ar ages and with a thermal gradient of 25-30 degrees C/km, the studied fault zones were located at a depth of 4-7 km. If they were active until now as observed in field, the exhumation was approximately 2.5-3.0 km for the last 12 Ma with a mean velocity of 0.4 mm/y. Comparison with available models on the recent Alpine evolution shows that the tectonic activity in the area relates to a continuum of the back-thrusting movements of the Canavese Line, and/or to several late-extensional phases of the Rhone-Simplon line. The Centovalli-Val Vigezzo zone therefore represents a major tectonic zone of the Central-Western Alps resulting from different interacting tectonic events. (C) 2011 Elsevier B.V. All rights reserved.

Computational and Experimental Analysis of Flow Field in the Diffusers of Centrifugal Compressors

Centrifugal compressors are widely used for example in process industry, oil and gas industry, in small gas turbines and turbochargers. In order to achieve lower consumption of energy and operation costs the efficiency of the compressor needs to be improve. In the present work different pinches and low solidity vaned diffusers were utilized in order to improve the efficiency of a medium size centrifugal compressor. In this study, pinch means the decrement of the diffuser flow passage height. First different geometries were analyzed using computational fluid dynamics. The flow solver Finflo was used to solve the flow field. Finflo is a Navier-Stokes solver. The solver is capable to solve compressible, incompressible, steady and unsteady flow fields. Chien's k-e turbulence model was used. One of the numerically investigated pinched diffuser and one low solidity vaned diffuser were studied experimentally. The overall performance of the compressor and the static pressure distribution before and after the diffuser were measured. The flow entering and leaving the diffuser was measured using a three-hole Cobra-probe and Kiel-probes. The pinch and the low solidity vaned diffuser increased the efficiency of the compressor. Highest isentropic efficiency increment obtained was 3\% of the design isentropic efficiency of the original geometry. It was noticed in the numerical results that the pinch made to the hub and the shroud wall was most beneficial to the operation of the compressor. Also the pinch made to the hub was better than the pinchmade to the shroud. The pinch did not affect the operation range of the compressor, but the low solidity vaned diffuser slightly decreased the operation range.The unsteady phenomena in the vaneless diffuser were studied experimentally andnumerically. The unsteady static pressure was measured at the diffuser inlet and outlet, and time-accurate numerical simulation was conducted. The unsteady static pressure showed that most of the pressure variations lay at the passing frequency of every second blade. The pressure variations did not vanish in the diffuser and were visible at the diffuser outlet. However, the amplitude of the pressure variations decreased in the diffuser. The time-accurate calculations showed quite a good agreement with the measured data. Agreement was very good at the design operation point, even though the computational grid was not dense enough inthe volute and in the exit cone. The time-accurate calculation over-predicted the amplitude of the pressure variations at high flow.