A Nonparametric Test for Differences in the Dispersion of Dependent Samples

A new statistic Δ to test the hypothesis of a difference in the dispersion of two dependent samples of ordinal data quality is proposed. It draws on the idea of rank assignment originally forwarded by Siegel and Tukey (1960). No exact probability levels can be given for this statistic for the time being, but it is shown that the statistic is linearly related to the so-called Hotelling-Pabst statistic D, and that one can use exact tables of the latter as a substitute in the statistical decision process with small samples. For larger samples, an approximation of Δ to the standard normal distribution is given. The problem of tied observations is not sufficiently solved yet. A conservative procedure of rank assignment is proposed as long as the exact distribution of Δ in the presence of ties is unknown.     

Nonparametric All‐Pairs Multiple Comparisons

Nonparametric all‐pairs multiple comparisons based on pairwise rankings can be performed in the one‐way design with the Steel‐Dwass procedure. To apply this test, Wilcoxon's rank sum statistic is calculated for all pairs of groups; the maximum of the rank sums is the test statistic. We provide exact calculations of the asymptotic critical values (and P‐values, respectively) even for unbalanced designs. We recommend this asymptotic method whenever large sample sizes are present. For small sample sizes we recommend the use of the new statistic according to Baumgartner , Weiss , and Schindler (1998, Biometrics 54 , 1129–1135) instead of Wilcoxon's rank sum for the multiple comparisons. We show that the resultant procedure can be less conservative and, according to simulation results, more powerful than the original Steel‐Dwass procedure. We illustrate the methods with a practical data set.     

Stratified Fisher's exact test and its sample size calculation

Chi‐squared test has been a popular approach to the analysis of a 2 × 2 table when the sample sizes for the four cells are large. When the large sample assumption does not hold, however, we need an exact testing method such as Fisher's test. When the study population is heterogeneous, we often partition the subjects into multiple strata, so that each stratum consists of homogeneous subjects and hence the stratified analysis has an improved testing power. While Mantel–Haenszel test has been widely used as an extension of the chi‐squared test to test on stratified 2 × 2 tables with a large‐sample approximation, we have been lacking an extension of Fisher's test for stratified exact testing. In this paper, we discuss an exact testing method for stratified 2 × 2 tables that is simplified to the standard Fisher's test in single 2 × 2 table cases, and propose its sample size calculation method that can be useful for designing a study with rare cell frequencies.     

How feeding competition determines female chimpanzee gregariousness and ranging in the Taï National Park, Côte d'Ivoire

Socioecological theory suggests that feeding competition shapes female social relationships. Chimpanzees (Pan troglodytes) live in fission–fusion societies that allow them to react flexibly to increased feeding competition by forming smaller foraging parties when food is scarce. In chimpanzees at Gombe and Kibale, female dominance rank can crucially influence feeding competition and reproductive success as high‐ranking females monopolize core areas of relatively high quality, are more gregarious, and have higher body mass and reproductive success than low‐ranking females. Chimpanzee females in Taï National Park do not monopolize core areas; they use the entire territory as do the males of their community and are highly gregarious. Although female chimpanzees in Taï generally exhibit a linear dominance hierarchy benefits of high rank are currently not well understood. We used a multivariate analysis of long‐term data from two Taï chimpanzee communities to test whether high‐ranking females (1) increase gregariousness and (2) minimize their travel costs. We found that high‐ranking females were more gregarious than low‐rankers only when food was scarce. During periods of food scarcity, high rank allowed females to enjoy benefits of gregariousness, while low‐ranking females strongly decreased their gregariousness. High‐ranking females traveled more than low‐ranking females, suggesting that low‐rankers might follow a strategy to minimize energy expenditure. Our results suggest that, in contrast to other chimpanzee populations and depending on the prevailing ecological conditions, female chimpanzees at Taï respond differently to varying levels of feeding competition. Care needs to be taken before generalizing results found in any one chimpanzee population to the species level. Am. J. Primatol. 73:305–313, 2011. © 2010 Wiley‐Liss, Inc.     

Confidence Interval for Rate Ratio in a 2 × 2 Table with Structural Zero: An Application in Assessing False-Negative Rate Ratio When Combining Two Diagnostic Tests

