Estimation in independent observer line transect surveys for clustered populations

This paper introduces a framework for animal abundance estimation in independent observer line transect surveys of clustered populations. The framework generalizes an approach given in Chen (1999, Environmental and Ecological Statistics 6, in press) to accommodate heterogeneity in detection caused by cluster size and other covariates. Both parametric and nonparametric estimators for the local effective search widths, given the covariates, can be derived from the framework. A nonparametric estimator based on conditional kernel density estimation is proposed and studied owing to its flexibility in modeling the detection functions. A real data set on harbor porpoise in the North Sea is analyzed.     

Comparison of parameters from rhythmometric models with multiple components on hybrid data

Population multiple components is a statistical tool useful for the analysis of time-dependent hybrid data. With a small number of parameters, it is possible to model and to predict the periodic behavior of a population. In this article, we propose two methods to compare among populations rhythmometric parameters obtained by multiple component analysis. The first is a parametric method based in the usual statistical techniques for comparison of mean vectors in multivariate normal populations. The method, through MANOVA analysis, allows comparison of the MESOR and amplitude-acrophase pair of each component among two or more populations. The second is a nonparametric method, based in bootstrap techniques, to compare parameters from two populations. This test allows one to compare the MESOR, the amplitude, and the acrophase of each fitted component, as well as the global amplitude, orthophase, and bathyphase estimated when all fitted components are harmonics of a fundamental period. The idea is to calculate a confidence interval for the difference of the parameters of interest. If this interval does not contain zero, it can be concluded that the parameters from the two models are different with high probability. An estimation of p-value for the corresponding test can also be calculated. Both methods are illustrated with an example, based on clinical data. The nonparametric test can also be applied to paired data, a special situation of great interest in practice. By the use of similar bootstrap techniques, we illustrate how to construct confidence intervals for any rhythmometric parameter estimated from population multiple components models, including the orthophase, bathyphase, and global amplitude. These tests for comparison of parameters among populations are a needed tool when modeling the nonsinusoidal rhythmic behavior of hybrid data by population multiple component analysis.     

Control of kernel weight and kernel water relations by post-flowering source-sink ratio in maize

The maize (Zea mays L.) kernel undergoes large changes in water content during its development. Whether such changes regulate the pattern of kernel development or are simply a consequence of it has not yet been established because other factors, such as assimilate supply, can also affect the rate and duration of kernel growth. This study was conducted to determine whether variation in kernel weight (KW) in response to source-sink treatments is mediated by a change in kernel water relations. Two hybrids were sown at three stand densities (one, eight and 18 plants m-2), and kernel numbers were restricted to control the post-flowering source-sink ratio within each stand density. Kernel development and water relations [water content, water potential (psiw), osmotic potential (psis) and turgor] were monitored throughout grain filling. Final KW varied from 253 to 372 mg per kernel in response to source-sink treatments. For both genotypes, variation in KW was a result of a change in kernel growth rate (r2 = 0.91; P < 0.001) and not in the duration of kernel filling. Final KW was closely correlated with maximum kernel water content (r2 = 0.94; P < 0.001) achieved during rapid dry matter accumulation. However, variation in KW was not reflected in kernel water status parameters (psiw, psis or turgor), which remained fairly stable across treatments. These results indicate that maximum water content provides an easily quantifiable measure of kernel sink capacity in maize. Kernel water status parameters may affect the duration of grain filling, but have no discernible impact on kernel growth rate.     

Semiparametric estimation of tag loss and reporting rates for tag-recovery experiments using exact time-at-liberty data

We present a semiparametric likelihood approach to estimating reporting rates and tag-loss rates from the tags returned from capture-recapture studies. Such studies are commonly used to estimate critical population parameters. Tag loss rates are estimated using double-tagged animals, while reporting rates are estimated using information from high-reward tags. A likelihood function is constructed based on the conditional distribution of the type of tag returned (low or high reward, single or double tag), given that a tag has been returned. This involves many sparse 5 x 1 tag-return contingency tables, and choosing a good functional form for the tag loss rate is difficult with such data. We model tag-loss rates using monotone-smoothing splines, and use these nonparametric estimates to diagnose the parametric form of the tag-loss rate. The nonparametric methods can also be used directly to model tag-loss rates.     

Estimating the accuracy of polymerase chain reaction-based tests using endpoint dilution

Polymerase chain reaction (PCR)-based tests for various microorganisms or target DNA sequences are generally acknowledged to be highly "sensitive," yet the concept of sensitivity is ill-defined in the literature on these tests. We propose that sensitivity should be expressed as a function of the number of target DNA molecules in the sample (or specificity, when the target number is 0). However, estimating this "sensitivity curve" is problematic, since it is difficult to construct samples with a fixed number of targets. Nonetheless, using serially diluted replicate aliquots of a known concentration of the target DNA sequence, we show that it is possible to disentangle random variations in the number of target DNA molecules from the underlying test sensitivity. We develop parametric, nonparametric, and semiparametric (spline-based) models for the sensitivity curve. The methods are compared on a new test for M. genitalium.     

Contrasting patterns of radiation in African and Australian Restionaceae

