A Diagnostic Test for the Mixing Distribution in a Generalised Linear Mixed Model

We introduce a diagnostic test for the mixing distribution in a generalised linear mixed model. The test is based on the difference between the marginal maximum likelihood and conditional maximum likelihood estimates of a subset of the fixed effects in the model. We derive the asymptotic variance of this difference, and propose a test statistic that has a limiting chi-square distribution under the null hypothesis that the mixing distribution is correctly specified. For the important special case of the logistic regression model with random intercepts, we evaluate via simulation the power of the test in finite samples under several alternative distributional forms for the mixing distribution. We illustrate the method by applying it to data from a clinical trial investigating the effects of hormonal contraceptives in women.

A null mutation in the inflammation-associated S100 protein S100A8 causes early resorption of the mouse embryo

S100A8 (also known as CP10 or MRP8) was the first member of the S100 family of calcium-binding proteins shown to be chemotactic for myeloid cells. The gene is expressed together with its dimerization partner S100A9 during myelopoiesis in the fetal liver and in adult bone marrow as well as in mature granulocytes. In this paper we show that S100A8 mRNA is expressed without S100A9 mRNA between 6.5 and 8.5 days postcoitum within fetal cells infiltrating the deciduum in the vicinity of the ectoplacental cone. Targeted disruption of the S100A8 gene caused rapid and synchronous embryo resorption by day 9.5 of development in 100% of homozygous null embryos. Until this point there was no evidence of developmental delay in S100A8(-/-) embryos and decidualization was normal. The results of PCR genotyping around 7.5-8.5 days postcoitum suggest that the null embryos are infiltrated with maternal cells before overt signs of resorption. This work is the first evidence for nonredundant function of a member of the S100 gene family and implies a role in prevention of maternal rejection of the implanting embryo. The S100A8 null provides a new model for studying fetal-maternal interactions during implantation.

A landscape-scale test of the predictive ability of a spatially explicit model for population viability analysis

1. Although population viability analysis (PVA) is widely employed, forecasts from PVA models are rarely tested. This study in a fragmented forest in southern Australia contrasted field data on patch occupancy and abundance for the arboreal marsupial greater glider Petauroides volans with predictions from a generic spatially explicit PVA model. This work represents one of the first landscape-scale tests of its type. 2. Initially we contrasted field data from a set of eucalypt forest patches totalling 437 ha with a naive null model in which forecasts of patch occupancy were made, assuming no fragmentation effects and based simply on remnant area and measured densities derived from nearby unfragmented forest. The naive null model predicted an average total of approximately 170 greater gliders, considerably greater than the true count (n = 81). 3. Congruence was examined between field data and predictions from PVA under several metapopulation modelling scenarios. The metapopulation models performed better than the naive null model. Logistic regression showed highly significant positive relationships between predicted and actual patch occupancy for the four scenarios (P = 0.001-0.006). When the model-derived probability of patch occupancy was high (0.50-0.75, 0.75-1.00), there was greater congruence between actual patch occupancy and the predicted probability of occupancy. 4. For many patches, probability distribution functions indicated that model predictions for animal abundance in a given patch were not outside those expected by chance. However, for some patches the model either substantially over-predicted or under-predicted actual abundance. Some important processes, such as inter-patch dispersal, that influence the distribution and abundance of the greater glider may not have been adequately modelled. 5. Additional landscape-scale tests of PVA models, on a wider range of species, are required to assess further predictions made using these tools. This will help determine those taxa for which predictions are and are not accurate and give insights for improving models for applied conservation management.

Repairs of databases with null values

Dissertação para obtenção do Grau de Mestre em Engenharia Informática

Increased therapeutic window for the R91W mutant form of Rpe65 compared to Rpe65 null backgroud

