A Functional-Based Distribution Diagnostic for a Linear Model with Correlated Outcomes: Technical Report

Despite the widespread popularity of linear models for correlated outcomes (e.g. linear mixed modesl and time series models), distribution diagnostic methodology remains relatively underdeveloped in this context. In this paper we present an easy-to-implement approach that lends itself to graphical displays of model fit. Our approach involves multiplying the estimated marginal residual vector by the Cholesky decomposition of the inverse of the estimated marginal variance matrix. Linear functions or the resulting "rotated" residuals are used to construct an empirical cumulative distribution function (ECDF), whose stochastic limit is characterized. We describe a resampling technique that serves as a computationally efficient parametric bootstrap for generating representatives of the stochastic limit of the ECDF. Through functionals, such representatives are used to construct global tests for the hypothesis of normal margional errors. In addition, we demonstrate that the ECDF of the predicted random effects, as described by Lange and Ryan (1989), can be formulated as a special case of our approach. Thus, our method supports both omnibus and directed tests. Our method works well in a variety of circumstances, including models having independent units of sampling (clustered data) and models for which all observations are correlated (e.g., a single time series).

The Importance of Scale for Spatial-confounding Bias and Precision of Spatial Regression Estimators

Increasingly, regression models are used when residuals are spatially correlated. Prominent examples include studies in environmental epidemiology to understand the chronic health effects of pollutants. I consider the effects of residual spatial structure on the bias and precision of regression coefficients, developing a simple framework in which to understand the key issues and derive informative analytic results. When the spatial residual is induced by an unmeasured confounder, regression models with spatial random effects and closely-related models such as kriging and penalized splines are biased, even when the residual variance components are known. Analytic and simulation results show how the bias depends on the spatial scales of the covariate and the residual; bias is reduced only when there is variation in the covariate at a scale smaller than the scale of the unmeasured confounding. I also discuss how the scales of the residual and the covariate affect efficiency and uncertainty estimation when the residuals can be considered independent of the covariate. In an application on the association between black carbon particulate matter air pollution and birth weight, controlling for large-scale spatial variation appears to reduce bias from unmeasured confounders, while increasing uncertainty in the estimated pollution effect.

IMPROVING STARTABILITY AND REDUCING EMISSIONS IN FLEXFUEL SPARK IGNITION DIRECT INJECTION VARIABLE CAM TIMING ENGINE

Experimental work and analysis was done to investigate engine startup robustness and emissions of a flex-fuel spark ignition (SI) direct injection (DI) engine. The vaporization and other characteristics of ethanol fuel blends present a challenge at engine startup. Strategies to reduce the enrichment requirements for the first engine startup cycle and emissions for the second and third fired cycle at 25°C ± 1°C engine and intake air temperature were investigated. Research work was conducted on a single cylinder SIDI engine with gasoline and E85 fuels, to study the effect on first fired cycle of engine startup. Piston configurations that included a compression ratio change (11 vs 15.5) and piston geometry change (flattop vs bowl) were tested, along with changes in intake cam timing (95,110,125) and fuel pressure (0.4 MPa vs 3 MPa). The goal was to replicate the engine speed, manifold pressure, fuel pressure and testing temperature from an engine startup trace for investigating the first fired cycle for the engine. Results showed bowl piston was able to enable lower equivalence ratio engine starts with gasoline fuel, while also showing lower IMEP at the same equivalence ratio compared to flat top piston. With E85, bowl piston showed reduced IMEP as compression ratio increased at the same equivalence ratio. A preference for constant intake valve timing across fuels seemed to indicate that flattop piston might be a good flex-fuel piston. Significant improvements were seen with higher CR bowl piston with high fuel pressure starts, but showed no improvement with low fuel pressures. Simulation work was conducted to analyze initial three cycles of engine startup in GT-POWER for the same set of hardware used in the experimentations. A steady state validated model was modified for startup conditions. The results of which allowed an understanding of the relative residual levels and IMEP at the test points in the cam phasing space. This allowed selecting additional test points that enable use of higher residual levels, eliminating those with smaller trapped mass incapable of producing required IMEP for proper engine turnover. The second phase of experimental testing results for 2nd and 3rd startup cycle revealed both E10 and E85 prefer the same SOI of 240°bTDC at second and third startup cycle for the flat top piston and high injection pressures. E85 fuel optimal cam timing for startup showed that it tolerates more residuals compared to E10 fuel. Higher internal residuals drives down the Ø requirement for both fuels up to their combustion stability limit, this is thought to be direct benefit to vaporization due to increased cycle start temperature. Benefits are shown for an advance IMOP and retarded EMOP strategy at engine startup. Overall the amount of residuals preferred by an engine for E10 fuel at startup is thought to be constant across engine speed, thus could enable easier selection of optimized cam positions across the startup speeds.

