Comparison of sample size formulae for 2x2 cross-over designs applied to bioequivalence studies

We consider the comparison of two formulations in terms of average bioequivalence using the 2 × 2 cross-over design. In a bioequivalence study, the primary outcome is a pharmacokinetic measure, such as the area under the plasma concentration by time curve, which is usually assumed to have a lognormal distribution. The criterion typically used for claiming bioequivalence is that the 90% confidence interval for the ratio of the means should lie within the interval (0.80, 1.25), or equivalently the 90% confidence interval for the differences in the means on the natural log scale should be within the interval (-0.2231, 0.2231). We compare the gold standard method for calculation of the sample size based on the non-central t distribution with those based on the central t and normal distributions. In practice, the differences between the various approaches are likely to be small. Further approximations to the power function are sometimes used to simplify the calculations. These approximations should be used with caution, because the sample size required for a desirable level of power might be under- or overestimated compared to the gold standard method. However, in some situations the approximate methods produce very similar sample sizes to the gold standard method. Copyright © 2005 John Wiley & Sons, Ltd.

Some general comparative points on Chao's and Zelterman's estimators of the population size

Two simple and frequently used capture–recapture estimates of the population size are compared: Chao's lower-bound estimate and Zelterman's estimate allowing for contaminated distributions. In the Poisson case it is shown that if there are only counts of ones and twos, the estimator of Zelterman is always bounded above by Chao's estimator. If counts larger than two exist, the estimator of Zelterman is becoming larger than that of Chao's, if only the ratio of the frequencies of counts of twos and ones is small enough. A similar analysis is provided for the binomial case. For a two-component mixture of Poisson distributions the asymptotic bias of both estimators is derived and it is shown that the Zelterman estimator can experience large overestimation bias. A modified Zelterman estimator is suggested and also the bias-corrected version of Chao's estimator is considered. All four estimators are compared in a simulation study.

Extending Zelterman’s approach for robust estimation of population size to zero-truncated clustered data

Estimation of population size with missing zero-class is an important problem that is encountered in epidemiological assessment studies. Fitting a Poisson model to the observed data by the method of maximum likelihood and estimation of the population size based on this fit is an approach that has been widely used for this purpose. In practice, however, the Poisson assumption is seldom satisfied. Zelterman (1988) has proposed a robust estimator for unclustered data that works well in a wide class of distributions applicable for count data. In the work presented here, we extend this estimator to clustered data. The estimator requires fitting a zero-truncated homogeneous Poisson model by maximum likelihood and thereby using a Horvitz-Thompson estimator of population size. This was found to work well, when the data follow the hypothesized homogeneous Poisson model. However, when the true distribution deviates from the hypothesized model, the population size was found to be underestimated. In the search of a more robust estimator, we focused on three models that use all clusters with exactly one case, those clusters with exactly two cases and those with exactly three cases to estimate the probability of the zero-class and thereby use data collected on all the clusters in the Horvitz-Thompson estimator of population size. Loss in efficiency associated with gain in robustness was examined based on a simulation study. As a trade-off between gain in robustness and loss in efficiency, the model that uses data collected on clusters with at most three cases to estimate the probability of the zero-class was found to be preferred in general. In applications, we recommend obtaining estimates from all three models and making a choice considering the estimates from the three models, robustness and the loss in efficiency. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Equivalence of truncated count mixture distributions and mixtures of truncated count distributions

This article is about modeling count data with zero truncation. A parametric count density family is considered. The truncated mixture of densities from this family is different from the mixture of truncated densities from the same family. Whereas the former model is more natural to formulate and to interpret, the latter model is theoretically easier to treat. It is shown that for any mixing distribution leading to a truncated mixture, a (usually different) mixing distribution can be found so. that the associated mixture of truncated densities equals the truncated mixture, and vice versa. This implies that the likelihood surfaces for both situations agree, and in this sense both models are equivalent. Zero-truncated count data models are used frequently in the capture-recapture setting to estimate population size, and it can be shown that the two Horvitz-Thompson estimators, associated with the two models, agree. In particular, it is possible to achieve strong results for mixtures of truncated Poisson densities, including reliable, global construction of the unique NPMLE (nonparametric maximum likelihood estimator) of the mixing distribution, implying a unique estimator for the population size. The benefit of these results lies in the fact that it is valid to work with the mixture of truncated count densities, which is less appealing for the practitioner but theoretically easier. Mixtures of truncated count densities form a convex linear model, for which a developed theory exists, including global maximum likelihood theory as well as algorithmic approaches. Once the problem has been solved in this class, it might readily be transformed back to the original problem by means of an explicitly given mapping. Applications of these ideas are given, particularly in the case of the truncated Poisson family.