Purpose:Given the advances of gene therapy studies to cure RPE65-derived Leber Congenital Amaurosis (LCA) (clinical trials phase I) and the heterogeneity of the targeted patients both genetically and phenotypically, it is of prime importance to examine the rescue efficiency of gene transfer in different mutant contexts. Indeed, half of these mutations are missense mutations, leading to potential residual RPE65 activity. Consequently, we wanted to evaluate the effect on retinal activity and cone survival of lentivirus-mediated gene therapy in the R91W knock-in mouse model expressing the mutant Rpe65R91W gene (Samardzija et al. 2008), a mutation found in LCA patients. Notably we investigated whether if the therapeutic window is prolonged in comparison to null mutations. Methods:An HIV-1-derived lentiviral vector (LV) expressing either the GFP or the mouse Rpe65 cDNA under the control of a 0.8 kb fragment of the human Rpe65 promoter (R0.8) was produced by transient transfection of 293T cells. LV-R0.8-RPE65 or GFP was injected into 5-days-old (P5) or 1 month-old R91W mice. Functional rescue was assessed by ERG (1 and 4 months post-injection) and pupillary light response (PLR) recordings and cone survival by histological analysis. Results:Increased light sensitivity was detected by scotopic ERG in animals injected with LV-R0.8-RPE65 at both P5 and 1 month compared to GFP-treated animals or untreated mice. PLR was also improved in some eyes and histological analysis of cone markers showed that the density of cones reached the wild type level in the region of wt RPE65 delivery after treatment at P5. However, the rescue effect of the injection at 1 month was limited and attained 60% of the wild type level, but still more cones were observed in the treated area than in 1 month-old untreated Rpe65R91W mice. Conclusions:We were able to show that lentivirus-mediated Rpe65 gene transfer not only increases retinal activity of the Rpe65R91W mouse and survival of cones after treatment at P5 but also after treatment at 1 month. However even if the treatment at 1 month is more limited (60% of the wild type level) than treatment at P5, the amount of cone markers is increased compared to the proportion found at 1 month of age in untreated animals. This results contrast with the lack of cone rescue by treatment at 1 month of age in Rpe65-/- (Bemelmans et al, 2006). Thus patient suffering from R91W mutation might benefit from a prolonged therapeutic window.

Chaotic systems with a null conditional Lyapunov exponent under nonlinear driving.

Control of a chaotic system by homogeneous nonlinear driving, when a conditional Lyapunov exponent is zero, may give rise to special and interesting synchronizationlike behaviors in which the response evolves in perfect correlation with the drive. Among them, there are the amplification of the drive attractor and the shift of it to a different region of phase space. In this paper, these synchronizationlike behaviors are discussed, and demonstrated by computer simulation of the Lorentz model [E. N. Lorenz, J. Atmos. Sci. 20 130 (1963)] and the double scroll [T. Matsumoto, L. O. Chua, and M. Komuro, IEEE Trans. CAS CAS-32, 798 (1985)].

Improvements to the accuracy of measurements of NO2 by zenith-sky visible spectrometers II: errors in zero using a more complete chemical model

Using a flexible chemical box model with full heterogeneous chemistry, intercepts of chemically modified Langley plots have been computed for the 5 years of zenith-sky NO2 data from Faraday in Antarctica (65°S). By using these intercepts as the effective amount in the reference spectrum, drifts in zero of total vertical NO2 were much reduced. The error in zero of total NO2 is ±0.03×1015 moleccm−2 from one year to another. This error is small enough to determine trends in midsummer and any variability in denoxification between midwinters. The technique also suggests a more sensitive method for determining N2O5 from zenith-sky NO2 data.

Morphology and myofiber composition of skeletal musculature of the forelimb in young and aged wild type and myostatin null mice

