Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

Numbers of antral follicles during follicular waves in cattle: Evidence for high variation among animals, very high repeatability in individuals, and an inverse association with serum follicle-stimulating hormone concentrations

The extent, causes, and physiological significance of the variation in number of follicles growing during ovarian follicular waves in human beings and cattle are unknown. Therefore, the present study examined the variability and repeatability in numbers of follicles 3 mm or greater in diameter during the follicular waves in bovine estrous cycles, and we determined if the variation in number of follicles during waves was associated with alterations in secretion of FSH, estradiol, inhibin, and insulin-like growth factor I (IGF-I). Dairy cattle were subjected to twice-daily ultrasound analysis to count total number of antral follicles 3 mm or greater in diameter throughout 138 different follicular waves. In another study, blood samples were taken at frequent intervals from cows that consistently had low or very high numbers of follicles during waves and were subjected to immunoassays. Results indicate the following: First, despite an approximately sevenfold variation in number of follicles during waves among animals and marked differences in age, stage of lactation, and season of the year, a very highly repeatable (0.95) number of follicles 3 mm or greater in diameter is maintained during the ovulatory and nonovulatory follicular waves of individuals. Second, variation in number of follicles 3 mm or greater in diameter during waves and the inverse association of number of follicles during waves with FSH are not directly explained by alterations in the patterns of secretion of estradiol, inhibin, or IGF-I. Third, ovarian ultrasound analysis can be used reliably by investigators to identify cattle that consistently have low or high numbers of follicles during waves, thus providing a novel experimental model to determine the causes and physiological significance of the high variation in antral follicle number during follicular waves among single-ovulating species, such as cattle or humans.

Pharmacological analysis of a dopamine D-2Short: G(alpha o) fusion protein expressed in Sf9 cells

A dopamine D-2Short receptor:G(alphao) fusion protein was expressed in Sf9 cells using the baculovirus expression system. [H-3]Spiperone bound to D-2Short:G(alphao) with a pK(d) approximate to 10. Dopamine stimulated the binding of [S-35]guanosine-5'-O-(3-thio)triphosphate (GTPgammaS) to D-2Short:G(alphao) expressed with Gbeta(1)gamma(2) (E-max > 460%; pEC(50) 5.43 +/- 0.06). Most of the putative D-2 antagonists behaved as inverse agonists (suppressing basal [S-35]GTPgammaS binding) at D-2Short:G(alphao)/Gbeta(1)gamma(2) although (-)-suipiride and ziprasidone were neutral antagonists. Competition of [H-3]spiperone binding by dopamine and 10,11-dihydroxy-N-n-propylnorapo-morphine revealed two, binding sites of different affinities, even in the presence of GTP (100 muM). The D-2Short:G(alphao) fusion protein is therefore a good model for characterising D-2 receptors. (C) 2003 Published by Elsevier Science B.V. on behalf of the Federation of European Biochemical Societies.

A study of the effect of heat source location in a ventilated room using multiple regression analysis

Multiple regression analysis is a statistical technique which allows to predict a dependent variable from m ore than one independent variable and also to determine influential independent variables. Using experimental data, in this study the multiple regression analysis is applied to predict the room mean velocity and determine the most influencing parameters on the velocity. More than 120 experiments for four different heat source locations were carried out in a test chamber with a high level wall mounted air supply terminal at air change rates 3-6 ach. The influence of the environmental parameters such as supply air momentum, room heat load, Archimedes number and local temperature ratio, were examined by two methods: a simple regression analysis incorporated into scatter matrix plots and multiple stepwise regression analysis. It is concluded that, when a heat source is located along the jet centre line, the supply momentum mainly influences the room mean velocity regardless of the plume strength. However, when the heat source is located outside the jet region, the local temperature ratio (the inverse of the local heat removal effectiveness) is a major influencing parameter.

Effects of beverage alcohol taxes and prices on drinking: a meta-analysis of 1003 estimates from 112 studies

Aims: We conducted a systematic review of studies examining relationships between measures of beverage alcohol tax or price levels and alcohol sales or self-reported drinking. A total of 112 studies of alcohol tax or price effects were found, containing 1003 estimates of the tax/price–consumption relationship. Design: Studies included analyses of alternative outcome measures, varying subgroups of the population, several statistical models, and using different units of analysis. Multiple estimates were coded from each study, along with numerous study characteristics. Using reported estimates, standard errors, t-ratios, sample sizes and other statistics, we calculated the partial correlation for the relationship between alcohol price or tax and sales or drinking measures for each major model or subgroup reported within each study. Random-effects models were used to combine studies for inverse variance weighted overall estimates of the magnitude and significance of the relationship between alcohol tax/price and drinking. Findings: Simple means of reported elasticities are -0.46 for beer, -0.69 for wine and -0.80 for spirits. Meta-analytical results document the highly significant relationships (P < 0.001) between alcohol tax or price measures and indices of sales or consumption of alcohol (aggregate-level r = -0.17 for beer, -0.30 for wine, -0.29 for spirits and -0.44 for total alcohol). Price/tax also affects heavy drinking significantly (mean reported elasticity = -0.28, individual-level r = -0.01, P < 0.01), but the magnitude of effect is smaller than effects on overall drinking. Conclusions: A large literature establishes that beverage alcohol prices and taxes are related inversely to drinking. Effects are large compared to other prevention policies and programs. Public policies that raise prices of alcohol are an effective means to reduce drinking.

