Meta-analysis: a unifying meta-likelihood approach framing unobserved heterogeneity, study covariates, publication bias, and study quality

OBJECTIVES: This contribution provides a unifying concept for meta-analysis integrating the handling of unobserved heterogeneity, study covariates, publication bias and study quality. It is important to consider these issues simultaneously to avoid the occurrence of artifacts, and a method for doing so is suggested here. METHODS: The approach is based upon the meta-likelihood in combination with a general linear nonparametric mixed model, which lays the ground for all inferential conclusions suggested here. RESULTS: The concept is illustrated at hand of a meta-analysis investigating the relationship of hormone replacement therapy and breast cancer. The phenomenon of interest has been investigated in many studies for a considerable time and different results were reported. In 1992 a meta-analysis by Sillero-Arenas et al. concluded a small, but significant overall effect of 1.06 on the relative risk scale. Using the meta-likelihood approach it is demonstrated here that this meta-analysis is due to considerable unobserved heterogeneity. Furthermore, it is shown that new methods are available to model this heterogeneity successfully. It is argued further to include available study covariates to explain this heterogeneity in the meta-analysis at hand. CONCLUSIONS: The topic of HRT and breast cancer has again very recently become an issue of public debate, when results of a large trial investigating the health effects of hormone replacement therapy were published indicating an increased risk for breast cancer (risk ratio of 1.26). Using an adequate regression model in the previously published meta-analysis an adjusted estimate of effect of 1.14 can be given which is considerably higher than the one published in the meta-analysis of Sillero-Arenas et al. In summary, it is hoped that the method suggested here contributes further to a good meta-analytic practice in public health and clinical disciplines.

A practical approach to a multi-level analysis with a sparse binary outcome within a large surgical trial

Objectives: To assess the potential source of variation that surgeon may add to patient outcome in a clinical trial of surgical procedures. Methods: Two large (n = 1380) parallel multicentre randomized surgical trials were undertaken to compare laparoscopically assisted hysterectomy with conventional methods of abdominal and vaginal hysterectomy; involving 43 surgeons. The primary end point of the trial was the occurrence of at least one major complication. Patients were nested within surgeons giving the data set a hierarchical structure. A total of 10% of patients had at least one major complication, that is, a sparse binary outcome variable. A linear mixed logistic regression model (with logit link function) was used to model the probability of a major complication, with surgeon fitted as a random effect. Models were fitted using the method of maximum likelihood in SAS((R)). Results: There were many convergence problems. These were resolved using a variety of approaches including; treating all effects as fixed for the initial model building; modelling the variance of a parameter on a logarithmic scale and centring of continuous covariates. The initial model building process indicated no significant 'type of operation' across surgeon interaction effect in either trial, the 'type of operation' term was highly significant in the abdominal trial, and the 'surgeon' term was not significant in either trial. Conclusions: The analysis did not find a surgeon effect but it is difficult to conclude that there was not a difference between surgeons. The statistical test may have lacked sufficient power, the variance estimates were small with large standard errors, indicating that the precision of the variance estimates may be questionable.

Design and analysis of climate model experiments for the efficient estimation of anthropogenic signals

Presented herein is an experimental design that allows the effects of several radiative forcing factors on climate to be estimated as precisely as possible from a limited suite of atmosphere-only general circulation model (GCM) integrations. The forcings include the combined effect of observed changes in sea surface temperatures, sea ice extent, stratospheric (volcanic) aerosols, and solar output, plus the individual effects of several anthropogenic forcings. A single linear statistical model is used to estimate the forcing effects, each of which is represented by its global mean radiative forcing. The strong colinearity in time between the various anthropogenic forcings provides a technical problem that is overcome through the design of the experiment. This design uses every combination of anthropogenic forcing rather than having a few highly replicated ensembles, which is more commonly used in climate studies. Not only is this design highly efficient for a given number of integrations, but it also allows the estimation of (nonadditive) interactions between pairs of anthropogenic forcings. The simulated land surface air temperature changes since 1871 have been analyzed. The changes in natural and oceanic forcing, which itself contains some forcing from anthropogenic and natural influences, have the most influence. For the global mean, increasing greenhouse gases and the indirect aerosol effect had the largest anthropogenic effects. It was also found that an interaction between these two anthropogenic effects in the atmosphere-only GCM exists. This interaction is similar in magnitude to the individual effects of changing tropospheric and stratospheric ozone concentrations or to the direct (sulfate) aerosol effect. Various diagnostics are used to evaluate the fit of the statistical model. For the global mean, this shows that the land temperature response is proportional to the global mean radiative forcing, reinforcing the use of radiative forcing as a measure of climate change. The diagnostic tests also show that the linear model was suitable for analyses of land surface air temperature at each GCM grid point. Therefore, the linear model provides precise estimates of the space time signals for all forcing factors under consideration. For simulated 50-hPa temperatures, results show that tropospheric ozone increases have contributed to stratospheric cooling over the twentieth century almost as much as changes in well-mixed greenhouse gases.