Summary . In this article, we consider problems with correlated data that can be summarized in a 2 × 2 table with structural zero in one of the off‐diagonal cells. Data of this kind sometimes appear in infectious disease studies and two‐step procedure studies. Lui (1998, Biometrics 54, 706–711) considered confidence interval estimation of rate ratio based on Fieller‐type, Wald‐type, and logarithmic transformation statistics. We reexamine the same problem under the context of confidence interval construction on false‐negative rate ratio in diagnostic performance when combining two diagnostic tests. We propose a score statistic for testing the null hypothesis of nonunity false‐negative rate ratio. Score test–based confidence interval construction for false‐negative rate ratio will also be discussed. Simulation studies are conducted to compare the performance of the new derived score test statistic and existing statistics for small to moderate sample sizes. In terms of confidence interval construction, our asymptotic score test–based confidence interval estimator possesses significantly shorter expected width with coverage probability being close to the anticipated confidence level. In terms of hypothesis testing, our asymptotic score test procedure has actual type I error rate close to the pre‐assigned nominal level. We illustrate our methodologies with real examples from a clinical laboratory study and a cancer study.     

Non-inferiority of new procedure to standard procedure in stratified matched-pair design

We consider the statistical testing for non-inferiority of a new treatment compared with the standard one under matched-pair setting in a stratified study or in several trials. A non-inferiority test based on the efficient scores and a Mantel-Haenszel (M-H) like procedure with restricted maximum likelihood estimators (RMLEs) of nuisance parameters and their corresponding sample size formulae are presented. We evaluate the above tests and the M-H type Wald test in level and power. The stratified score test is conservative and provides the best power. The M-H like procedure with RMLEs gives an accurate level. However, the Wald test is anti-conservative and we suggest caution when it is used. The unstratified score test is not biased but it is less powerful than the stratified score test when base-line probabilities related to strata are not the same. This investigation shows that the stratified score test possesses optimum statistical properties in testing non-inferiority. A common difference between two proportions across strata is the basic assumption of the stratified tests, we present appropriate tests to validate the assumption and related remarks.     

Incorporation of clustering effects for the Wilcoxon rank sum test: a large-sample approach

The Wilcoxon rank sum test is frequently used in statistical practice for the comparison of measures of location when the underlying distributions are far from normal or not known in advance. An assumption of the ordinary rank sum test is that individual sampling units are independent. In many ophthalmologic clinical trials, the Early Treatment for Diabetic Retinopathy Scale (ETDRS) is a principal endpoint used for measuring the level of diabetic retinopathy. This is an ordinal scale, and it is natural to consider the Wilcoxon rank sum test for the comparison of the level of diabetic retinopathy between treatment groups. However, under this design, unlike the usual Wilcoxon rank sum test, the subject is the unit of randomization, but the eye is the unit of analysis. Furthermore, a person will tend to have different, but correlated, ETDRS scores for fellow eyes. Thus, we propose a correction to the variance of the Wilcoxon rank sum statistic that accounts for clustering effects and that can be used for both balanced (same number of subunits per cluster) or unbalanced (different number of subunits per cluster) data, both in the presence or absence of ties, with p-value adjusted accordingly. In this article, we present large-sample theory and simulation results for this test procedure and apply it to diabetic retinopathy data from type I diabetics in the Sorbinil Retinopathy Trial.     

A Test for Homogeneity of Odds Ratios in Ordered 2 × 2 Tables

Consider K ordered 2 × 2 contingency tables. A new test of the null hypothesis that the odds ratios of these tables are equal vs the alternative hypothesis that the odds ratios are nondecreasing, is recommended. The test is exact (non‐asymptotic), is easily carried out (software is available), and has other favorable properties. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Social and menstrual cycle phase influences on the behavior of group-housed Cebus apella

This study examined sexual and other social behavior in relation to menstrual cycle phase in four mixed-sex social groups of Cebus apella. Groups consisted of two adult males and either six or nine adult females. Menstrual cycles of high (rank 1–2) and low (rank 3–7) dominance-rank females from each group were monitored via vaginal swabs, and correlated with data collected from ongoing behavioral observations. Only cycles bounded by positive detection of menstrual blood were included in this analysis (n = 15 females, 182 cycles; mean ± s.d., cycle length = 20.8 ± 1.2 days). Rates of copulation and female solicitation of males varied significantly with cycle phase, with highest rates at midcycle. While total rates of solicitation and copulation did not vary with female dominance rank, copulation rates with the dominant male were significantly greater for high ranking females than for low ranking ones. Variance observed in affiliative and agonistic behaviors, including those with males, was attributable to female rank rather than cycle phase. Females and males were also observed attempting to interfere with copulations of lower ranking same-sex individuals. Although further study, particularly of wild populations, is needed, these results indicate that female-female reproductive competition should be included as a component of the capuchin breeding system. As in other primate species, both social and hormone-related factors influence the sexual behavior of Cebus apella. © 1995 Wiley-Liss, Inc.     

