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.

Cutaneous Fibroma in a Captive Common Snapping Turtle (Chelydra serpentina)

An adult female common snapping turtle (Chelydra serpentina) had a mass on the plantar surface of the right forelimb that was removed surgically. Microscopical examination revealed many spindle cells with mild anisocytosis and anisokaryosis and a surrounding collagenous stroma. There were no mitoses. Immunohistochemistry showed that the spindle cells expressed vimentin, but not desmin. A diagnosis of cutaneous fibroma was made. Tumours are reported uncommonly in chelonian species. Cutaneous fibroma has been diagnosed in an alligator snapping turtle (Macrochelys temminckii), but not previously in a common snapping turtle. Crown Copyright (C) 2012 Published by Elsevier Ltd. All rights reserved.

A modified Primal-Dual Logarithmic-Barrier Method for solving the Optimal Power Flow problem with discrete and continuous control variables

The aim of solving the Optimal Power Flow problem is to determine the optimal state of an electric power transmission system, that is, the voltage magnitude and phase angles and the tap ratios of the transformers that optimize the performance of a given system, while satisfying its physical and operating constraints. The Optimal Power Flow problem is modeled as a large-scale mixed-discrete nonlinear programming problem. This paper proposes a method for handling the discrete variables of the Optimal Power Flow problem. A penalty function is presented. Due to the inclusion of the penalty function into the objective function, a sequence of nonlinear programming problems with only continuous variables is obtained and the solutions of these problems converge to a solution of the mixed problem. The obtained nonlinear programming problems are solved by a Primal-Dual Logarithmic-Barrier Method. Numerical tests using the IEEE 14, 30, 118 and 300-Bus test systems indicate that the method is efficient. (C) 2012 Elsevier B.V. All rights reserved.

A Chaotic Approach of Differential Evolution Optimization Applied to Loudspeaker Design Problem

Over the past few years, the field of global optimization has been very active, producing different kinds of deterministic and stochastic algorithms for optimization in the continuous domain. These days, the use of evolutionary algorithms (EAs) to solve optimization problems is a common practice due to their competitive performance on complex search spaces. EAs are well known for their ability to deal with nonlinear and complex optimization problems. Differential evolution (DE) algorithms are a family of evolutionary optimization techniques that use a rather greedy and less stochastic approach to problem solving, when compared to classical evolutionary algorithms. The main idea is to construct, at each generation, for each element of the population a mutant vector, which is constructed through a specific mutation operation based on adding differences between randomly selected elements of the population to another element. Due to its simple implementation, minimum mathematical processing and good optimization capability, DE has attracted attention. This paper proposes a new approach to solve electromagnetic design problems that combines the DE algorithm with a generator of chaos sequences. This approach is tested on the design of a loudspeaker model with 17 degrees of freedom, for showing its applicability to electromagnetic problems. The results show that the DE algorithm with chaotic sequences presents better, or at least similar, results when compared to the standard DE algorithm and other evolutionary algorithms available in the literature.

Changes in spatial and temporal gene expression during incompatible interaction between common bean and anthracnose pathogen

Common bean, one of the most important legumes for human consumption, may have drastic reduction in yield due to anthracnose, a disease caused by the fungus Colletotrichum lindemuthianum. Rapid induction of the plant defense mechanisms is essential to establish an incompatible interaction with this pathogenic fungus. In this study, we evaluated spatial (leaves, epicotyls and hypocotyls) and temporal (24, 48, 72 and 96 hours after inoculation [HAI]) relative expression (RE) of 12 defense-related transcripts selected from previously developed ESTs libraries, during incompatible interaction between the resistant common bean genotype SEL 1308 and the avirulent anthracnose pathogen race 73, using real time quantitative RT-PCR (RT-qPCR) analysis. All selected transcripts, including the ones coding for pathogenesis-related (PR) proteins (PR1a, PR1b, PR2, and PR16a and PR16b) were differentially regulated upon pathogen inoculation. The expression levels of these transcripts were dependent on the tissue and time post inoculation. This study contributes to a better understanding of the kinetics of induced defenses against a fungal pathogen of common bean and may be used as a base line to study defenses against a broad range of pathogens including bacteria as well as non-host resistance. (C) 2012 Elsevier GmbH. All rights reserved.

A new evolutionary clustering search for a no-wait flow shop problem with set-up times

