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.

Robust modeling using the generalized epsilon-skew-t distribution

In this paper, an alternative skew Student-t family of distributions is studied. It is obtained as an extension of the generalized Student-t (GS-t) family introduced by McDonald and Newey [10]. The extension that is obtained can be seen as a reparametrization of the skewed GS-t distribution considered by Theodossiou [14]. A key element in the construction of such an extension is that it can be stochastically represented as a mixture of an epsilon-skew-power-exponential distribution [1] and a generalized-gamma distribution. From this representation, we can readily derive theoretical properties and easy-to-implement simulation schemes. Furthermore, we study some of its main properties including stochastic representation, moments and asymmetry and kurtosis coefficients. We also derive the Fisher information matrix, which is shown to be nonsingular for some special cases such as when the asymmetry parameter is null, that is, at the vicinity of symmetry, and discuss maximum-likelihood estimation. Simulation studies for some particular cases and real data analysis are also reported, illustrating the usefulness of the extension considered.

The role of asymmetries in rock-paper-scissors biodiversity models.

The maintenance of biodiversity is a long standing puzzle in ecology. It is a classical result that if the interactions of the species in an ecosystem are chosen in a random way, then complex ecosystems can't sustain themselves, meaning that the structure of the interactions between the species must be a central component on the preservation of biodiversity and on the stability of ecosystems. The rock-paper-scissors model is one of the paradigmatic models that study how biodiversity is maintained. In this model 3 species dominate each other in a cyclic way (mimicking a trophic cycle), that is, rock dominates scissors, that dominates paper, that dominates rock. In the original version of this model, this dominance obeys a 'Z IND 3' symmetry, in the sense that the strength of dominance is always the same. In this work, we break this symmetry, studying the effects of the addition of an asymmetry parameter. In the usual model, in a two dimensional lattice, the species distribute themselves according to spiral patterns, that can be explained by the complex Landau-Guinzburg equation. With the addition of asymmetry, new spatial patterns appear during the transient and the system either ends in a state with spirals, similar to the ones of the original model, or in a state where unstable spatial patterns dominate or in a state where only one species survives (and biodiversity is lost).

Leptogenesis via hypermagnetic fields and baryon asymmetry

We study baryon asymmetry generation originated from the leptogenesis in the presence of hypermagnetic fields in the early Universe plasma before the electroweak phase I ransition (EWPT). For the simplest Chern-Simons (CS) wave configuration of hypermagnetic field we find the baryon asymmetry growth when the hypermagnetic field value changes due to alpha(2)-dynamo and the lepton asymmetry rises due to the Abelian anomaly. We solve the corresponding integro-differential equations for the lepton asymmetries describing such selfconsistent dynamics for lepto- and baryogenesis in the two scenarios: (i) when a primordial lepton asymmetry sits in right electrons e(R); and (ii) when, in addition to e(R), a left lepton asyninwtty for e(L) and v(eL) at due to chirality flip reactions provided by in Iiigg,s decays at the temperatures, T < T-RL similar to 10 TeV. We find that the baryon asymmetry of the Universe (BAU) rises very fast through such leptogenesis, especially, in strong hypermagnetic fields. Varying (decreasing) the CS wave number parameter k(0) < 10(-7) T-EW one can recover the observable value of BAU, eta(B) similar to 10(-9), where k(0) = 10(-7) T-EW corresponds to the ataxinittat value for CS wave number surviving ohmic dissipation of hypermagnetic field. In the scenario (ii) one predicts the essential difference of the lepton numbers of right- and left electrons at EWPT time, L-eR - L-eL similar to (mu(eR) / mu(eL))/T-EW = Delta mu/T-EW similar or equal to 10(-5) that can be used as an initial condition for chiral asymmetry after EWPT.

Reducing the parameter space for unparticle-inspired models using white dwarf masses

Based on astrophysical constraints derived from Chandrasekhar's mass limit for white dwarfs, we study the effects of the model on the parameters of unparticle-inspired gravity, on scales Lambda(U) > 1 TeV and d(U) approximate to 1.

Charcot foot: Skin temperature as a good clinical parameter for predicting disease outcome

Twenty-eight diabetics presenting with acute Charcot foot were immobilized and the temperature difference between limbs measured at each month. All patients had monthly follow-up visits for a year and the relapse rate was zero. We found that skin temperature is a good parameter to ensure safe immobilization withdrawal. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Transformer Operation at Deep Saturation: Model and Parameter Determination

This paper proposes a model that adequately describes the operation of the transformer at deep saturation, suitable for power-electronics applications, and a method for determining its parameters. Simulation and experimental results are presented to confirm the validity of the model and the method.

Linear parameter varying control design for rotating systems supported by journal bearings

Linear parameter varying (LPV) control is a model-based control technique that takes into account time-varying parameters of the plant. In the case of rotating systems supported by lubricated bearings, the dynamic characteristics of the bearings change in time as a function of the rotating speed. Hence, LPV control can tackle the problem of run-up and run-down operational conditions when dynamic characteristics of the rotating system change significantly in time due to the bearings and high vibration levels occur. In this work, the LPV control design for a flexible shaft supported by plain journal bearings is presented. The model used in the LPV control design is updated from unbalance response experimental results and dynamic coefficients for the entire range of rotating speeds are obtained by numerical optimization. Experimental implementation of the designed LPV control resulted in strong reduction of vibration amplitudes when crossing the critical speed, without affecting system behavior in sub- or supercritical speeds. (C) 2012 Elsevier Ltd. All rights reserved.

Is early manual preference in infants defined by intermanual performance asymmetry in reaching?

