Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

A Computerized Method for the Hydrologic Design of Culverts, February 1974

Nationwide, about five cents of each highway construction dollar is spent on culverts. In Iowa, average annual construction costs on the interstate, primary, and federal-aid secondary systems are about $120,000,000. Assuming the national figure applies to Iowa, about $6,000,000 are spent on culvert construction annually. For each one percent reduction in overall culvert costs, annual construction costs would be reduced by $60,000. One area of potential cost reduction lies in the sizing of the culvert. Determining the flow area and hydraulic capacity is accomplished in the initial design of the culvert. The normal design sequence is accomplished in two parts. The hydrologic portion consists of the determination of a design discharge in cubic feet per second using one of several available methods. This discharge is then used directly in the hydraulic portion of the design to determine the proper type, size, and shape of culvert to be used, based on various site and design restrictions. More refined hydrologic analyses, including rainfall-runoff analysis, flood hydrograph development, and streamflow routing techniques, are not pursued in the existing design procedure used by most county and state highway engineers.

Method to estimate soil macroporosity and microporosity based on sand content and bulk density

Macroporosity is often used in the determination of soil compaction. Reduced macroporosity can lead to poor drainage, low root aeration and soil degradation. The aim of this study was to develop and test different models to estimate macro and microporosity efficiently, using multiple regression. Ten soils were selected within a large range of textures: sand (Sa) 0.07-0.84; silt 0.03-0.24; clay 0.13-0.78 kg kg-1 and subjected to three compaction levels (three bulk densities, BD). Two models with similar accuracy were selected, with a mean error of about 0.02 m³ m-3 (2 %). The model y = a + b.BD + c.Sa, named model 2, was selected for its simplicity to estimate Macro (Ma), Micro (Mi) or total porosity (TP): Ma = 0.693 - 0.465 BD + 0.212 Sa; Mi = 0.337 + 0.120 BD - 0.294 Sa; TP = 1.030 - 0.345 BD 0.082 Sa; porosity values were expressed in m³ m-3; BD in kg dm-3; and Sa in kg kg-1. The model was tested with 76 datum set of several other authors. An error of about 0.04 m³ m-3 (4 %) was observed. Simulations of variations in BD as a function of Sa are presented for Ma = 0 and Ma = 0.10 (10 %). The macroporosity equation was remodeled to obtain other compaction indexes: a) to simulate maximum bulk density (MBD) as a function of Sa (Equation 11), in agreement with literature data; b) to simulate relative bulk density (RBD) as a function of BD and Sa (Equation 13); c) another model to simulate RBD as a function of Ma and Sa (Equation 16), confirming the independence of this variable in relation to Sa for a fixed value of macroporosity and, also, proving the hypothesis of Hakansson & Lipiec that RBD = 0.87 corresponds approximately to 10 % macroporosity (Ma = 0.10 m³ m-3).

Hysteresis in a system driven by either generalized force or displacement variables: martensitic phase transition in single-crystalline Cu-Zn-Al

We report on experiments aimed at comparing the hysteretic response of a Cu-Zn-Al single crystal undergoing a martensitic transition under strain-driven and stress-driven conditions. Strain-driven experiments were performed using a conventional tensile machine while a special device was designed to perform stress-driven experiments. Significant differences in the hysteresis loops were found. The strain-driven curves show reentrant behavior yield point which is not observed in the stress-driven case. The dissipated energy in the stress-driven curves is larger than in the strain-driven ones. Results from recently proposed models qualitatively agree with experiments.

A specific HPLC method for determination of beta-blockers topically used in ophthalmological diseases.

A HPLC method is presented for the identification and quantification in plasma and urine of beta-adrenergic receptor antagonists (betaxolol, carteolol, metipranolol, and timolol) commonly prescribed in ophthalmology. An extraction method is described using pindolol as an internal standard. An RSIL 10 micron column was used. The lower detection limits of the beta-blockers were found to be 4-27 ng/ml. This method is simple, rapid and sensitive; moreover, it allows the determination of 8 other beta-blockers.

Bargmann-Wigner method and (6S + 1)-component wave equations

A modified Bargmann-Wigner method is used to derive (6s + 1)-component wave equations. The relation between different forms of these equations is shown.

A method for using random parameters in analyzing permanent sample plots

Summary

Dating human skeletal remains using a radiometric method: Biogenic versus diagenetic (90)Sr and (210)Pb in vertebrae.

