Formulating Mixed Models for Experiments, Including Longitudinal Experiments

Mixed models have become important in analyzing the results of experiments, particularly those that require more complicated models (e.g., those that involve longitudinal data). This article describes a method for deriving the terms in a mixed model. Our approach extends an earlier method by Brien and Bailey to explicitly identify terms for which autocorrelation and smooth trend arising from longitudinal observations need to be incorporated in the model. At the same time we retain the principle that the model used should include, at least, all the terms that are justified by the randomization. This is done by dividing the factors into sets, called tiers, based on the randomization and determining the crossing and nesting relationships between factors. The method is applied to formulate mixed models for a wide range of examples. We also describe the mixed model analysis of data from a three-phase experiment to investigate the effect of time of refinement on Eucalyptus pulp from four different sources. Cubic smoothing splines are used to describe differences in the trend over time and unstructured covariance matrices between times are found to be necessary.

Private valuation of carbon sequestration in forest plantations

Approval of the Clean Development Mechanism, provided for in the Kyoto Protocol, enables countries with afforested land to trade in carbon emissions reduction certificates related to carbon dioxide equivalent quantities (CO(2-e)) stored within a certain forest area. Potential CO(2-e) above base line sequestration was determined for two forest sites on commercial eucalyptus plantations in northern Brazil (Bahia). Compensation values for silvicultural regimes involving rotation lengths greater than economically optimal were computed using the Faustmann formula. Mean values obtained were US$8.16 (MgCO(2-e))(-1) and US $7.19 (MgCO(2-e))(-1) for average and high site indexes, respectively. Results show that carbon supply is more cost-efficient in highly productive sites. Annuities of US$18.8 Mg C(-1) and US$35.1 Mg C(-1) and yearly payments of US$4.4 m(-3) and US$8.2 m(-3) due for each marginal cubic meter produced were computed for high and average sites, respectively. The estimated value of the tonne of carbon defines minimum values to be paid to forest owners, in order to induce a change in silvicultural management regimes. A reduction of carbon supply could be expected as a result of an increase in wood prices, although it would not respond in a regular manner. For both sites, price elasticity of supply was found to be inelastic and increased as rotation length moved further away from economically optimal: 0.24 and 0.27 for age 11 years in average- and high-productivity sites, respectively. This would be due to biomass production potential as a limiting factor; beyond a certain threshold value. an increase in price does not sustain a proportional change in carbon storage supply. The environmental service valuation model proposed might be adequate for assessing potential supply in plantation forestry, from a private landowner perspective, with an economic opportunity cost. The model is not applicable to low commercial value forest plantations. (C) 2009 Elsevier B.V. All rights reserved.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

Milk yield, milk composition, eating behavior, and lamb performance of ewes fed diets containing soybean hulls replacing coastcross (Cynodon species) hay

The objective of this experiment was to evaluate the effects of replacing coastcross hay NDF by soybean hull (SH) NDF on the lactation performance and eating behavior of ewes and also on the performance of their lambs. Fifty-six Santa Ines lactating ewes (56.1 +/- 6.8 kg of initial BW; mean +/- SD) were penned individually and used in a randomized complete block design with 14 blocks and 4 treatments. Diets were formulated to provide similar concentrations of NDF (56%) and CP (16%). The SH NDF replaced 33 (SH33), 67 (SH67), or 100% (SH100) of the NDF contributed by coastcross hay in a 70% forage-based diet (SH0), resulting in SH inclusion rates of 0, 25, 54, and 85% of the dietary DM. Once a week, from the second to the eighth week of lactation (weaning time), ewes were separated from their lambs, stimulated by a 6-IU i.v. oxytocin injection, and hand milked to empty the udder. After 3 h, milk production was obtained after the same procedure. Quadratic effect for milk production (142.4, 179.8, 212.6, and 202.9 g/3 h) and cubic effect for DMI (2.27, 2.69, 3.25, and 3.00 kg/d) were observed as SH inclusion increased from 0 to 85% of the dietary DM. Milk fat (7.59, 7.86, 7.59, and 7.74%), protein (4.53, 4.43, 4.40, and 4.55%), and total solids (18.24, 18.54, 18.39, and 18.64%) did not differ among the 70% forage-based diet and diets with SH NDF replacing 33, 67, or 100% of the NDF. A linear increase in lactose concentration was observed with SH inclusion. Ewe BW gain during the trial showed a cubic response (0.37, 0.03, 4.80, and 2.80 kg) with SH inclusion. The preweaning ADG of lambs increased linearly, and ADG of lambs after weaning decreased linearly with SH inclusion. Final BW of lambs (2 wk after weaning) did not differ among treatments. Eating behavior observations were conducted with 44 ewes. The same facilities, experimental design, dietary treatments, and feeding management were used. Observations were visually recorded every 5 min for a 24-h period when ewes were 46 +/- 6.8 d in milk. Eating time (min/d, min/g of DMI, and min/g of NDF intake) and time expended in rumination and chewing activities (min/g of DMI and min/g of NDF intake) decreased linearly with the addition of SH in the diets. The inclusion of SH improved DMI and milk production, also reflecting on the BW of lambs at weaning. Milk performance was not affected when SH NDF replaced 100% of hay NDF.

