Fitting model-based psychometric functions to simultaneity and temporal-order judgment data: MATLAB and R routines

Research on temporal-order perception uses temporal-order judgment (TOJ) tasks or synchrony judgment (SJ) tasks in their binary SJ2 or ternary SJ3 variants. In all cases, two stimuli are presented with some temporal delay, and observers judge the order of presentation. Arbitrary psychometric functions are typically fitted to obtain performance measures such as sensitivity or the point of subjective simultaneity, but the parameters of these functions are uninterpretable. We describe routines in MATLAB and R that fit model-based functions whose parameters are interpretable in terms of the processes underlying temporal-order and simultaneity judgments and responses. These functions arise from an independent-channels model assuming arrival latencies with exponential distributions and a trichotomous decision space. Different routines fit data separately for SJ2, SJ3, and TOJ tasks, jointly for any two tasks, or also jointly for the three tasks (for common cases in which two or even the three tasks were used with the same stimuli and participants). Additional routines provide bootstrap p-values and confidence intervals for estimated parameters. A further routine is included that obtains performance measures from the fitted functions. An R package for Windows and source code of the MATLAB and R routines are available as Supplementary Files.

3D ellipsoid fitting for multi-view gait recognition

Gait recognition approaches continue to struggle with challenges including view-invariance, low-resolution data, robustness to unconstrained environments, and fluctuating gait patterns due to subjects carrying goods or wearing different clothes. Although computationally expensive, model based techniques offer promise over appearance based techniques for these challenges as they gather gait features and interpret gait dynamics in skeleton form. In this paper, we propose a fast 3D ellipsoidal-based gait recognition algorithm using a 3D voxel model derived from multi-view silhouette images. This approach directly solves the limitations of view dependency and self-occlusion in existing ellipse fitting model-based approaches. Voxel models are segmented into four components (left and right legs, above and below the knee), and ellipsoids are fitted to each region using eigenvalue decomposition. Features derived from the ellipsoid parameters are modeled using a Fourier representation to retain the temporal dynamic pattern for classification. We demonstrate the proposed approach using the CMU MoBo database and show that an improvement of 15-20% can be achieved over a 2D ellipse fitting baseline.

A parametric duration model of the reaction times of drivers distracted by mobile phone conversations

The use of mobile phones while driving is more prevalent among young drivers—a less experienced cohort with elevated crash risk. The objective of this study was to examine and better understand the reaction times of young drivers to a traffic event originating in their peripheral vision whilst engaged in a mobile phone conversation. The CARRS-Q Advanced Driving Simulator was used to test a sample of young drivers on various simulated driving tasks, including an event that originated within the driver’s peripheral vision, whereby a pedestrian enters a zebra crossing from a sidewalk. Thirty-two licensed drivers drove the simulator in three phone conditions: baseline (no phone conversation), hands-free and handheld. In addition to driving the simulator each participant completed questionnaires related to driver demographics, driving history, usage of mobile phones while driving, and general mobile phone usage history. The participants were 21 to 26 years old and split evenly by gender. Drivers’ reaction times to a pedestrian in the zebra crossing were modelled using a parametric accelerated failure time (AFT) duration model with a Weibull distribution. Also tested where two different model specifications to account for the structured heterogeneity arising from the repeated measures experimental design. The Weibull AFT model with gamma heterogeneity was found to be the best fitting model and identified four significant variables influencing the reaction times, including phone condition, driver’s age, license type (Provisional license holder or not), and self-reported frequency of usage of handheld phones while driving. The reaction times of drivers were more than 40% longer in the distracted condition compared to baseline (not distracted). Moreover, the impairment of reaction times due to mobile phone conversations was almost double for provisional compared to open license holders. A reduction in the ability to detect traffic events in the periphery whilst distracted presents a significant and measurable safety concern that will undoubtedly persist unless mitigated.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Multi-model inference ranking and applications to physics at the Large Hadron Collider

In the measurement of the Higgs Boson decaying into two photons the parametrization of an appropriate background model is essential for fitting the Higgs signal mass peak over a continuous background. This diphoton background modeling is crucial in the statistical process of calculating exclusion limits and the significance of observations in comparison to a background-only hypothesis. It is therefore ideal to obtain knowledge of the physical shape for the background mass distribution as the use of an improper function can lead to biases in the observed limits. Using an Information-Theoretic (I-T) approach for valid inference we apply Akaike Information Criterion (AIC) as a measure of the separation for a fitting model from the data. We then implement a multi-model inference ranking method to build a fit-model that closest represents the Standard Model background in 2013 diphoton data recorded by the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC). Potential applications and extensions of this model-selection technique are discussed with reference to CMS detector performance measurements as well as in potential physics analyses at future detectors.

Estimating the size of the population of persons who inject drugs in the island of Montréal, Canada, using a six-source capture-recapture model