Modelling and control of Hammerstein system using B-spline approximation and the inverse of De Boor algorithm

In this article a simple and effective controller design is introduced for the Hammerstein systems that are identified based on observational input/output data. The nonlinear static function in the Hammerstein system is modelled using a B-spline neural network. The controller is composed by computing the inverse of the B-spline approximated nonlinear static function, and a linear pole assignment controller. The contribution of this article is the inverse of De Boor algorithm that computes the inverse efficiently. Mathematical analysis is provided to prove the convergence of the proposed algorithm. Numerical examples are utilised to demonstrate the efficacy of the proposed approach.

An inverse method for determining source characteristics for emergency response applications

Following a malicious or accidental atmospheric release in an outdoor environment it is essential for first responders to ensure safety by identifying areas where human life may be in danger. For this to happen quickly, reliable information is needed on the source strength and location, and the type of chemical agent released. We present here an inverse modelling technique that estimates the source strength and location of such a release, together with the uncertainty in those estimates, using a limited number of measurements of concentration from a network of chemical sensors considering a single, steady, ground-level source. The technique is evaluated using data from a set of dispersion experiments conducted in a meteorological wind tunnel, where simultaneous measurements of concentration time series were obtained in the plume from a ground-level point-source emission of a passive tracer. In particular, we analyze the response to the number of sensors deployed and their arrangement, and to sampling and model errors. We find that the inverse algorithm can generate acceptable estimates of the source characteristics with as few as four sensors, providing these are well-placed and that the sampling error is controlled. Configurations with at least three sensors in a profile across the plume were found to be superior to other arrangements examined. Analysis of the influence of sampling error due to the use of short averaging times showed that the uncertainty in the source estimates grew as the sampling time decreased. This demonstrated that averaging times greater than about 5min (full scale time) lead to acceptable accuracy.

Long-term cultivation-independent microbial diversity analysis demonstrates that bacterial communities infecting the adult cystic fibrosis lung show stability and resilience.

Whilst not true in all cases, the microbial communities that chronically infect the airways of patients with CF can vary little over a year despite antibiotic perturbation. The species present tended to vary more between than within subjects, suggesting that each CF airway infection is unique, with relatively stable and resilient bacterial communities. The inverse relationship between community richness and disease severity is similar to findings reported in other mucosal infections.

Parametric preference functionals under risk in the gain domain: a Bayesian analysis

The performance of rank dependent preference functionals under risk is comprehensively evaluated using Bayesian model averaging. Model comparisons are made at three levels of heterogeneity plus three ways of linking deterministic and stochastic models: the differences in utilities, the differences in certainty equivalents and contextualutility. Overall, the"bestmodel", which is conditional on the form of heterogeneity is a form of Rank Dependent Utility or Prospect Theory that cap tures the majority of behaviour at both the representative agent and individual level. However, the curvature of the probability weighting function for many individuals is S-shaped, or ostensibly concave or convex rather than the inverse S-shape commonly employed. Also contextual utility is broadly supported across all levels of heterogeneity. Finally, the Priority Heuristic model, previously examined within a deterministic setting, is estimated within a stochastic framework, and allowing for endogenous thresholds does improve model performance although it does not compete well with the other specications considered.

Kraichnan–Leith–Batchelor similarity theory and two-dimensional inverse cascades

We study the scaling properties and Kraichnan–Leith–Batchelor (KLB) theory of forced inverse cascades in generalized two-dimensional (2D) fluids (α-turbulence models) simulated at resolution 8192x8192. We consider α=1 (surface quasigeostrophic flow), α=2 (2D Euler flow) and α=3. The forcing scale is well resolved, a direct cascade is present and there is no large-scale dissipation. Coherent vortices spanning a range of sizes, most larger than the forcing scale, are present for both α=1 and α=2. The active scalar field for α=3 contains comparatively few and small vortices. The energy spectral slopes in the inverse cascade are steeper than the KLB prediction −(7−α)/3 in all three systems. Since we stop the simulations well before the cascades have reached the domain scale, vortex formation and spectral steepening are not due to condensation effects; nor are they caused by large-scale dissipation, which is absent. One- and two-point p.d.f.s, hyperflatness factors and structure functions indicate that the inverse cascades are intermittent and non-Gaussian over much of the inertial range for α=1 and α=2, while the α=3 inverse cascade is much closer to Gaussian and non-intermittent. For α=3 the steep spectrum is close to that associated with enstrophy equipartition. Continuous wavelet analysis shows approximate KLB scaling ℰ(k)∝k−2 (α=1) and ℰ(k)∝k−5/3 (α=2) in the interstitial regions between the coherent vortices. Our results demonstrate that coherent vortex formation (α=1 and α=2) and non-realizability (α=3) cause 2D inverse cascades to deviate from the KLB predictions, but that the flow between the vortices exhibits KLB scaling and non-intermittent statistics for α=1 and α=2.