In forensic science, there is a strong interest in determining the post-mortem interval (PMI) of human skeletal remains up to 50 years after death. Currently, there are no reliable methods to resolve PMI, the determination of which relies almost exclusively on the experience of the investigating expert. Here we measured (90)Sr and (210)Pb ((210)Po) incorporated into bones through a biogenic process as indicators of the time elapsed since death. We hypothesised that the activity of radionuclides incorporated into trabecular bone will more accurately match the activity in the environment and the food chain at the time of death than the activity in cortical bone because of a higher remodelling rate. We found that determining (90)Sr can yield reliable PMI estimates as long as a calibration curve exists for (90)Sr covering the studied area and the last 50 years. We also found that adding the activity of (210)Po, a proxy for naturally occurring (210)Pb incorporated through ingestion, to the (90)Sr dating increases the reliability of the PMI value. Our results also show that trabecular bone is subject to both (90)Sr and (210)Po diagenesis. Accordingly, we used a solubility profile method to determine the biogenic radionuclide only, and we are proposing a new method of bone decontamination to be used prior to (90)Sr and (210)Pb dating.

An improved method for measuring soil microbial activity by gas phase flow injection analysis

The rate of carbon dioxide production is commonly used as a measure of microbial activity in the soil. The traditional method of CO2 determination involves trapping CO2 in an alkali solution and then determining CO2 concentration indirectly by titration of the remaining alkali in the solution. This method is still commonly employed in laboratories throughout the world due to its relative simplicity and the fact that it does not require expensive, specific equipment. However, there are several drawbacks: the method is time-consuming, requires large amounts of chemicals and the consistency of results depends on the operator's skills. With this in mind, an improved method was developed to analyze CO2 captured in alkali traps, which is cheap and relatively simple, with a substantially shorter sample handling time and reproducibility equivalent to the traditional titration method. A comparison of the concentration values determined by gas phase flow injection analysis (GPFIA) and titration showed no significant difference (p > 0.05), but GPFIA has the advantage that only a tenth of the sample volume of the titration method is required. The GPFIA system does not require the purchase of new, costly equipment but the device was constructed from items commonly found in laboratories, with suggestions for alternative configurations for other detection units. Furthermore, GPFIA for CO2 analysis can be equally applied to samples obtained from either the headspace of microcosms or from a sampling chamber that allows CO2 to be released from alkali trapping solutions. The optimised GPFIA method was applied to analyse CO2 released from degrading hydrocarbons from a site contaminated by diesel spillage.

Soil CO2 flux: a method comparison of closed static chambers in a sugarcane field

A large variety of techniques have been used to measure soil CO2 released from the soil surface, and much of the variability observed between locations must be attributed to the different methods used by the investigators. Therefore, a minimum protocol of measurement procedures should be established. The objectives of this study were (a) to compare different absorption areas, concentrations and volumes of the alkali trapping solution used in closed static chambers (CSC), and (b) to compare both, the optimized alkali trapping solution and the soda-lime trapping using CSC to measure soil respiration in sugarcane areas. Three CO2 absorption areas were evaluated (7; 15 and 20 % of the soil emission area or chamber); two volumes of NaOH (40 and 80 mL) at three concentrations (0.1, 0.25 and 0.5 mol L-1). Three different types of alkaline traps were tested: (a), 80 mL of 0.5 mol L-1 NaOH in glass containers, absorption area 15 % (V0.5); (b) 40 mL of 2 mol L-1 NaOH retained in a sponge, absorption area 80 % (S2) and (c) 40 g soda lime, absorption area 15 % (SL). NaOH concentrations of 0.5 mol L-1 or lower underestimated the soil CO2-C flux or CO2 flux. The lower limit of the alkali trap absorption area should be a minimum of 20 % of the area covered by the chamber. The 2 mol L-1 NaOH solution trap (S2) was the most efficient (highest accuracy and highest CO2 fluxes) in measuring soil respiration.

A simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Radial displacement of a fluid annulus in a rotating Hele¿Shaw cell

The radial displacement of a fluid annulus in a rotating circular Hele-Shaw cell has been investigated experimentally. It has been found that the flow depends sensitively on the wetting conditions at the outer interface. Displacements in a prewet cell are well described by Darcy's law in a wide range of experimental parameters, with little influence of capillary effects. In a dry cell, however, a more careful analysis of the interface motion is required; the interplay between a gradual loss of fluid at the inner interface, and the dependence of capillary forces at the outer interface on interfacial velocity and dynamic contact angle, result in a constant velocity for the interfaces. The experimental results in this case correlate in the form of an empirical scaling relation between the capillary number Ca and a dimensionless group, related to the ratio of centrifugal to capillary forces, which spans about three orders of magnitude in both quantities. Finally, the relative thickness of the coating film left by the inner interface, alpha i, is obtained as a function of Ca.

Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.