Enantioselective Arylations Catalyzed by Carbohydrate-Based Chiral Amino Alcohols

The application of carbohydrate-derived amino alcohols in the asymmetric arylation of aldehydes by using arylboronic acids as the source of transferable aryl groups is described. The best ligand is derived from the readily available sugar D-xylose and it mediates the addition of a range of arylboronic acids to various aromatic aldehydes in excellent yields and high enantiomeric excesses.

Optimization of preservative system in ophthalmic suspension with dexamethasone and polymyxin B

A simplex-lattice statistical project was employed to study an optimization method for a preservative system in an ophthalmic suspension of dexametasone and polymyxin B. The assay matrix generated 17 formulas which were differentiated by the preservatives and EDTA (disodium ethylene diamine-tetraacetate), being the independent variable: X-1 = chlorhexidine digluconate (0.010 % w/v); X-2 = phenylethanol (0.500 % w/v); X-3 = EDTA (0.100 % w/v). The dependent variable was the Dvalue obtained from the microbial challenge of the formulas and calculated when the microbial killing process was modeled by an exponential function. The analysis of the dependent variable, performed using the software Design Expert/W, originated cubic equations with terms derived from stepwise adjustment method for the challenging microorganisms: Pseudomonas aeruginosa, Burkholderia cepacia, Staphylococcus aureus, Candida albicans and Aspergillus niger. Besides the mathematical expressions, the response surfaces and the contour graphics were obtained for each assay. The contour graphs obtained were overlaid in order to permit the identification of a region containing the most adequate formulas (graphic strategy), having as representatives: X-1 = 0.10 ( 0.001 % w/v); X-2 = 0.80 (0.400 % w/v); X-3 = 0.10 (0.010 % w/v). Additionally, in order to minimize responses (Dvalue), a numerical strategy corresponding to the use of the desirability function was used, which resulted in the following independent variables combinations: X-1 = 0.25 (0.0025 % w/v); X-2 = 0.75 (0.375 % w/v); X-3 = 0. These formulas, derived from the two strategies (graphic and numerical), were submitted to microbial challenge, and the experimental Dvalue obtained was compared to the theoretical Dvalue calculated from the cubic equation. Both Dvalues were similar to all the assays except that related to Staphylococcus aureus. This microorganism, as well as Pseudomonas aeruginosa, presented intense susceptibility to the formulas independently from the preservative and EDTA concentrations. Both formulas derived from graphic and numerical strategies attained the recommended criteria adopted by the official method. It was concluded that the model proposed allowed the optimization of the formulas in their preservation aspect.

Synergism on antioxidant activity between natural compounds optimized by response surface methodology

Despite the increase in the use of natural compounds in place of synthetic derivatives as antioxidants in food products, the extent of this substitution is limited by cost constraints. Thus, the objective of this study was to explore the synergism on the antioxidant activity of natural compounds, for further application in food products. Three hydrosoluble compounds (x(1) = caffeic acid, x(2) = carnosic acid, and x(3) = glutathione) and three liposoluble compounds (x(1) = quercetin, x(2) = rutin, and x(3) = genistein) were mixed according to a ""centroid simplex design"". The antioxidant activity of the mixtures was analyzed by the ferric reducing antioxidant power (FRAP) and oxygen radical absorbance capacity (ORAL) methodologies, and activity was also evaluated in an oxidized mixed micelle prepared with linoleic acid (LAOX). Cubic polynomial models with predictive capacity were obtained when the mixtures were submitted to the LAOX methodology ((y) over cap = 0.56 x(1) + 0.59 x(2) + 0.04 x(3) + 0.41 x(1)x(2) - 0.41 x(1)x(3) - 1.12 x(2)x(3) - 4.01 x(1)x(2)x(3)) for the hydrosoluble compounds, and to FRAP methodology ((y) over cap = 3.26 x(1) + 2.39 x(2) + 0.04 x(3) + 1.51 x(1)x(2) + 1.03 x(1)x(3) + 0.29 x(1)x(3) + 3.20 x(1)x(2)x(3)) for the liposoluble compounds. Optimization of the models suggested that a mixture containing 47% caffeic acid + 53% carnosic acid and a mixture containing 67% quercetin + 33% rutin were potential synergistic combinations for further evaluation using a food matrix.