In the present study we evaluated the relationship between manual preference and intermanual performance asymmetry in reaching of 5-month-old infants. Manual preference was assessed through frequency of reaches toward toys presented at midline, left or right in egocentric coordinates. Intermanual performance asymmetry was evaluated through kinematic analysis. Results showed that performance was predominantly symmetric between hands. Lateral toy positions induced predominance of ipsilateral reaching, while the midline position led to equivalent distribution between right and left handed reaches. No significant correlation between manual preference and intermanual performance asymmetry was observed. These results converge against the notion that manual preference derives from a genetically determined advantage of movement control favoring the right hand. (C) 2012 Elsevier Inc. All rights reserved.

Analysis of slow- and fast-alpha band asymmetry during performance of a saccadic eye movement task: Dissociation between memory- and attention-driven systems

This study aimed at analyzing the relationship between slow- and fast-alpha asymmetry within frontal cortex and the planning, execution and voluntary control of saccadic eye movements (SEM), and quantitative electroencephalography (qEEG) was recorded using a 20-channel EEG system in 12 healthy participants performing a fixed (i.e., memory-driven) and a random SEM (i.e., stimulus-driven) condition. We find main effects for SEM condition in slow- and fast-alpha asymmetry at electrodes F3-F4, which are located over premotor cortex, specifically a negative asymmetry between conditions. When analyzing electrodes F7-F8, which are located over prefrontal cortex, we found a main effect for condition in slow-alpha asymmetry, particularly a positive asymmetry between conditions. In conclusion, the present approach supports the association of slow- and fast-alpha bands with the planning and preparation of SEM, and the specific role of these sub-bands for both, the attention network and the coordination and integration of sensory information with a (oculo)-motor response. (C) 2011 Elsevier B.V. All rights reserved.

Low resolution structural characterization of the Hsp70-interacting protein - Hip - from Leishmania braziliensis emphasizes its high asymmetry

The Hsp70 is an essential molecular chaperone in protein metabolism since it acts as a pivot with other molecular chaperone families. Several co-chaperones act as regulators of the Hsp70 action cycle, as for instance Hip (Hsp70-interacting protein). Hip is a tetratricopeptide repeat protein (TPR) that interacts with the ATPase domain in the Hsp70-ADP state, stabilizing it and preventing substrate dissociation. Molecular chaperones from protozoans, which can cause some neglected diseases, are poorly studied in terms of structure and function. Here, we investigated the structural features of Hip from the protozoa Leishmania braziliensis (LbHip), one of the causative agents of the leishmaniasis disease. LbHip was heterologously expressed and purified in the folded state, as attested by circular dichroism and intrinsic fluorescence emission techniques. LbHip forms an elongated dimer, as observed by analytical gel filtration chromatography, analytical ultracentrifugation and small angle X-ray scattering (SAXS). With the SAXS data a low resolution model was reconstructed, which shed light on the structure of this protein, emphasizing its elongated shape and suggesting its domain organization. We also investigated the chemical-induced unfolding behavior of LbHip and two transitions were observed. The first transition was related to the unfolding of the TPR domain of each protomer and the second transition of the dimer dissociation. Altogether. LbHip presents a similar structure to mammalian Hip, despite their low level of conservation, suggesting that this class of eukaryotic protein may use a similar mechanism of action. (C) 2012 Elsevier Inc. All rights reserved.

Describing fatigue crack growth and load ratio effects in Al 2524 T3 alloy with an enhanced exponential model

The fatigue crack behavior in metals and alloys under constant amplitude test conditions is usually described by relationships between the crack growth rate da/dN and the stress intensity factor range Delta K. In the present work, an enhanced two-parameter exponential equation of fatigue crack growth was introduced in order to describe sub-critical crack propagation behavior of Al 2524-T3 alloy, commonly used in aircraft engineering applications. It was demonstrated that besides adequately correlating the load ratio effects, the exponential model also accounts for the slight deviations from linearity shown by the experimental curves. A comparison with Elber, Kujawski and "Unified Approach" models allowed for verifying the better performance, when confronted to the other tested models, presented by the exponential model. (C) 2012 Elsevier Ltd. All rights reserved.

A DEVELOPMENTAL PERSPECTIVE OF INTERMANUAL PERFORMANCE ASYMMETRY IN AIMING

This study compared intermanual performance asymmetry between adults and children of different ages in timing components of an aiming task. Participants manipulated a computer mouse with either hand, aiming to rapidly reach targets of different widths with a pointer on a monitor. Results showed longer deceleration movements with the left hand in children, leading to increased intermanual performance asymmetry when aiming at narrower targets. Increased intermanual performance asymmetry in children is conceptualized to derive from planning to minimize trajectory variability when controlling the nondominant hand.

Exploration of the parameter space in an agent-based model of tuberculosis spread: emergence of drug resistance in developing vs developed countries

In this work we present an agent-based model for the spread of tuberculosis where the individuals can be infected with either drug-susceptible or drug-resistant strains and can also receive a treatment. The dynamics of the model and the role of each one of the parameters are explained. The whole set of parameters is explored to check their importance in the numerical simulation results. The model captures the beneficial impact of the adequate treatment on the prevalence of tuberculosis. Nevertheless, depending on the treatment parameters range, it also captures the emergence of drug resistance. Drug resistance emergence is particularly likely to occur for parameter values corresponding to less efficacious treatment, as usually found in developing countries.

THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

Be stars possess gaseous circumstellar decretion disks, which are well described using standard alpha-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter alpha = 1.0 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass-loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.