Analysis of dietary pattern impact on weight status for personalised nutrition through on-line advice: the Food4Me Spanish cohort

Obesity prevalence is increasing. The management of this condition requires a detailed analysis of the global risk factors in order to develop personalised advice. This study is aimed to identify current dietary patterns and habits in Spanish population interested in personalised nutrition and investigate associations with weight status. Self-reported dietary and anthropometrical data from the Spanish participants in the Food4Me study, were used in a multidimensional exploratory analysis to define specific dietary profiles. Two opposing factors were obtained according to food groups’ intake: Factor 1 characterised by a more frequent consumption of traditionally considered unhealthy foods; and Factor 2, where the consumption of “Mediterranean diet” foods was prevalent. Factor 1 showed a direct relationship with BMI (β = 0.226; r2 = 0.259; p < 0.001), while the association with Factor 2 was inverse (β = −0.037; r2 = 0.230; p = 0.348). A total of four categories were defined (Prudent, Healthy, Western, and Compensatory) through classification of the sample in higher or lower adherence to each factor and combining the possibilities. Western and Compensatory dietary patterns, which were characterized by high-density foods consumption, showed positive associations with overweight prevalence. Further analysis showed that prevention of overweight must focus on limiting the intake of known deleterious foods rather than exclusively enhance healthy products.

The exponentiated generalized inverse Gaussian distribution

The modeling and analysis of lifetime data is an important aspect of statistical work in a wide variety of scientific and technological fields. Good (1953) introduced a probability distribution which is commonly used in the analysis of lifetime data. For the first time, based on this distribution, we propose the so-called exponentiated generalized inverse Gaussian distribution, which extends the exponentiated standard gamma distribution (Nadarajah and Kotz, 2006). Various structural properties of the new distribution are derived, including expansions for its moments, moment generating function, moments of the order statistics, and so forth. We discuss maximum likelihood estimation of the model parameters. The usefulness of the new model is illustrated by means of a real data set. (c) 2010 Elsevier B.V. All rights reserved.

Material processing by plasma shock wave generated in an inverse Z-pinch plasma expander

The applicability of plasma shock wave for material processing was investigated using modified inverse Z-pinch device. Shock wave expanding speed and plasma spectral analysis were studied using an internal magnetic,probe and spatially collimated light spectroscopy. The material processing capability of the device was shown by many different surface analysis techniques such as AES, IRS, EPM and SEM. The interactions between a plasma shock wave of similar to4x10(6) cm/s speed with a Si substrate surface shows some ion implantation capability using a nitrogen plasma and thin film formation using a methane plasma.

Computer-assisted analysis of p53 and PCNA expression in oral lesions infected with human papillomavirus

OBJECTIVE: To carry out a retrospective study to determine whether human papillomavirus (HPV) infection and immunohistochemical expression of p53 and proliferating cell nuclear antigen (PCNA) are related to the risk of oral cancer. STUDY DESIGN: Fifty-seven oral biopsies, consisting of 30 oral squamous papillomas (OSPs) and 27 oral squamous cell carcinomas (OSCCs) were tested for the presence of HPV 6/11 and 16/18 by in situ hybridization using catalyzed signal amplification and in situ hybridization. p53 And PCNA expression was analyzed by immunohistochemistry and evaluated quantitatively by image analysis. RESULTS: Nineteen of the 57 oral lesions (33.3%) were positive for HPV. HPV 6/11 was found in 6 of 30 (20%) OSPs and 1 of 27 (3.7%) OSCCs. HPV 16/18 was found in 10 of 27 (37%) OSCCs and 2 of 30 (6.7%) OSPs. Sixteen of the 19 HPV-positive cases (84.2%) were p53 negative; 5 (9%) were HPV 6/11 and 11 (19%) HPV 16/18, with an inverse correlation between the presence of HPV DNA and p53 expression (P=.017, P < .05). PCNA expression appeared in 18 (94.7%) of HPV positive cases, showing that HPV 16/18 was associated with intensity of PCNA expression and with OSCCs (P=.037, P < .05). CONCLUSION: Quantitative evaluation of p53 by image analysis showed an inverse correlation between p53 expression and HPV presence, suggesting protein degradation. Image analysis also demonstrated that PCNA expression was more intense in HPV DNA 16/18 OSCCs. These findings suggest involvement of high-risk HPV types in oral carcinogenesis.