Designs for generalized linear models with several variables and model uncertainty

Standard factorial designs sometimes may be inadequate for experiments that aim to estimate a generalized linear model, for example, for describing a binary response in terms of several variables. A method is proposed for finding exact designs for such experiments that uses a criterion allowing for uncertainty in the link function, the linear predictor, or the model parameters, together with a design search. Designs are assessed and compared by simulation of the distribution of efficiencies relative to locally optimal designs over a space of possible models. Exact designs are investigated for two applications, and their advantages over factorial and central composite designs are demonstrated.

Effects of phytase supplementation of diets with two tiers of nutrient specifications on growth performance and protein efficiency ratios of broiler chickens

in two feeding experiments male and mixed-sex broiler chicks were offered diets based on sorghum and a wheat-sorghum blend with two tiers of nutrient specifications, without and with microbial phytase (600 and 800 FTU/kg), from 7-25 and 1-42 days post-hatch, respectively. The nutrient specifications for protein, amino acids, energy density and phosphorus (P) of standard diets were reduced to formulate the modified diets on a least-cost basis. Calculated differences in nutrient specifications between standard and modified diets ranged from 14.3 to 17.1 g/kg crude protein, 0.24 to 0.40 MJ/kg apparent metabolisable energy (AME) and 1.06 to 1.20 g/kg available P. In both experiments, reduced nutrient specifications had a negative impact on growth rates and feed efficiency and phytase supplementation had a positive influence on growth performance and protein efficiency ratios (PER). Phytase addition to the less expensive, modified diets either partially or entirely compensated for reduced growth performance and, consequently, feed costs per kg of live weight gain were reduced. In Experiment 1, phytase increased (p<0.001) nitrogen-corrected AME (AMEn) from 15.39 to 15.89 MJ/kg dry matter. For nitrogen (N) retention there was an interaction (p<0.05) between diet type and phytase as the effects of phytase on N retention were more pronounced in the modified diets, with an increase from 0.512 to 0.561. These results demonstrate the positive effects of phytase on protein and energy utilisation, in addition to its established liberation of phytate-bound P and illustrate the feasibility of assigning nutrient replacement values to the feed enzyme for consideration in least-cost ration formulations. Further work is, however, required to define the most appropriate reductions in nutrient specifications in association with phytase supplementation.

Qubit channels that achieve capacity with two states

This paper considers a class of qubit channels for which three states are always sufficient to achieve the Holevo capacity. For these channels, it is known that there are cases where two orthogonal states are sufficient, two nonorthogonal states are required, or three states are necessary. Here a systematic theory is given which provides criteria to distinguish cases where two states are sufficient, and determine whether these two states should be orthogonal or nonorthogonal. In addition, we prove a theorem on the form of the optimal ensemble when three states are required, and present efficient methods of calculating the Holevo capacity.

Extinction times for a birth-death process with two phases

Many populations have a negative impact on their habitat or upon other species in the environment if their numbers become too large. For this reason they are often subjected to some form of control. One common control regime is the reduction regime: when the population reaches a certain threshold it is controlled (for example culled) until it falls below a lower predefined level. The natural model for such a controlled population is a birth-death process with two phases, the phase determining which of two distinct sets of birth and death rates governs the process. We present formulae for the probability of extinction and the expected time to extinction, and discuss several applications. (c) 2006 Elsevier Inc. All rights reserved.

Phosphorus nutrition of Caustis blakei grown with two phosphorus sources of different solubility in two soils of differing phosphorus adsorption capacity

A barrier to the domestication of the phosphorus (P) sensitive Australian species Caustis blakei (Cyperaceae) is the standard production systems used commercially which invariably result in problems associated either with P deficiency or P toxicity. This paper reports on the growth responses of Caustis blakei cv. M63 to applications of fertiliser P as either monocalcium phosphate (MCP) or granulated Guano Gold (R) rock phosphate (RP) in two soils with different capacities to adsorb P. The Caustis M63 plants grown in the two soils did not show P toxicity symptoms when fertilised with RP, but shoot dry weight was 30-60% lower than the control in both soils at the highest rate of MCP-P application (156 kg ha(-1), 184 g m(-3)) and this was associated with visible symptoms of drying of the tips of the ultimate branchlets, in the Mt Cotton soil only. The greatest shoot and root dry weights were achieved by plants grown in the higher P adsorbing Palmwoods soil fertilised with RP at P rates of 30-184 g m(-3). Caustis plants grown in the Palmwoods soil had 2.3 times greater root dry weights than plants grown in the Mt Cotton soil irrespective of the P fertiliser type used. Caustis plants growing in Mt Cotton soil which did not receive P showed significantly lower shoot and root dry weight when compared to plants in the Palmwoods soil, probably due to the low initial bicarbonate-extractable P and the high buffering capacity of the Mt Cotton soil. The P concentration in shoots of Caustis fertilised with MCP at 184 g m(-3) was higher when grown in Mt Cotton soil (0.22%) than in the Palmwoods soil (0.15%). The P concentration was lower in the terminal ultimate branchlets (TUB); 0.15% for the Mt Cotton soil and 0.10% for the Palmwoods soil, suggesting that shoots would provide a more useful indicator of P toxicity than the TUB. It is interesting to speculate as to why plants in the Palmwoods soil showed greater root growth and fewer symptoms of P toxicity. This could be because the Palmwoods soil had the greater P adsorption capacity. These results indicate in ground production of Caustis cut foliage will require careful management of P nutrition and understanding of the complex soil/plant interactions associated with the acquisition of P. (c) 2006 Elsevier B.V. All rights reserved.

