Modeling data with structural and temporal correlation using lower level and higher level multilevel models

Novel imaging techniques are playing an increasingly important role in drug development, providing insight into the mechanism of action of new chemical entities. The data sets obtained by these methods can be large with complex inter-relationships, but the most appropriate statistical analysis for handling this data is often uncertain - precisely because of the exploratory nature of the way the data are collected. We present an example from a clinical trial using magnetic resonance imaging to assess changes in atherosclerotic plaques following treatment with a tool compound with established clinical benefit. We compared two specific approaches to handle the correlations due to physical location and repeated measurements: two-level and four-level multilevel models. The two methods identified similar structural variables, but higher level multilevel models had the advantage of explaining a greater proportion of variation, and the modeling assumptions appeared to be better satisfied.

Development of a semantic graph-based scene description

Software representations of scenes, i.e. the modelling of objects in space, are used in many application domains. Current modelling and scene description standards focus on visualisation dimensions, and are intrinsically limited by their dependence upon their semantic interpretation and contextual application by humans. In this paper we propose the need for an open, extensible and semantically rich modelling language, which facilitates a machine-readable semantic structure. We critically review existing standards and techniques, and highlight a need for a semantically focussed scene description language. Based on this defined need we propose a preliminary solution, based on hypergraph theory, and reflect on application domains.

RAMI4PILPS: An intercomparison of formulations for the partitioning of solar radiation in land surface models

[1] Remotely sensed, multiannual data sets of shortwave radiative surface fluxes are now available for assimilation into land surface schemes (LSSs) of climate and/or numerical weather prediction models. The RAMI4PILPS suite of virtual experiments assesses the accuracy and consistency of the radiative transfer formulations that provide the magnitudes of absorbed, reflected, and transmitted shortwave radiative fluxes in LSSs. RAMI4PILPS evaluates models under perfectly controlled experimental conditions in order to eliminate uncertainties arising from an incomplete or erroneous knowledge of the structural, spectral and illumination related canopy characteristics typical for model comparison with in situ observations. More specifically, the shortwave radiation is separated into a visible and near-infrared spectral region, and the quality of the simulated radiative fluxes is evaluated by direct comparison with a 3-D Monte Carlo reference model identified during the third phase of the Radiation transfer Model Intercomparison (RAMI) exercise. The RAMI4PILPS setup thus allows to focus in particular on the numerical accuracy of shortwave radiative transfer formulations and to pinpoint to areas where future model improvements should concentrate. The impact of increasing degrees of structural and spectral subgrid variability on the simulated fluxes is documented and the relevance of any thus emerging biases with respect to gross primary production estimates and shortwave radiative forcings due to snow and fire events are investigated.

Genetic variability in partitioning to the yield component of cacao (Theobroma cacao L.)

Biomass partitioning of cacao (Theobroma cacao L.) was studied in seven clones and five hybrids in a replicated experiment in Bahia, Brazil. Over an eighteen month period, a seven- fold difference in dry bean yield was demonstrated between genotypes, ranging from the equivalent of 200 to 1389 kg.ha-1. During the same interval, the increase in trunk cross-sectional area ranged from 11.1 cm2 for clone EEG-29 to 27.6 cm2 for hybrid PA-150 * MA-15. Yield efficiency increment (the ratio of cumulative yield to the increase in trunk circumference), which indicated partitioning between the vegetative and reproductive components, ranged from 0.008 kg.cm-2 for clone CP-82 to 0.08 kg.cm-2 for clone EEG-29. An examination of biomass partitioning within the pod of the seven clones revealed that the beans accounted for between 32.0% (CP-82) and 44.5% (ICS-9) of the pod biomass. The study demonstrated the potential for yield improvement in cacao by selectively breeding for more efficient partitioning to the yield component.

Complexation of lanthanides, actinides and transition metal cations with a 6-(1,2,4-triazin-3-yl)-2,2′:6′,2′′-terpyridine ligand: implications for actinide(III)/lanthanide(III) partitioning

