Hommel's procedure in linear time

Hommel's and Hochberg's procedures for familywise error control are both derived as shortcuts in a closed testing procedure with the Simes local test. Hommel's shortcut is exact but takes quadratic time in the number of hypotheses. Hochberg's shortcut takes only linear time after the P‐values are sorted, but is conservative. In this paper, we present an exact shortcut in linear time on sorted P‐values, combining the strengths of both procedures. The novel shortcut also applies to a robust variant of Hommel's procedure that does not require the assumption of the Simes inequality.     

On Method of Moments Estimation in Linear Mixed Effects Models with Measurement Error on Covariates and Response with Application to a Longitudinal Study of Gene-Environment Interaction

We study a linear mixed effects model for longitudinal data, where the response variable and covariates with fixed effects are subject to measurement error. We propose a method of moment estimation that does not require any assumption on the functional forms of the distributions of random effects and other random errors in the model. For a classical measurement error model we apply the instrumental variable approach to ensure identifiability of the parameters. Our methodology, without instrumental variables, can be applied to Berkson measurement errors. Using simulation studies, we investigate the finite sample performances of the estimators and show the impact of measurement error on the covariates and the response on the estimation procedure. The results show that our method performs quite satisfactory, especially for the fixed effects with measurement error (even under misspecification of measurement error model). This method is applied to a real data example of a large birth and child cohort study.     

Statistical tests for identifying differentially expressed genes in time-course microarray experiments

MOTIVATION: Microarray technology allows the monitoring of expression levels for thousands of genes simultaneously. In time-course experiments in which gene expression is monitored over time, we are interested in testing gene expression profiles for different experimental groups. However, no sophisticated analytic methods have yet been proposed to handle time-course experiment data. RESULTS: We propose a statistical test procedure based on the ANOVA model to identify genes that have different gene expression profiles among experimental groups in time-course experiments. Especially, we propose a permutation test which does not require the normality assumption. For this test, we use residuals from the ANOVA model only with time-effects. Using this test, we detect genes that have different gene expression profiles among experimental groups. The proposed model is illustrated using cDNA microarrays of 3840 genes obtained in an experiment to search for changes in gene expression profiles during neuronal differentiation of cortical stem cells.     

Area under the free-response ROC curve (FROC) and a related summary index

Summary .  Free-response assessment of diagnostic systems continues to gain acceptance in areas related to the detection, localization, and classification of one or more "abnormalities" within a subject. A free-response receiver operating characteristic (FROC) curve is a tool for characterizing the performance of a free-response system at all decision thresholds simultaneously. Although the importance of a single index summarizing the entire curve over all decision thresholds is well recognized in ROC analysis (e.g., area under the ROC curve), currently there is no widely accepted summary of a system being evaluated under the FROC paradigm. In this article, we propose a new index of the free-response performance at all decision thresholds simultaneously, and develop a nonparametric method for its analysis. Algebraically, the proposed summary index is the area under the empirical FROC curve penalized for the number of erroneous marks, rewarded for the fraction of detected abnormalities, and adjusted for the effect of the target size (or "acceptance radius"). Geometrically, the proposed index can be interpreted as a measure of average performance superiority over an artificial "guessing" free-response process and it represents an analogy to the area between the ROC curve and the "guessing" or diagonal line. We derive the ideal bootstrap estimator of the variance, which can be used for a resampling-free construction of asymptotic bootstrap confidence intervals and for sample size estimation using standard expressions. The proposed procedure is free from any parametric assumptions and does not require an assumption of independence of observations within a subject. We provide an example with a dataset sampled from a diagnostic imaging study and conduct simulations that demonstrate the appropriateness of the developed procedure for the considered sample sizes and ranges of parameters.     

Asymptotic behaviour of solutions of the diffusive Lotka-Volterra equations

A spatially discrete version of the diffusive Lotka-Volterra equations is considered. Asymptotical spatial homogeneity of solutions of the equations with equilibrium, periodic or zero flux boundary conditions is proved without regard to crowding effects. The proof does not require the assumption of equal diffusion coefficients and the restrictions on the dimension of space and on the initial data, which are necessary in the spatially continuous model.     

Simultaneous confidence intervals from randomization tests with application in testing bioequivalence with multiple endpoints

