Detecting outlying studies in meta-regression models using a forward search algorithm.

When considering data from many trials, it is likely that some of them present a markedly different intervention effect or exert an undue influence on the summary results. We develop a forward search algorithm for identifying outlying and influential studies in meta-analysis models. The forward search algorithm starts by fitting the hypothesized model to a small subset of likely outlier-free studies and proceeds by adding studies into the set one-by-one that are determined to be closest to the fitted model of the existing set. As each study is added to the set, plots of estimated parameters and measures of fit are monitored to identify outliers by sharp changes in the forward plots. We apply the proposed outlier detection method to two real data sets; a meta-analysis of 26 studies that examines the effect of writing-to-learn interventions on academic achievement adjusting for three possible effect modifiers, and a meta-analysis of 70 studies that compares a fluoride toothpaste treatment to placebo for preventing dental caries in children. A simple simulated example is used to illustrate the steps of the proposed methodology, and a small-scale simulation study is conducted to evaluate the performance of the proposed method. Copyright © 2016 John Wiley & Sons, Ltd.

Improved VMAT planning for head and neck tumors with an advanced optimization algorithm

OBJECTIVE In this study, the "Progressive Resolution Optimizer PRO3" (Varian Medical Systems) is compared to the previous version "PRO2" with respect to its potential to improve dose sparing to the organs at risk (OAR) and dose coverage of the PTV for head and neck cancer patients. MATERIALS AND METHODS For eight head and neck cancer patients, volumetric modulated arc therapy (VMAT) treatment plans were generated in this study. All cases have 2-3 phases and the total prescribed dose (PD) was 60-72Gy in the PTV. The study is mainly focused on the phase 1 plans, which all have an identical PD of 54Gy, and complex PTV structures with an overlap to the parotids. Optimization was performed based on planning objectives for the PTV according to ICRU83, and with minimal dose to spinal cord, and parotids outside PTV. In order to assess the quality of the optimization algorithms, an identical set of constraints was used for both, PRO2 and PRO3. The resulting treatment plans were investigated with respect to dose distribution based on the analysis of the dose volume histograms. RESULTS For the phase 1 plans (PD=54Gy) the near maximum dose D2% of the spinal cord, could be minimized to 22±5 Gy with PRO3, as compared to 32±12Gy with PRO2, averaged for all patients. The mean dose to the parotids was also lower in PRO3 plans compared to PRO2, but the differences were less pronounced. A PTV coverage of V95%=97±1% could be reached with PRO3, as compared to 86±5% with PRO2. In clinical routine, these PRO2 plans would require modifications to obtain better PTV coverage at the cost of higher OAR doses. CONCLUSION A comparison between PRO3 and PRO2 optimization algorithms was performed for eight head and neck cancer patients. In general, the quality of VMAT plans for head and neck patients are improved with PRO3 as compared to PRO2. The dose to OARs can be reduced significantly, especially for the spinal cord. These reductions are achieved with better PTV coverage as compared to PRO2. The improved spinal cord sparing offers new opportunities for all types of paraspinal tumors and for re-irradiation of recurrent tumors or second malignancies.

A 22,000 14C year BP sediment and pollen record of climate change from Laguna Miscanti (23°S), northern Chile

