Joint Confidence Sets in Multiple Dilution Assays

Let {Q₁, …, Q_k} be the potencies of k substances relative to a standard in a multiple dilution assay. Joint confidence bounds for these are given with confidence coefficient at least 1-α. These bounds are easily interpreted; they appeal to available tables; they improve Scheffe's bounds; they are based on applicable probability inequalities together with extensions of Fieller's theorem; and they are genuinely nonparametric. The procedures are illustrated using data from parallel-line and slope-ratio assays.     

An Unbalanced Nested Model with Random Effects Having Unequal Variances

Consider the model Y_ijk=μ + a_i + b_ij + e_ijk (i=1, 2,…, t; j=1,2,…, B_i; k=1,2…,n_ij), where μ is a constant and a¹,b_ij and e_ijk are distributed independently and normally with zero means and variances σ²_ad_ij and σ², respectively, where it is assumed that the d_i's and d_ij's are known. In this paper procedures for estimating the variance components (σ², σ²_a and σ²_b) and for testing the hypothesis σ²_b = 0 and σ²_a = 0 are presented. In the last section the mixed model y_ijk, where x_ijkkm are known constants and β_m's are unknown fixed effects (m = 1, 2,…,p), is transformed to a fixed effect model with equal variances so that least squares theory can be used to draw inferences about the β_m's.     

Characterization and Distribution of a Cloned Rat μ-Opioid Receptor

Abstract: We have cloned and expressed a rat brain cDNA, TS11, that encodes a μ-opioid receptor based on pharmacological, physiological, and anatomical criteria. Membranes were prepared from COS-7 cells transiently expressing TS11 bound [³H]diprenorphine with high affinity (K_D = 0.23 ± 0.04 nM). The rank order potency of drugs competing with [³H]diprenorphine was as follows: levorphanol (K_i = 0.6 ± 0.2 nM) ≈β-endorphin (K_i = 0.7 ± 0.5 nM) ≈ morphine (K_i = 0.8 ± 0.5 nM) ≈ [d -Ala², N-Me-Phe⁴,Gly-ol⁵]-enkephalin (DAMGO; K_i = 1.6 ± 0.5 nM) ? U50,488 (K_i = 910 ± 0.78 nM) > [d -Pen^2,5]-enkephalin (K_i = 3,170 ± 98 nM) > dextrorphan (K_i = 4,100 ± 68 nM). The rank order potencies of these ligands, the stereospecificity of levorphanol, and morphine's subnanomolar K_i are consistent with a μ-opioid binding site. Two additional experiments provided evidence that this opioid-binding site is functionally coupled to G proteins: (a) In COS-7 cells 50 µM 5′-guanylylimidodiphosphate shifted a fraction of receptors with high affinity for DAMGO (IC₅₀ = 3.4 ± 0.5 nM) to a lower-affinity state (IC₅₀ = 89.0 ± 19.0 nM), and (b) exposure of Chinese hamster ovary cells stably expressing the cloned μ-opioid receptor to DAMGO resulted in a dose-dependent, naloxone-sensitive inhibition of forskolin-stimulated cyclic AMP production. The distribution of mRNA corresponding to the μ-opioid receptor encoded by TS11 was determined by in situ hybridization to brain sections prepared from adult female rats. The highest levels of μ-receptor mRNA were detected in the thalamus, medial habenula, and the caudate putamen; however, significant hybridization was also observed in many other brain regions, including the hypothalamus.     

Effect of Radiotransmitters on Northern Bobwhite Annual Survival