Structural and equilibrium studies on mixed ligand complexes of Co(II), Ni(II), Cu(II) and Zn(II) with N-(2-benzimidazolyl)methyliminodiacetic acid and typical N,N donor ligands

Mixed ligand complexes: [Co(L)(bipy)] (.) 3H(2)O (1), [Ni(L)(phen)] (.) H2O (2), [Cu(L)(phen)] (.) 3H(2)O (3) and [Zn(L)(bipy)] (.) 3H(2)O (4), where L2- = two -COOH deprotonated dianion of N-(2-benzimidazolyl)methyliminodiacetic acid (H(2)bzimida, hereafter, H,L), bipy = 2,2' bipyridine and phen = 1,10-phenanthroline have been isolated and characterized by elemental analysis, spectral and magnetic measurements and thermal studies. Single crystal X-ray diffraction studies show octahedral geometry for 1, 2 and 4 and square pyramidal geometry for 3. Equilibrium studies in aqueous solution (ionic strength I = 10(-1) mol dm(-3) (NaNO3), at 25 +/- 1 degrees C) using different molar proportions of M(II):H2L:B, where M = Co, Ni, Cu and Zn and B = phen, bipy and en (ethylene diamine), however, provides evidence of formation of mononuclear and binuclear binary and mixed ligand complexes: M(L), M(H-1L)(-), M(B)(2+), M(L)(B), M(H-1L)(B)(-), M-2(H-1L)(OH), (B)M(H-1L)M(B)(+), where H-1L3- represents two -COOH and the benzimidazole NI-H deprotonated quadridentate (O-, N, O-, N), or, quinquedentate (O-, N, O-, N, N-) function of the coordinated ligand H,L. Binuclear mixed ligand complex formation equilibria: M(L)(B) + M(B)(2+) = (B)M(H-1L)M(B)(+) + H+ is favoured with higher pi-acidity of the B ligands. For Co(II), Ni(II) and Cu(II), these equilibria are accompanied by blue shift of the electronic absorption maxima of M(II) ions, as a negatively charged bridging benzimidazolate moiety provides stronger ligand field than a neutral one. Solution stability of the mixed ligand complexes are in the expected order: Co(II) < Ni(II) < Cu(II) > Zn(II). The Delta logK(M) values are less negetive than their statistical values, indicating favoured formation of the mixed ligand complexes over the binary ones. (c) 2005 Elsevier B.V. All rights reserved.

Construction of coordination polymers of cadmium(II) with mixed hexamethylenetetramine and terephthalate or thiocyanate ligands

Three new polymeric complexes [Cd(hmt)(SCN)(2)(H2O)(2)](n) (1), [Cd-3(mu(2)-hmt)(2)(SCN)(6)(H2O)(2)](n) (2), and [Cd-2(hmt)(2)(tP)(2)(H2O)(6)](n) (3) [hmt = hexamethylenetetramine, tp = terephthalate] have been synthesized and characterized by single crystal X-ray diffraction. Both the compounds 1 and 2 are 1-D polymers where Cd units are linked by double end-to-end thiocyanate bridges but in 2 the chain is wider containing three cadmium atoms instead of one as found in 1. In both compounds the coordination environment around cadmium atom is distorted octahedral. Compound 3 is a three-dimensional polymer having water-filled microporous channels. Both tp and brut are mu(2)-bridged. One of the acid groups of tp is coordinated in chelating bidentate and the other in monodentate fashion. Half of its Cd atoms are hexa-coordinated and the rest are hepta-coordinated. Thermogravimetric analysis and X-ray diffraction study of 3 show that its framework remains intact upon removal of water molecules. The flexibility of coordination number around cadmium atoms (six or seven) probably plays an important role in establishing the rigidity of the framework. (C) 2003 Elsevier Ltd. All rights reserved.