Activation of pre-synaptic M-type K+ channels inhibits [3H]d-aspartate release by reducing Ca2+ entry through P/Q-type voltage-gated Ca2+channels

In this study, the functional consequences of the pharmacological modulation of the M‐current (I_KM) on cytoplasmic Ca²⁺ intracellular Ca²⁺concentration ([Ca²⁺]_i) changes and excitatory neurotransmitter release triggered by various stimuli from isolated rat cortical synaptosomes have been investigated. K_v7.2 immunoreactivity was identified in pre‐synaptic elements in cortical slices and isolated glutamatergic cortical synaptosomes. In cerebrocortical synaptosomes exposed to 20 mM [K⁺]_e, the I_KM activator retigabine (RT, 10 μM) inhibited [³H]d ‐aspartate ([³H]d ‐Asp) release and caused membrane hyperpolarization; both these effects were prevented by the I_KM blocker XE‐991 (20 μM). The I_KM activators RT (0.1–30 μM), flupirtine (10 μM) and BMS‐204352 (10 μM) inhibited 20 mM [K⁺]_e‐induced synaptosomal [Ca²⁺]_i increases; XE‐991 (20 μM) abolished RT‐induced inhibition of depolarization‐triggered [Ca²⁺]_i transients. The P/Q‐type voltage‐sensitive Ca²⁺channel (VSCC) blocker ω‐agatoxin IVA prevented RT‐induced inhibition of depolarization‐induced [Ca²⁺]_i increase and [³H]d ‐Asp release, whereas the N‐type blocker ω‐conotoxin GVIA failed to do so. Finally, 10 μM RT did not modify the increase of [Ca²⁺]_i and the resulting enhancement of [³H]d ‐Asp release induced by [Ca²⁺]_i mobilization from intracellular stores, or by store‐operated Ca²⁺channel activation. Collectively, the present data reveal that the pharmacological activation of I_KM regulates depolarization‐induced [³H]d ‐Asp release from cerebrocortical synaptosomes by selectively controlling the changes of [Ca²⁺]_i occurring through P/Q‐type VSCCs.     

A New Nonparametric Approach to the Comparison of K Independent Samples of Response Curves

Two simultaneous distribution free test procedures to test the homogeneity of K samples of response curves measured at T points of time each are proposed which react upon the shape of the response curve. One procedure is based on the comparison of K T-tupels of FRIEDMAN-rank sums, the other on the comparison of K T × T-incidence tables, which first were introduced by KANNEMANN (1976). The application of the two procedures is illustrated by a numerical example. In the discussion the tests are compared with each other and with concurring tests.     

Characterization of 2,4-Dinitrotoluene Dioxygenase from Recombinant Esherichia coli Strain PFJS39: Its Direct Interaction with Vitreoscilla Hemoglobin

Burkholderia dinitrotoluene (DNT) dioxygenase in this study (from recombinant Esherichia coli strain PFJS39) is probably a multicomponent enzyme system that oxidizes 2,4-dinitrotoluene (DNT) to 4-methyl-5-nitrocatechol (MNC). DNT dioxygenase was purified by a four-step procedure that utilized consecutive gel filtration chromatography and a nondenaturing gel system. The purified enzyme oxidized DNT only in the presence of NADH and its yield increased by lipase pretreatment of crude cytosol. An estimated molecular weight of 100,000 was obtained by gel filtration. Polyacrylamide gel electrophoresis (PAGE) in the presence of sodium dodecyl sulfate (SDS) revealed the presence of three subunits for the samples from consecutive gel filtration chromatography and nondenaturing PAGE. Their molecular weights were 52,000–71,000, 23,000–25,500, and 12,000–16,500. These results suggest that DNT dioxygenase exists as a heterotrimer. The K _M of DNT dioxygenase for O₂ is 50 μ M, consistent with inhibition results of DNT dioxygenase by Vitreoscilla hemoglobin (its K _M for O₂ is 7 μ M). The K _M for DNT is 180 μ M. The purified enzyme is relatively stable below 40°C, retains activity over a broad pH range, and is stimulated by several cofactors in addition to NADH.     

Order‐restricted Scores Test for the Evaluation of Population‐based Case–control Studies when the Genetic Model is Unknown