Abstract: Radiotelemetry has been widely used in northern bobwhite (Colinus virginianus) research to estimate survival rates and other demographic parameters. Biologists have used this knowledge to study bobwhite ecology, develop management theory, and base management actions. We tested the assumption that radiotransmitters do not bias survival rates of bobwhites by comparing survival rates of banded bobwhites with and without radiotransmitters on Tall Timbers Research Station (TTRS) from 1999 to 2004. We used Burnham's model in Program MARK and model-selection procedures to determine relative importance of year, gender, and radiotagged status on annual survival rates and recovery processes. Three plausible models (relative quasi-likelihood Akaike's Information Criterion [ΔQAIC_c] < 3) included year dependence in survival and an additive effect of gender but no radiotransmitter effect. Models including a radiotransmitter effect in survival were >8 ΔQAIC_c from the top models, had low Akaike model weights (w_i < 0.007), and low importance weight (S̀w_i(radio) = 0.01). We also compared band—recapture survival estimates from the QAIC_c minimizing model to staggered entry Kaplan—Meier (KM) survival estimates from 2000 to 2003. Annual KM survival estimates of male and female bobwhites were within the 95% confidence interval of band—recapture estimates in 7 of 8 comparisons. We conclude that radiotelemetry is a reliable technique for determining bobwhite survival. Managers should view information from properly conducted telemetry research as reliable and useful for management.     

Unconditional Non-Asymptotic One-Sided Tests for Independent Binomial Proportions When the Interest Lies in Showing Non-Inferiority and/or Superiority

For two independent binomial proportions Barnard (1947) has introduced a method to construct a non-asymptotic unconditional test by maximisation of the probabilities over the ‘classical’ null hypothesis H₀= {(θ₁, θ₂) ∈ [0, 1]²: θ₁ = θ₂}. It is shown that this method is also useful when studying test problems for different null hypotheses such as, for example, shifted null hypotheses of the form H₀ = {(θ₁, θ₂) ∈ [0, 1]²: θ₂ ≤ θ₁ ± Δ } for non-inferiority and 1-sided superiority problems (including the classical null hypothesis with a 1-sided alternative hypothesis). We will derive some results for the more general ‘shifted’ null hypotheses of the form H₀ = {(θ₁, θ₂) ∈ [0, 1]²: θ₂ ≤ g(θ₁ )} where g is a non decreasing curvilinear function of θ₁. Two examples for such null hypotheses in the regulatory setting are given. It is shown that the usual asymptotic approximations by the normal distribution may be quite unreliable. Non-asymptotic unconditional tests (and the corresponding p-values) may, therefore, be an alternative, particularly because the effort to compute non-asymptotic unconditional p-values for such more complex situations does not increase as compared to the classical situation. For ‘classical’ null hypotheses it is known that the number of possible p-values derived by the unconditional method is very large, albeit finite, and the same is true for the null hypotheses studied in this paper. In most of the situations investigated it becomes obvious that Barnard's CSM test (1947) when adapted to the respective null space is again a very powerful test. A theorem is provided which in addition to allowing fast algorithms to compute unconditional non-asymptotical p-values fills a methodological gap in the calculation of exact unconditional p-values as it is implemented, for example, in Stat Xact 3 for Windows (1995).     

Fat mass compared to four body condition scoring systems in the Asian elephant (Elephas maximus)

Captive elephant populations are not self‐sustaining due to health concerns possibly related to obesity. Categorizing obesity relies on qualitative analyses like body condition scores (BCS). However, elephant indices have not been validated against measured body composition. The objective was to compare BCS systems to body composition determined by deuterium dilution in 28 zoo‐kept Asian elephants. Elephants were weighed and given deuterated water orally (0.05 ml/kg). Blood was collected at ~0, 24, 120, 240, 360, and 480 hr after dosing. Photographs were taken to score the elephant based on four BCS systems (BCS_Wemmer [0 to 11 scoring], BCS_Morfeld [1 to 5 scoring], BCS_Fernando [0 to 10 scoring], BCS_Wijeyamohan [1 to 10 scoring]). Based on regression analysis, relative fat ranged from ?305 kg to 515 kg, where negative values indicate less and positive values indicate more fat than expected for the elephant's mass in this population. BCS_Fernando was associated with relative fat (p = .020, R² = 0.194). Relative fat, adjusted for sex and age in the statistical model, was associated with BCS_Wemmer (p = .027, R² = 0.389), BCS_Fernando (p = .002, R² = 0.502), and BCS_Wijeyamohan (p = .011, R² = 0.426). Inclusion of zoo and familial relatedness resulted in all BCS systems associated with relative fat (p ≤ .015). Only BCS_Fernando predicted relative fat, unadjusted, suggesting it is the most capable system for practical use. Compared to absolute fat, relative fat may be more biologically relevant as greater fat relative to body mass is more likely to lead to health issues.     