BACKGROUND: To plan and implement services to persons who inject drugs (PWID), knowing their number is essential. For the island of Montréal, Canada, the only estimate, of 11,700 PWID, was obtained in 1996 through a capture-recapture method. Thirteen years later, this study was undertaken to produce a new estimate. METHODS: PWID were defined as individuals aged 14-65 years, having injected recently and living on the island of Montréal. The study period was 07/01/2009 to 06/30/2010. An estimate was produced using a six-source capture-recapture log-linear regression method. The data sources were two epidemiological studies and four drug dependence treatment centres. Model selection was conducted in two steps, the first focusing on interactions between sources and the second, on age group and gender as covariates and as modulators of interactions. RESULTS: A total of 1480 PWID were identified in the six capture sources. They corresponded to 1132 different individuals. Based on the best-fitting model, which included age group and sex as covariates and six two-source interactions (some modulated by age), the estimated population was 3910 PWID (95% confidence intervals (CI): 3180-4900) which represents a prevalence of 2.8 (95% CI: 2.3-3.5) PWID per 1000 persons aged 14-65 years. CONCLUSIONS: The 2009-2010 estimate represents a two-third reduction compared to the one for 1996. The multisource capture-recapture method is useful to produce estimates of the size of the PWID population. It is of particular interest when conducted at regular intervals thus allowing for close monitoring of the injection phenomenon.

Reliable knowledge discovery with a minimal causal model inducer

This paper presents a Minimal Causal Model Inducer that can be used for the reliable knowledge discovery. The minimal-model semantics of causal discovery is an essential concept for the identification of a best fitting model in the sense of satisfactory consistent with the given data and be the simpler, less expressive model. Consistency is one of major measures of reliability in knowledge discovery. Therefore to develop an algorithm being able to derive a minimal model is an interesting topic in the are of reliable knowledge discovery. various causal induction algorithms and tools developed so far can not guarantee that the derived model is minimal and consistent. It was proved the MML induction approach introduced by Wallace, Keven and Honghua Dai is a minimal causal model learner. In this paper, we further prove that the developed minimal causal model learner is reliable in the sense of satisfactory consistency. The experimental results obtained from the tests on a number of both artificial and real models provided in this paper confirm this theoretical result.

A genetic two-factor model of the covariation among a subset of Multidimensional Aptitude Battery and Wechsler Adult Intelligence Scale - Revised subtests

The phenotypic and genetic factor structure of performance on five Multidimensional Aptitude Battery (MAB) subtests and one Wechsler Adult Intelligence Scale-Revised (WAIS-R) subtest was explored in 390 adolescent twin pairs (184 monozygotic [MZ]; 206 dizygotic (DZ)). The temporal stability of these measures was derived from a subsample of 49 twin pairs, with test-retest correlations ranging from .67 to .85. A phenotypic factor model, in which performance and verbal factors were correlated, provided a good fit to the data. Genetic modeling was based on the phenotypic factor structure, but also took into account the additive genetic (A), common environmental (C), and unique environmental (E) parameters derived from a fully saturated ACE model. The best fitting model was characterized by a genetic correlated two-factor structure with specific effects, a general common environmental factor, and overlapping unique environmental effects. Results are compared to multivariate genetic models reported in children and adults, with the most notable difference being the growing importance of common genes influencing diverse abilities in adolescence. (C) 2003 Elsevier Inc. All rights reserved.

Adapting to the work-life interface : the influence of individual differences, work and family on well-being, mental health and work engagement