The Cochran–Armitage (CA) linear trend test for proportions is often used for genotype‐based analysis of candidate gene association. Depending on the underlying genetic mode of inheritance, the use of model‐specific scores maximises the power. Commonly, the underlying genetic model, i.e. additive, dominant or recessive mode of inheritance, is a priori unknown. Association studies are commonly analysed using permutation tests, where both inference and identification of the underlying mode of inheritance are important. Especially interesting are tests for case–control studies, defined by a maximum over a series of standardised CA tests, because such a procedure has power under all three genetic models. We reformulate the test problem and propose a conditional maximum test of scores‐specific linear‐by‐linear association tests. For maximum‐type, sum and quadratic test statistics the asymptotic expectation and covariance can be derived in a closed form and the limiting distribution is known. Both the limiting distribution and approximations of the exact conditional distribution can easily be computed using standard software packages. In addition to these technical advances, we extend the area of application to stratified designs, studies involving more than two groups and the simultaneous analysis of multiple loci by means of multiplicity‐adjusted p‐values for the underlying multiple CA trend tests. The new test is applied to reanalyse a study investigating genetic components of different subtypes of psoriasis. A new and flexible inference tool for association studies is available both theoretically as well as practically since already available software packages can be easily used to implement the suggested test procedures.     

Tolerance of Australian tropical and subtropical Acacias to acid soil

Two glasshouse experiments were conducted to refine screening procedures and to rank 36 native Acacia species for acid-soil tolerance. The first experiment showed that nitrogen fixation was more sensitive to soil acidity than plant growth per se. This suggests that, in screening experiments where acid soil tolerant rhizobia are not available, inorganic nitrogen should be supplied in determining the species' intrinsic ability to tolerate soil acidity and to avoid confounding sensitivity to acid soil with ineffectiveness of Rhizobium/Bradyrhizobium. The second experiment enabled us to rank 36 species into 4 classes of differing tolerance to acid soil. This experiment also indicated that the ranking order of species differed with the soil type, suggesting the need to use more than one soil type in screening experiments in order to obtain a better and more reliable ranking of species for tolerance to acid soil.     

Rank Tests for Complete Block Designs

In this article a general univariate K-sample rank test for complete block designs with proportional cell frequencies is derived. It is shown that the test statistic has under H₀ and for arbitrary scores asymptotically a X²-distribution with K — 1 degrees of freedom. Special cases of this test are the Kruskal-Wallis test and the Friedman test. The test is compared with the Benard-van-Elteren test, the Mack-Skillings test and a test proposed by Downton. Finally the application of the test is illustrated by two examples.     

Atp1a3‐deficient heterozygous mice show lower rank in the hierarchy and altered social behavior

Atp1a3 is the Na‐pump alpha3 subunit gene expressed mainly in neurons of the brain. Atp1a3‐deficient heterozygous mice (Atp1a3^+/?) show altered neurotransmission and deficits of motor function after stress loading. To understand the function of Atp1a3 in a social hierarchy, we evaluated social behaviors (social interaction, aggression, social approach and social dominance) of Atp1a3^+/? and compared the rank and hierarchy structure between Atp1a3^+/? and wild‐type mice within a housing cage using the round‐robin tube test and barbering observations. Formation of a hierarchy decreases social conflict and promote social stability within the group. The hierarchical rank is a reflection of social dominance within a cage, which is heritable and can be regulated by specific genes in mice. Here we report: (1) The degree of social interaction but not aggression was lower in Atp1a3^+/? than wild‐type mice, and Atp1a3^+/? approached Atp1a3^+/? mice more frequently than wild type. (2) The frequency of barbering was lower in the Atp1a3^+/? group than in the wild‐type group, while no difference was observed in the mixed‐genotype housing condition. (3) Hierarchy formation was not different between Atp1a3^+/? and wild type. (4) Atp1a3^+/? showed a lower rank in the mixed‐genotype housing condition than that in the wild type, indicating that Atp1a3 regulates social dominance. In sum, Atp1a3^+/? showed unique social behavior characteristics of lower social interaction and preference to approach the same genotype mice and a lower ranking in the hierarchy.     

Role of Human Nucleoside Transporters in the Uptake and Cytotoxicity of Azacitidine and Decitabine