Point and Interval Estimation in the Combination of Bioassay Results with Unequal Variances

Symmetric parallel‐line biological assays involve the estimation of (log) relative potencies. The class of p(≥ 2) combination of symmetric parallel line bioassays are considered in this study. A large sample test for the equality of the several potencies is developed. An estimator and a confidence interval are proposed for the common relative potency parameter. The asymptotic distribution of the proposed test‐statistic under the null hypothesis as well as under a contagious hypothesis is derived.     

SPERM COMPETITIVE ABILITY AND INDICES OF LIFETIME REPRODUCTIVE SUCCESS

Understanding the selection pressures shaping components of male reproductive success is essential for assessing the role of sexual selection on phenotypic evolution. A male's competitive reproductive success is often measured in sequential mating tests by recording P1 (first mating male) and P2 (second mating male) paternity scores. How each of these scores relates to a male's overall fitness, for example, lifetime reproductive success is, however, not known. This information is needed to determine whether males benefit from maximizing both P1 and P2 or by trading off P1 against P2 ability. We measured P1, P2, and an index of lifetime reproductive success (LRS_i, a male's competitive reproductive success measured over 12 days) for individual male Drosophila melanogaster. We found no evidence for phenotypic correlations between P1 and P2. In addition, whereas both P1 and P2 were associated with relative LRS_i, only P2 predicted absolute LRS_i. The results suggest that P2 was most closely linked to LRS_i in the wild‐type population studied, a finding which may be common to species with strong second male sperm precedence. The study illustrates how P1 and P2 can have differing relationships with a male's overall reproductive success, and highlights the importance of understanding commonly used measures of sperm competition in the currency of fitness.     

Intensity of nitrification in soil as a function of time and soil moisture

Proceeding from three previously derived expressions for the intensity of nitrification in soil as a function of time (logΣN=K.logt+q), as a function of incubation moisture (logΣN=A.pF _i+B), as a function of initial moisture (logΣN=C.pF _v+D), it was shown that the nitrification intensity as a function of time and of moisture can be expressed by the bilinear function log ΣN=a.pF _i.logT+b.pF _i+c.logt+d; as a function of time and of initial moisture by the bilinear function logΣ=N=a.pF _v.logt+b.pF _v+c.logt+d; as a function of initial and incubation moisture by the bilinear function log ΣN=a.pF _ipF _v+b.pF _i+c.pF _v+d. The intensity of nitrification as a function of time, incubation moisture and initial moisture may be expressed by the multilinear function log ΣN=a.pF _i.pF _v.logt+b.pF _i.pF _v+c.pF _i.logt+d.pF _v.logt+e .pF _i+f.pF _v=g.logt+h. This function is valid for all the incubation moistures lying between pF _i 3.0 and 4.0 and for all initial moistures between 3.5 and 5.9 provided that the incubation temperature remains constant.     

An Unbalanced Two-way Model With Random Effects Having Unequal Variances

Consider the model y_ijk=u ± a_i ± b_i ± c_ij ± e_ijk i=1, 2,…, t; j=1, 2,…b; k=1, 2,…,n_ij where μ is a constant and a_i, b_i, c_ij are distributed independently and normally with zero means and variances Δ₂ Δ2/bd_ij and δ₂ respectively. It is assumed that d_i's, and d_ij's are known (positive) constants (for all i and j). In this paper procedures for estimating the variance components (Δ₂, Δ²_b and Δ²_a) and for testing the hypothesis H_oc:Δ²_c/Δ² = y₃ and H_oa:Δ²_b/Δ² = y₄ (where y₂, y₃, and y₄, are specified constants) are presented. A generalization for the mixed model case is discussed in the last section.     

A Note on the Mann-Whitney Statistic for Lehmann Alternatives

For the case of the LEHMANN alternatives the paper presents some new facts on the MANN -WHITNEY statistic and, in particular, its variance V(p, m, n), where p = P(x_i<y_i). Explicit formulas for U and V are used to prove, among other things, the following propositions: For any m, n, V is a one-hump function of p, and the hump always lies in the interval (1/2(3 - √5), 1/2(√5 - 1)). If no restrictions are imposed on p the boundaries of this interval are sharp. Given s = m + n, V(1/2, s/2,s/2) is maximal among all values V(p, m, n). The formulas allow, moreover, the improvement of the known bounds for the variance of p? = U/mn.     

Multivariate Homogeneity Testing Using an Extended Concept of Nearest Neighbors

Given independent multivariate random samples {X_ij: j = 1, …, n_i} from F_i, for i = 1,2, a test is desired for H₀: F₁ = F₂ against general alternatives. Consider the k · (n₁ + n₂) possible ways of choosing one observation from the combined samples and then one of its k nearest neighbors, and let S_k be the proportion of these choices in which the point and neighbor are in the same sample. Schilling (1986) proposed S_k as a test statistic, but did not indicate how to determine k. We suggest as test statistic W = N Σ kS_k, which we show is equivalent to a sum of N Wilcoxon rank sums, and also to a sum of two two-sample U-statistics of degrees (1, 2) and (2, 1). Simulation with multivariate normal data suggests that our test is generally more powerful than Schilling's test using k = 1, 2, or 3. We illustrate its use with Fisher's iris data.     

Hirschsprung's disease prevalence in Europe: A register based study

Background: Hirschsprung's disease is a congenital gut motility disorder, characterised by the absence of the enteric ganglion cells along the distal gut. The aim of this study was to describe the epidemiology of Hirschsprung's disease, including additional congenital anomalies, total prevalence, trends, and association with maternal age. Methods: Cases of Hirschsprung's disease delivered during 1980 to 2009 notified to 31 European Surveillance of Congenital Anomaly registers formed the population‐based case‐series. Prevalence rates and 95% confidence intervals were calculated as the number of cases per 10,000 births. Multilevel Poisson regression was performed to investigate trends in prevalence, geographical variation and the association with maternal age. Results: There were 1,322 cases of Hirschsprung's disease among 12,146,210 births. The total prevalence was 1.09 (95% confidence interval, 1.03–1.15) per 10,000 births and there was a small but significant increase in prevalence over time (relative risk = 1.01; 95% credible interval, 1.00–1.02; p = 0.004). There was evidence of geographical heterogeneity in prevalence (p < 0.001). Excluding 146 (11.0%) cases with chromosomal anomalies or genetic syndromes, there were 1,176 cases (prevalence = 0.97; 95% confidence interval, 0.91–1.03 per 10,000 births), of which 137 (11.6%) had major structural anomalies. There was no evidence of a significant increased risk of Hirschsprung's disease in cases born to women aged ≥35 years compared with those aged 25 to 29 (relative risk = 1.09; 95% credible interval, 0.91–1.31; p = 0.355). Conclusion: This large population‐based study found evidence of a small increasing trend in Hirschsprung's disease and differences in prevalence by geographic location. There was also no evidence of an association with maternal age. Birth Defects Research (Part A), 100:695–702, 2014. © 2014 Wiley Periodicals, Inc.     

The Estimation of a Multivariate Fitness Function from Several Samples Taken from a Population

The situation is considered where the multivariate distribution of certain variables X₁, X₂, …, X_p is changing with time in a population because natural selection related to the X's is taking place. It is assumed that random samples taken from the population at times t₁, t₂, …, t_s are available and it is desirable to estimate the fitness function w_t(x₁, x₂,…,x_p) which shows how the number of individuals with X_i = x_i, i = 1, 2, …, p at time t is related to the number of individuals with the same X values at time zero. Tests for population changes are discussed and indices of the selection on the population dispersion and the population mean are proposed. The situation with a multivariate normal distribution is considered as a special case. A maximum likelihood method that can be applied with any form of population distribution is proposed for estimating w_t. The methods discussed in the paper are illustrated with data on four dimensions of male Egyptian skulls covering a time span from about 4500 B.C. to about 300 A.D. In this case there seems to have been very little selection on the population dispersion but considerable selection on means.     

Kinetic Analysis of the Noncompetitive Inhibition of the Lignin-Peroxidase-catalyzed Reaction by Oxalic Acid

Oxalic acid was found to inhibit noncompetitively the Cα-Cβ bond cleavage of veratrylglycerol catalyzed by a lignin peroxidase (LiP) isozyme of the white-rot fungus P. chrysosporium. With greater amounts of oxalic acid in the LiP system, the substrate was not converted to veratraldehyde but was almost all recovered. Oxalic acid was shown to be decomposed to CO₂ during the enzymatic reaction. The results clearly indicate that oxalic acid reduced the cation radical intermediate formed in the reaction back to the substrate to block the production of veratraldehyde. A novel equation has been derived to explain the mechanism for this unique non-competitive inhibition that is different from the classical noncompetitive one. The inhibition constant K_i obtained here, which is different from the classical inhibition constant K_i, is defined as the ratio of the rate constant (k_p) for product formation to the rate constant (k_i) for the reduction of the cation radical to the substrate.     

On Confidence Regions in Multiple Assays

In multiple assays a ‘Standard’ preparation is compared in a single assay with two or more ‘Test’ preparations (Finney, 1978, ch. 11). In addition to the estimation of the relative potency of the ‘Test’ preparations the question arises as to simultaneous confidence limits, or else a joint confidence region for these relative potencies. This paper compares a joint confidence region using a LR statistic with separate limits based on Fieller's theorem. These results are presented in an Example on a multiple assay of two tuberculin preparations A, B compared with a Standard one (data of Wadley, Finney, 1978, Table 11.1).     

A Note on Lancaster's Procedure for the Combination of Independent Events

Lancaster (1961) generalized Fisher's (1932) nonparametric procedure for combining independent p-values by transforming P_i from the i-th experiment to a chi-squared random variable with d_i degrees of freedom, with d_i not necessarily equal to 2. We explore the relationship between Lancaster's procedure and a weighted Lipták procedure (Koziol and Tuckwell, 1994) under which P_i is transformed to the standard normal scale. We investigate approximations to the null distribution of Lancaster's test procedure, chi-squared with d degrees of freedom. We find that the Cornish-Fisher (1960) expansions and the Lugannani-Rice (1980) saddlepoint approximations are quite accurate, for non-integral values of d, and for values of d as low as 20.     

Purification and characterization of a trypsin-like enzyme with fibrinolytic activity present in the abdomen of horn fly, Haematobia irritans irritans (Diptera: Muscidae)

This work describes the purification and characterization of a trypsin-like enzyme with fibrinolytic activity present in the abdomen of Haematobia irritans irritans (Diptera: Muscidae). The enzyme was purified using a one-step process, consisting of affinity chromatography on SBTI-Sepharose. The purified protease showed one major active proteinase band on reverse zymography with 0.15% gelatin, corresponding to a molecular mass of 25.5 kDa, with maximum activity at pH 9.0. The purified trypsin-like enzyme preferentially hydrolyzed synthetic substrates with arginine residue at the P1 position. The K _m values determined for three different substrates were 1.88 × 10^–4, 1.28 × 10^–4, and 1.40 × 10^–4 M for H--benzoyl-Ile-Glu-Gly-Arg-p-nitroanilide (S2222), dl-Ile-Pro-Arg-p-nitroanilide (S2288), and DL-Phe-Pip-Arg-p-nitroanilide (S2238), respectively. The enzyme was strongly inhibited by typical serine proteinase inhibitors such as SBTI (soybean trypsin inhibitor, K _i = 0.19 nM) and BuXI (Bauhinia ungulata factor Xa inhibitor, K _i = 0.48 nM), and less inhibited by LDTI (leech-derived tryptase inhibitor, K _i = 1.5 nM) and its variants LDTI 2T and 5T (0.8 and 1.5 nM, respectively). The most effective inhibitor for this protease was r-aprotinin (r-BPTI) with a K _i value of 39 pM. Synthetic serine protease inhibitors presented only weak inhibition, e.g., benzamidine with K _i = 3.0 × 10^–4 M and phenylmethylsulfonyl fluoride (PMSF) showed traces of inhibition. The purified trypsin-like enzyme also digested natural substrates such as fibrinogen and fibrin net. The protease showed higher activity against fibrinogen and fibrin than did bovine trypsin. These data suggest that the proteolytic enzyme of H. irritans irritans is more specific to proteins from blood than are the vertebrate digestive enzymes. This enzyme's characteristics may be an adaptation resulting from the feeding behavior of this hematophagous insect.     

Cerebral Cortex Gsα Protein Levels and Forskolin-Stimulated Cyclic AMP Formation Are Increased in Bipolar Affective Disorder

Abstract: Experimental animal and peripheral blood cell studies point to guanine nucleotide regulatory (G) protein disturbances in bipolar affective disorder. We have previously reported elevated prefrontal cortex G_sα protein in bipolar affective disorder and have now extended these preliminary observations in a larger number of subjects, assessing the brain regional specificity of these changes in greater detail, determining the functional biochemical correlates of such changes, and evaluating their diagnostic specificity. Membrane G protein (G_sα, G_iα, G_oα, and Gβ) immunoreactivities were estimated by western blotting in postmortem brain regions obtained from 10 patients with a DSMIII-R diagnosis of bipolar affective disorder and 10 nonpsychiatric controls matched on the basis of age, postmortem delay, and brain pH. To examine whether there were functional correlates to the observed elevated G_sα levels, basal and GTPγS-and forskolin-stimulated cyclic AMP production was determined in the same brain regions. Compared with controls, G_sα (52-kDa species) immunoreactivity was significantly (p < 0.05) elevated in prefrontal (+36%), temporal (+65%), and occipital (+96%) cortex but not in hippocampus (+28%), thalamus (-23%), or cerebellum (+21%). In contrast, no significant differences were found in the other G protein subunits (G_iα, G_oα, Gβ) measured in these regions. Forskolin-stimulated cyclic AMP production was significantly increased in temporal (+31%) and occipital (+96%) cortex but not in other regions. No significant differences were apparent in basal or GTPγS-stimulated cyclic AMP production. A significant correlation (r= 0.60, p < 0.001) was observed between forskolin-stimulated cyclic AMP formation and G_sα (52 kDa) immunoreactivity when examined across these cortical regions. The observed increase in G_sα may be specific to bipolar disorders as no significant differences were detected in G_sα levels in temporal cortex from patients with either schizophrenia (n = 7) or Alzheimer's disease (n = 7). In summary, the present study confirms and extends our earlier findings and supports the notion that increased G_sα levels and possibly G_sα-adenylyl cyclase-mediated signal transduction are relevant to the pathophysiology of bipolar affective disorder.     

Association between television in bedroom and adiposity throughout adolescence

Objective: The objective was to determine if having a television (TV) in the bedroom is associated with physical activity (PA), TV/video viewing, and adiposity throughout adolescence. Research Methods and Procedures: Longitudinal data (September 2002 through June 2005) were analyzed of 379 initially 12‐year‐old French adolescents participating as controls in the Intervention Centered on Adolescents’ Physical activity and Sedentary behavior (ICAPS). Presence of a TV set in the bedroom (TV_bedroom) and leisure activities were obtained by questionnaire. There was annual assessment of BMI, waist circumference, and body fat by bioimpedance. Results: In boys but not girls, baseline TV_bedroom was associated with higher TV/video viewing over time [odds ratio (OR) of high TV/video = 1.87; 95% confidence interval, 1.2 to 2.8] and less no‐sport club participation (OR = 0.59; 95% confidence interval, 0.35 to 1.0). Both boys and girls with baseline TV_bedroom had lower reading time (p < 0.0001 in boys; p = 0.04 in girls), while PA did not differ according to TV_bedroom for boys or for girls. For boys only, baseline TV_bedroom was associated with higher BMI (mean BMI over time 20.5 ± 0.5 vs. 19.0 ± 0.5 kg/m²; p = 0.001), waist circumference (70.9 ± 0.9 vs. 67.2 ± 0.8 cm; p < 0.001), and body fat (15.9 ± 0.9% vs. 13.5 ± 0.9%; p < 0.001), without interaction with time. These relationships remained significant after adjustment for socioeconomic status. TV/video viewing explained 26%, 42%, and 36% of the relationships of TV_bedroom with BMI, waist circumference, and body fat, respectively, while addition of other leisure activities in the models only marginally reduced the effects. Discussion: These results suggest the importance of keeping TV out of an adolescent's bedroom from an obesity prevention perspective but show gender differences.