Steiner triple systems with two disjoint subsystems

It is shown that there exists a triangle decomposition of the graph obtained from the complete graph of order v by removing the edges of two vertex disjoint complete subgraphs of orders u and w if and only if u, w, and v are odd, ((v)(2)) - ((u)(2)) - ((w)(2)) equivalent to 0 (mod 3), and v >= w + u + max {u, w}. Such decompositions are equivalent to group divisible designs with block size 3, one group of size u, one group of size w, and v - u - w groups of size 1. This result settles the existence problem for Steiner triple systems having two disjoint specified subsystems, thereby generalizing the well-known theorem of Doyen and Wilson on the existence of Steiner triple systems with a single specified subsystem. (c) 2005 Wiley Periodicals, Inc.

Drying process of microcrystalline cellulose studied by attenuated total reflection IR spectroscopy with two-dimensional correlation spectroscopy and principal component analysis

Molecular interactions between microcrystalline cellulose (MCC) and water were investigated by attenuated total reflection infrared (ATR/IR) spectroscopy. Moisture-content-dependent IR spectra during a drying process of wet MCC were measured. In order to distinguish overlapping O–H stretching bands arising from both cellulose and water, principal component analysis (PCA) and, generalized two-dimensional correlation spectroscopy (2DCOS) and second derivative analysis were applied to the obtained spectra. Four typical drying stages were clearly separated by PCA, and spectral variations in each stage were analyzed by 2DCOS. In the drying time range of 0–41 min, a decrease in the broad band around 3390 cm−1 was observed, indicating that bulk water was evaporated. In the drying time range of 49–195 min, decreases in the bands at 3412, 3344 and 3286 cm−1 assigned to the O6H6cdots, three dots, centeredO3′ interchain hydrogen bonds (H-bonds), the O3H3cdots, three dots, centeredO5 intrachain H-bonds and the H-bonds in Iβ phase in MCC, respectively, were observed. The result of the second derivative analysis suggests that water molecules mainly interact with the O6H6cdots, three dots, centeredO3′ interchain H-bonds. Thus, the H-bonding network in MCC is stabilized by H-bonds between OH groups constructing O6H6cdots, three dots, centeredO3′ interchain H-bonds and water, and the removal of the water molecules induces changes in the H-bonding network in MCC.

A new analytical solution for water table fluctuations in coastal aquifers with sloping beaches

The Boussinesq equation appears as the zeroth-order term in the shallow water flow expansion of the non-linear equation describing the flow of fluid in an unconfined aquifer. One-dimensional models based on the Boussinesq equation have been used to analyse tide-induced water table fluctuations in coastal aquifers. Previous analytical solutions for a sloping beach are based on the perturbation parameter, epsilon(N) = alphaepsilon cot beta (in which beta is the beach slope, alpha is the amplitude parameter and epsilon is the shallow water parameter) and are limited to tan(-1) (alphaepsilon) much less than beta less than or equal to pi/2. In this paper, a new higher-order solution to the non-linear boundary value problem is derived. The results demonstrate the significant influence of the higher-order components and beach slope on the water table fluctuations. The relative difference between the linear solution and the present solution increases as 6 and a increase, and reaches 7% of the linear solution. (C) 2003 Elsevier Ltd. All rights reserved.

Comparative proteomic analysis of GS-NSO murine myeloma cell lines with varying recombinant monoclonal antibody production rate