Chemical, physical and electrical properties of aged dodecylbenzene: thermal ageing of mixed isomers in air

A commercial dodecylbenzene (DDB) cable oil was aged at temperatures between 90 and 135 degrees C in air and was analyzed using various analytical techniques including optical and infra-red spectroscopy and dielectric analysis. On ageing, the oil darkened, significant oxidation features were found by infra-red spectroscopy and the acid number, water content and dielectric loss all increased. Ageing in the presence of paper or aluminum did not affect the ageing process, whereas ageing was significantly modified by the presence of copper. An absorption at 680 nm ("red absorbers") was detected by ultra-violet/visible spectroscopy followed by the production of an opaque precipitate. A reaction between copper and the acid generated on ageing is thought to produce copper carboxylates, and X-ray fluorescence confirmed that copper was indeed present in both the aged oil and the precipitate. Significantly, once red absorbers were detected, the dielectric loss increased to catastrophically high values and, therefore, the appearance of these compounds may serve as a useful diagnostic indicator. The development of acidity on ageing appears to be key in initiating the destructive copper conversion reaction and hence the control of oil acidity may be key to prolonging the life of DDB cable oils.

Chemistry of mixed-ligand methoxy bonded oxidovanadium(V) complexes with a family of hydrazone ligands containing VO3+ core and their substituent controlled methoxy-bridged dimeric forms

[(VO)-O-IV(acac)(2)] reacts with an equimolar amount of benzoyl hydrazones of 2-hydroxyacetophenone (H2L1), 2-hydroxy-5-methylacetophenone (H2L2) and 5-chloro-2-hydroxyacetophenone (H2L4) in methanol to afford the penta-coordinated mixed-ligand methoxy bonded oxidovanadium(V) complexes [(VO)-O-V(L-1)-(OCHA(3))](1). [(VO)-O-V(L-2)(OCH3)](2), and [(VO)-O-V(L-4)(OCH3)](4), respectively, whereas, the similar reaction with the benzoyl hydrazone of 2-hydroxy-5-methoxyacetophenone (H2L3) producing only the hexa-coordinated dimethoxy-bridged dimeric complex [(VO)-O-V(L-3)(OCH3)](2) (3A). Similar type of hexa-coordinated dimeric analogue of 1 i.e., [(VO)-O-V(L-1)(OCH3)](2) (1A) was obtained from the reaction of [(VO)-O-IV(acac)(2)] with the equimolar amount of H2L1 in presence of half equivalent 4,4'-bipyridine in methanol while the decomposition of [(VO)-O-IV(L-2)(bipy)] complex in methanol afforded the dimeric analogue of 2 i.e., [(VO)-O-V(L-2)(OCH3)](2) (2A). All these dimeric complexes 1A-3A react with an excess amount of imidazole in methanol producing the respective monomeric complex. The X-ray structural analysis of 1-3 and their dimeric analogues 1A-3A indicates that the geometry around the vanadium center in the monomeric form is distorted square-pyramidal while that of their respective dimeric forms is distorted octahedral, where the ligands are bonded to vanadium meridionally in their fully deprotonated enol forms. Due to the formation of bridge, the V-O(methoxy) bond in the dimeric complexes is lengthened to such an extent that it becomes equal in length with the V-O(phenolate) bond in 3A and even longer in 1A and 2A, which is unprecedented. The H-1 NMR spectra of the complexes 1A-3A in CDCl3 solution, indicates that these dimeric complexes are converted appreciably into their respective monomeric form. Complexes are electro-active displaying one quasi-reversible reduction peak near +0.25 V versus SCE in CH2Cl2 solution. The E-1/2 values of the complexes show linear relationship with the Hammett parameter (sigma) of the substituents. All these VO3+-complexes are converted to the corresponding complexes with V2O34+ motif simply on refluxing them in acetone and to the complexes with VO2+ motif on reaction with 2 KOH in methanol. (C) 2008 Elsevier Ltd. All rights reserved.