Lake sediments and pollen, spores and algae from the high-elevation endorheic Laguna Miscanti (22°45′S, 67°45′W, 4140 m a.s.l., 13.5 km2 water surface, 10 m deep) in the Atacama Desert of northern Chile provide information about abrupt and high amplitude changes in effective moisture. Although the lack of terrestrial organic macrofossils and the presence of a significant 14C reservoir effect make radiocarbon dating of lake sediments very difficult, we propose the following palaeoenvironmental history. An initial shallow freshwater lake (ca. 22,000 14C years BP) disappeared during the extremely dry conditions of the Last Glacial Maximum (LGM; 18,000 14C years BP). That section is devoid of pollen. The late-glacial lake transgression started around 12,000 14C years BP, peaked in two phases between ca. 11,000 and <9000 14C years BP, and terminated around 8000 14C years BP. Effective moisture increased more than three times compared to modern conditions (∼200 mm precipitation), and a relatively dense terrestrial vegetation was established. Very shallow hypersaline lacustrine conditions prevailed during the mid-Holocene until ca. 3600 14C years BP. However, numerous drying and wetting cycles suggest frequent changes in moisture, maybe even individual storms during the mid-Holocene. After several humid spells, modern conditions were reached at ca. 3000 14C years BP. Comparison between limnogeological data and pollen of terrestrial plants suggest century-scale response lags. Relatively constant concentrations of long-distance transported pollen from lowlands east of the Andes suggest similar atmospheric circulation patterns (mainly tropical summer rainfall) throughout the entire period of time. These findings compare favorably with other regional paleoenvironmental data.

Treatment allocation in hepatocellular carcinoma: Assessment of the BCLC algorithm

BACKGROUND AND AIMS The Barcelona Clinic Liver Cancer (BCLC) staging system is the algorithm most widely used to manage patients with hepatocellular carcinoma (HCC). We aimed to investigate the extent to which the BCLC recommendations effectively guide clinical practice and assess the reasons for any deviation from the recommendations. MATERIAL AND METHODS The first-line treatments assigned to patients included in the prospective Bern HCC cohort were analyzed. RESULTS Among 223 patients included in the cohort, 116 were not treated according to the BCLC algorithm. Eighty percent of the patients in BCLC stage 0 (very early HCC) and 60% of the patients in BCLC stage A (early HCC) received recommended curative treatment. Only 29% of the BCLC stage B patients (intermediate HCC) and 33% of the BCLC stage C patients (advanced HCC) were treated according to the algorithm. Eighty-nine percent of the BCLC stage D patients (terminal HCC) were treated with best supportive care, as recommended. In 98 patients (44%) the performance status was disregarded in the stage assignment. CONCLUSION The management of HCC in clinical practice frequently deviates from the BCLC recommendations. Most of the curative therapy options, which have well-defined selection criteria, were allocated according to the recommendations, while the majority of the palliative therapy options were assigned to patients with tumor stages not aligned with the recommendations. The only parameter which is subjective in the algorithm, the performance status, is also the least respected.

Metal Artifact Reduction in Pelvic Computed Tomography With Hip Prostheses: Comparison of Virtual Monoenergetic Extrapolations From Dual-Energy Computed Tomography and an Iterative Metal Artifact Reduction Algorithm in a Phantom Study.

OBJECTIVE The aim of this study was to directly compare metal artifact reduction (MAR) of virtual monoenergetic extrapolations (VMEs) from dual-energy computed tomography (CT) with iterative MAR (iMAR) from single energy in pelvic CT with hip prostheses. MATERIALS AND METHODS A human pelvis phantom with unilateral or bilateral metal inserts of different material (steel and titanium) was scanned with third-generation dual-source CT using single (120 kVp) and dual-energy (100/150 kVp) at similar radiation dose (CT dose index, 7.15 mGy). Three image series for each phantom configuration were reconstructed: uncorrected, VME, and iMAR. Two independent, blinded radiologists assessed image quality quantitatively (noise and attenuation) and subjectively (5-point Likert scale). Intraclass correlation coefficients (ICCs) and Cohen κ were calculated to evaluate interreader agreements. Repeated measures analysis of variance and Friedman test were used to compare quantitative and qualitative image quality. Post hoc testing was performed using a corrected (Bonferroni) P < 0.017. RESULTS Agreements between readers were high for noise (all, ICC ≥ 0.975) and attenuation (all, ICC ≥ 0.986); agreements for qualitative assessment were good to perfect (all, κ ≥ 0.678). Compared with uncorrected images, VME showed significant noise reduction in the phantom with titanium only (P < 0.017), and iMAR showed significantly lower noise in all regions and phantom configurations (all, P < 0.017). In all phantom configurations, deviations of attenuation were smallest in images reconstructed with iMAR. For VME, there was a tendency toward higher subjective image quality in phantoms with titanium compared with uncorrected images, however, without reaching statistical significance (P > 0.017). Subjective image quality was rated significantly higher for images reconstructed with iMAR than for uncorrected images in all phantom configurations (all, P < 0.017). CONCLUSIONS Iterative MAR showed better MAR capabilities than VME in settings with bilateral hip prosthesis or unilateral steel prosthesis. In settings with unilateral hip prosthesis made of titanium, VME and iMAR performed similarly well.