Most current research into therapeutic approaches to muscle diseases involves the use of the mouse as an experimental model. Furthermore, a major strategy to alleviate myopathic symptoms through enhancing muscle growth and regeneration is to inhibit the action of myostatin (Mstn), a transforming growth factor-beta (TGF-beta) family member that inhibits muscle growth. Presently, however, no study has expanded the morphological analysis of mouse skeletal muscle beyond a few individual muscles of the distal hindlimb, through which broad conclusions have been based. Therefore, we have initially undertaken an expansive analysis of the skeletal musculature of the mouse forelimb and highlighted the species-specific differences between equivalent muscles of the rat, another prominently used experimental model. Subsequently, we examined the musculature of the forelimb in both young and old adult wild-type (mstn(+/+)) and myostatin null (mstn(-/-)) mice and assessed the potential beneficial and detrimental effects of myostatin deletion on muscle morphology and composition during the aging process. We showed that: (1) the forelimb muscles of the mouse display a more glycolytic phenotype than those of the rat; (2) in the absence of myostatin, the induced myofiber hyperplasia, hypertrophy, and glycolytic conversion all occur in a muscle-specific manner; and, importantly, (3) the loss of myostatin significantly alters the dynamics of postnatal muscle growth and impairs age-related oxidative myofiber conversion.

Characterisation of connective tissue from the hypertrophic skeletal muscle of myostatin null mice

Myostatin is a potent inhibitor of muscle development. Genetic deletion of myostatin in mice results in muscle mass increase, with muscles often weighing three times their normal values. Contracting muscle transfers tension to skeletal elements through an elaborate connective tissue network. Therefore, the connective tissue of skeletal muscle is an integral component of the contractile apparatus. Here we examine the connective tissue architecture in myostatin null muscle. We show that the hypertrophic muscle has decreased connective tissue content compared with wild-type muscle. Secondly, we show that the hypertrophic muscle fails to show the normal increase in muscle connective tissue content during ageing. Therefore, genetic deletion of myostatin results in an increase in contractile elements but a decrease in connective tissue content. We propose a model based on the contractile profile of muscle fibres that reconciles this apparent incompatible tissue composition phenotype.

Food restriction reverses the hyper-muscular phenotype and force generation capacity deficit of the myostatin null mouse

Food restriction has a great impact on skeletal muscle mass by inducing muscle protein breakdown to provide substrates for energy production through gluconeogenesis. Genetic models of hyper-muscularity interfere with the normal balance between protein synthesis and breakdown which eventually results in extreme muscle growth. Mutations or deletions in the myostatin gene result in extreme muscle mass. Here we evaluated the impact of food restriction for a period of 5 weeks on skeletal muscle size (i.e., fibre cross-sectional area), fibre type composition and contractile properties (i.e., tetanic and specific force) in myostatin null mice. We found that this hyper-muscular model was more susceptible to catabolic processes than wild type mice. The mechanism of skeletal muscle mass loss was examined and our data shows that the myostatin null mice placed on a low calorie diet maintained the activity of molecules involved in protein synthesis and did not up-regulate the expression of genes pivotal in ubiquitin-mediated protein degradation. However, we did find an increase in the expression of genes associated with autophagy. Surprisingly, the reduction on muscle size was followed by improved tetanic and specific force in the null mice compared to wild type mice. These data provide evidence that food restriction may revert the hyper-muscular phenotype of the myostatin null mouse restoring muscle function.

The variant rpoS allele of S-enteritidis strain 27655R does not affect virulence in a chick model nor constitutive curliation but does generate a cold-sensitive phenotype

Adherence of pathogenic Escherichia coli and Salmonella spp. to host cells is in part mediated by curli fimbriae which, along with other virulence determinants, are positively regulated by RpoS. Interested in the role and regulation of curli (SEF17) fimbriae of Salmonella enteritidis in poultry infection, we tested the virulence of naturally occurring S. enteritidis PT4 strains 27655R and 27655S which displayed constitutive and null expression of curli (SEF17) fimbriae, respectively, in a chick invasion assay and analysed their rpoS alleles. Both strains were shown to be equally invasive and as invasive as a wild-type phage type 4 strain and an isogenic derivative defective for the elaboration of curli. We showed that the rpoS allele of 27655S was intact even though this strain was non-curliated and we confirmed that a S. enteritidis rpoS::str(r) null mutant was unable to express curli, as anticipated. Strain 27655R, constitutively curliated, possessed a frameshift mutation at position 697 of the rpoS coding sequence which resulted in a truncated product and remained curliated even when transduced to rpoS::str(r). Additionally, rpoS mutants are known to be cold-sensitive, a phenotype confirmed for strain 27655R. Collectively, these data indicated that curliation was not a significant factor for pathogenesis of S. enteritidis in this model and that curliation of strains 27655R and 27655S was independent of RpoS. Significantly, strain 27655R possessed a defective rpoS allele and remained virulent. Here was evidence that supported the concept that different naturally occurring rpoS alleles may generate varying virulence phenotypic traits. (C) 1998 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