Analysis of red and yellow ochre samples from Clearwell Caves and Catalhoyuk by vibrational spectroscopy and other techniques

Ochre samples excavated from the neolithic site at Qatalhoyuk, Turkey have been compared with "native" ochres from Clearwell Caves, UK using infrared spectroscopy backed up by Raman spectroscopy, scanning electron microscopy (with energy-dispersive X-rays (EDX) analysis), powder X-ray diffraction, diffuse reflection UV-Vis and atomic absorption spectroscopies. For the Clearwell Caves ochres, which range in colour from yellow-orange to red-brown, it is shown that the colour is related to the nature of the chromophore present and not to any differences in particle size. The darker red ochres contain predominantly haematite while the yellow ochre contains only goethite. The ochres from Qatalhoyuk contain only about one-twentieth of the levels of iron found in the Clearwell Caves ochres. The iron oxide pigment (haematite in all cases studied here) has been mixed with a soft lime plaster which also contains calcite and silicate (clay) minerals. (C) 2003 Elsevier B.V. All rights reserved.

Validating a new clinical subtyping scheme for delirium with electronic motion analysis

The usefulness of motor subtypes of delirium is unclear due to inconsistency in sub-typing methods and a lack of validation with objective measures of activity. The activity of 40 patients was measured with 24 h accelerometry monitoring. Patients with Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) delirium (n = 30) were allocated into hyperactive, hypoactive and mixed motor subtypes. Delirium subtypes differed in relation to overall amount of activity, including movement in both sagittal and transverse planes. Differences were greater in the daytime and during the early evening ‘sundowning’ period. Frequency of postural changes was the most discriminating measure examined. Clinical subtypes of delirium defined by observed motor behaviour on the ward differ in electronically measured activity levels.

The structure of the mixed OH+H2O overlayer on Pt{111}

The structure of the mixed p(3x3)-(3OH+3H(2)O) phase on Pt{111} has been investigated by low-energy electron diffraction-IV structure analysis. The OH+H2O overlayer consists of hexagonal rings of coplanar oxygen atoms interlinked by hydrogen bonds. Lateral shifts of the O atoms away from atop sites result in different O-O separations and hexagons with only large separations (2.81 and 3.02 angstrom) linked by hexagons with alternating separations of 2.49 and 2.81/3.02 A. This unusual pattern is consistent with a hydrogen-bonded network in which water is adsorbed in cyclic rings separated by OH in a p(3x3) structure. The topmost two layers of the Pt atoms relax inwards with respect to the clean surface and both show vertical buckling of up to 0.06 angstrom. In addition, significant shifts away from the lateral bulk positions have been found for the second layer of Pt atoms. (C) 2005 American Institute of Physics.

A refined LEED analysis of water on Ru{0001}: an experimental test of the partial dissociation model

Despite a number of earlier studies which seemed to confirm molecular adsorption of water on close-packed surfaces of late transition metals, new controversy has arisen over a recent theoretical work by Feibelman, according to which partial dissociation occurs on the Ru{0001} surface leading to a mixed (H2O + OH + H) superstructure. Here, we present a refined LEED-IV analysis of the (root3 x root3)R30degrees-D2O-Ru{0001} structure, testing explicitly this new model by Feibelman. Our results favour the model proposed earlier by Held and Menzel assuming intact water molecules with almost coplanar oxygen atoms and out-of-plane hydrogen atoms atop the slightly higher oxygen atoms. The partially dissociated model with an almost identical arrangement of oxygen atoms can, however, not unambiguously be excluded, especially when the single hydrogen atoms are not present in the surface unit cell. In contrast to the earlier LEED-IV analysis, we can, however, clearly exclude a buckled geometry of oxygen atoms.

Geostatistical analysis of rainfall

Rainfall can be modeled as a spatially correlated random field superimposed on a background mean value; therefore, geostatistical methods are appropriate for the analysis of rain gauge data. Nevertheless, there are certain typical features of these data that must be taken into account to produce useful results, including the generally non-Gaussian mixed distribution, the inhomogeneity and low density of observations, and the temporal and spatial variability of spatial correlation patterns. Many studies show that rigorous geostatistical analysis performs better than other available interpolation techniques for rain gauge data. Important elements are the use of climatological variograms and the appropriate treatment of rainy and nonrainy areas. Benefits of geostatistical analysis for rainfall include ease of estimating areal averages, estimation of uncertainties, and the possibility of using secondary information (e.g., topography). Geostatistical analysis also facilitates the generation of ensembles of rainfall fields that are consistent with a given set of observations, allowing for a more realistic exploration of errors and their propagation in downstream models, such as those used for agricultural or hydrological forecasting. This article provides a review of geostatistical methods used for kriging, exemplified where appropriate by daily rain gauge data from Ethiopia.

