New Algorithms for M-Estimation of Multivariate Scatter and Location

We present new algorithms for M-estimators of multivariate scatter and location and for symmetrized M-estimators of multivariate scatter. The new algorithms are considerably faster than currently used fixed-point and related algorithms. The main idea is to utilize a second order Taylor expansion of the target functional and to devise a partial Newton-Raphson procedure. In connection with symmetrized M-estimators we work with incomplete U-statistics to accelerate our procedures initially.

Development of a low energy ion source for ROSINA ion mode calibration

The European Rosetta mission on its way to comet 67P/Churyumov-Gerasimenko will remain for more than a year in the close vicinity (1 km) of the comet. The two ROSINA mass spectrometers on board Rosetta are designed to analyze the neutral and ionized volatile components of the cometary coma. However, the relative velocity between the comet and the spacecraft will be minimal and also the velocity of the outgassing particles is below 1km∕s. This combination leads to very low ion energies in the surrounding plasma of the comet, typically below 20eV. Additionally, the spacecraft may charge up to a few volts in this environment. In order to simulate such plasma and to calibrate the mass spectrometers, a source for ions with very low energies had to be developed for the use in the laboratory together with the different gases expected at the comet. In this paper we present the design of this ion source and we discuss the physical parameters of the ion beam like sensitivity, energy distribution, and beam shape. Finally, we show the first ion measurements that have been performed together with one of the two mass spectrometers.

Introduction to the Maïdo Lidar Calibration Campaign dedicated to the validation of upper air meteorological parameters

The first operations at the new High-altitude Maïdo Observatory at La Réunion began in 2013. The Maïdo Lidar Calibration Campaign (MALICCA) was organized there in April 2013 and has focused on the validation of the thermodynamic parameters (temperature, water vapor, and wind) measured with many instruments including the new very large lidar for water vapor and temperature profiles. The aim of this publication consists of providing an overview of the different instruments deployed during this campaign and their status, some of the targeted scientific questions and associated instrumental issues. Some specific detailed studies for some individual techniques were addressed elsewhere. This study shows that temperature profiles were obtained from the ground to the mesopause (80 km) thanks to the lidar and regular meteorological balloon-borne sondes with an overlap range showing good agreement. Water vapor is also monitored from the ground to the mesopause by using the Raman lidar and microwave techniques. Both techniques need to be pushed to their limit to reduce the missing range in the lower stratosphere. Total columns obtained from global positioning system or spectrometers are valuable for checking the calibration and ensuring vertical continuity. The lidar can also provide the vertical cloud structure that is a valuable complementary piece of information when investigating the water vapor cycle. Finally, wind vertical profiles, which were obtained from sondes, are now also retrieved at Maïdo from the newly implemented microwave technique and the lidar. Stable calibrations as well as a small-scale dynamical structure are required to monitor the thermodynamic state of the middle atmosphere, ensure validation of satellite sensors, study the transport of water vapor in the vicinity of the tropical tropopause and study their link with cirrus clouds and cyclones and the impact of small-scale dynamics (gravity waves) and their link with the mean state of the mesosphere.

Sensitivity and Fragmentation Calibration of the ROSINA Reflectron-type Time-Of-Flight Mass Spectrometer

The goal of this work was to increase the performance and to calibrate one of the ROSINA sensors, the Reflectron-type Time-Of-Flight mass spectrometer, currently flying aboard the ESA Rosetta spacecraft. Different optimization techniques were applied to both the lab and space models, and a static calibration was performed using different gas species expected to be detected in the vicinity of comet 67P/Churyumov-Gerasimenko. The database thus created was successfully applied to space data, giving consistent results with the other ROSINA sensors.

Influence of prosthesis type and retention mechanism on complications with fixed implant-supported prostheses: a systematic review applying multivariate analyses