The floras of the Mediterranean-climate areas of southern Africa and southwestern Australia are remarkably species rich. Because the two areas are at similar latitudes and in similar positions on their respective continents, they have probably had similar Cenozoic climatic histories. Here we test the prediction that the evolution of the species richness in the two areas followed a similar temporal progression by comparing the rates of lineage accumulation for African and Australian Restionaceae. Restionaceae (Poales) are typical and often dominant elements in the fynbos vegetation of the Cape Floristic Region of southern Africa and the kwongan vegetation of the Southwestern Floristic Province of Western Australia. The phylogeny of the family was estimated from combined datasets for rbcL and trnL-F sequences and a large morphological dataset; these datasets are largely congruent. The monophyly of Restionaceae is supported and a basal division into an African clade (approximately 350 species) and an Australian clade (146 species) corroborated. There is also support for a futher subdivision of these two large sister-clades, but the terminal resolution within the African clade is very weak. Fossil pollen records provided a minimum age of the common ancestor of Australian and African Restionaceae as 64-71 million years ago, and this date was used to calibrate a molecular clock. A molecular clock was rejected by a likelihood ratio test; therefore, rate changes between the lineages were smoothed using nonparametric rate smoothing. The rate-corrected ages were used to construct a plot of lineages through time. During the Palaeogene the Australian lineage diversity increased consistent with the predictions of the constant birthrate model, while the African lineage diversity showed a dramatic increase in diversification rate in the Miocene. Incomplete sampling obscures the patterns in the Neogene, but extending the trends to the modern extant diversity suggests that this acceleration in the speciation rate continued in the African clade, whereas the Australian clade retained a constant diversification rate. The substantial morphological and anatomical similarity between the African and Australian Restionaceae appear to preclude morphological innovations as possible explanations for the intercontinental differences. Most likely these differences are due to the greater geographical extent and ecological variation in temperate Australia than temperate Africa, which might have provided refugia for basal Restionaceae lineages, whereas the more mountainous terrain of southern Africa might have provided the selective regimes for a more rapid, recent speciation.     

Statistical inference for simultaneous clustering of gene expression data

Current methods for analysis of gene expression data are mostly based on clustering and classification of either genes or samples. We offer support for the idea that more complex patterns can be identified in the data if genes and samples are considered simultaneously. We formalize the approach and propose a statistical framework for two-way clustering. A simultaneous clustering parameter is defined as a function theta=Phi(P) of the true data generating distribution P, and an estimate is obtained by applying this function to the empirical distribution P(n). We illustrate that a wide range of clustering procedures, including generalized hierarchical methods, can be defined as parameters which are compositions of individual mappings for clustering patients and genes. This framework allows one to assess classical properties of clustering methods, such as consistency, and to formally study statistical inference regarding the clustering parameter. We present results of simulations designed to assess the asymptotic validity of different bootstrap methods for estimating the distribution of Phi(P(n)). The method is illustrated on a publicly available data set.     

Membrane oxidation assay--a novel lipoxygenase-based evaluation of membrane oxidizability

A membrane oxidation assay is presented which uses isolated erythrocyte membranes ("ghosts") and lipoxygenase as a selective catalyst for the transfer of oxygen to cis - cis -1,4-pentadiene-moieties. The latter are, for instance, present in linoleic and arachidonic acids, both of which are integral parts of membranes. These non-conjugated double bonds represent energetically favorable sites for oxidative attack and therefore, may be rearranged and partially consumed during oxidative stress. Consequently, the measurement of oxygen consumption in the course of the lipoxygenase-mediated oxidation provides a tool for the quick and reliable determination of such double bonds. Significant inter-individual differences have been noted in 11 subjects, which also correlate to the total radical antioxidant parameter (TRAP) values obtained. This assay will be helpful in the assessment of oxidizable structures in erythrocyte membranes that may be diminished as a consequence of oxidative damage suffered by an individual. In conclusion, a simple and rapid assay for the assessment of the oxidizability of erythrocyte membranes is presented complementing the TRAP assay for plasma antioxidative status.     

Three Types of Membrane Excitations in the Marine Diatom Coscinodiscus wailesii

Three types of electrical excitation have been investigated in the marine diatom Coscinodiscus wailesii. I: Depolarization-triggered, transient Cl⁻ conductance, G _Cl (t), followed by a transient, voltage-gated K⁺ conductance, G _K , with an active state a and two inactive states i ₁ and i ₂ in series (a-i ₁-i ₂). II: Similar G _Cl (t) as in Type-I but triggered by hyperpolarization; a subsequent increase of G _K in this type is indicated but not analyzed in detail. III: Hyperpolarization-induced transient of a voltage-gated activity of an electrogenic pump (i ₂-a-i ₂), followed by G _Cl (t) as in Type-II excitations. Type-III with pump gating is novel as such. G _Cl (t) in all types seems to reflect the mechanism of InsP⁻ ₃ and Ca²⁺-mediated G _Cl (t) in the action potential in Chara (Biskup et al., 1999). The nonlinear current-voltage-time relationships of Type-I and Type-III excitations have been recorded under voltage-clamp using single saw-tooth command voltages (voltage range: −200 to +50 mV, typical slope: ±1 Vs⁻¹). Fits of the corresponding models to the experimental data provided numerical values of the model parameters. The statistical significance of these solutions is investigated. We suggest that the original function of electrical excitability of biological membranes is related to osmoregulation which has persisted through evolution in plants, whereas the familiar and osmotically neutral action potentials in animals have evolved later towards the novel function of rapid transmission of information over long distances. Received: 2 December 1999/Revised: 3 March 2000