We propose a method to construct simultaneous confidence intervals for a parameter vector from inverting a series of randomization tests (RT). The randomization tests are facilitated by an efficient multivariate Robbins–Monro procedure that takes the correlation information of all components into account. The estimation method does not require any distributional assumption of the population other than the existence of the second moments. The resulting simultaneous confidence intervals are not necessarily symmetric about the point estimate of the parameter vector but possess the property of equal tails in all dimensions. In particular, we present the constructing the mean vector of one population and the difference between two mean vectors of two populations. Extensive simulation is conducted to show numerical comparison with four methods. We illustrate the application of the proposed method to test bioequivalence with multiple endpoints on some real data.     

Joint models for a primary endpoint and multiple longitudinal covariate processes

Summary .   Joint models are formulated to investigate the association between a primary endpoint and features of multiple longitudinal processes. In particular, the subject-specific random effects in a multivariate linear random-effects model for multiple longitudinal processes are predictors in a generalized linear model for primary endpoints. Li, Zhang, and Davidian (2004, Biometrics 60 , 1–7) proposed an estimation procedure that makes no distributional assumption on the random effects but assumes independent within-subject measurement errors in the longitudinal covariate process. Based on an asymptotic bias analysis, we found that their estimators can be biased when random effects do not fully explain the within-subject correlations among longitudinal covariate measurements. Specifically, the existing procedure is fairly sensitive to the independent measurement error assumption. To overcome this limitation, we propose new estimation procedures that require neither a distributional or covariance structural assumption on covariate random effects nor an independence assumption on within-subject measurement errors. These new procedures are more flexible, readily cover scenarios that have multivariate longitudinal covariate processes, and can be implemented using available software. Through simulations and an analysis of data from a hypertension study, we evaluate and illustrate the numerical performances of the new estimators.     

Isolation of hen's egg yolk very low density lipoproteins by DEAE-cellulose chromatography

A procedure has been developed for the isolation of very low density lipoproteins from hen's egg yolk plasma using DEAE-cellulose chromatography. This procedure is rapid and does not require ultracentrifugation and should, therefore, serve as a useful procedure for use in laboratories where this facility does not exist.     

Random sampling does not exclude spatial dependence: The importance of neutral models for ecological hypothesis testing

Lájer (2007) notes that, to investigate phytosociological and ecological relationships, many authors apply traditional inferential tests to sets of relevés obtained by non-random methods. Unfortunately, this procedure does not provide reliable support for hypothesis testing because non-random sampling violates the assumptions of independence required by many parametric inferential tests. Instead, a random sampling scheme is recommended. Nonetheless, random sampling will not eliminate spatial autocorrelation. For instance, a classical law of geography holds that everything in a piece of (biotic) space is interrelated, but near objects are more related than distant ones. Because most ecological processes that shape community structure and species coexistence are spatially explicit, spatial autocorrelation is a vital part of almost all ecological data. This means that, independently from the underlying sampling design, ecological data are generally spatially autocorrelated, violating the assumption of independence that is generally required by traditional inferential tests. To overcome this drawback, randomization tests may be used. Such tests evaluate statistical significance based on empirical distributions generated from the sample and do not necessarily require data independence. However, as concerns hypothesis testing, randomization tests are not the universal remedy for ecologists, because the choice of inadequate null models can have significant effects on the ecological hypotheses tested. In this paper, I emphasize the need of developing null models for which the statistical assumptions match the underlying biological mechanisms.     

Dynamic ion-ion and water-ion interactions in ion channels.

