Monte Carlo methods for matrix computations on the grid

Many scientific and engineering applications involve inverting large matrices or solving systems of linear algebraic equations. Solving these problems with proven algorithms for direct methods can take very long to compute, as they depend on the size of the matrix. The computational complexity of the stochastic Monte Carlo methods depends only on the number of chains and the length of those chains. The computing power needed by inherently parallel Monte Carlo methods can be satisfied very efficiently by distributed computing technologies such as Grid computing. In this paper we show how a load balanced Monte Carlo method for computing the inverse of a dense matrix can be constructed, show how the method can be implemented on the Grid, and demonstrate how efficiently the method scales on multiple processors. (C) 2007 Elsevier B.V. All rights reserved.

Agonist-dependent internalization of D2 receptors: imaging quantification by confocal microscopy

In positron emission tomography and single photon emission computed tomography studies using D2 dopamine (DA) receptor radiotracers, a decrease in radiotracer binding potential (BP) is usually interpreted in terms of increased competition with synaptic DA. However, some data suggest that this signal may also reflect agonist (DA)-induced increases in D2 receptor (D2R) internalization, a process which would presumably also decrease the population of receptors available for binding to hydrophilic radioligands. To advance interpretation of alterations in D2 radiotracer BP, direct methods of assessment of D2R internalization are required. Here, we describe a confocal microscopy-based approach for the quantification of agonist-dependent receptor internalization. The method relies upon double-labeling of the receptors with antibodies directed against intracellular as well as extracellular epitopes. Following agonist stimulation, DA D2R internalization was quantified by differentiating, in optical cell sections, the signal due to the staining of the extracellular from intracellular epitopes of D2Rs. Receptor internalization was increased in the presence of the D2 agonists DA and bromocriptine, but not the D1 agonist SKF38393. Pretreatment with either the D2 antagonist sulpiride, or inhibitors of internalization (phenylarsine oxide and high molarity sucrose), blocked D2-agonist induced receptor internalization, thus validating this method in vitro. This approach therefore provides a direct and streamlined methodology for investigating the pharmacological and mechanistic aspects of D2R internalization, and should inform the interpretation of results from in vivo receptor imaging studies.

Dung and nest surveys: estimating decay rates

1. Wildlife managers often require estimates of abundance. Direct methods of estimation are often impractical, especially in closed-forest environments, so indirect methods such as dung or nest surveys are increasingly popular. 2. Dung and nest surveys typically have three elements: surveys to estimate abundance of the dung or nests; experiments to estimate the production (defecation or nest construction) rate; and experiments to estimate the decay or disappearance rate. The last of these is usually the most problematic, and was the subject of this study. 3. The design of experiments to allow robust estimation of mean time to decay was addressed. In most studies to date, dung or nests have been monitored until they disappear. Instead, we advocate that fresh dung or nests are located, with a single follow-up visit to establish whether the dung or nest is still present or has decayed. 4. Logistic regression was used to estimate probability of decay as a function of time, and possibly of other covariates. Mean time to decay was estimated from this function. 5. Synthesis and applications. Effective management of mammal populations usually requires reliable abundance estimates. The difficulty in estimating abundance of mammals in forest environments has increasingly led to the use of indirect survey methods, in which abundance of sign, usually dung (e.g. deer, antelope and elephants) or nests (e.g. apes), is estimated. Given estimated rates of sign production and decay, sign abundance estimates can be converted to estimates of animal abundance. Decay rates typically vary according to season, weather, habitat, diet and many other factors, making reliable estimation of mean time to decay of signs present at the time of the survey problematic. We emphasize the need for retrospective rather than prospective rates, propose a strategy for survey design, and provide analysis methods for estimating retrospective rates.

Agonist-dependent internalization of D2 receptors: Imaging quantification by confocal microscopy

