Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Response plateau models fitting via isotonic regression

Within the nutritional context, the supplementation of microminerals in bird food is often made in quantities exceeding those required in the attempt to ensure the proper performance of the animals. The experiments of type dosage x response are very common in the determination of levels of nutrients in optimal food balance and include the use of regression models to achieve this objective. Nevertheless, the regression analysis routine, generally, uses a priori information about a possible relationship between the response variable. The isotonic regression is a method of estimation by least squares that generates estimates which preserves data ordering. In the theory of isotonic regression this information is essential and it is expected to increase fitting efficiency. The objective of this work was to use an isotonic regression methodology, as an alternative way of analyzing data of Zn deposition in tibia of male birds of Hubbard lineage. We considered the models of plateau response of polynomial quadratic and linear exponential forms. In addition to these models, we also proposed the fitting of a logarithmic model to the data and the efficiency of the methodology was evaluated by Monte Carlo simulations, considering different scenarios for the parametric values. The isotonization of the data yielded an improvement in all the fitting quality parameters evaluated. Among the models used, the logarithmic presented estimates of the parameters more consistent with the values reported in literature.

Predicting the behaviour of structures under impact loads using geometrically distorted scaled models

When a scaled structure (model or replica) is used to predict the response of a full-size compound (prototype), the model geometric dimensions should relate to the corresponding prototype dimensions by a single scaling factor. However, owing to manufacturing technical restrictions, this condition cannot be accomplished for some of the dimensions in real structures. Accordingly, the distorted geometry will not comply with the overall geometric scaling factor, infringing the Pi theorem requirements for complete dynamic similarity. In the present study, a method which takes geometrical distortions into account is introduced, leading to a model similar to the prototype. As a means to infer the performance of this method, three analytical problems of structures subjected to dynamic loads are analysed. It is shown that the replica developed applying this technique is able to accurately predict the full-size structure behaviour even when the studied models have some of their dimensions severely distorted. (C) 2012 Elsevier Ltd. All rights reserved.

PD controller synthesis from open-loop response measurements of rotating system

In this work, a method of computing PD stabilising gains for rotating systems is presented based on the D-decomposition technique, which requires the sole knowledge of frequency response functions. By applying this method to a rotating system with electromagnetic actuators, it is demonstrated that the stability boundary locus in the plane of feedback gains can be easily plotted, and the most suitable gains can be found to minimise the resonant peak of the system. Experimental results for a Laval rotor show the feasibility of not only controlling lateral shaft vibration and assuring stability, but also helps in predicting the final vibration level achieved by the closed-loop system. These results are obtained based solely on the input-output response information of the system as a whole.

Design-based random permutation models with auxiliary information

We extend the random permutation model to obtain the best linear unbiased estimator of a finite population mean accounting for auxiliary variables under simple random sampling without replacement (SRS) or stratified SRS. The proposed method provides a systematic design-based justification for well-known results involving common estimators derived under minimal assumptions that do not require specification of a functional relationship between the response and the auxiliary variables.

Plasma corticosterone levels in mouse models of pain

Background: Pain markedly activates the hypothalamic-pituitary-adrenal (HPA) axis and increases plasma corticosterone release interfering significantly with nociceptive behaviour as well as the mechanism of action of analgesic drugs. Aims/Methods: In the present study, we monitored the time course of circulating corticosterone in two mouse strains (C57Bl/6 and Balb/C) under different pain models. In addition, the stress response was investigated following animal handling, intrathecal (i.t.) manipulation and habituation to environmental conditions commonly used in nociceptive experimental assays. We also examined the influence of within-cage order of testing on plasma corticosterone. Results: Subcutaneous injection of capsaicin precipitated a prompt stress response whereas carrageenan and complete Freund's adjuvant induced an increased corticosterone release around the third hour post-injection. However, carrageenan induced a longer increased corticosterone in C57Bl/6 mice. In partial sciatic nerve ligation, neuropathic pain model corticosterone increased only in the first days whereas mechanical hypersensitivity remained much longer. Animal handling also represents an important stressor whereas the i.t. injection per se does not exacerbate the handling-induced stress response. Moreover, the order of testing animals from the same cage does not interfere with plasma corticosterone levels in the intrathecal procedure. Animal habituation to the testing apparatus also does not reduce the immediate corticosterone increase as compared with non-habituated mice. Conclusion: Our data indicate that HPA axis activation in acute and chronic pain models is time dependent and may be dissociated from evoked hyperalgesia. Therefore, HPA-axis activation represents an important variable to be considered when designing experimental assays of persistent pain as well as for interpretation of data.

Assessment of variance components in nonlinear mixed-effects elliptical models