In situ analysis and structural elucidation of sainfoin (Onobrychis viciifolia) tannins for high-throughput germplasm screening

A rapid thiolytic degradation and cleanup procedure was developed for analyzing tannins directly in chlorophyll-containing sainfoin (Onobrychis viciifolia) plants. The technique proved suitable for complex tannin mixtures containing catechin, epicatechin, gallocatechin, and epigallocatechin flavan-3-ol units. The reaction time was standardized at 60 min to minimize the loss of structural information as a result of epimerization and degradation of terminal flavan-3-ol units. The results were evaluated by separate analysis of extractable and unextractable tannins, which accounted for 63.6−113.7% of the in situ plant tannins. It is of note that 70% aqueous acetone extracted tannins with a lower mean degree of polymerization (mDP) than was found for tannins analyzed in situ. Extractable tannins had between 4 and 29 lower mDP values. The method was validated by comparing results from individual and mixed sample sets. The tannin composition of different sainfoin accessions covered a range of mDP values from 16 to 83, procyanidin/prodelphinidin (PC/PD) ratios from 19.2/80.8 to 45.6/54.4, and cis/trans ratios from 74.1/25.9 to 88.0/12.0. This is the first high-throughput screening method that is suitable for analyzing condensed tannin contents and structural composition directly in green plant tissue.

Impacts of pollution and climate change on ombrotrophic Sphagnum species in the UK: analysis of uncertainties in two empirical niche models

A significant challenge in the prediction of climate change impacts on ecosystems and biodiversity is quantifying the sources of uncertainty that emerge within and between different models. Statistical species niche models have grown in popularity, yet no single best technique has been identified reflecting differing performance in different situations. Our aim was to quantify uncertainties associated with the application of 2 complimentary modelling techniques. Generalised linear mixed models (GLMM) and generalised additive mixed models (GAMM) were used to model the realised niche of ombrotrophic Sphagnum species in British peatlands. These models were then used to predict changes in Sphagnum cover between 2020 and 2050 based on projections of climate change and atmospheric deposition of nitrogen and sulphur. Over 90% of the variation in the GLMM predictions was due to niche model parameter uncertainty, dropping to 14% for the GAMM. After having covaried out other factors, average variation in predicted values of Sphagnum cover across UK peatlands was the next largest source of variation (8% for the GLMM and 86% for the GAMM). The better performance of the GAMM needs to be weighed against its tendency to overfit the training data. While our niche models are only a first approximation, we used them to undertake a preliminary evaluation of the relative importance of climate change and nitrogen and sulphur deposition and the geographic locations of the largest expected changes in Sphagnum cover. Predicted changes in cover were all small (generally <1% in an average 4 m2 unit area) but also highly uncertain. Peatlands expected to be most affected by climate change in combination with atmospheric pollution were Dartmoor, Brecon Beacons and the western Lake District.

Drivers of fund performance: a panel data analysis

The principle aim of this research is to elucidate the factors driving the total rate of return of non-listed funds using a panel data analytical framework. In line with previous results, we find that core funds exhibit lower yet more stable returns than value-added and, in particular, opportunistic funds, both cross-sectionally and over time. After taking into account overall market exposure, as measured by weighted market returns, the excess returns of value-added and opportunity funds are likely to stem from: high leverage, high exposure to development, active asset management and investment in specialized property sectors. A random effects estimation of the panel data model largely confirms the findings obtained from the fixed effects model. Again, the country and sector property effect shows the strongest significance in explaining total returns. The stock market variable is negative which hints at switching effects between competing asset classes. For opportunity funds, on average, the returns attributable to gearing are three times higher than those for value added funds and over five times higher than for core funds. Overall, there is relatively strong evidence indicating that country and sector allocation, style, gearing and fund size combinations impact on the performance of unlisted real estate funds.