In positron emission tomography and single photon emission computed tomography studies using D2 dopamine (DA) receptor radiotracers, a decrease in radiotracer binding potential (BP) is usually interpreted in terms of increased competition with synaptic DA. However, some data suggest that this signal may also reflect agonist (DA)-induced increases in D2 receptor (D2R) internalization, a process which would presumably also decrease the population of receptors available for binding to hydrophilic radioligands. To advance interpretation of alterations in D2 radiotracer BP, direct methods of assessment of D2R internalization are required. Here, we describe a confocal microscopy-based approach for the quantification of agonist-dependent receptor internalization. The method relies upon double-labeling of the receptors with antibodies directed against intracellular as well as extracellular epitopes. Following agonist stimulation, DA D2R internalization was quantified by differentiating, in optical cell sections, the signal due to the staining of the extracellular from intracellular epitopes of D2Rs. Receptor internalization was increased in the presence of the D2 agonists DA and bromocriptine, but not the D1 agonist SKF38393. Pretreatment with either the D2 antagonist sulpiride, or inhibitors of internalization (phenylarsine oxide and high molarity sucrose), blocked D2-agonist induced receptor internalization, thus validating this method in vitro. This approach therefore provides a direct and streamlined methodology for investigating the pharmacological and mechanistic aspects of D2R internalization, and should inform the interpretation of results from in vivo receptor imaging studies.

A differential thermal expansion approach to crystal structure determination from powder diffraction data

Differential thermal expansion over the range 90-210 K has been applied successfully to determine the crystal structure of chlorothiazide from synchrotron powder diffraction data using direct methods. Key to the success of the approach is the use of a multi-data-set Pawley refinement to extract a set of reflection intensities that is more 'single-crystal-like' than those extracted from a single data set. The improvement in reflection intensity estimates is quantified by comparison with reference single-crystal intensities. (C) 2008 International Union of Crystallography Printed in Singapore - all rights reserved

On computational algorithms for pole assignment

A number of computationally reliable direct methods for pole assignment by feedback have recently been developed. These direct procedures do not necessarily produce robust solutions to the problem, however, in the sense that the assigned poles are insensitive to perturbalions in the closed-loop system. This difficulty is illustrated here with results from a recent algorithm presented in this TRANSACTIONS and its causes are examined. A measure of robustness is described, and techniques for testing and improving robustness are indicated.

Approximate Gauss-Newton methods for nonlinear least squares problems

The Gauss–Newton algorithm is an iterative method regularly used for solving nonlinear least squares problems. It is particularly well suited to the treatment of very large scale variational data assimilation problems that arise in atmosphere and ocean forecasting. The procedure consists of a sequence of linear least squares approximations to the nonlinear problem, each of which is solved by an “inner” direct or iterative process. In comparison with Newton’s method and its variants, the algorithm is attractive because it does not require the evaluation of second-order derivatives in the Hessian of the objective function. In practice the exact Gauss–Newton method is too expensive to apply operationally in meteorological forecasting, and various approximations are made in order to reduce computational costs and to solve the problems in real time. Here we investigate the effects on the convergence of the Gauss–Newton method of two types of approximation used commonly in data assimilation. First, we examine “truncated” Gauss–Newton methods where the inner linear least squares problem is not solved exactly, and second, we examine “perturbed” Gauss–Newton methods where the true linearized inner problem is approximated by a simplified, or perturbed, linear least squares problem. We give conditions ensuring that the truncated and perturbed Gauss–Newton methods converge and also derive rates of convergence for the iterations. The results are illustrated by a simple numerical example. A practical application to the problem of data assimilation in a typical meteorological system is presented.

Calibrating flood inundation models: methods for assessing model performance

