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.

Manipulation effects of prior exercise intensity feedback by the Borg scale during open-loop cycling

Objective To verify the effects of exercise intensity deception by the Borg scale on the ratings of perceived exertion (RPE), heart rate (HR) and performance responses during a constant power output open-loop exercise. Methods Eight healthy men underwent a maximal incremental test on a cycle ergometer to identify the peak power output (PPO) and heart rate deflection point (HRDP). Subsequently, they performed a constant power output trial to exhaustion set at the HRDP intensity, in deception (DEC) and informed (INF) conditions: DEC-subjects were told that they would be cycling at an intensity corresponding to two categories below the RPE quantified at the HRDP; INF-subjects were told that they would cycle at the exact intensity corresponding to the RPE quantified at the HRDP. Results The PPO and power output at the HRDP obtained in maximal incremental tests were 247.5 +/- 32.1 W and 208.1 +/- 27.1 W, respectively. No significant difference in the time to exhaustion was found between DEC (525 +/- 244 s) or INF (499 +/- 224 s) trials. The slope and the first and second measurements of the RPE and HR parameters showed no significant difference between trials. Conclusions Psychophysiological variables such as RPE and HR as well as performance were not affected when exercise intensity was deceptively manipulated via RPE scores. This may suggest that unaltered RPE during exercise is a regulator of performance in this open-loop exercise.

Prior history of sexually transmitted diseases in women living with AIDS in Sao Paulo, Brazil

Objectives: To describe the epidemiological profile, risk behaviors, and the prior history of sexually transmitted diseases (STDs) in women living with acquired immunodeficiency syndrome (AIDS). Methods: Cross-sectional study, performed at the Centro de Referencia e Treinamento em DST/AIDS of Sao Paulo. The social, demographic, behavioral, and clinical data such as age, schooling, marital status, age at first sexual intercourse, number of sexual partners, parity, use of drugs, time of HIV diagnosis, CD4 count, and viral load determination were abstracted from the medical records of women living with AIDS who had gynecological consultation scheduled in the period from June 2008 to May 2009. Results: Out of 710 women who were scheduled to a gynecological consultation during the period of the study, 598 were included. Previous STD was documented for 364 (60.9%; 95% CI: 56.9%-64.8%) women. The associated factors with previous STDs and their respective risks were: human development index (HDI) <0.50 (ORaj = 5.5; 95% CI: 2.8-11.0); non-white race (ORaj = 5.2; 95% CI: 2.5-11.0); first sexual intercourse at or before 15 years of age (ORaj = 4.4; 95% CI: 2.3-8.3); HIV infection follow-up time of nine years or more (ORaj = 4.2; 95% CI: 2.3-7.8)]; number of sexual partners during the entire life between three and five partners (ORaj = 2.2; 95% CI: 1.1-4.6), and six or more sexual partners (ORaj = 3.9; 95% CI: 1.9-8.0%); being a sex worker (ORaj = 1.9; 95% CI: 1.1-3.1). Conclusions: A high prevalence of a prior history of STDs in the studied population was found. It is essential to find better ways to access HIV infection prevention, so that effective interventions can be more widely implemented. (C) 2012 Elsevier Editora Ltda. All rights reserved.

Mating system parameters at hierarchical levels of fruits, individuals and populations in the Brazilian insect-pollinated tropical tree, Tabebuia roseo-alba (Bignoniaceae)

For many tree species, mating system analyses have indicated potential variations in the selfing rate and paternity correlation among fruits within individuals, among individuals within populations, among populations, and from one flowering event to another. In this study, we used eight microsatellite markers to investigate mating systems at two hierarchical levels (fruits within individuals and individuals within populations) for the insect pollinated Neotropical tree Tabebuia roseo-alba. We found that T. roseo-alba has a mixed mating system with predominantly outcrossed mating. The outcrossing rates at the population level were similar across two T. roseo-alba populations; however, the rates varied considerably among individuals within populations. The correlated paternity results at different hierarchical levels showed that there is a high probability of shared paternal parentage when comparing seeds within fruits and among fruits within plants and full-sibs occur in much higher proportion within fruits than among fruits. Significant levels of fixation index were found in both populations and biparental inbreeding is believed to be the main cause of the observed inbreeding. The number of pollen donors contributing to mating was low. Furthermore, open-pollinated seeds varied according to relatedness, including half-sibs, full-sibs, self-sibs and self- half-sibs. In both populations, the effective population size within a family (seed-tree and its offspring) was lower than expected for panmictic populations. Thus, seeds for ex situ conservation genetics, progeny tests and reforestation must be collected from a large number of seed-trees to guarantee an adequate effective population in the sample.