Influence of environmental stress on distributions of times to first division in Escherichia coli populations, as determined by digital-image analysis of individual cells

The distributions of times to first cell division were determined for populations of Escherichia coli stationary-phase cells inoculated onto agar media. This was accomplished by using automated analysis of digital images of individual cells growing on agar and calculation of the "box area ratio." Using approximately 300 cells per experiment, the mean time to first division and standard deviation for cells grown in liquid medium at 37 degrees C and inoculated on agar and incubated at 20 degrees C were determined as 3.0 h and 0.7 h, respectively. Distributions were observed to tail toward the higher values, but no definitive model distribution was identified. Both preinoculation stress by heating cultures at 50 degrees C and postinoculation stress by growth in the presence of higher concentrations of NaCl increased mean times to first division. Both stresses also resulted in an increase in the spread of the distributions that was proportional to the mean division time, the coefficient of variation being constant at approximately 0.2 in all cases. The "relative division time," which is the time to first division for individual cells expressed in terms of the cell size doubling time, was used as measure of the "work to be done" to prepare for cell division. Relative division times were greater for heat-stressed cells than for those growing under osmotic stress.

A bagging-based correction for the mixture model estimator of population size

Estimation of a population size by means of capture-recapture techniques is an important problem occurring in many areas of life and social sciences. We consider the frequencies of frequencies situation, where a count variable is used to summarize how often a unit has been identified in the target population of interest. The distribution of this count variable is zero-truncated since zero identifications do not occur in the sample. As an application we consider the surveillance of scrapie in Great Britain. In this case study holdings with scrapie that are not identified (zero counts) do not enter the surveillance database. The count variable of interest is the number of scrapie cases per holding. For count distributions a common model is the Poisson distribution and, to adjust for potential heterogeneity, a discrete mixture of Poisson distributions is used. Mixtures of Poissons usually provide an excellent fit as will be demonstrated in the application of interest. However, as it has been recently demonstrated, mixtures also suffer under the so-called boundary problem, resulting in overestimation of population size. It is suggested here to select the mixture model on the basis of the Bayesian Information Criterion. This strategy is further refined by employing a bagging procedure leading to a series of estimates of population size. Using the median of this series, highly influential size estimates are avoided. In limited simulation studies it is shown that the procedure leads to estimates with remarkable small bias.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.

Synergetic effect of carbon nanopore size and surface oxidation on CO2 capture from CO2/CH4 mixtures

We have studied the synergetic effect of confinement (carbon nanopore size) and surface chemistry (the number of carbonyl groups) on CO2 capture from its mixtures with CH4 at typical operating conditions for industrial adsorptive separation (298 K and compressed CO2CH4 mixtures). Although both confinement and surface oxidation have an impact on the efficiency of CO2/CH4 adsorptive separation at thermodynamics equilibrium, we show that surface functionalization is the most important factor in designing an efficient adsorbent for CO2 capture. Systematic Monte Carlo simulations revealed that adsorption of CH4 either pure or mixed with CO2 on oxidized nanoporous carbons is only slightly increased by the presence of functional groups (surface dipoles). In contrast, adsorption of CO2 is very sensitive to the number of carbonyl groups, which can be examined by a strong electric quadrupolar moment of CO2. Interestingly, the adsorbed amount of CH4 is strongly affected by the presence of the co-adsorbed CO2. In contrast, the CO2 uptake does not depend on the molar ratio of CH4 in the bulk mixture. The optimal carbonaceous porous adsorbent used for CO2 capture near ambient conditions should consist of narrow carbon nanopores with oxidized pore walls. Furthermore, the equilibrium separation factor was the greatest for CO2/CH4 mixtures with a low CO2 concentration. The maximum equilibrium separation factor of CO2 over CH4 of ∼18–20 is theoretically predicted for strongly oxidized nanoporous carbons. Our findings call for a review of the standard uncharged model of carbonaceous materials used for the modeling of the adsorption separation processes of gas mixtures containing CO2 (and other molecules with strong electric quadrupolar moment or dipole moment).