PURPOSE To identify the influence of fixed prosthesis type on biologic and technical complication rates in the context of screw versus cement retention. Furthermore, a multivariate analysis was conducted to determine which factors, when considered together, influence the complication and failure rates of fixed implant-supported prostheses. MATERIALS AND METHODS Electronic searches of MEDLINE (PubMed), EMBASE, and the Cochrane Library were conducted. Selected inclusion and exclusion criteria were used to limit the search. Data were analyzed statistically with simple and multivariate random-effects Poisson regressions. RESULTS Seventy-three articles qualified for inclusion in the study. Screw-retained prostheses showed a tendency toward and significantly more technical complications than cemented prostheses with single crowns and fixed partial prostheses, respectively. Resin chipping and ceramic veneer chipping had high mean event rates, at 10.04 and 8.95 per 100 years, respectively, for full-arch screwed prostheses. For "all fixed prostheses" (prosthesis type not reported or not known), significantly fewer biologic and technical complications were seen with screw retention. Multivariate analysis revealed a significantly greater incidence of technical complications with cemented prostheses. Full-arch prostheses, cantilevered prostheses, and "all fixed prostheses" had significantly higher complication rates than single crowns. A significantly greater incidence of technical and biologic complications was seen with cemented prostheses. CONCLUSION Screw-retained fixed partial prostheses demonstrated a significantly higher rate of technical complications and screw-retained full-arch prostheses demonstrated a notably high rate of veneer chipping. When "all fixed prostheses" were considered, significantly higher rates of technical and biologic complications were seen for cement-retained prostheses. Multivariate Poisson regression analysis failed to show a significant difference between screw- and cement-retained prostheses with respect to the incidence of failure but demonstrated a higher rate of technical and biologic complications for cement-retained prostheses. The incidence of technical complications was more dependent upon prosthesis and retention type than prosthesis or abutment material.

Sensitivity and fragmentation calibration of the ROSINA double focusing mass spectrometer

The goal of this work has been to calibrate sensitivities and fragmentation pattern of various molecules as well as further characterize the lab model of the ROSINA Double Focusing Mass Spectrometer (DFMS) on board ESA’s Rosetta spacecraft bound to comet 67P/Churyumov-Gerasimenko. The detailed calibration and characterization of the instrument is key to understand and interpret the results in the coma of the comet. A static calibration was performed for the following species: Ne, Ar, Kr, Xe, H2O, N2, CO2, CH4, C2H6, C3H8, C4H10, and C2H4. The purpose of the calibration was to obtain sensitivities for all detectors and emissions, the fragmentation behavior of the ion source and to show the capabilities to measure isotopic ratios at the comet. The calibration included the recording of different correction factors to evaluate the data, including a detailed investigation of the detector gain. The quality of the calibration that could be tested for different gas mixtures including the calibration of the density inside the ion source when calibration gas from the gas calibration unit is introduced. In conclusion the calibration shows that DFMS meets the design requirements and that DFMS will be able to measure the D/H at the comet and help shed more light on the puzzle about the origin of water on Earth.

Calibration of cosmogenic noble gas production based on ³⁶Cl-³⁶Ar ages. Part 2. The ⁸¹Kr-Kr dating technique