Bronfenbrenner.s Bioecological Model, expressed as the developmental equation, D f PPCT, is the theoretical framework for two studies that bring together diverse strands of psychology to study the work-life interface of working adults. Occupational and organizational psychology is focused on the demands and resources of work and family, without emphasising the individual in detail. Health and personality psychology examine the individual but without emphasis on the individual.s work and family roles. The current research used Bronfenbrenner.s theoretical framework to combine individual differences, work and family to understand how these factors influence the working adult.s psychological functioning. Competent development has been defined as high well-being (measured as life satisfaction and psychological well-being) and high work engagement (as work vigour, work dedication and absorption in work) and as the absence of mental illness (as depression, anxiety and stress) and the absence of burnout (as emotional exhaustion, cynicism and professional efficacy). Study 1 and 2 were linked, with Study 1 as a cross-sectional survey and Study 2, a prospective panel study that followed on from the data used in Study1. Participants were recruited from a university and from a large public hospital to take part in a 3-wave, online study where they completed identical surveys at 3-4 month intervals (N = 470 at Time 1 and N = 198 at Time 3). In Study 1, hierarchical multiple regressions were used to assess the effects of individual differences (Block 1, e.g. dispositional optimism, coping self-efficacy, perceived control of time, humour), work and family variables (Block 2, e.g. affective commitment, skill discretion, work hours, children, marital status, family demands) and the work-life interface (Block 3, e.g. direction and quality of spillover between roles, work-life balance) on the outcomes. There were a mosaic of predictors of the outcomes with a group of seven that were the most frequent significant predictors and which represented the individual (dispositional optimism and coping self-efficacy), the workplace (skill discretion, affective commitment and job autonomy) and the work-life interface (negative work-to-family spillover and negative family-to-work spillover). Interestingly, gender and working hours were not important predictors. The effects of job social support, generally and for work-life issues, perceived control of time and egalitarian gender roles on the outcomes were mediated by negative work-to-family spillover, particularly for emotional exhaustion. Further, the effect of negative spillover on depression, anxiety and work engagement was moderated by the individual.s personal and workplace resources. Study 2 modelled the longitudinal relationships between the group of the seven most frequent predictors and the outcomes. Using a set of non-nested models, the relative influences of concurrent functioning, stability and change over time were assessed. The modelling began with models at Time 1, which formed the basis for confirmatory factor analysis (CFA) to establish the underlying relationships between the variables and calculate the composite variables for the longitudinal models. The CFAs were well fitting with few modifications to ensure good fit. However, using burnout and work engagement together required additional analyses to resolve poor fit, with one factor (representing a continuum from burnout to work engagement) being the only acceptable solution. Five different longitudinal models were investigated as the Well-Being, Mental Distress, Well-Being-Mental Health, Work Engagement and Integrated models using differing combinations of the outcomes. The best fitting model for each was a reciprocal model that was trimmed of trivial paths. The strongest paths were the synchronous correlations and the paths within variables over time. The reciprocal paths were more variable with weak to mild effects. There was evidence of gain and loss spirals between the variables over time, with a slight net gain in resources that may provide the mechanism for the accumulation of psychological advantage over a lifetime. The longitudinal models also showed that there are leverage points at which personal, psychological and managerial interventions can be targeted to bolster the individual and provide supportive workplace conditions that also minimise negative spillover. Bronfenbrenner.s developmental equation has been a useful framework for the current research, showing the importance of the person as central to the individual.s experience of the work-life interface. By taking control of their own life, the individual can craft a life path that is most suited to their own needs. Competent developmental outcomes were most likely where the person was optimistic and had high self-efficacy, worked in a job that they were attached to and which allowed them to use their talents and without too much negative spillover between their work and family domains. In this way, individuals had greater well-being, better mental health and greater work engagement at any one time and across time.

Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis

Objectives - It has long been suspected that susceptibility to ankylosing spondylitis (AS) is influenced by genes lying distant to the major histocompatibility complex. This study compares genetic models of AS to assess the most likely mode of inheritance, using recurrence risk ratios in relatives of affected subjects. Methods - Recurrence risk ratios in different degrees of relatives were determined using published data from studies specifically designed to address the question. The methods of Risch were used to determine the expected recurrence risk ratios in different degrees of relatives, assuming equal first degree relative recurrence risk between models. Goodness of fit was determined by χ2 comparison of the expected number of affected subjects with the observed number, given equal numbers of each type of relative studied. Results - The recurrence risks in different degrees of relatives were: monozygotic (MZ) twins 63% (17/27), first degree relatives 8.2% (441/5390), second degree relatives 1.0% (8/834), and third degree relatives 0.7% (7/997). Parent-child recurrence risk (7.9%, 37/466) was not significantly different from the sibling recurrence risk (8.2%, 404/4924), excluding a significant dominance genetic component to susceptibility. Poor fitting models included single gene, genetic heterogeneity, additive, two locus multiplicative, and one locus and residual polygenes (χ2 > 32 (two degrees of freedom), p < 10-6 for all models). The best fitting model studied was a five locus model with multiplicative interaction between loci (χ2 = 1.4 (two degrees of freedom), p = 0.5). Oligogenic multiplicative models were the best fitting over a range of population prevalences and first degree recurrence risk rates. Conclusions - This study suggests that of the genetic models tested, the most likely model operating in AS is an oligogenic model with predominantly multiplicative interaction between loci.

Estimating the Threshold Surface Density of Gp120-CCR5 Complexes Necessary for HIV-1 Envelope-Mediated Cell-Cell Fusion

