Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

The effect of heat transfer mode on heart rate responses and hysteresis during heating and cooling in the estuarine crocodile Crocodylus porosus

The effect of heating and cooling on heart rate in the estuarine crocodile Crocodylus porosus was studied in response to different heat transfer mechanisms and heat loads. Three heating treatments were investigated. C. porosus were: (1) exposed to a radiant heat source under dry conditions; (2) heated via radiant energy while half-submerged in flowing water at 23degreesC and (3) heated via convective transfer by increasing water temperature from 23degreesC to 35degreesC. Cooling was achieved in all treatments by removing the heat source and with C. porosus half-submerged in flowing water at 23degreesC. In all treatments, the heart rate of C. porosus increased markedly in response to heating and decreased rapidly with the removal of the heat source. Heart rate during heating was significantly faster than during cooling at any given body temperature, i.e. there was a significant heart rate hysteresis. There were two identifiable responses to heating and cooling. During the initial stages of applying or removing the heat source, there was a dramatic increase or decrease in heart rate ('rapid response'), respectively, indicating a possible cardiac reflex. This rapid change in heart rate with only a small change or no change in body temperature (

Heat transfer and flow behind a step in high enthalpy superorbital flow

Heat transfer levels have been investigated behind a rearward-facing step in a superorbital expansion tube. The heat transfer was measured along a flat plate and behind 2 and 3mm steps with the same length to step height ratio. Results were obtained with air as the test gas at speeds of 6.76kms(-1) and 9-60kms(-1) corresponding to stagnation enthalpies of 26MJ/kg and 48MJ/kg respectively. A laminar boundary layer was established on the flat plate and measured heat transfer levels were consistent with classical empirical correlations. In the case of flow behind a step, the measurements showed a gradual rise in heat transfer from the rear of the step to a plateau several step heights downstream for both flow conditions. Reattachment distance was estimated to be approximately 1.6 step heights downstream of the 2mm step at the low enthalpy condition through the use of flow visualisation.

Influence of ageing, heat shock treatment and in vivo total antioxidant status on gene-expression profile and protein synthesis in human peripheral lymphocytes

Ageing results in a progressive, intrinsic and generalised imbalance of the control of regulatory systems. A key manifestation of this complex biological process includes the attenuation of the universal stress response. Here we provide the first global assessment of the ageing process as it affects the heat shock response, utilising human peripheral lymphocytes and cDNA microarray analysis. The genomic approach employed in our preliminary study was supplemented with a proteomic approach. In addition, the current study correlates the in vivo total antioxidant status with the age-related differential gene expression as well as the translational kinetics of heat shock proteins (hsps). Most of the genes encoding stress response proteins on the 4224 element microarray used in this study were significantly elevated after heat shock treatment of lymphocytes obtained from both young and old individuals albeit to a greater extent in the young. Cell signaling and signal transduction genes as well as some oxidoreductases showed varied response. Results from translational kinetics of induction of major hsps, from 0 to 24 It recovery period were broadly consistent with the differential expression of HSC 70 and HSP 40 genes. Total antioxidant levels in plasma from old individuals were found to be significantly lower by comparison with young, in agreement with the widely acknowledged role of oxidant homeostasis in the ageing process. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Telehealth responses to bio-terrorism and emerging infections

Emerging infectious diseases, such as severe acute respiratory syndrome (SARS), are of huge economic importance. They are difficult to predict. The World Health Organization has a Global Outbreak Alert and Response Network, which was involved at an early stage in the SARS outbreak in 2003. Three major lessons were learned as a result of the SARS epidemic in 2003, involving communication, evidence-based action and global partnerships. It is proposed that a series of broadband global response networks should be developed. At a technical level the networks are essentially in place, such as the Internet2 global network. Suitable peripheral devices also exist. What has not yet been created is the appropriate software to allow the use of these networks, although a number of commercial products are in the process of development.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Link invariants associated with gauge equivalent solutions of the Yang-Baxter equation: The one-parameter family of minimal typical representations of Uq[gl(2|1)]

In this paper we investigate the construction of state models for link invariants using representations of the braid group obtained from various gauge choices for a solution of the trigonometric Yang-Baxter equation. Our results show that it is possible to obtain invariants of regular isotopy (as defined by Kauffman) which may not be ambient isotopic. We illustrate our results with explicit computations using solutions of the trigonometric Yang-Baxter equation associated with the one-parameter family of minimal typical representations of the quantum superalgebra U-q,[gl(2/1)]. We have implemented MATHEMATICA code to evaluate the invariants for all prime knots up to 10 crossings.

Heat resistance of Bacillus spores when adhered to stainless steel and its relationship to spore hydrophobicity

Twenty-one strains of Bacillus (10 B. stearothermophilus, 3 B. cereus, and 8 B. licheniformis strains) were assayed for spore surface hydrophobicity on the basis of three measures: contact angle measurement (CAM), microbial adhesion to hydrocarbons (MATH), and hydrophobic interaction chromatography (HIC). On the basis of the spore surface characteristics obtained from these assays, along with data on the heat resistance of these spores in water, eight strains of Bacillus (three B. stearothermophilus, three B. cereus, and two B. licheniformis strains) either suspended in water or adhering to stainless steel were exposed to sublethal heat treatments at 90 to 110degreesC to determine heat resistance (D-value). Significant increases in heat resistance (ranging from 3 to 400%) were observed for the eight strains adhering to stainless steel. No significant correlation was found between these heat resistance increases and spore surface characteristics as determined by the three hydrophobicity assays. There was a significant positive correlation between the hydrophobicity data obtained by the MATH assay and those obtained by the HIC assay, but these data did not correlate with those obtained by the CAM assay.

Heat-induced changes in UHT milks - Part 1

Raw milk was stored for 0, 2 and 4 days and processed in a UHT pilot plant by either direct or indirect heating. The unstored raw milk was also pasteurised. The thermally induced changes resulting from these treatments were investigated by examining a number of indices of heat damage. Lactulose, furosine, total and free hydroxymethylfurfural (HMF) and acid-soluble beta-lactoglobulin were analysed by high performance liquid chromatography (HPLC) while soluble tryptophan was examined by fluorescence spectroscopy. The directly heated UHT milk showed less heat damage than the indirectly heated milk, while the pasteurised milk displayed the least heat damage. During storage of the UHT milk for 12 weeks at similar to20degreesC, the levels of lactulose remained constant, while the furosine concentration increased. Both the total HMF and undenatured beta-lactoglobulin contents showed a general decrease during storage; however free HMF values initially rose but then decreased after four weeks' storage. As the age of the milk at the time of UHT processing increased, the levels of some of the indicators decreased. It is concluded that lactulose is the most reliable index of heat treatment, as it is virtually unaffected by refrigerated storage of the milk before or ambient storage after UHT processing. Reliance on other indicators may give misleading information on the heat load that UHT milk has received during processing.

Effects of heat treatment on the nutritional value of raw soybean selected for low trypsin inhibitor activity

1. Two broiler experiments and a layer experiments were conducted on Kunitz trypsin inhibitor (Kti) soybeans (SB) of low trypsin inhibitor (TI) activity to determine their nutritive value when included as mash in least-cost poultry diets. 2. Experiment 1 compared chick performance on the Kti or raw SB using a commercial full-fat SB meal (FFSBM) and a solvent extracted SB meal (SBM) as controls during a 20 d experimental period. Broiler experiment 2 compared Kti and raw SB, non-steamed, or steam-pelleted with and without DL-methionine supplementation added to every treatment containing 170 g SB/kg. For each broiler experiment the levels of each SB were 70, 120 and 170 g/kg with the control birds fed only 170 g SB/kg. 3. The layer experiment, compared steam-pelleted Kti and raw SB against a non-steamed Kti and raw SB each fed at two levels (70 and 110 g/kg) x 30 replicates from 29 weeks of age for 19 weeks in a completely randomised design. Production parameters were measured when diets were formulated to contain minimum required specifications and calculated apparent metabolisable energy (AME). At the completion of each trial, 2 broiler birds from each cage and 5 layer birds per treatment were killed, weighed, and their liver and pancreas weighed. 4. Both broiler experiments indicated that production parameters on the Kti SB treatments were significantly lower (P < 0.05) than on the two commercial control SB treatments. However, the Kti treatments were superior to the raw SB treatments. 5. Pancreas weight increased with increasing inclusion of both raw and Kti SB, suggesting that a TI was causing the depression in performance. The AME of the Kti SB was similar to that of commercial FFSB meal. After steam conditioning, the raw SB meal AME value of 9.5 MJ/kg dry matter (DM) was improved to 14.1 MJ/kg DM by reduced TI activity, but this AME improvement with TI activity reduction, plus the supplementation with DL-methionine on birds fed the raw SB had no effect (P > 0.05) on any parameter evaluated in experiment 2. 6. The layer experiment showed that hens on the Kti SB treatments had significantly greater live weight gain (LWG), egg weight and daily egg mass than birds given raw SB. A reduced food intake (FI) was observed in the Kti treatments but egg mass was generally similar to that on the FFSB control diet, indicating that Kti SB supported excellent egg production at an inclusion of 110 g/kg. The depressed performance observed for broiler chicks suggest that younger birds are more susceptible to the effects of SB TI.

Dryout phenomena in a three-phase fixed-bed reactor

Understanding the mechanism of liquid-phase evaporation in a three-phase fixed-bed reactor is of practical importance, because the reaction heat is usually 7-10 times the vaporization heat of the liquid components. Evaporation, especially the liquid dryout, can largely influence the reactor performance and even safety. To predict the vanishing condition of the liquid phase, Raoult's law was applied as a preliminary approach, with the liquid vanishing temperature defined based on a liquid flow rate of zero. While providing correct trends, Raoult's law exhibits some limitation in explaining the temperature profile in the reactor. To comprehensively understand the whole process of liquid evaporation, a set of experiments on inlet temperature, catalyst activity, liquid flow rate, gas flow rate, and operation pressure were carried out. A liquid-region length-predicting equation is suggested based on these experiments and the principle of heat balance.

Comparação de desempenho entre a formulação direta do Método dos Elementos de Contorno com Funções Radiais e o Método dos Elementos Finitos em problemas de Poisson e Helmholtz

O presente trabalho objetiva avaliar o desempenho do MECID (Método dos Elementos de Contorno com Interpolação Direta) para resolver o termo integral referente à inércia na Equação de Helmholtz e, deste modo, permitir a modelagem do Problema de Autovalor assim como calcular as frequências naturais, comparando-o com os resultados obtidos pelo MEF (Método dos Elementos Finitos), gerado pela Formulação Clássica de Galerkin. Em primeira instância, serão abordados alguns problemas governados pela equação de Poisson, possibilitando iniciar a comparação de desempenho entre os métodos numéricos aqui abordados. Os problemas resolvidos se aplicam em diferentes e importantes áreas da engenharia, como na transmissão de calor, no eletromagnetismo e em problemas elásticos particulares. Em termos numéricos, sabe-se das dificuldades existentes na aproximação precisa de distribuições mais complexas de cargas, fontes ou sorvedouros no interior do domínio para qualquer técnica de contorno. No entanto, este trabalho mostra que, apesar de tais dificuldades, o desempenho do Método dos Elementos de Contorno é superior, tanto no cálculo da variável básica, quanto na sua derivada. Para tanto, são resolvidos problemas bidimensionais referentes a membranas elásticas, esforços em barras devido ao peso próprio e problemas de determinação de frequências naturais em problemas acústicos em domínios fechados, dentre outros apresentados, utilizando malhas com diferentes graus de refinamento, além de elementos lineares com funções de bases radiais para o MECID e funções base de interpolação polinomial de grau (um) para o MEF. São geradas curvas de desempenho através do cálculo do erro médio percentual para cada malha, demonstrando a convergência e a precisão de cada método. Os resultados também são comparados com as soluções analíticas, quando disponíveis, para cada exemplo resolvido neste trabalho.