SOMS: SurrOgate MultiStart algorithm for use with nonlinear programming for global optimization

SOMS is a general surrogate-based multistart algorithm, which is used in combination with any local optimizer to find global optima for computationally expensive functions with multiple local minima. SOMS differs from previous multistart methods in that a surrogate approximation is used by the multistart algorithm to help reduce the number of function evaluations necessary to identify the most promising points from which to start each nonlinear programming local search. SOMS’s numerical results are compared with four well-known methods, namely, Multi-Level Single Linkage (MLSL), MATLAB’s MultiStart, MATLAB’s GlobalSearch, and GLOBAL. In addition, we propose a class of wavy test functions that mimic the wavy nature of objective functions arising in many black-box simulations. Extensive comparisons of algorithms on the wavy testfunctions and on earlier standard global-optimization test functions are done for a total of 19 different test problems. The numerical results indicate that SOMS performs favorably in comparison to alternative methods and does especially well on wavy functions when the number of function evaluations allowed is limited.

Paleolimnological evidence for increased landslide activity due to forest clearing and land-use since 3600 cal BP in the western Swiss Alps

Schwarzsee is located in the western Swiss Alps, in a region that has been affected by numerous landslides during the Holocene, as evidenced by geological surveys. Lacustrine sediments were cored to a depth of 13 m. The vegetation history of the lake's catchment was reconstructed and investigated to identify possible impacts on slope stability. The pollen analyses record development of forest cover during the middle and late Holocene, and provide strong evidence for regional anthropogenic influence such as forest clearing and agricultural activity. Vegetation change is characterized by continuous landscape denudation that begins at ca. 4300 cal. yrs BP, with five distinct pulses of increased deforestation, at 3650, 2700, 1500, 900, and 450 cal. yrs BP. Each pulse can be attributed to increased human impact, recorded by the appearance or increase of specific anthropogenic indicator plant taxa. These periods of intensified deforestation also appear to be correlated with increased landslide activity in the lake's catchment and increased turbidite frequency in the sediment record. Therefore, this study gives new evidence for a strong influence of vegetation changes on slope stability during the middle and late Holocene in the western Swiss Alps, and may be used as a case study for anthropogenically induced landslide activity.

CHEOPS Performance for Exomoons: The Detectability of Exomoons by Using Optimal Decision Algorithm

Many attempts have already been made to detect exomoons around transiting exoplanets, but the first confirmed discovery is still pending. The experiences that have been gathered so far allow us to better optimize future space telescopes for this challenge already during the development phase. In this paper we focus on the forthcoming CHaraterising ExOPlanet Satellite (CHEOPS), describing an optimized decision algorithm with step-by-step evaluation, and calculating the number of required transits for an exomoon detection for various planet moon configurations that can be observable by CHEOPS. We explore the most efficient way for such an observation to minimize the cost in observing time. Our study is based on PTV observations (photocentric transit timing variation) in simulated CHEOPS data, but the recipe does not depend on the actual detection method, and it can be substituted with, e.g., the photodynamical method for later applications. Using the current state-of-the-art level simulation of CHEOPS data we analyzed transit observation sets for different star planet moon configurations and performed a bootstrap analysis to determine their detection statistics. We have found that the detection limit is around an Earth-sized moon. In the case of favorable spatial configurations, systems with at least a large moon and a Neptune-sized planet, an 80% detection chance requires at least 5-6 transit observations on average. There is also a nonzero chance in the case of smaller moons, but the detection statistics deteriorate rapidly, while the necessary transit measurements increase quickly. After the CoRoT and Kepler spacecrafts, CHEOPS will be the next dedicated space telescope that will observe exoplanetary transits and characterize systems with known Doppler-planets. Although it has a smaller aperture than Kepler (the ratio of the mirror diameters is about 1/3) and is mounted with a CCD that is similar to Kepler's, it will observe brighter stars and operate with larger sampling rate; therefore, the detection limit for an exomoon can be the same as or better, which will make CHEOPS a competitive instruments in the quest for exomoons.

