Bayesian analysis of an inverse Gaussian correlated frailty model

In survival analysis frailty is often used to model heterogeneity between individuals or correlation within clusters. Typically frailty is taken to be a continuous random effect, yielding a continuous mixture distribution for survival times. A Bayesian analysis of a correlated frailty model is discussed in the context of inverse Gaussian frailty. An MCMC approach is adopted and the deviance information criterion is used to compare models. As an illustration of the approach a bivariate data set of corneal graft survival times is analysed. (C) 2006 Elsevier B.V. All rights reserved.

A novel method for measuring lag times in division of individual bacterial cells using image analysis

A method is presented for determining the time to first division of individual bacterial cells growing on agar media. Bacteria were inoculated onto agar-coated slides and viewed by phase-contrast microscopy. Digital images of the growing bacteria were captured at intervals and the time to first division estimated by calculating the "box area ratio". This is the area of the smallest rectangle that can be drawn around an object, divided by the area of the object itself. The box area ratios of cells were found to increase suddenly during growth at a time that correlated with cell division as estimated by visual inspection of the digital images. This was caused by a change in the orientation of the two daughter cells that occurred when sufficient flexibility arose at their point of attachment. This method was used successfully to generate lag time distributions for populations of Escherichia coli, Listeria monocytogenes and Pseudomonas aeruginosa, but did not work with the coccoid organism Staphylococcus aureus. This method provides an objective measure of the time to first cell division, whilst automation of the data processing allows a large number of cells to be examined per experiment. (c) 2005 Elsevier B.V. All rights reserved.

The effect of phase-correlated returns and spatial smoothing on the accuracy of radar refractivity retrievals

Radar refractivity retrievals have the potential to accurately capture near-surface humidity fields from the phase change of ground clutter returns. In practice, phase changes are very noisy and the required smoothing will diminish large radial phase change gradients, leading to severe underestimates of large refractivity changes (ΔN). To mitigate this, the mean refractivity change over the field (ΔNfield) must be subtracted prior to smoothing. However, both observations and simulations indicate that highly correlated returns (e.g., when single targets straddle neighboring gates) result in underestimates of ΔNfield when pulse-pair processing is used. This may contribute to reported differences of up to 30 N units between surface observations and retrievals. This effect can be avoided if ΔNfield is estimated using a linear least squares fit to azimuthally averaged phase changes. Nevertheless, subsequent smoothing of the phase changes will still tend to diminish the all-important spatial perturbations in retrieved refractivity relative to ΔNfield; an iterative estimation approach may be required. The uncertainty in the target location within the range gate leads to additional phase noise proportional to ΔN, pulse length, and radar frequency. The use of short pulse lengths is recommended, not only to reduce this noise but to increase both the maximum detectable refractivity change and the number of suitable targets. Retrievals of refractivity fields must allow for large ΔN relative to an earlier reference field. This should be achievable for short pulses at S band, but phase noise due to target motion may prevent this at C band, while at X band even the retrieval of ΔN over shorter periods may at times be impossible.

Correlated observation errors in data assimilation

Data assimilation provides techniques for combining observations and prior model forecasts to create initial conditions for numerical weather prediction (NWP). The relative weighting assigned to each observation in the analysis is determined by its associated error. Remote sensing data usually has correlated errors, but the correlations are typically ignored in NWP. Here, we describe three approaches to the treatment of observation error correlations. For an idealized data set, the information content under each simplified assumption is compared with that under correct correlation specification. Treating the errors as uncorrelated results in a significant loss of information. However, retention of an approximated correlation gives clear benefits.

Relative humidity changes in a warmer climate

Key climate feedbacks due to water vapor and clouds rest largely on how relative humidity R changes in a warmer climate, yet this has not been extensively analyzed in models. General circulation models (GCMs) from the CMIP3 archive and several higher resolution atmospheric GCMs examined here generally predict a characteristic pattern of R trend with global temperature that has been reported previously in individual models, including increase around the tropopause, decrease in the tropical upper troposphere, and decrease in midlatitudes. This pattern is very similar to that previously reported for cloud cover in the same GCMs, confirming the role of R in controlling changes in simulated cloud. Comparing different models, the trend in each part of the troposphere is approximately proportional to the upward and/or poleward gradient of R in the present climate. While this suggests that the changes simply reflect a shift of the R pattern upward with the tropopause and poleward with the zonal jets, the drying trend in the subtropics is roughly three times too large to be attributable to shifts of subtropical features, and the subtropical R minima deepen in most models. R trends are correlated with horizontal model resolution, especially outside the tropics, where they show signs of convergence and latitudinal gradients become close to available observations for GCM resolutions near T85 and higher. We argue that much of the systematic change in R can be explained by the local specific humidity having been set (by condensation) in remote regions with different temperature changes, hence the gradients and trends each depend on a model’s ability to resolve moisture transport. Finally, subtropical drying trends predicted from the warming alone fall well short of those observed in recent decades. While this discrepancy supports previous reports of GCMs underestimating Hadley Cell expansion, our results imply that shifts alone are not a sufficient interpretation of changes.

Chlorophyll-a in the rivers of eastern England