This paper addresses the m-machine no-wait flow shop problem where the set-up time of a job is separated from its processing time. The performance measure considered is the total flowtime. A new hybrid metaheuristic Genetic Algorithm-Cluster Search is proposed to solve the scheduling problem. The performance of the proposed method is evaluated and the results are compared with the best method reported in the literature. Experimental tests show superiority of the new method for the test problems set, regarding the solution quality. (c) 2012 Elsevier Ltd. All rights reserved.

Biochemical and metabolic profiles of Trichoderma strains isolated from common bean crops in the Brazilian Cerrado, and potential antagonism against Sclerotinia sclerotiorum

Some species of Trichoderma have successfully been used in the commercial biological control of fungal pathogens, e.g., Sclerotinia sclerotiorum, an economically important pathogen of common beans (Phaseolus vulgaris L.). The objectives of the present study were (1) to provide molecular characterization of Trichoderma strains isolated from the Brazilian Cerrado; (2) to assess the metabolic profile of each strain by means of Biolog FF Microplates; and (3) to evaluate the ability of each strain to antagonize S. sclerotiorum via the production of cell wall-degrading enzymes (CWDEs), volatile antibiotics, and dual-culture tests. Among 21 isolates, we identified 42.86 % as Trichoderma asperellum, 33.33 % as Trichoderma harzianum, 14.29 % as Trichoderma tomentosum, 4.76 % as Trichoderma koningiopsis, and 4.76 % as Trichoderma erinaceum. Trichoderma asperellum showed the highest CWDE activity. However, no species secreted a specific group of CWDEs. Trichoderma asperellum 364/01, T. asperellum 483/02, and T. asperellum 356/02 exhibited high and medium specific activities for key enzymes in the mycoparasitic process, but a low capacity for antagonism. We observed no significant correlation between CWDE and antagonism, or between metabolic profile and antagonism. The diversity of Trichoderma species, and in particular of T. harzianum, was clearly reflected in their metabolic profiles. Our findings indicate that the selection of Trichoderma candidates for biological control should be based primarily on the environmental fitness of competitive isolates and the target pathogen. (C) 2012 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

Geometrical and kinematic properties of interfacial waves in stratified oil-water flow in inclined pipe

The stratified oil-water flow pattern is common in the petroleum industry, especially in offshore directional wells and pipelines. Previous studies have shown that the phenomenon of flow pattern transition in stratified flow can be related to the interfacial wave structure (problem of hydrodynamic instability). The study of the wavy stratified flow pattern requires the characterization of the interfacial wave properties, i.e., average shape, celerity and geometric properties (amplitude and wavelength) as a function of holdup, inclination angle and phases' relative velocity. However, the data available in the literature on wavy stratified flow is scanty, especially in inclined pipes and when oil is viscous. This paper presents new geometric and kinematic interfacial wave properties as a function of a proposed two-phase Froude number in the wavy-stratified liquid-liquid flow. The experimental work was conducted in a glass test line of 12 m and 0.026 m id., oil (density and viscosity of 828 kg/m(3) and 0.3 Pa s at 20 degrees C, respectively) and water as the working fluids at several inclinations from horizontal (-20 degrees, -10 degrees, 0 degrees, 10 degrees, 20 degrees). The results suggest a physical relation between wave shape and the hydrodynamic stability of the stratified liquid-liquid flow pattern. (C) 2011 Elsevier Inc. All rights reserved.

Modeling and measuring the spatial relation "along": Regions, contours and fuzzy sets

The analysis of spatial relations among objects in an image is an important vision problem that involves both shape analysis and structural pattern recognition. In this paper, we propose a new approach to characterize the spatial relation along, an important feature of spatial configurations in space that has been overlooked in the literature up to now. We propose a mathematical definition of the degree to which an object A is along an object B, based on the region between A and B and a degree of elongatedness of this region. In order to better fit the perceptual meaning of the relation, distance information is included as well. In order to cover a more wide range of potential applications, both the crisp and fuzzy cases are considered. In the crisp case, the objects are represented in terms of 2D regions or ID contours, and the definition of the alongness between them is derived from a visibility notion and from the region between the objects. However, the computational complexity of this approach leads us to the proposition of a new model to calculate the between region using the convex hull of the contours. On the fuzzy side, the region-based approach is extended. Experimental results obtained using synthetic shapes and brain structures in medical imaging corroborate the proposed model and the derived measures of alongness, thus showing that they agree with the common sense. (C) 2011 Elsevier Ltd. All rights reserved.