The issue of assessing variance components is essential in deciding on the inclusion of random effects in the context of mixed models. In this work we discuss this problem by supposing nonlinear elliptical models for correlated data by using the score-type test proposed in Silvapulle and Silvapulle (1995). Being asymptotically equivalent to the likelihood ratio test and only requiring the estimation under the null hypothesis, this test provides a fairly easy computable alternative for assessing one-sided hypotheses in the context of the marginal model. Taking into account the possible non-normal distribution, we assume that the joint distribution of the response variable and the random effects lies in the elliptical class, which includes light-tailed and heavy-tailed distributions such as Student-t, power exponential, logistic, generalized Student-t, generalized logistic, contaminated normal, and the normal itself, among others. We compare the sensitivity of the score-type test under normal, Student-t and power exponential models for the kinetics data set discussed in Vonesh and Carter (1992) and fitted using the model presented in Russo et al. (2009). Also, a simulation study is performed to analyze the consequences of the kurtosis misspecification.

Hydrothermal pretreatment of sugarcane bagasse using response surface methodology improves digestibility and ethanol production by SSF

Sugarcane bagasse was characterized as a feedstock for the production of ethanol using hydrothermal pretreatment. Reaction temperature and time were varied between 160 and 200A degrees C and 5-20 min, respectively, using a response surface experimental design. The liquid fraction was analyzed for soluble carbohydrates and furan aldehydes. The solid fraction was analyzed for structural carbohydrates and Klason lignin. Pretreatment conditions were evaluated based on enzymatic extraction of glucose and xylose and conversion to ethanol using a simultaneous saccharification and fermentation scheme. SSF experiments were conducted with the washed pretreated biomass. The severity of the pretreatment should be sufficient to drive enzymatic digestion and ethanol yields, however, sugars losses and especially sugar conversion into furans needs to be minimized. As expected, furfural production increased with pretreatment severity and specifically xylose release. However, provided that the severity was kept below a general severity factor of 4.0, production of furfural was below an inhibitory concentration and carbohydrate contents were preserved in the pretreated whole hydrolysate. There were significant interactions between time and temperature for all the responses except cellulose digestion. The models were highly predictive for cellulose digestibility (R (2) = 0.8861) and for ethanol production (R (2) = 0.9581), but less so for xylose extraction. Both cellulose digestion and ethanol production increased with severity, however, high levels of furfural generated under more severe pretreatment conditions favor lower severity pretreatments. The optimal pretreatment condition that gave the highest conversion yield of ethanol, while minimizing furfural production, was judged to be 190A degrees C and 17.2 min. The whole hydrolysate was also converted to ethanol using SSF. To reduce the concentration of inhibitors, the liquid fraction was conditioned prior to fermentation by removing inhibitory chemicals using the fungus Coniochaeta ligniaria.

Changes in intrapopulation resource use patterns of an endangered raptor in response to a disease-mediated crash in prey abundance

1. A long-standing question in ecology is how natural populations respond to a changing environment. Emergent optimal foraging theory-based models for individual variation go beyond the population level and predict how its individuals would respond to disturbances that produce changes in resource availability. 2. Evaluating variations in resource use patterns at the intrapopulation level in wild populations under changing environmental conditions would allow to further advance in the research on foraging ecology and evolution by gaining a better idea of the underlying mechanisms explaining trophic diversity. 3. In this study, we use a large spatio-temporal scale data set (western continental Europe, 19682006) on the diet of Bonellis Eagle Aquila fasciata breeding pairs to analyse the predator trophic responses at the intrapopulation level to a prey population crash. In particular, we borrow metrics from studies on network structure and intrapopulation variation to understand how an emerging infectious disease [the rabbit haemorrhagic disease (RHD)] that caused the density of the eagles primary prey (rabbit Oryctolagus cuniculus) to dramatically drop across Europe impacted on resource use patterns of this endangered raptor. 4. Following the major RHD outbreak, substantial changes in Bonellis Eagles diet diversity and organisation patterns at the intrapopulation level took place. Dietary variation among breeding pairs was larger after than before the outbreak. Before RHD, there were no clusters of pairs with similar diets, but significant clustering emerged after RHD. Moreover, diets at the pair level presented a nested pattern before RHD, but not after. 5. Here, we reveal how intrapopulation patterns of resource use can quantitatively and qualitatively vary, given drastic changes in resource availability. 6. For the first time, we show that a pathogen of a prey species can indirectly impact the intrapopulation patterns of resource use of an endangered predator.

Comparative analysis of the pathological events involved in immune and non-immune TRALI models