The quadridentate N-heterocyclic ligand 6-(5,5,8,8-tetramethyl-5,6,7,8-tetrahydro-1,2,4-benzotriazin- 3-yl)-2,2′ : 6′,2′′-terpyridine (CyMe4-hemi-BTBP) has been synthesized and its interactions with Am(III),U(VI), Ln(III) and some transition metal cations have been evaluated by X-ray crystallographic analysis, Am(III)/Eu(III) solvent extraction experiments, UVabsorption spectrophotometry, NMR studies and ESI-MS. Structures of 1 : 1 complexes with Eu(III), Ce(III) and the linear uranyl (UO2 2+) ion were obtained by X-ray crystallographic analysis, and they showed similar coordination behavior to related BTBP complexes. In methanol, the stability constants of the Ln(III) complexes are slightly lower than those of the analogous quadridentate bis-triazine BTBP ligands, while the stability constant for the Yb(III)complex is higher. 1H NMR titrations and ESI-MS with lanthanide nitrates showed that the ligand forms only 1 : 1 complexes with Eu(III), Ce(III) and Yb(III), while both 1 : 1 and 1 : 2 complexes were formed with La(III) and Y(III) in acetonitrile. A mixture of isomeric chiral 2 : 2 helical complexes was formed with Cu(I), with a slight preference (1.4 : 1) for a single directional isomer. In contrast, a 1 : 1 complex was observed with the larger Ag(I) ion. The ligand was unable to extract Am(III) or Eu(III) from nitric acid solutions into 1-octanol, except in the presence of a synergist at low acidity. The results show that the presence of two outer 1,2,4-triazine rings is required for the efficient extraction and separation of An(III)from Ln(III) by quadridentate N-donor ligands.

Association mapping of partitioning loci in barley

Background: Association mapping, initially developed in human disease genetics, is now being applied to plant species. The model species Arabidopsis provided some of the first examples of association mapping in plants, identifying previously cloned flowering time genes, despite high population sub-structure. More recently, association genetics has been applied to barley, where breeding activity has resulted in a high degree of population sub-structure. A major genotypic division within barley is that between winter- and spring-sown varieties, which differ in their requirement for vernalization to promote subsequent flowering. To date, all attempts to validate association genetics in barley by identifying major flowering time loci that control vernalization requirement (VRN-H1 and VRN-H2) have failed. Here, we validate the use of association genetics in barley by identifying VRN-H1 and VRN-H2, despite their prominent role in determining population sub-structure. Results: By taking barley as a typical inbreeding crop, and seasonal growth habit as a major partitioning phenotype, we develop an association mapping approach which successfully identifies VRN-H1 and VRN-H2, the underlying loci largely responsible for this agronomic division. We find a combination of Structured Association followed by Genomic Control to correct for population structure and inflation of the test statistic, resolved significant associations only with VRN-H1 and the VRN-H2 candidate genes, as well as two genes closely linked to VRN-H1 (HvCSFs1 and HvPHYC). Conclusion: We show that, after employing appropriate statistical methods to correct for population sub-structure, the genome-wide partitioning effect of allelic status at VRN-H1 and VRN-H2 does not result in the high levels of spurious association expected to occur in highly structured samples. Furthermore, we demonstrate that both VRN-H1 and the candidate VRN-H2 genes can be identified using association mapping. Discrimination between intragenic VRN-H1 markers was achieved, indicating that candidate causative polymorphisms may be discerned and prioritised within a larger set of positive associations. This proof of concept study demonstrates the feasibility of association mapping in barley, even within highly structured populations. A major advantage of this method is that it does not require large numbers of genome-wide markers, and is therefore suitable for fine mapping and candidate gene evaluation, especially in species for which large numbers of genetic markers are either unavailable or too costly.

BTBPs versus BTPhens: some reasons for their differences in properties concerning the partitioning of minor actinides and the advantages of BTPhens

Two members of the tetradentate N-donor ligand families 6,6′-bis(1,2,4-triazin-3-yl)-2,2′-bipyridine (BTBP) and 2,9-bis(1,2,4-triazin-3-yl)-1,10-phenanthroline (BTPhen) currently being developed for separating actinides from lanthanides have been studied. It has been confirmed that CyMe4-BTPhen 2 has faster complexation kinetics than CyMe4-BTBP 1. The values for the HOMO−LUMO gap of 2 are comparable with those of CyMe4-BTBP 1 for which the HOMO−LUMO gap was previously calculated to be 2.13 eV. The displacement of BTBP from its bis-lanthanum(III) complex by BTPhen was observed by NMR, and constitutes the only direct evidence for the greater thermodynamic stability of the complexes of BTPhen. NMR competition experiments suggest the following order of bis-complex stability: 1:2 bis-BTPhen complex ≥ heteroleptic BTBP/BTPhen 1:2 bis-complex > 1:2 bis-BTBP complex. Kinetics studies on some bis-triazine N-donor ligands using the stopped-flow technique showed a clear relationship between the rates of metal ion complexation and the degree to which the ligand is preorganized for metal binding. The BTBPs must overcome a significant (ca. 12 kcal mol−1) energy barrier to rotation about the central biaryl C−C axis in order to achieve the cis−cis conformation that is required to form a complex, whereas the cis−cis conformation is fixed in the BTPhens. Complexation thermodynamics and kinetics studies in acetonitrile show subtle differences between the thermodynamic stabilities of the complexes formed, with similar stability constants being found for both ligands. The first crystal structure of a 1:1 complex of CyMe4-BTPhen 2 with Y(NO3)3 is also reported. The metal ion is 10- coordinate being bonded to the tetradentate ligand 2 and three bidentate nitrate ions. The tetradentate ligand is nearly planar with angles between consecutive rings of 16.4(2)°, 6.4(2)°, 9.7(2)°, respectively.