Relationship between Rheological Properties and One-Step W/O/W Multiple Emulsion Formation

Formation of a normal (not temporary) W/O/W multiple emulsion via the one-step method as a result of the simultaneous occurrence of catastrophic and transitional phase inversion processes has been recently reported. Critical features of this process include the emulsification temperature (corresponding to the ultralow surface tension point), the use of a specific nonionic surfactant blend and the surfactant blend/oil phase ratio, and the addition of the surfactant blend to the oil phase. The purpose of this study was to investigate physicochemical properties in an effort to gain a mechanistic understanding of the formation of these emulsions. Bulk, surface, and interfacial theological properties of adsorbed nonionic surfactant (CremophorRH40 and Span80) films were investigated under conditions known to affect W/O/W emulsion formation. Bulk viscosity results demonstrated that CremophorRH40 has a higher mobility in oil compared than in water, explaining the significance of the solvent phase. In addition, the bulk viscosity profile of aqueous solutions containing CremophorRH40 indicated a phase transition at around 78 +/- 2 degrees C, which is in agreement with cubic phase formation in the Winsor III region. The similarity in the interfacial elasticity values of CremophorRH40 and Span80 indicated that canola oil has a major effect on surface activity, showing the significance of vegetable oil. The highest interfacial shear elasticity and viscosity were observed when both surfactants were added to the oil phase, indicating the importance of the microstructural arrangement. CremophorRH40/Span80 complexes tended to desorb from the solution/solution interface with increasing temperature, indicating surfactant phase formation as is theoretically predicted in the Winsor III region. Together these interfacial and bulk rheology data demonstrate that one-step W/O/W emulsions form as a result of the simultaneous occurrence of phase-transition processes in the Winsor III region and explain the critical formulation and processing parameters necessary to achieve the formation of these normal W/O/W emulsions.

Oxidation mechanism of Si3N4-bonded SiC ceramics by CO, CO2 and steam

This paper presents a theoretical and experimental investigation into the oxidation reactions of Si3N4-bonded SiC ceramics. Such ceramics which contain a small amount of silicon offer increased oxidation and wear resistance and are widely used as lining refractories in blast furnaces. The thermodynamics of oxidation reactions were studied using the JANAF tables. The weight gain was measured using a thermogravimetric analysis technique to study the kinetics. The temperature range of oxidation measurements is from 1073 to 1573 K and the oxidation atmosphere is water vapour, pure CO and CO-CO2 gas mixtures with various CO-to-CO2 ratios. Thermodynamic simulations showed that the oxidation mechanism of Si3N4-bonded SiC ceramics is passive oxidation and all components contribute to the formation of a silica film. The activated energies of the reactions follow the sequence Si3N4>SiC>Si. The kinetic study revealed that the oxidation of Si3N4-bonded SiC ceramics occurred in a mixed regime controlled by both interface reaction and diffusion through the silica film. Under the atmosphere conditions prevailing in the blast furnace, this ceramic is predicted to be passively oxidized with the chemical reaction rate becoming more dominant as the CO concentration increases. (C) 1998 Chapman & Hall.

Three-dimensional quantum solitons with parametric coupling

We consider the quantum field theory of two bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium, this corresponds to the process of second-harmonic generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. The quantum solitons or energy eigenstates (bound-state solutions) are obtained exactly in the simplest case of two-particle binding, in one, two, and three space dimensions. We also investigate three-particle binding in one space dimension. The results indicate that the exact quantum solitons of this field theory have a singular, pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. To estimate the physically accessible radii and binding energies of the bound states, we impose a momentum cutoff on the nonlinear couplings. In the case of nonlinear optical interactions, the resulting radii and binding energies of these photonic particlelike excitations in highly nonlinear parametric media appear to be close to physically observable values.