A first chironomid-based summer temperature reconstruction (13–5 ka BP) around 49°N in inland Europe compared with local lake development

Temperature reconstructions for the end of the Pleistocene and the first half of the Holocene based on biotic proxies are rare for inland Europe around 49°N. We analysed a 7 m long sequence of lake deposits in the Vihorlat Mts in eastern Slovakia (820 m a.s.l.). Chironomid head capsules were used to reconstruct mean July temperature (TJuly), other proxies (diatoms, green algae, pollen, geochemistry) were used to reconstruct local environmental changes that might have affected the climate reconstruction, such as epilimnetic total phosphorus concentrations (TP), lake level changes and development of surrounding vegetation. During the Younger Dryas (YD), temperature fluctuated between 7 and 11 °C, with distinct, decadal to centennial scale variations, that agree with other palaeoclimate records in Europe such as δ18O content in stalagmites or Greenland ice cores. The results indicate that the site was somewhat colder than expected from the general south-to-north YD temperature gradient within Europe, possibly because of north-facing exposition. The warmer phases of the YD were characterised by low water level or even complete desiccation of the lake (12,200-12,400 cal yr BP). At the Late-Glacial/Holocene transition TJuly steeply increased from from 11 to 15.5 °C (11,700-11,400 cal yr BP) - the highest TJuly for entire sequence. This rapid climate change was reflected by all proxies as a compositional change and increasing species diversity. The open woodlands of Pinus, Betula, Larix and Picea were replaced by broad-leaved temperate forests dominated by Betula, later by Ulmus and finally by Corylus (ca 9700 cal yr BP). At the same time, input of eroded coarse-grained material into the lake decreased and organic matter (LOI) and biogenic silica increased. The Early-Holocene climate was rather stable till 8700 cal yr BP, with temporary decrease in TJuly around 11,200 cal yr BP. The lake was productive with a well-developed littoral, as indicated by both diatoms and chironomids. A distinct decline of TJuly to 10 °C between 8700 and 8000 cal yr BP was associated with decreasing chironomid diversity and increasing climate moistening indicated by pollen. Tychoplanktonic and phosphorus-demanding diatoms increased which might be explained by hydrological and land-cover changes. Later, a gradual warming started after 7000 cal yr BP and representation of macrophytes, periphytic diatoms and littoral chironomids increased. Our results suggest that the Holocene thermal maximum was taking place unusually early in the Holocene at our study site, but its timing might be affected by topography and mesoclimate. We further demonstrated that temperature changes had coincided with variations in local hydrology

History of Atmospheric Lead Deposition Since 12,370 14C yr BP from a Peat Bog, Jura Mountains, Switzerland

A continuous record of atmospheric lead since 12,370 carbon-14 years before the present (14C yr BP) is preserved in a Swiss peat bog. Enhanced fluxes caused by climate changes reached their maxima 10,590 14C yr BP (Younger Dryas) and 823014C yr BP. Soil erosion caused by forest clearing and agricultural tillage increased lead deposition after 532014C yr BP. Increasing lead/scandium and decreasing lead-206/lead-207 beginning 3000 14C yr BP indicate the beginning of lead pollution from mining and smelting, and anthropogenic sources have dominated lead emissions ever since. The greatest lead flux (15.7 milligrams per square meter per year in A.D. 1979) was 1570 times the natural, background value (0.01 milligram per square meter per year from 8030 to 5320 14C yr BP).