Graph-based algorithms for comparison and prediction of household-level energy use profiles

We present an efficient graph-based algorithm for quantifying the similarity of household-level energy use profiles, using a notion of similarity that allows for small time–shifts when comparing profiles. Experimental results on a real smart meter data set demonstrate that in cases of practical interest our technique is far faster than the existing method for computing the same similarity measure. Having a fast algorithm for measuring profile similarity improves the efficiency of tasks such as clustering of customers and cross-validation of forecasting methods using historical data. Furthermore, we apply a generalisation of our algorithm to produce substantially better household-level energy use forecasts from historical smart meter data.

Characterization of energy flux partitioning in urban environments: links with surface seasonal properties

A better understanding of links between the properties of the urban environment and the exchange to the atmosphere is central to a wide range of applications. The numerous measurements of surface energy balance data in urban areas enable intercomparison of observed fluxes from distinct environments. This study analyzes a large database in two new ways. First, instead of normalizing fluxes using net all-wave radiation only the incoming radiative fluxes are used, to remove the surface attributes from the denominator. Second, because data are now available year-round, indices are developed to characterize the fraction of the surface (built; vegetation) actively engaged in energy exchanges. These account for shading patterns within city streets and seasonal changes in vegetation phenology; their impact on the partitioning of the incoming radiation is analyzed. Data from 19 sites in North America, Europe, Africa, and Asia (including 6-yr-long observation campaigns) are used to derive generalized surface–flux relations. The midday-period outgoing radiative fraction decreases with an increasing total active surface index, the stored energy fraction increases with an active built index, and the latent heat fraction increases with an active vegetated index. Parameterizations of these energy exchange ratios as a function of the surface indices [i.e., the Flux Ratio–Active Index Surface Exchange (FRAISE) scheme] are developed. These are used to define four urban zones that characterize energy partitioning on the basis of their active surface indices. An independent evaluation of FRAISE, using three additional sites from the Basel Urban Boundary Layer Experiment (BUBBLE), yields accurate predictions of the midday flux partitioning at each location.

Utilizing electronic effects in the modulation of BTPhen ligands with respect to the partitioning of minor actinides from lanthanides

Effects of bromine substitution at the 5 and 5,6-positions of the 1,10-phenanthroline nucleus of BTPhen ligand on their extraction properties for Ln(III) andAn(III) cations have been studied. Compared to C5-BTPhen, electronic modulation in BrC5-BTPhen and Br2C5-BTPhen enabled these ligands to be fine-tuned in order to enhance the separation selectivity of Am(III) from Eu(III)

Mg/Ca partitioning between aqueous solution and aragonite mineral: a molecular dynamics study

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO3 in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

Improved evaporative flux partitioning and carbon flux in the land surface model JULES: impact on the simulation of land surface processes in temperate Europe

The primary role of land surface models embedded in climate models is to partition surface available energy into upwards, radiative, sensible and latent heat fluxes. Partitioning of evapotranspiration, ET, is of fundamental importance: as a major component of the total surface latent heat flux, ET affects the simulated surface water balance, and related energy balance, and consequently the feedbacks with the atmosphere. In this context it is also crucial to credibly represent the CO2 exchange between ecosystems and their environment. In this study, JULES, the land surface model used in UK weather and climate models, has been evaluated for temperate Europe. Compared to eddy covariance flux measurements, the CO2 uptake by the ecosystem is underestimated and the ET overestimated. In addition, the contribution to ET from soil and intercepted water evaporation far outweighs the contribution of plant transpiration. To alleviate these biases, adaptations have been implemented in JULES, based on key literature references. These adaptations have improved the simulation of the spatio-temporal variability of the fluxes and the accuracy of the simulated GPP and ET, including its partitioning. This resulted in a shift of the seasonal soil moisture cycle. These adaptations are expected to increase the fidelity of climate simulations over Europe. Finally, the extreme summer of 2003 was used as evaluation benchmark for the use of the model in climate change studies. The improved model captures the impact of the 2003 drought on the carbon assimilation and the water use efficiency of the plants. It, however, underestimates the 2003 GPP anomalies. The simulations showed that a reduction of evaporation from the interception and soil reservoirs, albeit not of transpiration, largely explained the good correlation between the carbon and the water fluxes anomalies that was observed during 2003. This demonstrates the importance of being able to discriminate the response of individual component of the ET flux to environmental forcing.