Background and Objectives Transfusion-related acute lung injury (TRALI) is characterized by leukocyte transmigration and alveolar capillary leakage shortly after transfusion. TRALI pathogenesis has not been fully elucidated. In some cases, the infusion of alloantibodies (immune model), whereas in others the combination of neutrophil priming by proinflammatory molecules with the subsequent infusion of biological response modifiers (BRMs) in the hemocomponent (non-immune model) have been implicated. Our aim was to compare the pathological events involved in TRALI induced by antibodies or BRMs using murine models. Materials and Methods In the immune model, human HNA-2+ neutrophils were incubated in vitro with a monoclonal antibody (anti-CD177, clone 7D8) directed against the HNA-2 antigen and injected i.v. in NOD/SCID mice. In the non-immune model, BALB/c mice were treated with low doses of lipopolysaccharide (LPS) followed by platelet-activating factor (PAF) infusion 2 h later. Forty minutes after PAF administration, or 6 h after neutrophil injection, lungs were isolated and histological analysis, determination of a variety of cytokines and chemokines including keratinocyte-derived chemokine (KC), MIP-2, the interleukins IL-1 beta, IL-6, IL-8 as well as TNFa, cell influx and alveolar capillary leakage were performed. Results In both models, characteristic histological findings of TRALI and an increase in KC and MIP-2 levels were detected. In contrast to the immune model, in the non-immune model, there was a dramatic increase in IL-1 beta and TNFa. However, capillary leakage was only detected if PAF was administrated. Conclusions Regardless of the triggering event(s), KC, MIP-2 and integrins participate in TRALI pathogenesis, whereas PAF is essential for capillary leakage when two events are involved.

Heterogeneous time-dependent response of adipose tissue during the development of cancer cachexia

Cancer cachexia induces loss of fat mass that accounts for a large part of the dramatic weight loss observed both in humans and in animal models; however, the literature does not provide consistent information regarding the set point of weight loss and how the different visceral adipose tissue depots contribute to this symptom. To evaluate that, 8-week-old male Wistar rats were subcutaneously inoculated with 1 ml (2 x 10(7)) of tumour cells (Walker 256). Samples of different visceral white adipose tissue (WAT) depots were collected at days 0, 4, 7 and 14 and stored at -80 degrees C (seven to ten animals/each day per group). Mesenteric and retroperitoneal depot mass was decreased to the greatest extent on day 14 compared with day 0. Gene and protein expression of PPAR gamma(2) (PPARG) fell significantly following tumour implantation in all three adipose tissue depots while C/EBP alpha (CEBPA) and SREBP-1c (SREBF1) expression decreased over time only in epididymal and retroperitoneal depots. Decreased adipogenic gene expression and morphological disruption of visceral WAT are further supported by the dramatic reduction in mRNA and protein levels of perilipin. Classical markers of inflammation and macrophage infiltration (f4/80, CD68 and MIF-1 alpha) in WAT were significantly increased in the later stage of cachexia (although showing a incremental pattern along the course of cachexia) and presented a depot-specific regulation. These results indicate that impairment in the lipid-storing function of adipose tissue occurs at different times and that the mesenteric adipose tissue is more resistant to the 'fat-reducing effect' than the other visceral depots during cancer cachexia progression. Journal of Endocrinology (2012) 215, 363-373

Low-Level Laser Therapy Influences Mouse Odontoblast-Like Cell Response in Vitro

Objective: The purpose of this study was to analyze the influence of two different irradiation times with 85mW/cm(2) 830nm laser on the behavior of mouse odontoblast-like cells. Background data: The use of low-level laser therapy (LLLT) to stimulate pulp tissue is a reality, but few reports relate odontoblastic responses to irradiation in in vitro models. Methods: Odontoblast-like cells (MDPC-23) were cultivated and divided into three groups: control/nonirradiated (group 1); or irradiated with 85mW/cm(2), 830nm laser for 10 sec (0.8 J/cm(2)) (group 2); or for 50 sec (4.2 J/cm(2)) (group 3) with a wavelength of 830 nm. After 3, 7, and 10 days, it was analyzed: growth curve and cell viability, total protein content, alkaline phosphatase (ALP) activity, calcified nodules detection and quantification, collagen immunolocalization, vascular endothelial growth factor (VEGF) expression, and real-time polymerase chain reaction (PCR) for DMP1 gene. Data were analyzed by Kruskall-Wallis test (alpha = 0.05). Results: Cell growth was smaller in group 2 (p < 0.01), whereas viability was similar in all groups and at all periods. Total protein content and ALP activity increased on the 10th day with 0.8 J/cm(2) (p < 0.01), as well as the detection and quantification of mineralization nodules (p < 0.05), collagen, and VEGF expression (p < 0.01). The expression of DMP1 increased in all groups (p < 0.05) compared with control at 3 days, except for 0.8 J/cm(2) at 3 days and control at 10 days. Conclusions: LLLT influenced the behavior of odontoblast-like cells; the shorter time/smallest energy density promoted the expression of odontoblastic phenotype in a more significant way.

Percutaneous sciatic nerve block with tramadol induces analgesia and motor blockade in two animal pain models