Influence of prior 810-nm-diode intracanal laser irradiation on hydrophilic resin-based sealer obturation

Dentin wall structural changes caused by 810-nm-diode laser irradiation can influence the sealing ability of endodontic sealers. The objective of this study was to evaluate the apical leakage of AH Plus and RealSeal resin-based sealers with and without prior diode laser irradiation. Fifty-two single-rooted mandibular premolars were prepared and divided into 4 groups, according to the endodontic sealer used and the use or non-use of laser irradiation. The protocol for laser irradiation was 2.5W, continuous wave in scanning mode, with 4 exposures per tooth. After sample preparation, apical leakage of 50% ammoniacal silver nitrate impregnation was analyzed. When the teeth were not exposed to irradiation, the Real Seal sealer achieved the highest scores, showing the least leakage, with significant differences at the 5% level (Kruskal-Wallis test, p = 0.0004), compared with AH Plus. When the teeth were exposed to the 810-nm-diode laser irradiation, the sealing ability of AH Plus sealer was improved (p = 0282). In the Real Seal groups, the intracanal laser irradiation did not interfere with the leakage index, showing similar results in the GRS and GRSd groups (p = 0.1009).

Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

A decaying factor accounts for contained activity in neuronal networks with no need of hierarchical or modular organization

The mechanisms responsible for containing activity in systems represented by networks are crucial in various phenomena, for example, in diseases such as epilepsy that affect the neuronal networks and for information dissemination in social networks. The first models to account for contained activity included triggering and inhibition processes, but they cannot be applied to social networks where inhibition is clearly absent. A recent model showed that contained activity can be achieved with no need of inhibition processes provided that the network is subdivided into modules (communities). In this paper, we introduce a new concept inspired in the Hebbian theory, through which containment of activity is achieved by incorporating a dynamics based on a decaying activity in a random walk mechanism preferential to the node activity. Upon selecting the decay coefficient within a proper range, we observed sustained activity in all the networks tested, namely, random, Barabasi-Albert and geographical networks. The generality of this finding was confirmed by showing that modularity is no longer needed if the dynamics based on the integrate-and-fire dynamics incorporated the decay factor. Taken together, these results provide a proof of principle that persistent, restrained network activation might occur in the absence of any particular topological structure. This may be the reason why neuronal activity does not spread out to the entire neuronal network, even when no special topological organization exists.

Why Are Women With Cervical Cancer Not Being Diagnosed in Preinvasive Phase? An Analysis of Risk Factors Using a Hierarchical Model

Objective: To assess the risk factors for delayed diagnosis of uterine cervical lesions. Materials and Methods: This is a case-control study that recruited 178 women at 2 Brazilian hospitals. The cases (n = 74) were composed of women with a late diagnosis of a lesion in the uterine cervix (invasive carcinoma in any stage). The controls (n = 104) were composed of women with cervical lesions diagnosed early on (low-or high-grade intraepithelial lesions). The analysis was performed by means of logistic regression model using a hierarchical model. The socioeconomic and demographic variables were included at level I (distal). Level II (intermediate) included the personal and family antecedents and knowledge about the Papanicolaou test and human papillomavirus. Level III (proximal) encompassed the variables relating to individuals' care for their own health, gynecologic symptoms, and variables relating to access to the health care system. Results: The risk factors for late diagnosis of uterine cervical lesions were age older than 40 years (odds ratio [OR] = 10.4; 95% confidence interval [CI], 2.3-48.4), not knowing the difference between the Papanicolaou test and gynecological pelvic examinations (OR, = 2.5; 95% CI, 1.3-4.9), not thinking that the Papanicolaou test was important (odds ratio [OR], 4.2; 95% CI, 1.3-13.4), and abnormal vaginal bleeding (OR, 15.0; 95% CI, 6.5-35.0). Previous treatment for sexually transmissible disease was a protective factor (OR, 0.3; 95% CI, 0.1-0.8) for delayed diagnosis. Conclusions: Deficiencies in cervical cancer prevention programs in developing countries are not simply a matter of better provision and coverage of Papanicolaou tests. The misconception about the Papanicolaou test is a serious educational problem, as demonstrated by the present study.