We calibrated the ⁸¹Kr-Kr dating system for ordinary chondrites of different sizes using independent shielding-corrected ³⁶Cl-³⁶Ar ages. Krypton concentrations and isotopic compositions were measured in bulk samples from 14 ordinary chondrites of high petrologic type and the cosmogenic Kr component was obtained by subtracting trapped Kr from phase Q. The thus-determined average cosmogenic ⁷⁸Kr/⁸³Kr, ⁸⁰Kr/⁸³Kr, ⁸²Kr/⁸³Kr, and ⁸4Kr/⁸³Kr ratiC(Lavielle and Marti 1988; Wieler 2002). The cosmogenic ⁷⁸Kr/⁸³Kr ratio is correlated with the cosmogenic 22Ne/21Ne ratio, confirming that ⁷⁸Kr/⁸³Kr is a reliable shielding indicator. Previously, ⁸¹Kr-Kr ages have been determined by assuming the cosmogenic production rate of ⁸¹Kr, P(⁸¹Kr)c, to be 0.95 times the average of the cosmogenic production rates of ⁸⁰Kr and ⁸²Kr; the factor Y = 0.95 therefore accounts for the unequal production of the various Kr isotopes (Marti 1967a). However, Y should be regarded as an empirical adjustment. For samples whose ⁸⁰Kr and ⁸²Kr concentrations may be affected by neutron-capture reactions, the shielding-dependent cosmogenic (⁷⁸Kr/⁸³Kr)c ratio has been used instead to calculate P(⁸¹Kr)/P(⁸³Kr), as for some lunar samples, this ratio has been shown to linearly increase with (⁷⁸Kr/⁸³Kr)c (Marti and Lugmair 1971). However, the ⁸¹Kr-Kr ages of our samples calculated with these methods are on average ~30% higher than their ³⁶Cl-³⁶Ar ages, indicating that most if not all the ⁸¹Kr-Kr ages determined so far are significantly too high. We therefore re-evaluated both methods to determine P(⁸¹Kr)c/P(⁸³Kr)c. Our new Y value of 0.70 ± 0.04 is more than 25% lower than the value of 0.95 used so far. Furthermore, together with literature data, our data indicate that for chondrites, P(⁸¹Kr)c/P(⁸³Kr)c is rather constant at 0.43 ± 0.02, at least for the shielding range covered by our samples ([⁷⁸Kr/⁸³Kr]c = 0.119–0.185; [22Ne/21Ne]c = 1.083–1.144), in contrast to the observations on lunar samples. As expected considering the method used, ⁸¹Kr-Kr ages calculated either directly with this new P(⁸¹Kr)c/P(⁸³Kr)c value or with our new Y value both agree with the corresponding ³⁶Cl-³⁶Ar ages. However, the average deviation of 2% indicates the accuracy of both new ⁸¹Kr-Kr dating methods and the precision of the new dating systems of ~10% is demonstrated by the low scatter in the data. Consequently, this study indicates that the ⁸¹Kr-Kr ages published so far are up to 30% too high.

Design and Implementation of a Laparoscope Calibration Method for Augmented Reality Navigation

Intraoperative laparoscopic calibration remains a challenging task. In this work we present a new method and instrumentation for intraoperative camera calibration. Contrary to conventional calibration methods, the proposed technique allows intraoperative laparoscope calibration from single perspective observations, resulting in a standardized scheme for calibrating in a clinical scenario. Results show an average displacement error of 0.52 ± 0.19 mm, indicating sufficient accuracy for clinical use. Additionally, the proposed method is validated clinically by performing a calibration during the surgery.

Stable isotope response to lake eutrophication: Calibration of a high-resolution lacustrine sequence from Baldeggersee, Switzerland

Stable isotope analyses of discrete seasonal layers from a 108-yr annually laminated freeze-core from Baldeg-gersee, a small, eutrophic lake in central Switzerland, provide information on the climatological and environmental factors, including lake eutrophication, that control oxygen and carbon isotopic composition of epilimnic biologically induced calcite precipitate. During the last 100 yr, Baldeggersee has undergone major increases in productivity and eutrophication in response to nutrient loading from agriculture and industrialization in the lake's watershed. Calibration of the isotopic signal in Baldeggersee to historical limnological data quantitatively links evidence of isotopic depletion in the sedimented calcite to trophic state of the lake. δ18O values from the spring/summer “light” sediment layers steadily diverged to more depleted values in response to historical eutrophication: measured δ18O values were up to 21.5‰ more negative than calculated equilibrium δ18O values. Evidence for 13C depletion in the calcite, relative to equilibrium values, is more difficult to ascertain because of an overall dominance of isotopic enrichment in the dissolved inorganic pool as productivity in Baldeggersee increases. A positive association exists between the degree of oxygen-18 depletion and the calcite crystal size. Thus, large amorphous calcite grains can be used as a proxy for recognizing apparent isotopic nonequilibrium in sediment sequences from highly productive lacustrine environments from all geologic time scales. In contrast to the light layers, the oxygen isotopic composition of the calcite in the late summer/fall “dark” sediment layers is unaffected by the apparent isotope nonequilibrium. Oxygen and carbon isotope values from the dark laminae in the Baldeggersee sediment therefore provide environmental and climatological proxies that can be calibrated with known environmental and regional climate data for the last century.