Assessing the impact of deforestation and climate change on the range size and environmental niche of bird species in the Atlantic forests, Brazil

Aim Habitat loss and climate change are two major drivers of biological diversity. Here we quantify how deforestation has already changed, and how future climate scenarios may change, environmental conditions within the highly disturbed Atlantic forests of Brazil. We also examine how environmental conditions have been altered within the range of selected bird species. Location Atlantic forests of south-eastern Brazil. Methods The historical distribution of 21 bird species was estimated using Maxent. After superimposing the present-day forest cover, we examined the environmental niches hypothesized to be occupied by these birds pre- and post-deforestation using environmental niche factor analysis (ENFA). ENFA was also used to compare conditions in the entire Atlantic forest ecosystem pre- and post-deforestation. The relative influence of land use and climate change on environmental conditions was examined using analysis of similarity and principal components analysis. Results Deforestation in the region has resulted in a decrease in suitable habitat of between 78% and 93% for the Atlantic forest birds included here. Further, Atlantic forest birds today experience generally wetter and less seasonal forest environments than they did historically. Models of future environmental conditions within forest remnants suggest generally warmer conditions and lower annual variation in rainfall due to greater precipitation in the driest quarter of the year. We found that deforestation resulted in a greater divergence of environmental conditions within Atlantic forests than that predicted by climate change. Main conclusions The changes in environmental conditions that have occurred with large-scale deforestation suggest that selective regimes may have shifted and, as a consequence, spatial patterns of intra-specific variation in morphology, behaviour and genes have probably been altered. Although the observed shifts in available environmental conditions resulting from deforestation are greater than those predicted by climate change, the latter will result in novel environments that exceed temperatures in any present-day climates and may lead to biotic attrition unless organisms can adapt to these warmer conditions. Conserving intra-specific diversity over the long term will require considering both how changes in the recent past have influenced contemporary populations and the impact of future environmental change.

Immobilization of glucose oxidase enzyme (GOD) in large pore ordered mesoporous cage-like FDU-1 silica

Large pore ordered mesoporous silica FDU-1 with three-dimensional (3D) face-centered cubic, Fm3m arrangement of rnesopores, was synthesized under strong acid media using B-50-6600 poly(ethylene oxide)-poly(butylene oxide)-poly(ethylene oxide) triblock copolymer (EO(39)BO(47)EO(39)), tetraethyl orthosilicate (TEOS) and trimethyl-benzene (TMB). Large pore FDU-1 silica was obtained by using the following gel composition 1TEOS:0.00735B50-6600:0.00735TMB:6HCl:155H(2)O. The pristine material exhibited a BET specific surface area of 684 m(2) g(-1), total pore volume of 0.89 cm(3) g(-1), external surface area of 49 m(2) g(-1) and microporous volume of 0.09 cm(3) g(-1). The enzyme activity was determined by the Flow Injection Analysis-Chemiluminescence (FIA-CL) method. For GOD immobilized on the FDU-1 silica, GOD supernatant and GOD solution, the FIA-CL results were 9.0, 18.6 and 34.0 U, respectively. The value obtained for the activity of the GOD solution with FIA-CL method is in agreement with the 35 U, obtained by spectrophotometry. (C) 2011 Elsevier B.V. All rights reserved.

Center of mass energy and system-size dependence of photon production at forward rapidity at RHIC