The implementation of the Medical Regulation Office and Mobile Emergency Attendance System and its impact on the gravity profile of non-traumatic afflictions treated in a University Hospital: a research study

Abstract Background The public health system of Brazil is structured by a network of increasing complexity, but the low resolution of emergency care at pre-hospital units and the lack of organization of patient flow overloaded the hospitals, mainly the ones of higher complexity. The knowledge of this phenomenon induced Ribeirão Preto to implement the Medical Regulation Office and the Mobile Emergency Attendance System. The objective of this study was to analyze the impact of these services on the gravity profile of non-traumatic afflictions in a University Hospital. Methods The study conducted a retrospective analysis of the medical records of 906 patients older than 13 years of age who entered the Emergency Care Unit of the Hospital of the University of São Paulo School of Medicine at Ribeirão Preto. All presented acute non-traumatic afflictions and were admitted to the Internal Medicine, Surgery or Neurology Departments during two study periods: May 1996 (prior to) and May 2001 (after the implementation of the Medical Regulation Office and Mobile Emergency Attendance System). Demographics and mortality risk levels calculated by Acute Physiology and Chronic Health Evaluation II (APACHE II) were determined. Results From 1996 to 2001, the mean age increased from 49 ± 0.9 to 52 ± 0.9 (P = 0.021), as did the percentage of co-morbidities, from 66.6 to 77.0 (P = 0.0001), the number of in-hospital complications from 260 to 284 (P = 0.0001), the mean calculated APACHE II mortality risk increased from 12.0 ± 0.5 to 14.8 ± 0.6 (P = 0.0008) and mortality rate from 6.1 to 12.2 (P = 0.002). The differences were more significant for patients admitted to the Internal Medicine Department. Conclusion The implementation of the Medical Regulation and Mobile Emergency Attendance System contributed to directing patients with higher gravity scores to the Emergency Care Unit, demonstrating the potential of these services for hierarchical structuring of pre-hospital networks and referrals.

Prior history of sexually transmitted diseases in women living with AIDS in São Paulo, Brazil

OBJECTIVES: To describe the epidemiological profile, risk behaviors, and the prior history of sexually transmitted diseases (STDs) in women living with acquired immunodeficiency syndrome (AIDS). METHODS: Cross-sectional study, performed at the Centro de Referência e Treinamento em DST/AIDS of São Paulo. The social, demographic, behavioral, and clinical data such as age, schooling, marital status, age at first sexual intercourse, number of sexual partners, parity, use of drugs, time of HIV diagnosis, CD4 count, and viral load determination were abstracted from the medical records of women living with AIDS who had gynecological consultation scheduled in the period from June 2008 to May 2009. RESULTS: Out of 710 women who were scheduled to a gynecological consultation during the period of the study, 598 were included. Previous STD was documented for 364 (60.9%; 95% CI: 56.9%-64.8%) women. The associated factors with previous STDs and their respective risks were: human development index (HDI) < 0.50 (ORaj = 5.5; 95% CI: 2.8-11.0); non-white race (ORaj = 5.2; 95% CI: 2.5-11.0); first sexual intercourse at or before 15 years of age (ORaj = 4.4; 95% CI: 2.3-8.3); HIV infection follow-up time of nine years or more (ORaj = 4.2; 95% CI: 2.3-7.8)]; number of sexual partners during the entire life between three and five partners (ORaj = 2.2; 95% CI: 1.1-4.6), and six or more sexual partners (ORaj = 3.9; 95% CI: 1.9-8.0%); being a sex worker (ORaj = 1.9; 95% CI: 1.1-3.1). CONCLUSIONS: A high prevalence of a prior history of STDs in the studied population was found. It is essential to find better ways to access HIV infection prevention, so that effective interventions can be more widely implemented.

Effect of Prior Photodegradation on the Biodegradation of Polypropylene/Poly(3-hydroxybutyrate) Blends

