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 rheology of flotation froths

A literature review has highlighted the need to measure flotation froth rheology in order to fully characterise the role of the froth in the flotation process. The initial investigation using a coaxial cylinder viscometer for froth rheology measurement led to the development of a new device employing a vane measuring head. The modified rheometer was used in industrial scale flotation tests at Mt. Isa Copper Concentrator. The measured froth rheograms show a non-Newtonian nature for the flotation froths (pseudoplastic flow). The evidence of the non-Newtonian flow has questioned the validity of application of the Laplace equation in froth motion modelling as used by a number of researchers, since the assumption of irrotational flow is violated. Correlations between the froth rheology and the froth retention time, water hold-up in the froth and concentrate grades have been found. These correlations are independent of air flow rate (test data at various air flow rates fall on one similar trend line). This implies that froth rheology may be used as a lumped parameter for other operating variables in flotation modelling and scale up. (C) 2003 Elsevier Science B.V. All rights reserved.

A modified Weibull distribution

A new lifetime distribution capable of modeling a bathtub-shaped hazard-rate function is proposed. The proposed model is derived as a limiting case of the Beta Integrated Model and has both the Weibull distribution and Type I extreme value distribution as special cases. The model can be considered as another useful 3-parameter generalization of the Weibull distribution. An advantage of the model is that the model parameters can be estimated easily based on a Weibull probability paper (WPP) plot that serves as a tool for model identification. Model characterization based on the WPP plot is studied. A numerical example is provided and comparison with another Weibull extension, the exponentiated Weibull, is also discussed. The proposed model compares well with other competing models to fit data that exhibits a bathtub-shaped hazard-rate function.

The blood volumes of the primary and secondary circulatory system in the Atlantic cod Gadus morhua L., using plasma bound Evans Blue and compartmental analysis

The volume of the primary (PCS) and secondary (SCS) circulatory system in the Atlantic cod Gadus morhua was determined using a modified dye dilution technique. Cod (N=10) were chronically cannulated in the second afferent branchial artery with PE-50 tubing. Evans Blue dye was bound to harvested fish plasma at a concentration of 1 mg dye ml(-1) plasma, and injected at a concentration of 1 mg kg(-1) body mass. Serial sampling from the cannula produced a dye dilution curve, which could be described by a double exponential decay equation. Curve analysis enabled the calculation of the primary circulatory and total distribution volume. The difference between these volumes is assumed to be the volume of the SCS. From the dilution curve, it was also possible to calculate flow rates between and within the systems. The results of these experiments suggest a plasma volume in the PCS of 3.42+/-0.89 ml 100 g(-1) body mass, and in the SCS of 1.68+/-0.35 ml 100 g(-1) body mass (mean +/- S.D.) or approximately 50% that of the PCS. Flow rates to the SCS were calculated as 2.7% of the resting cardiac output. There was an allometric relationship between body mass and blood volumes. Increasing condition factor showed a tendency towards smaller blood volumes of the PCS, expressed as percentage body mass, but this was not evident for the volume of the SCS.

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.

Purification of post-translationally modified proteins from bacteria: homologous expression and purification of histidine-tagged pilin from Neisseria meningitidis

Until recently, glycosylation of proteins in prokaryotes was regarded as uncommon and thought to be limited to special cases such as S-layer proteins and some archeal outer membrane proteins. Now, there are an increasing number of reports of bacterial proteins that are glycosylated. Pilin of pathogenic Neisseria is one of the best characterised post-translation ally modified bacterial proteins, with four different types of modifications reported, including a novel glycosylation. Pilin monomers assemble to form pilus fibres, which are long protein filaments that protrude from the surface of bacterial cells and are key virulence factors. To aid in the investigation of these modifications, pure pilin is required. A number of pilin purification methods have been published, but none are appropriate for the routine purification of pilin from many different isolates. This study describes a novel, rapid, and simple method of pilin purification from Neisseria meningitidis C311#3, which facilitates the production of consistent quantities of pure, native pilin. A 6 x histidine tag was fused to the C-terminus of the pilin subunit structural gene, pilE, via homologous recombination placing the 6 x histidine-tagged allele in the chromosome of N. meningitidis C311#3. Pilin was purified under non-denaturing conditions via a two-step process using immobilised metal affinity chromatography (IMAC), followed by dye affinity chromatography. Analysis of the purified pilin confirmed that it retained both of the post-translational modifications examined. This novel approach may prove to be a generally applicable method for purification and analysis of post-translationally modified proteins in bacteria. (C) 2003 Elsevier Science (USA). All rights reserved.

Effect of pore-network connectivity on multicomponent adsorption of large molecules

The effect of pore-network connectivity on binary liquid-phase adsorption equilibria using the ideal adsorbed solution theory (LAST) was studied. The liquid-phase binary adsorption experiments used ethyl propionate, ethyl butyrate, and ethyl isovalerate as the adsorbates and commercial activated carbons Filtrasorb-400 and Norit ROW 0.8 as adsorbents. As the single-component isotherm, a modified Dubinin-Radushkevich equation was used. A comparison with experimental data shows that incorporating the connectivity of the pore network and considering percolation processes associated with different molecular sizes of the adsorptives in the mixture, as well as their different corresponding accessibility, can improve the prediction of binary adsorption equilibria using the LAST Selectivity of adsorption for the larger molecule in binary systems increases with an increase in the pore-network coordination number, as well with an increase in the mean pore width and in the spread of the pore-size distribution.

Modified logit life table system: Principles, empirical validation, and application

Despite its widespread use, the Coale-Demeny model life table system does not capture the extensive variation in age-specific mortality patterns observed in contemporary populations, particularly those of the countries of Eastern Europe and populations affected by HIV/AIDS. Although relational mortality models, such as the Brass logit system, can identify these variations, these models show systematic bias in their predictive ability as mortality levels depart from the standard. We propose a modification of the two-parameter Brass relational model. The modified model incorporates two additional age-specific correction factors (gamma(x), and theta(x)) based on mortality levels among children and adults, relative to the standard. Tests of predictive validity show deviations in age-specific mortality rates predicted by the proposed system to be 30-50 per cent lower than those predicted by the Coale-Demeny system and 15-40 per cent lower than those predicted using the original Brass system. The modified logit system is a two-parameter system, parameterized using values of l(5) and l(60).