We have employed an inverse engineering strategy based on quantitative proteome analysis to identify changes in intracellular protein abundance that correlate with increased specific recombinant monoclonal antibody production (qMab) by engineered murine myeloma (NSO) cells. Four homogeneous NSO cell lines differing in qMab were isolated from a pool of primary transfectants. The proteome of each stably transfected cell line was analyzed at mid-exponential growth phase by two-dimensional gel electrophoresis (2D-PAGE) and individual protein spot volume data derived from digitized gel images were compared statistically. To identify changes in protein abundance associated with qMab clatasets were screened for proteins that exhibited either a linear correlation with cell line qMab or a conserved change in abundance specific only to the cell line with highest qMab. Several proteins with altered abundance were identified by mass spectrometry. Proteins exhibiting a significant increase in abundance with increasing qMab included molecular chaperones known to interact directly with nascent immunoglobulins during their folding and assembly (e.g., BiP, endoplasmin, protein disulfide isomerase). 2D-PAGE analysis showed that in all cell lines Mab light chain was more abundant than heavy chain, indicating that this is a likely prerequisite for efficient Mab production. In summary, these data reveal both the adaptive responses and molecular mechanisms enabling mammalian cells in culture to achieve high-level recombinant monoclonal antibody production. (C) 2004 Wiley Periodicals, Inc.

Conditional production of superpositions of coherent states with inefficient photon detection

It is shown that a linear superposition of two macroscopically distinguishable optical coherent states can be generated using a single photon source and simple all-optical operations. Weak squeezing on a single photon, beam mixing with an auxiliary coherent state, and photon detecting with imperfect threshold detectors are enough to generate a coherent state superposition in a free propagating optical field with a large coherent amplitude (alpha>2) and high fidelity (F>0.99). In contrast to all previous schemes to generate such a state, our scheme does not need photon number resolving measurements nor Kerr-type nonlinear interactions. Furthermore, it is robust to detection inefficiency and exhibits some resilience to photon production inefficiency.

Constriction flows of monodisperse linear entangled polymers: Multiscale modeling and flow visualization

We explore both the rheology and complex flow behavior of monodisperse polymer melts. Adequate quantities of monodisperse polymer were synthesized in order that both the materials rheology and microprocessing behavior could be established. In parallel, we employ a molecular theory for the polymer rheology that is suitable for comparison with experimental rheometric data and numerical simulation for microprocessing flows. The model is capable of matching both shear and extensional data with minimal parameter fitting. Experimental data for the processing behavior of monodisperse polymers are presented for the first time as flow birefringence and pressure difference data obtained using a Multipass Rheometer with an 11:1 constriction entry and exit flow. Matching of experimental processing data was obtained using the constitutive equation with the Lagrangian numerical solver, FLOWSOLVE. The results show the direct coupling between molecular constitutive response and macroscopic processing behavior, and differentiate flow effects that arise separately from orientation and stretch. (c) 2005 The Society of Rheology.

A new repair-replace strategy for items sold with a two-dimensional warranty

For repairable items, the manufacturer has the option to either repair or replace a failed item that is returned under warranty. In this paper, we look at a new warranty servicing strategy for items sold with two-dimensional warranty where the failed item is replaced by a new one when it fails for the first time in a specified region of the warranty and all other failures are repaired minimally. The region is characterised by two parameters and we derive the optimal values for these to minimise the total expected warranty servicing cost. We compare the results with other repair-replace strategies reported in the literature. (C) 2003 Elsevier Ltd. All rights reserved.

Estimating variable returns to scale production frontiers with alternative stochastic assumptions

Two stochastic production frontier models are formulated within the generalized production function framework popularized by Zellner and Revankar (Rev. Econ. Stud. 36 (1969) 241) and Zellner and Ryu (J. Appl. Econometrics 13 (1998) 101). This framework is convenient for parsimonious modeling of a production function with returns to scale specified as a function of output. Two alternatives for introducing the stochastic inefficiency term and the stochastic error are considered. In the first the errors are added to an equation of the form h(log y, theta) = log f (x, beta) where y denotes output, x is a vector of inputs and (theta, beta) are parameters. In the second the equation h(log y,theta) = log f(x, beta) is solved for log y to yield a solution of the form log y = g[theta, log f(x, beta)] and the errors are added to this equation. The latter alternative is novel, but it is needed to preserve the usual definition of firm efficiency. The two alternative stochastic assumptions are considered in conjunction with two returns to scale functions, making a total of four models that are considered. A Bayesian framework for estimating all four models is described. The techniques are applied to USDA state-level data on agricultural output and four inputs. Posterior distributions for all parameters, for firm efficiencies and for the efficiency rankings of firms are obtained. The sensitivity of the results to the returns to scale specification and to the stochastic specification is examined. (c) 2004 Elsevier B.V. All rights reserved.