Satellite observed data for flood events have been used to calibrate and validate flood inundation models, providing valuable information on the spatial extent of the flood. Improvements in the resolution of this satellite imagery have enabled indirect remote sensing of water levels by using an underlying LiDAR DEM to extract the water surface elevation at the flood margin. Further to comparison of the spatial extent, this now allows for direct comparison between modelled and observed water surface elevations. Using a 12.5m ERS-1 image of a flood event in 2006 on the River Dee, North Wales, UK, both of these data types are extracted and each assessed for their value in the calibration of flood inundation models. A LiDAR guided snake algorithm is used to extract an outline of the flood from the satellite image. From the extracted outline a binary grid of wet / dry cells is created at the same resolution as the model, using this the spatial extent of the modelled and observed flood can be compared using a measure of fit between the two binary patterns of flooding. Water heights are extracted using points at intervals of approximately 100m along the extracted outline, and the students T-test is used to compare modelled and observed water surface elevations. A LISFLOOD-FP model of the catchment is set up using LiDAR topographic data resampled to the 12.5m resolution of the satellite image, and calibration of the friction parameter in the model is undertaken using each of the two approaches. Comparison between the two approaches highlights the sensitivity of the spatial measure of fit to uncertainty in the observed data and the potential drawbacks of using the spatial extent when parts of the flood are contained by the topography.

The 'direct costs' of livestock disease: The development of a system of models for the analysis of 30 endemic livestock diseases in Great Britain

The 'direct costs' attributable to 30 different endemic diseases of farm animals in Great Britain are estimated using a standardised method to construct a simple model for each disease that includes consideration of disease prevention and treatment costs. The models so far developed provide a basis for further analyses including cost-benefit analyses for the economic assessment of disease control options. The approach used reflects the inherent livestock disease information constraints, which limit the application of other economic analytical methods. It is a practical and transparent approach that is relatively easily communicated to veterinary scientists and policy makers. The next step is to develop the approach by incorporating wider economic considerations into the analyses in a way that will demonstrate to policy makers and others the importance of an economic perspective to livestock disease issues.

Participatory Farm Management methods for assessing the suitability of potential innovations. A case study on green manuring options for tomato producers in Ghana

This paper examines the potential of using Participatory Farm Management methods to examine the suitability of a technology with farmers prior to on-farm trials. A study examining the suitability of green manuring as a technology for use with wet season tomato producers in Ghana is described. Findings from this case-study demonstrate that Participatory Budgeting can be used by farmers and researchers to analyse current cultivation practices, identify the options for including green manures into the system and explore the direct and wider resource implications of the technology. Scored-Causal Diagrams can be used to identify farmers' perceptions of the relative importance of the problem that the technology seeks to address. The use of the methods in this examine evaluation process appears to have the potential to improve the effectiveness and efficiency of the adaptive research process. This ensures that technologies subsequently examined in trials ate relevant to farmers' interests, existing systems and resources, thereby increasing the chances of farmer adoption. It is concluded that this process has potential for use-with other technologies and in other farming systems. (C) 2002 Elsevier Science Ltd. All rights reserved.

The development of indigo reduction methods and pre-reduced indigo products

The continuing importance of blue denim maintains indigo as an important vat dye industrially. In this review, we examine the various methods that have been used in the past and are currently used to reduce and dissolve indigo for dyeing. We discuss recent insights into the bacterial fermentation technology, the advantages and disadvantages of the direct chemical methods that have predominated for the last century and potentially cleaner technologies of catalytic hydrogenation and electrochemistry, which are becoming increasingly important. With considerations of environmental impact high on the dyeing industry's agenda, we also discuss the developments that have led to the production of pre-reduced indigo.

Direct detection of dimethylstannylene and tetramethyldistannene in solution and the gas phase by laser flash photolysis of 1,1-dimethylstannacyclopent-3-enes