Direct MR arthrography at 1.5 and 3.0 T: signal dependence on gadolinium and iodine concentrations--phantom study

PURPOSE: To prospectively quantify in vitro the influence of gadopentetate dimeglumine and ioversol on the magnetic resonance (MR) imaging signal observed with a variety of musculoskeletal pulse sequences to predict optimum gadolinium concentrations for direct MR arthrography at 1.5 and 3.0 T. MATERIALS AND METHODS: In an in vitro study, T1 and T2 relaxation times of three dilution series of gadopentetate dimeglumine (concentration, 0-20.0 mmol gadolinium per liter) at ioversol concentrations with iodine concentration of 0, 236.4, and 1182 mmol iodine per liter (corresponding to 0, 30, and 150 mg of iodine per milliliter) were measured at 1.5 and 3.0 T. The relaxation rate dependence on concentrations of gadolinium and iodine was analytically modeled, and continuous profiles of signal versus gadolinium concentration were calculated for 10 pulse sequences used in current musculoskeletal imaging. After fitting to experimental discrete profiles, maximum signal-to-noise ratio (SNR), gadolinium concentration with maximum SNR, and range of gadolinium concentration with 90% of maximum SNR were derived. The overall influence of field strength and iodine concentration on these parameters was assessed by using t tests. The deviation of simulated from experimental signal-response profiles was assessed with the autocorrelation of the residuals. RESULTS: The model reproduced relaxation rates of 0.37-38.24 sec(-1), with a mean error of 4.5%. Calculated SNR profiles matched the discrete experimental profiles, with autocorrelation of the residuals divided by the mean of less than 5.0. Admixture of ioversol consistently reduced T1 and T2, narrowed optimum gadolinium concentration ranges (P = .004-.006), and reduced maximum SNR (P < .001 to not significant). Optimum gadolinium concentration was 0.7-3.4 mmol/L at both field strengths. At 3.0 T, maximum SNR was up to 75% higher than at 1.5 T. CONCLUSION: Admixture of ioversol to gadopentetate dimeglumine solutions results in a consistent additional relaxation enhancement, which can be analytically modeled to allow a near-quantitative a priori optimized match of contrast media concentrations and imaging protocol for a broad variety of pulse sequences.

Differences between dual X-ray absorptiometry using pencil beam and fan beam modes and their determinants in vivo and in vitro

The aim of this study was to determine the influence of individual factors on differences in bone mineral density (BMD) using dual X-ray absorptiometry pencil beam (PB) and fan beam (FB) modes in vivo and in vitro. PB.BMD and FB.BMD of 63 normal Caucasian females ages 21-80 yr were measured at the lumbar spine and hip. Residuals of the FB/PB regression were used to assess the impact of height, weight, adiposity index (AI) (= weight/height(3/2)), back tissue thickness, and PB.BMD, respectively, on FB/PB difference. The Hologic Anthropomorphic Spine Phantom (ASP) was measured using the PB and FB modes at two different levels to assess the impact of scanning mode and focus distance. The European Spine Phantom (ESP) prototype, a geometrically well-defined phantom with known vertebral densities, was measured using PB and FB modes and analyzed manually to determine the impact of bone density on FB/PB difference and automatically to determine the impact of edge detection on FB/PB difference. Population BMD results were perfectly correlated, but significantly overestimated by 1.5% at the lumbar spine and underestimated by 0.7% at the neck, 1.8% at the trochanter, and 2.0% at the total hip, respectively, when using the FB compared with PB mode. At the lumbar spine, the FB/PB residual correlated negatively with height (r = 0.34, p < 0.01) and PB.BMD (r = 0.48, p <: 0. 0001) and positively with AI (r = 0.26, p < 0.05). At the hip, residual of trochanter correlated positively with weight (r = 0.36, p < 0.01) and AI (r = 0.36, p < 0.01). The FB mode significantly increased ASP BMD by 0.7% compared with PB. Using the FB mode, increasing focus distance significantly (p < 0.001) decreased area and bone mineral content, but not BMD. By contrast, increasing focus distance significantly decreased PB.BMD by 0.7%. With the ESP, the PB mode supplied accurate projected are of the bone (AREA) results but significant underestimation of specified BMD in the manual analysis. The FB mode significantly underestimated PB. AREA by 2.9% but fitted specified BMD quite well. FB/PB overestimation was larger for the low-density (+8.7%) than for the high-density vertebra (+4. 9%). The automated analysis resulted in more than 14% underestimation of PB. AREA (low-density vertebra) and an almost 13% overestimation of PB.BMD (high-density vertebra) using FB. In conclusion, FB and PB measurements are highly correlated at the lumbar spine and hip with small but significant BMD differences related to height, adiposity, and BMD. In clinical practice, it can be erroneous to switch from one method to another, especially in women with low bone density.

