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.

Estimation of body fatness from body mass index and bioelectrical impedance: comparison of New Zealand European, Maori and Pacific Island children

Objective: To compare percentage body fat (%BF) for a given body mass index (BMI) among New Zealand European, Maori and Pacific Island children. To develop prediction equations based on bioimpedance measurements for the estimation of fat-free mass (FFM) appropriate to children in these three ethnic groups. Design: Cross-sectional study. Purposive sampling of schoolchildren aimed at recruiting three children of each sex and ethnicity for each year of age. Double cross-validation of FFM prediction equations developed by multiple regression. Setting: Local schools in Auckland. Subjects: Healthy European, Maori and Pacific Island children (n = 172, 83 M, 89 F, mean age 9.4 +/- 2.8(s. d.), range 5 - 14 y). Measurements: Height, weight, age, sex and ethnicity were recorded. FFM was derived from measurements of total body water by deuterium dilution and resistance and reactance were measured by bioimpedance analysis. Results: For fixed BMI, the Maori and Pacific Island girls averaged 3.7% lower % BF than European girls. For boys a similar relation was not found since BMI did not significantly influence % BF of European boys ( P = 0.18). Based on bioimpedance measurements a single prediction equation was developed for all children: FFM (kg) = 0.622 height (cm)(2)/ resistance +0.234 weight (kg)+1.166, R-2 = 0.96, s. e. e. = 2.44 kg. Ethnicity, age and sex were not significant predictors. Conclusions: A robust equation for estimation of FFM in New Zealand European, Maori and Pacific Island children in the 5 - 14 y age range that is more suitable than BMI for the determination of body fatness in field studies has been developed. Sponsorship: Maurice and Phyllis Paykel Trust, Auckland University of Technology Contestable Grants Fund and the Ministry of Health.

Assessment of bioelectrical impedance analysis for the prediction of total body water in cystic fibrosis

The aim of this study was to compare the measurement of total body water (TBW) by deuterium ((H2O)-H-2) dilution and bioelectrical impedance analysis (BIA) in patients with cystic fibrosis (CF) and healthy controls. Thirty-six clinically stable patients with CF (age 25.4 +/- 5.6 yrs) and 42 healthy controls (age 25.4 +/- 4.8) were recruited into this study. TBW was measured by (H2O)-H-2 dilution and predicted by BIA in patients and controls. The TBW predicted from BIA was significantly different from TBW as measured using (H2O)-H-2 in patients (P

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.

Measured and predicted total body water in children with myelomeningocele

Objective: Children with myelomeningocele (MMC) have an altered body composition and an atypical distribution of total body water (TBW). The aim of the present study was to determine the accuracy of current predictive equations, based on bioelectrical impedance analysis (BIA), in determining TBW when compared with measured TBW using deuterium dilution. Methods: Fourteen children with MMC were measured for whole body BIA and TBW (using deuterium dilution and the Plateau method). Total body water was predicted using equations based on the resistance and characteristic frequency from BIA measurements and heights of subjects. Results: The mean measured TBW was 15.46 +/- 8.28 L and the mean predictions for TBW using equations based on the resistance and characteristic frequency from BIA measurements and heights of subjects were 18.29 +/- 8.41 L, 17.72 +/- 11.42 L and 12.51 +/- 7.59 L, respectively. The best correlation was found using characteristic frequency. The limits of agreement between measured and predicted TBW values using Bland-Altman analysis were large. Conclusions: The present study suggests that the prediction of TBW in children with MMC can be made accurately using the equation of Cornish et al . based on BIA measurements of characteristic frequency.

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.

Prediction of mortality in patients without angina: Use of an exercise score and exercise echocardiography

Background Exercise testing has limited efficacy for identifying coronary artery disease (CAD) in the absence of anginal. symptoms. Exercise echocardiography is more accurate than standard exercise testing, but its efficacy in this situation has not been defined. We sought to identify whether the Duke treadmill. score or exercise echocardiography (ExE) could be used to identify risk in patients without anginal symptoms. Methods We studied 1859 patients without typical or atypical angina, heart failure, or a history or ECG evidence of infarction or CAD, who were referred for ExE, of whom 1832 (age 51 15 years, 944 men) were followed for up to 10 years. The presence and extent of ischaemia and scar were interpreted by expert reviewers at the time of the original study. Results Exercise provoked significant (>0.1 mV) ST segment depression in 215 patients (12%), and wall motion abnormalities in 137 (8%). Seventy-eight patients (4%) died before revascularization, only 17 from known cardiac causes. The independent predictors of death were age (RR 1.1, p

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

Solubility of selected esters in supercritical carbon dioxide

The solubility of ethyl propionate, ethyl butyrate, and ethyl isovalerate in supercritical carbon dioxide was measured at temperature ranging from 308.15 to 333.15 K and pressure ranging from 85 to 195 bar. At the same temperature, the solubility of these compounds increases with pressure. The crossover pressure region was also observed in this study. The experimental data were correlated by the semi-empirical Chrastil equation and Peng-Robinson equation of state (EOS) using several mixing rules. The Peng-Robinson EOS gives better solubility prediction than the empirical Chrastil equation. (C) 2002 Elsevier Science B.V. 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.

Prediction of multilayer adsorption and capillary condensation phenomena in cylindrical mesopores

MCM-41 periodic mesoporous silicates with a high degree of structural ordering are synthesized and used as model adsorbents to study the isotherm prediction of nitrogen adsorption. The nitrogen adsorption isotherm at 77 K for a macroporous silica is measured and used in high-resolution alpha(s)-plot comparative analysis to determine the external surface area, total surface area and primary mesopore volume of the MCM-41 materials. Adsorption equilibrium data of nitrogen on the different pore size MCM-41 samples (pore diameters from 2.40 to 4.92 nm) are also obtained. Based on the Broekhoff and de Boer' thermodynamic analysis, the nitrogen adsorption isotherms for the different pore size MCM-41 samples are interpreted using a novel strategy, in which the parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting only the multilayer region prior to capillary condensation for C-16 MCM-41. Subsequently the entire isotherm, including the phase transition, is predicted for all the different pore size MCM-41 samples without any fitting. The results show that the prediction of multilayer adsorption and total adsorbed amount are in good agreement with the experimental isotherms. The predictions of the relative pressure corresponding to capillary equilibrium (coexistence) transition agree remarkably with experimental data on the adsorption branch even for hysteretic isotherms, confirming that this is the branch appropriate for pore size distribution analysis. The impact of pore radius on the adsorption film thickness and capillary coexistence pressure is also investigated, and found to agree with the experimental data. (C) 2003 Elsevier Inc. All rights reserved.

Fatigue Prediction on Fibrin Poly-ε-caprolactone Macroporous Scaffolds

Tissue engineering applications rely on scaffolds that during its service life, either for in-vivo or in vitro applications, are under mechanical solicitations. The variation of the mechanical condition of the scaffold is strongly relevant for cell culture and has been scarcely addressed. Fatigue life cycle of poly-ε-caprolactone, PCL, scaffolds with and without fibrin as filler of the pore structure were characterized both dry and immersed in liquid water. It is observed that the there is a strong increase from 100 to 500 in the number of loading cycles before collapse in the samples tested in immersed conditions due to the more uniform stress distributions within the samples, the fibrin loading playing a minor role in the mechanical performance of the scaffolds