The tangential variation of a localized flux-type eigenvalue problem

In this work the differentiability of the principal eigenvalue lambda = lambda(1)(Gamma) to the localized Steklov problem -Delta u + qu = 0 in Omega, partial derivative u/partial derivative nu = lambda chi(Gamma)(x)u on partial derivative Omega, where Gamma subset of partial derivative Omega is a smooth subdomain of partial derivative Omega and chi(Gamma) is its characteristic function relative to partial derivative Omega, is shown. As a key point, the flux subdomain Gamma is regarded here as the variable with respect to which such differentiation is performed. An explicit formula for the derivative of lambda(1) (Gamma) with respect to Gamma is obtained. The lack of regularity up to the boundary of the first derivative of the principal eigenfunctions is a further intrinsic feature of the problem. Therefore, the whole analysis must be done in the weak sense of H(1)(Omega). The study is of interest in mathematical models in morphogenesis. (C) 2011 Elsevier Inc. All rights reserved.

Retinal projections and neurochemical characterization of the pregeniculate nucleus of the common marmoset (Callithrix jacchus)

In mammals, the suprachiasmatic nucleus (SCN) and the intergeniculate leaflet (IGL) are the main components of the circadian timing system. The SCN is the site of the endogenous biological clock that generates rhythms and synchronizes them to environmental cues. The IGL is a key structure that modulates SCN activity and is responsible for the transmission of non-photic information to the SCN, thus participating in the integration between photic and non-photic stimuli. Both the SCN and IGL receive projections of retinal ganglion cells and the IGL is connected to the SCN through the geniculohypothalamic tract. Little is known about these structures in the primate brain and the pregeniculate nucleus (PGN) has been suggested to be the primate equivalent of the rodent IGL. The aim of this study was to characterize the PGN of a primate, the common marmoset (Callithrix jacchus), and to analyze its retinal afferents. Here, the marmoset PGN was found to be organized into three subsectors based on neuronal size, pattern of retinal projections, and the distribution of neuropeptide Y-, GAD-, serotonin-, enkephalin- and substance P-labeled terminals. This pattern indicates that the marmoset PGN is equivalent to the IGL. This detailed description contributes to the understanding of the circadian timing system in this primate species considering the importance of the IGL within the context of circadian regulation. (C) 2012 Elsevier B.V. All rights reserved.

BayGO: Bayesian analysis of ontology term enrichment in microarray data

Abstract Background The search for enriched (aka over-represented or enhanced) ontology terms in a list of genes obtained from microarray experiments is becoming a standard procedure for a system-level analysis. This procedure tries to summarize the information focussing on classification designs such as Gene Ontology, KEGG pathways, and so on, instead of focussing on individual genes. Although it is well known in statistics that association and significance are distinct concepts, only the former approach has been used to deal with the ontology term enrichment problem. Results BayGO implements a Bayesian approach to search for enriched terms from microarray data. The R source-code is freely available at http://blasto.iq.usp.br/~tkoide/BayGO in three versions: Linux, which can be easily incorporated into pre-existent pipelines; Windows, to be controlled interactively; and as a web-tool. The software was validated using a bacterial heat shock response dataset, since this stress triggers known system-level responses. Conclusion The Bayesian model accounts for the fact that, eventually, not all the genes from a given category are observable in microarray data due to low intensity signal, quality filters, genes that were not spotted and so on. Moreover, BayGO allows one to measure the statistical association between generic ontology terms and differential expression, instead of working only with the common significance analysis.

Decreasing the number of false positives in sequence classification