Local anesthetic efficacy of tramadol has been reported following intradermal application. Our aim was to investigate the effect of perineural tramadol as the sole analgesic in two pain models. Male Wistar rats (280-380 g; N = 5/group) were used in these experiments. A neurostimulation-guided sciatic nerve block was performed and 2% lidocaine or tramadol (1.25 and 5 mg) was perineurally injected in two different animal pain models. In the flinching behavior test, the number of flinches was evaluated and in the plantar incision model, mechanical and heat thresholds were measured. Motor effects of lidocaine and tramadol were quantified and a motor block score elaborated. Tramadol, 1.25 mg, completely blocked the first and reduced the second phase of the flinching behavior test. In the plantar incision model, tramadol (1.25 mg) increased both paw withdrawal latency in response to radiant heat (8.3 +/- 1.1, 12.7 +/- 1.8, 8.4 +/- 0.8, and 11.1 +/- 3.3 s) and mechanical threshold in response to von Frey filaments (459 +/- 82.8, 447.5 +/- 91.7, 320.1 +/- 120, 126.43 +/- 92.8 mN) at 5, 15, 30, and 60 min, respectively. Sham block or contralateral sciatic nerve block did not differ from perineural saline injection throughout the study in either model. The effect of tramadol was not antagonized by intraperitoneal naloxone. High dose tramadol (5 mg) blocked motor function as well as 2% lidocaine. In conclusion, tramadol blocks nociception and motor function in vivo similar to local anesthetics.

Maternal immune activation in late gestation enhances locomotor response to acute but not chronic amphetamine treatment in male mice offspring: Role of the D1 receptor

Exposure to elevated levels of maternal cytokines can lead to functional abnormalities of the dopaminergic system in the adult offspring, including enhanced amphetamine (AMPH)-induced locomotion. Therefore, it seems reasonable to consider that offspring of challenged mothers would behave differently in models of addictive behavior, such as behavioral sensitization. Thus, we sought to evaluate the effects of prenatal exposure to lipopolysaccharide (LPS) on the locomotor response to acute and chronic AMPH treatment in male mice offspring. For this purpose, LPS (Escherichia coli 0127:B8; 120 mu g/kg) was administered intraperitoneally to pregnant Swiss mice on gestational day 17. At adulthood, male offspring were studied under one of the following conditions: (1) locomotor response to acute AMPH treatment (2.5 or 5.0 mg/kg) in an open field test; (2) behavioral sensitization paradigm, which consists of a daily injection of AMPH (1.0 mg/kg) for 10 days and observation of locomotion in the open field on days 1, 5, 10 (development phase), 15 and 17 (expression phase). The LPS stimulated offspring showed enhancement of the locomotor-stimulant effect after an acute AMPH challenge in comparison to baseline and saline pre-treated mice. They also showed development of behavioral sensitization earlier than the saline pre-treated group, although no changes between saline and LPS pre-treated groups were observed on development or expression of locomotor behavioral sensitization to AMPH. Furthermore, there was up-regulation of D1 receptor protein level within striatum in the LPS-stimulated offspring which was strongly correlated with increased grooming behavior. Taken together, our results indicate that motor and dopaminergic alterations caused by maternal immune activation are restricted to the acute AMPH challenge, mostly due to up-regulation of the D1 receptor within the mesolimbic and nigrostriatal pathways, but no locomotor differences were observed for behavioral sensitization to AMPH. (C) 2012 Elsevier B.V. All rights reserved.

Low-Resolution Molecular Models Reveal the Oligomeric State of the PPAR and the Conformational Organization of Its Domains in Solution

The peroxisome proliferator-activated receptors (PPARs) regulate genes involved in lipid and carbohydrate metabolism, and are targets of drugs approved for human use. Whereas the crystallographic structure of the complex of full length PPAR gamma and RXR alpha is known, structural alterations induced by heterodimer formation and DNA contacts are not well understood. Herein, we report a small-angle X-ray scattering analysis of the oligomeric state of hPPAR gamma alone and in the presence of retinoid X receptor (RXR). The results reveal that, in contrast with other studied nuclear receptors, which predominantly form dimers in solution, hPPAR gamma remains in the monomeric form by itself but forms heterodimers with hRXR alpha. The low-resolution models of hPPAR gamma/RXR alpha complexes predict significant changes in opening angle between heterodimerization partners (LBD) and extended and asymmetric shape of the dimer (LBD-DBD) as compared with X-ray structure of the full-length receptor bound to DNA. These differences between our SAXS models and the high-resolution crystallographic structure might suggest that there are different conformations of functional heterodimer complex in solution. Accordingly, hydrogen/deuterium exchange experiments reveal that the heterodimer binding to DNA promotes more compact and less solvent-accessible conformation of the receptor complex.