The nucleoside analogs 5-azacytidine (azacitidine) and 5-aza-2′-deoxycytidine (decitabine) are active against acute myeloid leukemia and myelodysplastic syndromes. Cellular transport across membranes is crucial for uptake of these highly polar hydrophilic molecules. We assessed the ability of azacitidine, decitabine, and, for comparison, gemcitabine, to interact with human nucleoside transporters (hNTs) in Saccharomyces cerevisiae cells (hENT1/2, hCNT1/2/3) or Xenopus laevis oocytes (hENT3/4). All three drugs inhibited hCNT1/3 potently (K _i values, 3–26 μM), hENT1/2 and hCNT2 weakly (K _i values, 0.5–3.1 mM), and hENT3/4 poorly if at all. Rates of transport of [³H]gemcitabine, [¹⁴C]azacitidine, and [³H]decitabine observed in Xenopus oocytes expressing individual recombinant hNTs differed substantially. Cytotoxicity of azacitidine and decitabine was assessed in hNT-expressing or hNT-deficient cultured human cell lines in the absence or presence of transport inhibitors where available. The rank order of cytotoxic sensitivities (IC ₅₀ values, μM) conferred by hNTs were hCNT1 (0.1) > hENT1 (0.3) ? hCNT2 (8.3), hENT2 (9.0) for azacitidine and hENT1 (0.3) > hCNT1 (0.8) ? hENT2, hCNT2 (>100) for decitabine. Protection against cytotoxicity was observed for both drugs in the presence of inhibitors of nucleoside transport, thus suggesting the importance of hNTs in manifestation of toxicity. In summary, all seven hNTs transported azacitidine, with hCNT3 showing the highest rates, whereas hENT1 and hENT2 showed modest transport and hCNT1 and hCNT3 poor transport of decitabine. Our results show for the first time that azacitidine and decitabine exhibit different human nucleoside transportability profiles and their cytotoxicities are dependent on the presence of hNTs, which could serve as potential biomarkers of clinical response.     

The Exact Distribution of the Anderson-Kannemann Statistic

The Anderson-Kannemann test is a rank test for treatment effects in a randomized block design with K treatments and N blocks. In this paper, an algorithm for computing the exact distribution of the Anderson-Kannemann test statistic under the null hypothesis is deduced. Then, the exact distribution is compared with the asymptotic χ²-distribution, and it is shown that the exact distribution is approximated fairly well by the asymptotic distribution. Tables of the exact distribution are given for K = 3, N = 3(1)15; K = 4, N = 3(1)11; K = 5, N = 3(1)7; and K = 6, N = 3(1)5.     

Comparative Effects of Hydrogen Sulfide-Releasing Compounds on [<Superscript>3</Superscript>H]D-Aspartate Release from Bovine Isolated Retinae

We investigated the pharmacological actions of a slow-releasing H₂S donor, GYY 4137; a substrate for the biosynthesis of H₂S, l-cysteine and its precursor, N-acetylcysteine on potassium (K⁺; 50 mM)-evoked [³H]D-aspartate release from bovine isolated retinae using the Superfusion Method. GYY 4137 (10 nM–10 µM), l-cysteine (100 nM–10 µM) and N-acetylcysteine (10 µM–1 mM) elicited a concentration-dependent decrease in K⁺-evoked [³H]D-aspartate release from isolated bovine retinae without affecting basal tritium efflux. At equimolar concentration of 10 µM, the rank order of activity was as follows: l-cysteine?>?GYY 4137?>?N-acetylcysteine. A dual inhibitor of the biosynthetic enzymes for H₂S, cystathionine β-synthase (CBS) and cystathionine γ-lyase (CSE), amino-oxyacetic acid (AOA; 3 mM) reversed the inhibitory responses caused by GYY 4137, l-cysteine and N-acetylcysteine on K⁺-evoked [³H]D-aspartate release. Glibenclamide (300 µM), an inhibitor of K_ATP channels blocked the inhibitory action of GYY 4137 and l-cysteine but not that elicited by N-acetylcysteine on K⁺-induced [³H]D-aspartate release. The inhibitory effect of GYY 4137 and l-cysteine on K⁺-evoked [³H]D-aspartate release was reversed by the non-specific inhibitor of nitric oxide synthase (NOS), l-NAME (300 µM). Furthermore, a specific inhibitor of inducible NOS (iNOS), aminoguanidine (10 µM) blocked the inhibitory action of l-cysteine on K⁺-evoked [³H]D-aspartate release. We conclude that both donors and substrates for H₂S production can inhibit amino acid neurotransmission in bovine isolated retinae, an effect that is dependent, at least in part, upon the intramural biosynthesis of this gas, and on the activity of K_ATP channels and NO synthase.     

An Efficient Alternative to Average Ranks for Testing with Incomplete Ranking Data

In this paper we consider the setting where a group of n judges are to independently rank a series of k objects, but the intended complete rankings are not realized and we are faced with analyzing randomly incomplete ranking vectors. In this paper we propose a new testing procedure for dealing with such data realizations. We concentrate on the problem of testing for no differences in the objects being ranked (i.e., they are indistinguishable) against general alternatives, but our approach could easily be extended to restricted (e.g., ordered or umbrella) alternatives. Using an improvement of a preliminary screening approach previously proposed by the authors, we present an algorithm for computation of the relevant Friedman‐type statistic in the general alternatives setting and present the results of an extensive simulation study comparing the new procedure with the standard approach of imputing average within‐judge ranks to the unranked objects.