Abstract Background A large number of probabilistic models used in sequence analysis assign non-zero probability values to most input sequences. To decide when a given probability is sufficient the most common way is bayesian binary classification, where the probability of the model characterizing the sequence family of interest is compared to that of an alternative probability model. We can use as alternative model a null model. This is the scoring technique used by sequence analysis tools such as HMMER, SAM and INFERNAL. The most prevalent null models are position-independent residue distributions that include: the uniform distribution, genomic distribution, family-specific distribution and the target sequence distribution. This paper presents a study to evaluate the impact of the choice of a null model in the final result of classifications. In particular, we are interested in minimizing the number of false predictions in a classification. This is a crucial issue to reduce costs of biological validation. Results For all the tests, the target null model presented the lowest number of false positives, when using random sequences as a test. The study was performed in DNA sequences using GC content as the measure of content bias, but the results should be valid also for protein sequences. To broaden the application of the results, the study was performed using randomly generated sequences. Previous studies were performed on aminoacid sequences, using only one probabilistic model (HMM) and on a specific benchmark, and lack more general conclusions about the performance of null models. Finally, a benchmark test with P. falciparum confirmed these results. Conclusions Of the evaluated models the best suited for classification are the uniform model and the target model. However, the use of the uniform model presents a GC bias that can cause more false positives for candidate sequences with extreme compositional bias, a characteristic not described in previous studies. In these cases the target model is more dependable for biological validation due to its higher specificity.

HEALTH PROBLEMS AWARENESS DURING TRAVEL AMONG FACULTY MEMBERS OF A LARGE UNIVERSITY IN LATIN AMERICA. PRELIMINARY REPORT

Health safety during trips is based on previous counseling, vaccination and prevention of infections, previous diseases or specific problems related to the destination. Our aim was to assess two aspects, incidence of health problems related to travel and the traveler’s awareness of health safety. To this end we phone-interviewed faculty members of a large public University, randomly selected from humanities, engineering and health schools. Out of 520 attempts, we were able to contact 67 (12.9%) and 46 (68.6%) agreed to participate in the study. There was a large male proportion (37/44, 84.1%), mature adults mostly in their forties and fifties (32/44, 72.7%), all of them with higher education, as you would expect of faculty members. Most described themselves as being sedentary or as taking occasional exercise, with only 15.9% (7/44) taking regular exercise. Preexisting diseases were reported by 15 travelers. Most trips lasted usually one week or less. Duration of the travel was related to the destination, with (12h) or longer trips being taken by 68.2% (30/44) of travelers, and the others taking shorter (3h) domestic trips. Most travelling was made by air (41/44) and only 31.8% (14/44) of the trips were motivated by leisure. Field research trips were not reported. Specific health counseling previous to travel was reported only by two (4.5%). Twenty seven of them (61.4%) reported updated immunization, but 11/30 reported unchecked immunizations. 30% (9/30) reported travel without any health insurance coverage. As a whole group, 6 (13.6%) travelers reported at least one health problem attributed to the trip. All of them were males travelling abroad. Five presented respiratory infections, such as influenza and common cold, one neurological, one orthopedic, one social and one hypertension. There were no gender differences regarding age groups, destination, type of transport, previous health counseling, leisure travel motivation or pre-existing diseases. Interestingly, the two cases of previous health counseling were made by domestic travelers. Our data clearly shows that despite a significant number of travel related health problems, these highly educated faculty members, had a low awareness of those risks, and a significant number of travels are made without prior counseling or health insurance. A counseling program conducted by a tourism and health professional must be implemented for faculty members in order to increase the awareness of travel related health problems.

Phosphite as phosphorus source to grain yield of common bean plants grown in soils under low or adequate phosphate availability

The effects of foliar and soil applied phosphite on grain yield in common bean (Phaseolus vulgaris L.) grown in a weathered soil under low and adequate phosphate availability were evaluated. In the first experiment, treatments were composed of a 2 x 7 + 2 factorial scheme, with 2 soil P levels supplied as phosphate (40 e 200 mg P dm-3 soil), 7 soil P levels supplied as phosphite (0-100 mg P dm-3 soil), and 2 additional treatments (without P supply in soil, and all P supplied as phosphite). In the second experiment, treatments were composed of a 2 x 3 x 2 factorial scheme, with 2 soil phosphate levels (40 e 200 mg P dm-3 soil), combined with 3 nutrient sources applied via foliar sprays (potassium phosphite, potassium phosphate, and potassium chloride as a control), and 2 foliar application numbers (single and two application). Additional treatments showed that phosphite is not P source for common bean nutrition. Phosphite supply in soil increased the P content in shoot (at full physiological maturity stage) and grains, but at the same time considerably decreased grain yield, regardless of the soil phosphate availability. Foliar sprays of phosphite decreased grain yield in plants grown under low soil phosphate availability, but no effect was observed in plants grown under adequate soil phosphate availability. In general, foliar sprays of phosphate did not satisfactorily improve grain yield of the common bean plants grown under low soil phosphate availability.