The photochemistry of 1,1-dimethyl- and 1,1,3,4-tetramethylstannacyclopent-3-ene (4a and 4b,respectively) has been studied in the gas phase and in hexane solution by steady-state and 193-nm laser flash photolysis methods. Photolysis of the two compounds results in the formation of 1,3-butadiene (from 4a) and 2,3-dimethyl-1,3-butadiene (from 4b) as the major products, suggesting that cycloreversion to yield dimethylstannylene (SnMe2) is the main photodecomposition pathway of these molecules. Indeed, the stannylene has been trapped as the Sn-H insertion product upon photolysis of 4a in hexane containing trimethylstannane. Flash photolysis of 4a in the gas phase affords a transient absorbing in the 450-520nm range that is assigned to SnMe2 by comparison of its spectrum and reactivity to those previously reported from other precursors. Flash photolysis of 4b in hexane solution affords results consistent with the initial formation of SnMe2 (lambda(max) approximate to 500 nm), which decays over similar to 10 mu s to form tetramethyldistannene (5b; lambda(max) approximate to 470 nm). The distannene decays over the next ca. 50 mu s to form at least two other longer-lived species, which are assigned to higher SnMe2 oligomers. Time-dependent DFT calculations support the spectral assignments for SnMe2 and Sn2Me4, and calculations examining the variation in bond dissociation energy with substituent (H, Me, and Ph) in disilenes, digermenes, and distannenes rule out the possibility that dimerization of SnMe2 proceeds reversibly. Addition of methanol leads to reversible reaction with SnMe2 to form a transient absorbing at lambda(max) approximate to 360 nm, which is assigned to the Lewis acid-base complex between SnMe2 and the alcohol.

Fecal water as a non-invasive biomarker in nutritional intervention: comparison of preparation methods and refinement of different endpoints

The assessment of cellular effects by the aqueous phase of human feces (fecal water, FW) is a useful biomarker approach to study cancer risks and protective activities of food. In order to refine and develop the biomarker, different protocols of preparing FW were compared. Fecal waters were prepared by 3 methods: (A) direct centrifugation; (B) extraction of feces in PBS before centrifugation; and (C) centrifugation of lyophilized and reconstituted feces. Genotoxicity was determined in colon cells using the Comet assay. Selected samples were investigated for additional parameters related to carcinogenesis. Two of 7 FWs obtained by methods A and B were similarly genotoxic. Method B, however, yielded higher volumes of FW, allowing sterile filtration for long-term culture experiments. Four of 7 samples were non-genotoxic when prepared according to all 3 methods. FW from lyophilized feces and from fresh samples were equally genotoxic. FWs modulated cytotoxicity, paracellular permeability, and invasion, independent of their genotoxicity. All 3 methods of FW preparation can be used to assess genotoxicity. The higher volumes of FWobtained by preparation method B greatly enhance the perspectives of measuring different types of biological parameters and using these to disclose activities related to cancer development.

Lighting energy saving due to the use of direct sunlight

This paper describe a simulation program, which uses Trengenza’s average room illuminance method in conjunction with hourly solar irradiance and luminous efficacy, to predict the potential lighting energy saving for a side-lit room. Two lighting control algorithms of photoelectric switching (on/off) and photoelectric dimming (top-up) have been coded in the program. A simulation for a typical UK office room has been conducted and the results show that energy saving due to the sunlight dependent on the various factors such as orientation, control methods, building depth, glazing area and shading types, etc. This simple tool can be used for estimating the potential lighting energy saving of the windows with various shading devices at the early design stage.

What have dendritic cells ever done for adjuvant design? Cellular and molecular methods for the rational development of vaccine adjuvants

Our new molecular understanding of immune priming states that dendritic cell activation is absolutely pivotal for expansion and differentiation of naïve T lymphocytes, and it follows that understanding DC activation is essential to understand and design vaccine adjuvants. This chapter describes how dendritic cells can be used as a core tool to provide detailed quantitative and predictive immunomics information about how adjuvants function. The role of distinct antigen, costimulation, and differentiation signals from activated DC in priming is explained. Four categories of input signals which control DC activation – direct pathogen detection, sensing of injury or cell death, indirect activation via endogenous proinflammatory mediators, and feedback from activated T cells – are compared and contrasted. Practical methods for studying adjuvants using DC are summarised and the importance of DC subset choice, simulating T cell feedback, and use of knockout cells is highlighted. Finally, five case studies are examined that illustrate the benefit of DC activation analysis for understanding vaccine adjuvant function.