Toward multifactorial null models of range contraction in terrestrial vertebrates

The contraction of a species’ distribution range, which results from the extirpation of local populations, generally precedes its extinction. Therefore, understanding drivers of range contraction is important for conservation and management. Although there are many processes that can potentially lead to local extirpation and range contraction, three main null models have been proposed: demographic, contagion, and refuge. The first two models postulate that the probability of local extirpation for a given area depends on its relative position within the range; but these models generate distinct spatial predictions because they assume either a ubiquitous (demographic) or a clinal (contagion) distribution of threats. The third model (refuge) postulates that extirpations are determined by the intensity of human impacts, leading to heterogeneous spatial predictions potentially compatible with those made by the other two null models. A few previous studies have explored the generality of some of these null models, but we present here the first comprehensive evaluation of all three models. Using descriptive indices and regression analyses we contrast the predictions made by each of the null models using empirical spatial data describing range contraction in 386 terrestrial vertebrates (mammals, birds, amphibians, and reptiles) distributed across the World. Observed contraction patterns do not consistently conform to the predictions of any of the three models, suggesting that these may not be adequate null models to evaluate range contraction dynamics among terrestrial vertebrates. Instead, our results support alternative null models that account for both relative position and intensity of human impacts. These new models provide a better multifactorial baseline to describe range contraction patterns in vertebrates. This general baseline can be used to explore how additional factors influence contraction, and ultimately extinction for particular areas or species as well as to predict future changes in light of current and new threats.

Size and power properties of some tests in the Birnbaum-Saunders regression model

The Birnbaum-Saunders distribution has been used quite effectively to model times to failure for materials subject to fatigue and for modeling lifetime data. In this paper we obtain asymptotic expansions, up to order n(-1/2) and under a sequence of Pitman alternatives, for the non-null distribution functions of the likelihood ratio, Wald, score and gradient test statistics in the Birnbaum-Saunders regression model. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the shape parameter. Monte Carlo simulation is presented in order to compare the finite-sample performance of these tests. We also present two empirical applications. (C) 2010 Elsevier B.V. All rights reserved.

Invariant tests in an instrumental variables model with unknown data generating process

In this work we focus on tests for the parameter of an endogenous variable in a weakly identi ed instrumental variable regressionmodel. We propose a new unbiasedness restriction for weighted average power (WAP) tests introduced by Moreira and Moreira (2013). This new boundary condition is motivated by the score e ciency under strong identi cation. It allows reducing computational costs of WAP tests by replacing the strongly unbiased condition. This latter restriction imposes, under the null hypothesis, the test to be uncorrelated to a given statistic with dimension given by the number of instruments. The new proposed boundary condition only imposes the test to be uncorrelated to a linear combination of the statistic. WAP tests under both restrictions to perform similarly numerically. We apply the di erent tests discussed to an empirical example. Using data from Yogo (2004), we assess the e ect of weak instruments on the estimation of the elasticity of inter-temporal substitution of a CCAPM model.

Hamiltonian Formulation of the Yang-Mills field on the null-plane

We have studied the null plane hamiltonian structure of the free Yang Mills fields. Following the Dirac's procedure for constrained systems we have performed a detailed analysis of the constraint structure of the model and we give the generalized Dirac brackets for the physical variables. Using the correspondence principle in the Dime's brackets we obtain the same commutators present in the literature and new ones.