Multidimensional parametric quantum solitons

We consider the parametric quantum field theory involving cubic and quartic couplings of two bosonic fields. This is exactly soluble for the two-particle energy eigenstates (or quantum solitons) in one, two, and three space dimensions. We estimate the binding energies and corresponding radii in the case of photonic fields in nonlinear optical materials, and Bose-Einstein condensates. [S1050-2947(98)51110-9].

Peak bone mineral accrual and age at menarche in adolescent girls: A 6-year longitudinal study

Background and objectives: The greatest increase in bone mineral content occurs during adolescence. The amount of bone accrued may significantly affect bone mineral status in later life. We carried out a longitudinal investigation of the magnitude and timing of peak bone mineral content velocity (PBMCV) in relation to peak height velocity (PHV) and the age at menarche in a group of adolescent girls over a 6-year period. Methods: The 53 girls in this study are a subset of the 115 girls (initially 8 to 16 years) in a g-year longitudinal study of bone mineral accretion. The ages at PBMCV and PHV were determined by using a cubic spline curve fitting procedure. Determinations were based on height (n = 12) and bone (n = 6) measurements over 6 years. Results: The timing of PBMCV and menarche were coincident, preceded approximately 1 year earlier by PHV. Correlation showed a negative relationship between age at menarche and both peak bone mineral accrual (r = -0.42, P

Novel solitons in parametric amplifiers and atom lasers

We review recent developments in quantum and classical soliton theory, leading to the possibility of observing both classical and quantum parametric solitons in higher-dimensional environments. In particular, we consider the theory of three bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium this corresponds to the process of sum frequency generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. Potential applications include an ultrafast photonic AND-gate. The simplest quantum solitons or energy eigenstates (bound-state solutions) of the interacting field Hamiltonian are obtained exactly in three space dimensions. They have a point-like structure-even though the corresponding classical theory is nonsingular. We show that the solutions can be regularized with the imposition of a momentum cut-off on the nonlinear couplings. The case of three-dimensional matter-wave solitons in coupled atomic/molecular Bose-Einstein condensates is discussed.

A six-year longitudinal study of the relationship of physical activity to bone mineral accrual in growing children: The University of Saskatchewan bone mineral accrual study

To investigate the influence of physical activity on bone mineral accrual during the adolescent years, we analyzed 6 years of data from 53 girls and 60 boys. Physical activity, dietary intakes, and anthropometry were measured every 6 months and dual-energy X-ray absorptiometry scans of the total body (TB), lumbar spine (LS), and proximal femur (Hologic 2000, array mode) were collected annually. Distance and velocity curves for height and bone mineral content (BMC) were fitted for each child at several skeletal sites using a cubic spline procedure, from which ages at peak height velocity (PHV) and peak BMC velocity (PBMCV) were identified. A mean age- and gender-specific standardized activity (Z) score was calculated for each subject based on multiple yearly activity assessments collected up until age of PHV. This score was used to identify active (top quartile), average (middle 2 quartiles), or inactive (bottom quartile) groups. Two-way analysis of covariance, with height and weight at PHV controlled for, demonstrated significant physical activity and gender main effects (but no interaction) for PBMCV, for BMC accrued for 2 years around peak velocity, and for BMC at 1 year post-PBMCV for the TB and femoral neck and for physical activity but not gender at the LS (all p < 0.05). Controlling for maturational and size differences between groups, we noted a 9% and 17% greater TB BMC for active boys and girls, respectively, over their inactive peers 1 year after the age of PBMCV. We also estimated that, on average, 26% of adult TB bone mineral was accrued during the 2 years around PBMCV.

Disagreement between correlations of quantum mechanics and stochastic electrodynamics in the damped parametric oscillator

Intracavity and external third order correlations in the damped nondegenerate parametric oscillator are calculated for quantum mechanics and stochastic electrodynamics (SED), a semiclassical theory. The two theories yield greatly different results, with the correlations of quantum mechanics being cubic in the system's nonlinear coupling constant and those of SED being linear in the same constant. In particular, differences between the two theories are present in at least a mesoscopic regime. They also exist when realistic damping is included. Such differences illustrate distinctions between quantum mechanics and a hidden variable theory for continuous variables.