We present the multiplicity and pseudorapidity distributions of photons produced in Au + Au and Cu + Cu collisions at root(s)NN = 62.4 and 200 GeV. The photons are measured in the region -3.7 < eta < -2.3 using the photon Multiplicity detector in the STAR experiment at RHIC. The number of photons produced per average number of participating nucleon pairs increases with the beam energy and is independent of (lie collision centrality. For collisions with similar average numbers of participating nucleons the photon multiplicities are observed to be similar for An + Au and Cu + Cu collisions at a given beam energy. The ratios of the number of charged particles to photons in the measured pseudorapidity range are found to be 1.4 +/- 0.1 and 1.2 +/- 0.1 for root(s)NN = 62.4 and 200 GeV, respectively. The energy dependence of this ratio could reflect varying contributions from baryons to charged particles, while mesons are the dominant contributors to photon production in the given kinematic region. The photon pseudorapidity distributions normalized by average number of participating nucleon pairs, when plotted as a function of eta-Y(beam), are found to follow a longitudinal scaling independent of centrality and colliding ion species at both beam energies. (C) 2009 Elsevier B.V. All rights reserved.

Energy and system size dependence of phi meson production in Cu plus Cu and Au plus Au collisions

We study the beam-energy and system-size dependence of phi meson production (using the hadronic decay mode phi -> K(+) K(-)) by comparing the new results from Cu + Cu collisions and previously reported Au + Au collisions at root s(NN) = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented in this Letter are from mid-rapidity (vertical bar y vertical bar < 0.5) for 0.4 < p(T) < 5 GeV/c. At a given beam energy, the transverse momentum distributions for phi mesons are observed to be similar in yield and shape for Cu + Cu and Au + Au colliding systems with similar average numbers of participating nucleons. The phi meson yields in nucleus-nucleus collisions, normalized by the average number of participating nucleons, are found to be enhanced relative to those from p + p collisions. The enhancement for phi mesons lies between strange hadrons having net strangeness = 1 (K(-) and <(A)over bar>) and net strangeness = 2 (Xi). The enhancement for phi mesons is observed to be higher at root s(NN) = 200 GeV compared to 62.4 GeV. These observations for the produced phi(s (s) over bar) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems. (C) 2009 Elsevier B.V. All rights reserved.

Size and power properties of some tests in the Birnbaum-Saunders regression model

The Birnbaum-Saunders distribution has been used quite effectively to model times to failure for materials subject to fatigue and for modeling lifetime data. In this paper we obtain asymptotic expansions, up to order n(-1/2) and under a sequence of Pitman alternatives, for the non-null distribution functions of the likelihood ratio, Wald, score and gradient test statistics in the Birnbaum-Saunders regression model. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the shape parameter. Monte Carlo simulation is presented in order to compare the finite-sample performance of these tests. We also present two empirical applications. (C) 2010 Elsevier B.V. All rights reserved.

Bias and size effects of price-comparison search engines : theory and experimental evidence

We analyze the impact on consumer prices of the size and bias of price comparison search engines. In the context of a model related to Burdett and Judd (1983) and Varian (1980), we develop and test experimentally several theoretical predictions. The experimental results confirm the model’s predictions regarding the impact of the number of firms, and the type of bias of the search engine, but reject the model’s predictions regarding changes in the size of the index. The explanatory power of an econometric model for the price distributions is significantly improved when variables accounting for risk attitudes are introduced.

Comparison of fracture surface and plane section analysis for ceramic grain size characterisation

A method has been developed to obtain quantitative information about grain size and shape from fractured surfaces of ceramic materials. One elaborated a routine to split intergranular and transgranular grains facets of ceramic fracture surfaces by digital image processing. A commercial ceramic (ALCOA A-16, Al2O3-1.5% of CrO) was used to test the proposed method. Microstructural measurements of grain shape and size taken from fracture surfaces have been compared through descriptive statistics of distributions, with the corresponding measurements from polished and etched surfaces. The agreement between results, with the expected bias on grain size values from fractures, obtained for both types of surfaces allowed to infer that this new technique can be used to extract the relevant microstructural information from fractured surfaces, thus minimising the time consuming steps of sample preparation. (C) 2003 Elsevier Ltd. All rights reserved.