Chlorophyll-a concentration variations are described for two major river basins in England, the Humber and the Thames and related to catchment characteristics and nutrient concentrations across a range of rural, agricultural and urban/industrial settings. For all the rivers there are strong seasonal variations, with concentrations peaking in the spring and summer time when biological activity is at its highest. However, there are large variations in the magnitude of the seasonal effects across the rivers. For the spring-summer low-flow periods, average concentrations of chlorophyll-a correlate with soluble reactive phosphor-us (SRP). Chlorophyll-a is also correlated with particulate nitrogen (PN), organic carbon (POC) and suspended sediments. However, the strongest relationships are with catchment area and flow, where two straight line relationships are observed. The results indicate the importance of residence times for determining planktonic growth within the rivers. This is also indicated by the lack of chlorophyll-a response to lowering of SRP concentrations in several of the rivers in the area due to phosphorus stripping of effluents at major sewage treatment works. A key control on chlorophyll-a concentration may be the input of canal and reservoir waters during the growing period: this too relates to issues of residence times. However, there may well be a complex series of factors influencing residence time across the catchments due to features such as inhomogeneous flow within the catchments, a fractal distribution of stream channels that leads to a distribution of residence times and differences in planktonic inoculation sources. Industrial pollution on the Aire and Calder seems to have affected the relationship of chlorophyll-a with PN and POC. The results are discussed in relation to the Water Framework Directive. (c) 2006 Elsevier B.V. All rights reserved.

Changing times, changing places: urban development and the politics of space - time

Much of the writing on urban regeneration in the UK has been focused on the types of urban spaces that have been created in city centres. Less has been written about the issue of when the benefits of regeneration could and should be delivered to a range of different interests, and the different time frames that exist in any development area. Different perceptions of time have been reflected in dominant development philosophies in the UK and elsewhere. The trickle-down agendas of the 1980s, for example, were criticised for their focus on the short-term time frames and needs of developers, often at the expense of those of local communities. The recent emergence of sustainability discourses, however, ostensibly changes the time focus of development and promotes a broader concern with new imagined futures. This paper draws on the example of development in Salford Quays, in the North West of England, to argue that more attention needs to be given to the politics of space-time in urban development processes. It begins by discussing the importance and relevance of this approach before turning to the case study and the ways in which the local politics of space-time has influenced development agendas and outcomes. The paper argues that such an approach harbours the potential for more progressive, far-reaching, and sustainable development agendas to be developed and implemented.

Remote estimation of thermal inertia and soil heat flux for bare soil

A method is presented which allows thermal inertia (the soil heat capacity times the square root of the soil thermal diffusivity, C(h)rootD(h)), to be estimated remotely from micrometeorological observations. The method uses the drop in surface temperature, T-s, between sunset and sunrise, and the average night-time net radiation during that period, for clear, still nights. A Fourier series analysis was applied to analyse the time series of T-s . The Fourier series constants, together with the remote estimate of thermal inertia, were used in an analytical expression to calculate diurnal estimates of the soil heat flux, G. These remote estimates of C(h)rootD(h) and G compared well with values derived from in situ sensors. The remote and in situ estimates of C(h)rootD(h) both correlated well with topsoil moisture content. This method potentially allows area-average estimates of thermal inertia and soil heat flux to be derived from remote sensing, e.g. METEOSAT Second Generation, where the area is determined by the sensor's height and viewing angle. (C) 2003 Elsevier B.V. All rights reserved.

Radiative forcing since preindustrial times due to ozone change in the troposphere and the lower stratosphere

Changes in atmospheric ozone have occurred since the preindustrial era as a result of increasing anthropogenic emissions. Within ACCENT, a European Network of Excellence, ozone changes between 1850 and 2000 are assessed for the troposphere and the lower stratosphere ( up to 30 km) by a variety of seven chemistry-climate models and three chemical transport models. The modeled ozone changes are taken as input for detailed calculations of radiative forcing. When only changes in chemistry are considered ( constant climate) the modeled global-mean tropospheric ozone column increase since preindustrial times ranges from 7.9 DU to 13.8 DU among the ten participating models, while the stratospheric column reduction lies between 14.1 DU and 28.6 DU in the models considering stratospheric chemistry. The resulting radiative forcing is strongly dependent on the location and altitude of the modeled ozone change and varies between 0.25 Wm(-2) and 0.45 Wm(-2) due to ozone change in the troposphere and - 0.123 Wm(-2) and + 0.066 Wm(-2) due to the stratospheric ozone change. Changes in ozone and other greenhouse gases since preindustrial times have altered climate. Six out of the ten participating models have performed an additional calculation taking into account both chemical and climate change. In most models the isolated effect of climate change is an enhancement of the tropospheric ozone column increase, while the stratospheric reduction becomes slightly less severe. In the three climate-chemistry models with detailed tropospheric and stratospheric chemistry the inclusion of climate change increases the resulting radiative forcing due to tropospheric ozone change by up to 0.10 Wm(-2), while the radiative forcing due to stratospheric ozone change is reduced by up to 0.034 Wm(-2). Considering tropospheric and stratospheric change combined, the total ozone column change is negative while the resulting net radiative forcing is positive.

Ambient solar wind's effect on ICME transit times

Most empirical and numerical models of Interplanetary Coronal Mass Ejection (ICME) propagation use the initial CME velocity as their primary, if not only, observational input. These models generally predict a wide spread of 1 AU transit times for ICMEs with the same initial velocity. We use a 3D coupled MHD model of the corona and heliosphere to determine the ambient solar wind's effect on the propagation of ICMEs from 30 solar radii to 1 AU. We quantitatively characterize this deceleration by the velocity of the upstream ambient solar wind. The effects of varying solar wind parameters on the ICME transit time are quantified and can explain the observed spread in transit times for ICMEs of the same initial velocity. We develop an adjustment formula that can be used in conjunction with other models to reduce the spread in predicted transit times of Earth-directed ICMEs.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

A parametric model for the interval censored survival times of Acacia mangium plantation in a spacing trial

Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.