A multivariate frailty model for disease recurrences and survival

A multivariate frailty hazard model is developed for joint-modeling of three correlated time-to-event outcomes: (1) local recurrence, (2) distant recurrence, and (3) overall survival. The term frailty is introduced to model population heterogeneity. The dependence is modeled by conditioning on a shared frailty that is included in the three hazard functions. Independent variables can be included in the model as covariates. The Markov chain Monte Carlo methods are used to estimate the posterior distributions of model parameters. The algorithm used in present application is the hybrid Metropolis-Hastings algorithm, which simultaneously updates all parameters with evaluations of gradient of log posterior density. The performance of this approach is examined based on simulation studies using Exponential and Weibull distributions. We apply the proposed methods to a study of patients with soft tissue sarcoma, which motivated this research. Our results indicate that patients with chemotherapy had better overall survival with hazard ratio of 0.242 (95% CI: 0.094 - 0.564) and lower risk of distant recurrence with hazard ratio of 0.636 (95% CI: 0.487 - 0.860), but not significantly better in local recurrence with hazard ratio of 0.799 (95% CI: 0.575 - 1.054). The advantages and limitations of the proposed models, and future research directions are discussed. ^

Multivariate methods for correlated effects sizes

Current statistical methods for estimation of parametric effect sizes from a series of experiments are generally restricted to univariate comparisons of standardized mean differences between two treatments. Multivariate methods are presented for the case in which effect size is a vector of standardized multivariate mean differences and the number of treatment groups is two or more. The proposed methods employ a vector of independent sample means for each response variable that leads to a covariance structure which depends only on correlations among the $p$ responses on each subject. Using weighted least squares theory and the assumption that the observations are from normally distributed populations, multivariate hypotheses analogous to common hypotheses used for testing effect sizes were formulated and tested for treatment effects which are correlated through a common control group, through multiple response variables observed on each subject, or both conditions.^ The asymptotic multivariate distribution for correlated effect sizes is obtained by extending univariate methods for estimating effect sizes which are correlated through common control groups. The joint distribution of vectors of effect sizes (from $p$ responses on each subject) from one treatment and one control group and from several treatment groups sharing a common control group are derived. Methods are given for estimation of linear combinations of effect sizes when certain homogeneity conditions are met, and for estimation of vectors of effect sizes and confidence intervals from $p$ responses on each subject. Computational illustrations are provided using data from studies of effects of electric field exposure on small laboratory animals. ^

THE USE OF MULTIVARIATE ANALYSIS OF COVARIANCE IN EVALUATING THE COMPARABILITY OF LABORATORY DETERMINATIONS IN COOPERATIVE STUDIES

The role of clinical chemistry has traditionally been to evaluate acutely ill or hospitalized patients. Traditional statistical methods have serious drawbacks in that they use univariate techniques. To demonstrate alternative methodology, a multivariate analysis of covariance model was developed and applied to the data from the Cooperative Study of Sickle Cell Disease.^ The purpose of developing the model for the laboratory data from the CSSCD was to evaluate the comparability of the results from the different clinics. Several variables were incorporated into the model in order to control for possible differences among the clinics that might confound any real laboratory differences.^ Differences for LDH, alkaline phosphatase and SGOT were identified which will necessitate adjustments by clinic whenever these data are used. In addition, aberrant clinic values for LDH, creatinine and BUN were also identified.^ The use of any statistical technique including multivariate analysis without thoughtful consideration may lead to spurious conclusions that may not be corrected for some time, if ever. However, the advantages of multivariate analysis far outweigh its potential problems. If its use increases as it should, the applicability to the analysis of laboratory data in prospective patient monitoring, quality control programs, and interpretation of data from cooperative studies could well have a major impact on the health and well being of a large number of individuals. ^