A systematic approach to analyse the sequelae of LCPD

The analysis and treatment of hips with healed Legg-Calvé-Perthes disease (LCPD) differs substantially from the treatment in the acute phase of the disease. More specifically, the treating orthopaedic surgeon is often faced with a complex three-dimensional pathomorphology of the hip that is difficult to understand and correct. To date, none of the current classification systems provide a useful decision-making algorithm with regards to the type of surgical intervention necessary to improve hip function in patients with sequelae of LCPD. The conceptual recognition of the femoroacetabular impingement (FAI) and the ability to safely dislocate the hip have revolutionised our diagnostic and therapeutic algorithm for joint-preserving surgery of hips with structural residuals of LCPD. We present a systematic approach to analyse femoral and acetabular pathomorphologic features. The resulting pathomechanisms and the surgical treatment options are presented.

The effects of thoracic epidurally administered drugs on urethral sphincter function in women: a pooled analysis

OBJECTIVES Thoracic epidural analgesia (TEA) has been shown to inhibit detrusor activity in patients undergoing open renal surgery, resulting in clinically relevant post-void residuals. However, the impact of different epidural drug mixtures on urethral sphincter function is not completely elucidated. DESIGN Pooled analysis of an open observational study and a double-blind randomized trial. SETTING Single tertiary centre. SUBJECTS Twenty-eight women without lower urinary tract symptoms and post-void residual <100 mL, who underwent open renal surgery with TEA. METHODS Pooling results in three groups with different epidural regimens (7 with bupivacaine 0.125%, 8 with bupivacaine 0.125% and fentanyl 2 μg/mL, and 13 with bupivacaine 0.1% plus fentanyl 2 μg/mL and epinephrine 2 μg/mL). All women underwent urethral pressure measurements before TEA and during TEA 2-3 days postoperatively. All patients received a TEA placed at the insertion site interspace T 8-9. RESULTS Maximum urethral closure pressure at rest decreased significantly during TEA with bupivacaine alone (median 70 cm H2 O [interquartile range 66-76] to 43 [43-65], P = 0.031) and with bupivacaine/fentanyl/epinephrine (75 cm H2 O [68-78] to 56 [52-75], P = 0.028), whereas with bupivacaine/fentanyl, no significant change could be detected (74 [51-88] vs 67 [46-70], P = 0.156). In all groups, functional profile length at rest was not influenced during TEA. CONCLUSION TEA with bupivacaine and the addition of fentanyl and epinephrine appears to decrease maximum urethral closure pressure at rest in women. The addition of fentanyl alone to bupivacaine may reduce this effect. Thus, the TEA effect on urethral sphincter function seems to depend on the drug mixture administered.

Joint Reprocessing of GPS, GLONASS and SLR

A joint reprocessing of GPS, GLONASS and SLR observations has been carried out at TU Dresden, TU Munich, AIUB and ETH Zurich. Common a priori models have been applied for the processing of all types of observation to ensure both consistent parameter estimates and the rigorous combination of microwave and optical measurements. Based on that reprocessing results, we evaluate the impact of adding GLONASS observations to the standard GPS data processing. In particular, changes in station position time series and day boundary overlaps of consecutive satellite arcs are analyzed. In addition, the GNSS orbits derived from microwave measurements are validated using independent SLR range measurements. Our SLR residuals indicate a significant improvement compared to previous results. Furthermore, we evaluate the performance of our high-rate (30s) combined GNSS satellite clocks and discuss associated zero-difference phase residuals.

Small aliquot and single grain IRSL and post-IR IRSL dating of fluvial and alluvial sediments from the Pativilca valley, Peru