The dynamic interactions among ions and water molecules in ion channels are treated based on an assumption that ions at binding sites can be knocked off by both transient entering ions and local water molecules. The theory, when applied to a single-site model K+ channel, provides solutions for super- and subsaturations, flux-ratio exponent (n') greater than 1, osmotic streaming current, activity-dependent reversal potentials, and anomalous mole-fraction behavior. The analysis predicts that: (a) the saturation may but, in general, does not follow the Michaelis-Menten relation; (b) streaming current results from imbalanced water-ion knock-off interactions; (c) n' greater than 1 even for single-site channels, but it is unlikely to exceed 1.4 unless the pore is occupied by one or more ion(s); (d) in the calculation involving two permeant ion species with similar radii, the heavier ions show higher affinity; the ion-ion knock-off dissociation from the site is more effective when two interacting ions are identical. Therefore, the "multi-ion behaviors" found in most ion channels are the consequences of dynamic ion-ion and water-ion interactions. The presence of these interactions does not require two or more binding sites in channels.     

On-line estimation of the maximum specific growth rate of nitrifiers in activated sludge systems.

The on-line estimation of the maximum specific growth rate of autotrophic biomass is addressed in this article. A general nitrification process model, which is valid for any realistic flow pattern, is used to develop the estimation algorithm. Depending on the measurements available, two estimation equations are derived. While both require measuring the nitrification activity of the activated sludge, one requires the additional measurement of the nitrifiable nitrogen concentrations at the two ends of the bioreactor, and the other requires the nitrate nitrogen concentrations at the same locations. The algorithm also requires some stoichiometric and kinetic parameters. However, sensitivity analysis shows that the estimate is insensitive to the parameters other than the autotrophic decay rate. Compared to the existing algorithms, the algorithm developed in this article does not rely on the assumption of ideal flow pattern in the plant and does not require an error-prone estimate of the autotrophic biomass concentration. Experimental and simulation studies show that the algorithm performs well and is robust to influent variations and accidental sludge losses.     

Too many mutants with multiple mutations

It has recently become clear that the classical notion of the random nature of mutation does not hold for the distribution of mutations among genes: most collections of mutants contain more isolates with two or more mutations than predicted by the mutant frequency on the assumption of a random distribution of mutations. Excesses of multiples are seen in a wide range of organisms, including riboviruses, DNA viruses, prokaryotes, yeasts, and higher eukaryotic cell lines and tissues. In addition, such excesses are produced by DNA polymerases in vitro. These "multiples" appear to be generated by transient, localized hypermutation rather than by heritable mutator mutations. The components of multiples are sometimes scattered at random and sometimes display an excess of smaller distances between mutations. As yet, almost nothing is known about the mechanisms that generate multiples, but such mutations have the capacity to accelerate those evolutionary pathways that require multiple mutations where the individual mutations are neutral or deleterious. Examples that impinge on human health may include carcinogenesis and the adaptation of microbial pathogens as they move between individual hosts.     

Non-linear benefits and the evolution of eusociality in the hymenoptera

A general model for examining the evolution of social behavior is developed which does not require that benefits received be linear functions of the number of social donors encountered. The subsocial route for the evolution of eusociality in haplodiploid organisms is then examined within the context of this model Non-linearities render conditions for frequency independent fixation or loss of sister helping alleles more stringent than expected from models based on the assumption of linear benefits. In particular, both stable polymorphisms and frequency dependent selective thresholds for sister helping behaviour may commonly obtain.     

Electron transfer between globular proteins. Dependence of the rate of transfer on distance

Based on the assumption that electron transfer between globular proteins occurs by a collective excitation of polaron type, the dependence of the rate of this process on the distance between the donor and acceptor centers with regard to their detailed electron structure was calculated. The electron structure of the heme was calculated by the quantum-chemical MNDO-PM3 method. The results were compared with experimental data on interprotein and intraglobular electron transfer. It is shown that, in the framework of this model, the electron transfer is not exponential and does not require a particular transfer pathway since the whole protein macromolecule is involved in the formation of the electron excited state.     

Preparation of Inositol Triphosphate from Brain: GLC of Trimethylsilyl Derivative

A simple method for the preparation of myo-inositol triphosphate from brain is described that does not require the prior isolation of phosphatidylinositol diphosphate (triphosphoinositide). A gas chromatographic procedure for detection of inositol triphosphate is reported.     

Disparity index: a simple statistic to measure and test the homogeneity of substitution patterns between molecular sequences.

A common assumption in comparative sequence analysis is that the sequences have evolved with the same pattern of nucleotide substitution (homogeneity of the evolutionary process). Violation of this assumption is known to adversely impact the accuracy of phylogenetic inference and tests of evolutionary hypotheses. Here we propose a disparity index, ID, which measures the observed difference in evolutionary patterns for a pair of sequences. On the basis of this index, we have developed a Monte Carlo procedure to test the homogeneity of the observed patterns. This test does not require a priori knowledge of the pattern of substitutions, extent of rate heterogeneity among sites, or the evolutionary relationship among sequences. Computer simulations show that the ID-test is more powerful than the commonly used chi2-test under a variety of biologically realistic models of sequence evolution. An application of this test in an analysis of 3789 pairs of orthologous human and mouse protein-coding genes reveals that the observed evolutionary patterns in neutral sites are not homogeneous in 41% of the genes, apparently due to shifts in G + C content. Thus, the proposed test can be used as a diagnostic tool to identify genes and lineages that have evolved with substantially different evolutionary processes as reflected in the observed patterns of change. Identification of such genes and lineages is an important early step in comparative genomics and molecular phylogenetic studies to discover evolutionary processes that have shaped organismal genomes.     

A rapid and inexpensive method for preparing E. coli plasmid-DNA

A simple, rapid and inexpensive scaled up miniprep procedure for preparing pure E. coli plasmid DNA is described. Bacterial cells were subjected to the boiling procedure and high molecular weight RNA was removed by LiCl-precipitation. Residual RNA and proteins were removed by subsequent treatment with RNase A and proteinase K/SDS respectively, followed by Sephadex G-50 and Sepharose 6B-Cl chromatography. The average yield from a 100 ml over-night bacterial suspension was 100 to 150 micrograms for pBR-322 DNA, and 250-500 micrograms for SP-6 derived recombinant plasmids. In addition, the described "scaled up" preparation does not require CsCl-ethidium bromide centrifugation; pure plasmid DNA can be prepared within 1 to 2 days.     

Healthcare and the hospital chaplain

Many chaplains and most chaplaincy programs in the United States--with encouragement from their accrediting organization, the Association for Clinical Pastoral Education (ACPE)--have begun to assume a more proactive stance toward patients, healthcare professionals, and healthcare facilities. Some chaplains and chaplaincy programs have begun to engage in activities that have ranged from initiating conversations with and perusing the medical records of patients who have not requested their services to proposing that they be permitted to do "spiritual assessments" on patients--in some instances whether these patients have been explicitly informed and have agreed to this beforehand. Moreover, many chaplains and chaplaincy programs have begun to assume that chaplains are full-fledged members of the healthcare team, complete with access to patients' medical records both to gather information and to make notations of their own. It would appear that such novel activities are being justified by a questionable set of claims and assumptions that includes: (1) the claim that chaplains have a spiritual--as opposed to purely religious--expertise that entitles them to interact with patients and/or significant others (even those who have not requested a chaplain)--presumably without in the least compromising patient autonomy or the confidentiality of the patient/healthcare professional relationship; (2) the assumption that the terms "spirituality" and "religiosity" mutually entail one another; (3) the claim that the Joint Commission on the Accreditation of Healthcare Organizations (JCAHO) mandates "spiritual assessments" (which it does not); (4) the assumption that chaplains are full-fledged members of the healthcare team; and (5) the claim that chaplains must, therefore, be permitted access to patients and patients' medical records both to gather information and to make notations of their own. We consider such claims and assumptions disquieting, and suggest that it is high time we revisit the terms "chaplaincy," "healthcare professional," and "member of the healthcare team" in reassessing what our professional commitments to respect and protect the bio-psycho-social integrity of patients require.     

Estimation of parity progression ratios from the truncated distribution of closed and open birth intervals.

A procedure to estimate parity progression ratios in a population from the truncated distribution of open and closed birth intervals is presented. The approach is quite simple in computation and data needs. It does not require any separate data on age at last birth to the women of completed fertility as in earlier methods. The procedure is illustrated with an observed set of data.     

Expected Error Rates for Logistic Regression Versus Normal Discriminant Analysis

The expected error rates associated with using the allocation rule based on logistic regression are derived in the context of two multivariate normal populations with a common covariance matrix and compared with the corresponding error rates of the classical rule based on this normality assumption. It is shown in terms of the actual sizes of the asymptotic expected error rates that the performance of the logistic procedure does not fall far short of the normality based method, even for widely separated populations. This latter result is not obvious from previously available work on the asymptotic relative efficiency of the logistic procedure.