Reduced expression of CCR5 on target CD4(+) cells lowers their susceptibility to infection by R5-tropic HIV-1, potentially preventing transmission of infection and delaying disease progression. Binding of the HIV-1 envelope (Env) protein gp120 with CCR5 is essential for the entry of R5 viruses into target cells. The threshold surface density of gp120-CCR5 complexes that enables HIV-1 entry remains poorly estimated. We constructed a mathematical model that mimics Env-mediated cell-cell fusion assays, where target CD4(+)CCR5(+) cells are exposed to effector cells expressing Env in the presence of a coreceptor antagonist and the fraction of target cells fused with effector cells is measured. Our model employs a reaction network-based approach to describe protein interactions that precede viral entry coupled with the ternary complex model to quantify the allosteric interactions of the coreceptor antagonist and predicts the fraction of target cells fused. By fitting model predictions to published data of cell-cell fusion in the presence of the CCR5 antagonist vicriviroc, we estimated the threshold surface density of gp120-CCR5 complexes for cell-cell fusion as similar to 20 mu m(-2). Model predictions with this threshold captured data from independent cell-cell fusion assays in the presence of vicriviroc and rapamycin, a drug that modulates CCR5 expression, as well as assays in the presence of maraviroc, another CCR5 antagonist, using sixteen different Env clones derived from transmitted or early founder viruses. Our estimate of the threshold surface density of gp120-CCR5 complexes necessary for HIV-1 entry thus appears robust and may have implications for optimizing treatment with coreceptor antagonists, understanding the non-pathogenic infection of non-human primates, and designing vaccines that suppress the availability of target CD4(+)CCR5(+) cells.

Performance Modeling based on Multidimensional Surface Learning for Performance Predictions of Parallel Applications in Non-Dedicated Environments

Modeling the performance behavior of parallel applications to predict the execution times of the applications for larger problem sizes and number of processors has been an active area of research for several years. The existing curve fitting strategies for performance modeling utilize data from experiments that are conducted under uniform loading conditions. Hence the accuracy of these models degrade when the load conditions on the machines and network change. In this paper, we analyze a curve fitting model that attempts to predict execution times for any load conditions that may exist on the systems during application execution. Based on the experiments conducted with the model for a parallel eigenvalue problem, we propose a multi-dimensional curve-fitting model based on rational polynomials for performance predictions of parallel applications in non-dedicated environments. We used the rational polynomial based model to predict execution times for 2 other parallel applications on systems with large load dynamics. In all the cases, the model gave good predictions of execution times with average percentage prediction errors of less than 20%

Identifying regions for conservation of sloth bears through occupancy modelling in north-eastern Karnataka, India

In the absence of information on species in decline with contracting ranges, management should emphasize remaining populations and protection of their habitats. Threatened by anthropogenic pressure including habitat degradation and loss, sloth bears (Melursus ursinus) in India have become limited in range, habitat, and population size. We identified ecological and anthropogenic determinants of occurrence within an occupancy framework to evaluate habitat suitability of non-protected regions (with sloth bears) in northeastern Karnataka, India. We employed a systematic sampling methodology to yield presence absence data to examine a priori hypotheses of determinants that affected occupancy. These covariates were broadly classified as habitat or anthropogenic factors. Mean number of termite mounds and trees positively influenced sloth bear occupancy, and grazing pressure expounded by mean number of livestock dung affected it negatively. Also, mean percentage of shrub coverage had no impact on bear inhabitance. The best fitting model further predicted habitats in Bukkasagara, Agoli, and Benakal reserved forests to have 38%, 75%, and 88%, respectively, of their sampled grid cells with high occupancies (>0.70) albeit little or no legal protection. We recommend a conservation strategy that includes protection of vegetation stand-structure, maintenance of soil moisture, and enrichment of habitat for the long-term welfare of this species.

Optymalizacja metod analizy Főrsterowskiego rezonansowego przeniesienia energii (FRET) pojedynczych cząsteczek

Wydział Chemii: Zakład Chemii Fizycznej

Proof of polar ejection from the close-binary core of the planetary nebula Abell 63

We present the first detailed kinematical analysis of the planetary nebula Abell 63, which is known to contain the eclipsing close-binary nucleus UU Sge. Abell 63 provides an important test case in investigating the role of close-binary central stars on the evolution of planetary nebulae. Longslit observations were obtained using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. The spectra reveal that the central bright rim of Abell 63 has a tube-like structure. A deep image shows collimated lobes extending from the nebula, which are shown to be high-velocity outflows. The kinematic ages of the nebular rim and the extended lobes are calculated to be 8400 +/- 500 and 12900 +/- 2800 yr, respectively, which suggests that the lobes were formed at an earlier stage than the nebular rim. This is consistent with expectations that disc-generated jets form immediately after the common envelope phase. A morphological-kinematical model of the central nebula is presented and the best-fitting model is found to have the same inclination as the orbital plane of the central binary system; this is the first proof that a close-binary system directly affects the shaping of its nebula. A Hubble-type flow is well-established in the morphological-kinematical modelling of the observed line profiles and imagery. Two possible formation models for the elongated lobes of Abell 63 are considered, (i) a low-density, pressure-driven jet excavates a cavity in the remnant asymptotic giant branch (AGB) envelope; (ii) high-density bullets form the lobes in a single ballistic ejection event.