Infrared stimulated luminescence (IRSL) and post-IR IRSL are applied to small aliquots and single grains to determine the equivalent dose (De) of eleven alluvial and fluvial sediment samples collected in the Pativilca valley, Central Peru at ca. 10°S latitude. Small aliquot De distributions are rather symmetric and display over-dispersion values between 15 and 46%. Small aliquot g-values range between 4 and 8% per decade for the IRSL and 1 and 2% per decade for the post-IR IRSL signal. The single grain De distributions are highly over-dispersed with some of them skewed to higher doses, implying partial bleaching; this is especially true for the post-IR IRSL. Measurements of a modern analog reveal that residuals due to partial bleaching are present in both the IRSL as well as the post-IR IRSL signal. The g-values of individual grains exhibit a wide range with high individual uncertainties and might contribute significantly to the spread of the single grain De values, at least for the IRSL data. Electron Microprobe Analysis performed on single grains reveal that a varying K-content can be excluded as the origin of over-dispersion. Final ages for the different approaches are calculated using the Central Age Model and the Minimum Age Model (MAM). The samples are grouped into well-beached, potentially well-bleached and partially bleached according to the evaluation of the single grain distributions and the agreement of age estimates between methods. The application of the MAM to the single grain data resulted in consistent age estimates for both the fading corrected IRSL and the post-IR IRSL ages, and suggests that both approaches are suitable for dating these samples. Keywords

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

Systematic analysis of the intravoxel incoherent motion threshold separating perfusion and diffusion effects: Proposal of a standardized algorithm

PURPOSE To systematically evaluate the dependence of intravoxel-incoherent-motion (IVIM) parameters on the b-value threshold separating the perfusion and diffusion compartment, and to implement and test an algorithm for the standardized computation of this threshold. METHODS Diffusion weighted images of the upper abdomen were acquired at 3 Tesla in eleven healthy male volunteers with 10 different b-values and in two healthy male volunteers with 16 different b-values. Region-of-interest IVIM analysis was applied to the abdominal organs and skeletal muscle with a systematic increase of the b-value threshold for computing pseudodiffusion D*, perfusion fraction Fp , diffusion coefficient D, and the sum of squared residuals to the bi-exponential IVIM-fit. RESULTS IVIM parameters strongly depended on the choice of the b-value threshold. The proposed algorithm successfully provided optimal b-value thresholds with the smallest residuals for all evaluated organs [s/mm2]: e.g., right liver lobe 20, spleen 20, right renal cortex 150, skeletal muscle 150. Mean D* [10(-3) mm(2) /s], Fp [%], and D [10(-3) mm(2) /s] values (±standard deviation) were: right liver lobe, 88.7 ± 42.5, 22.6 ± 7.4, 0.73 ± 0.12; right renal cortex: 11.5 ± 1.8, 18.3 ± 2.9, 1.68 ± 0.05; spleen: 41.9 ± 57.9, 8.2 ± 3.4, 0.69 ± 0.07; skeletal muscle: 21.7 ± 19.0; 7.4 ± 3.0; 1.36 ± 0.04. CONCLUSION IVIM parameters strongly depend upon the choice of the b-value threshold used for computation. The proposed algorithm may be used as a robust approach for IVIM analysis without organ-specific adaptation. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 % ; in the absence of the polar data gap, SLR residuals agree by ≈96 % . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.

Galileo Orbit and Clock Quality of the IGS Multi-GNSS Experiment

The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) aims at the data collection and analysis of all available satellite navigation systems. In particular the new global and regional satellite navigation systems are of interest, i.e., the European Galileo, the Chinese BeiDou, the Japanese QZSS as well as satellite based augmentation systems. This article analyzes the orbit and clock quality of the Galileo products of four MGEX analysis centers for a common time period of 20 weeks. Orbit comparisons of the individual analysis centers have a consistency at the 5–30 cm level. Day boundary discontinuities range from 4 to 28 cm whereas 2-day orbit fit RMS values vary between 1 and 7 cm. The accuracy evaluated by satellite laser ranging residuals is on the one decimeter level with a systematic bias of about −5 cm for all analysis centers. In addition, systematic errors on the decimeter level related to solar radiation pressure mismodeling are present in all orbit products. Due to the correlation of radial orbit errors with the clock parameters, these errors are also visible as a bump in the Allan deviation of the Galileo satellite clocks at the orbital frequency.