In this work, the effect of blend composition and previous photodegradation on the biodegradation of polypropylene/ poly(3-hydroxybutyrate) (PP/PHB) blends was studied. The individual polymers and blends with or without the addition of poly(ethylene-co-methyl acrylate- co-glycidyl methacrylate) [P(E-MA-GMA)] as a compatibilizer (in the case of 80/20 blend) were exposed to UV light for 4 weeks and their biodegradation was evaluated. The biodegradation of PHB phase within the blends was hindered as PHB was the dispersed phase and PP fibrous particles were observed at the surface of the blend samples after biodegradation. Previous photodegradation lessened PHB biodegradation but enhanced the biodegradation of PP and the blends within the biodegradation time studied. Photodegradation resulted in cracks at the surface of PP and the blends, which probably facilitated the biotic reactions due to an easier access of the enzymes to deeper polymer layers. It also resulted in a decrease of molecular weight of PP phase and formation of carbonyl and hydroxyl groups which were consumed during biodegradation. Size exclusion chromatography analysis revealed that only the short chains of PP were consumed during biodegradation.

Annealed Ising model on hierarchical structures: a transfer matrix approach.

Spin systems in the presence of disorder are described by two sets of degrees of freedom, associated with orientational (spin) and disorder variables, which may be characterized by two distinct relaxation times. Disordered spin models have been mostly investigated in the quenched regime, which is the usual situation in solid state physics, and in which the relaxation time of the disorder variables is much larger than the typical measurement times. In this quenched regime, disorder variables are fixed, and only the orientational variables are duly thermalized. Recent studies in the context of lattice statistical models for the phase diagrams of nematic liquid-crystalline systems have stimulated the interest of going beyond the quenched regime. The phase diagrams predicted by these calculations for a simple Maier-Saupe model turn out to be qualitative different from the quenched case if the two sets of degrees of freedom are allowed to reach thermal equilibrium during the experimental time, which is known as the fully annealed regime. In this work, we develop a transfer matrix formalism to investigate annealed disordered Ising models on two hierarchical structures, the diamond hierarchical lattice (DHL) and the Apollonian network (AN). The calculations follow the same steps used for the analysis of simple uniform systems, which amounts to deriving proper recurrence maps for the thermodynamic and magnetic variables in terms of the generations of the construction of the hierarchical structures. In this context, we may consider different kinds of disorder, and different types of ferromagnetic and anti-ferromagnetic interactions. In the present work, we analyze the effects of dilution, which are produced by the removal of some magnetic ions. The system is treated in a “grand canonical" ensemble. The introduction of two extra fields, related to the concentration of two different types of particles, leads to higher-rank transfer matrices as compared with the formalism for the usual uniform models. Preliminary calculations on a DHL indicate that there is a phase transition for a wide range of dilution concentrations. Ising spin systems on the AN are known to be ferromagnetically ordered at all temperatures; in the presence of dilution, however, there are indications of a disordered (paramagnetic) phase at low concentrations of magnetic ions.

Bias and learning in temporal binding: intervals between actions and outcomes are compressed by prior bias

It has consistently been shown that agents judge the intervals between their actions and outcomes as compressed in time, an effect named intentional binding. In the present work, we investigated whether this effect is result of prior bias volunteers have about the timing of the consequences of their actions, or if it is due to learning that occurs during the experimental session. Volunteers made temporal estimates of the interval between their action and target onset (Action conditions), or between two events (No-Action conditions). Our results show that temporal estimates become shorter throughout each experimental block in both conditions. Moreover, we found that observers judged intervals between action and outcomes as shorter even in very early trials of each block. To quantify the decrease of temporal judgments in experimental blocks, exponential functions were fitted to participants’ temporal judgments. The fitted parameters suggest that observers had different prior biases as to intervals between events in which action was involved. These findings suggest that prior bias might play a more important role in this effect than calibration-type learning processes.

Hierarchical multi-label classification using local neural networks

Hierarchical multi-label classification is a complex classification task where the classes involved in the problem are hierarchically structured and each example may simultaneously belong to more than one class in each hierarchical level. In this paper, we extend our previous works, where we investigated a new local-based classification method that incrementally trains a multi-layer perceptron for each level of the classification hierarchy. Predictions made by a neural network in a given level are used as inputs to the neural network responsible for the prediction in the next level. We compare the proposed method with one state-of-the-art decision-tree induction method and two decision-tree induction methods, using several hierarchical multi-label classification datasets. We perform a thorough experimental analysis, showing that our method obtains competitive results to a robust global method regarding both precision and recall evaluation measures.