abctools: an R package for tuning approximate Bayesian computation analyses

Approximate Bayesian computation (ABC) is a popular family of algorithms which perform approximate parameter inference when numerical evaluation of the likelihood function is not possible but data can be simulated from the model. They return a sample of parameter values which produce simulations close to the observed dataset. A standard approach is to reduce the simulated and observed datasets to vectors of summary statistics and accept when the difference between these is below a specified threshold. ABC can also be adapted to perform model choice. In this article, we present a new software package for R, abctools which provides methods for tuning ABC algorithms. This includes recent dimension reduction algorithms to tune the choice of summary statistics, and coverage methods to tune the choice of threshold. We provide several illustrations of these routines on applications taken from the ABC literature.

Closed-loop feedback computation model of dynamical reputation based on the local trust evaluation in business-to-consumer e-commerce

Trust and reputation are important factors that influence the success of both traditional transactions in physical social networks and modern e-commerce in virtual Internet environments. It is difficult to define the concept of trust and quantify it because trust has both subjective and objective characteristics at the same time. A well-reported issue with reputation management system in business-to-consumer (BtoC) e-commerce is the “all good reputation” problem. In order to deal with the confusion, a new computational model of reputation is proposed in this paper. The ratings of each customer are set as basic trust score events. In addition, the time series of massive ratings are aggregated to formulate the sellers’ local temporal trust scores by Beta distribution. A logical model of trust and reputation is established based on the analysis of the dynamical relationship between trust and reputation. As for single goods with repeat transactions, an iterative mathematical model of trust and reputation is established with a closed-loop feedback mechanism. Numerical experiments on repeated transactions recorded over a period of 24 months are performed. The experimental results show that the proposed method plays guiding roles for both theoretical research into trust and reputation and the practical design of reputation systems in BtoC e-commerce.

Dynamical trust and reputation computation model for B2C E-Commerce

Trust is one of the most important factors that influence the successful application of network service environments, such as e-commerce, wireless sensor networks, and online social networks. Computation models associated with trust and reputation have been paid special attention in both computer societies and service science in recent years. In this paper, a dynamical computation model of reputation for B2C e-commerce is proposed. Firstly, conceptions associated with trust and reputation are introduced, and the mathematical formula of trust for B2C e-commerce is given. Then a dynamical computation model of reputation is further proposed based on the conception of trust and the relationship between trust and reputation. In the proposed model, classical varying processes of reputation of B2C e-commerce are discussed. Furthermore, the iterative trust and reputation computation models are formulated via a set of difference equations based on the closed-loop feedback mechanism. Finally, a group of numerical simulation experiments are performed to illustrate the proposed model of trust and reputation. Experimental results show that the proposed model is effective in simulating the dynamical processes of trust and reputation for B2C e-commerce.

Scientific output trends in oral health in Brazil

Objective. To describe dental research trends in Brazil (especially population-based oral health) in the early Twenty-first Century. Methods. The abstracts of studies presented at meetings of the Brazilian Society for Dental Medicine Research (Sociedade Brasileira de Pesquisa Odontologica) from 2001-2006 were assessed in terms of methodological design (aggregate or population-based and individual-based studies, observational and intervention studies, cross-sectional and longitudinal studies); general type (literature review, studies with human subjects, and laboratory studies); and classification into one of the 19 specialty categories recognized by the Brazilian Federal Dentistry Council. Of the 10 406 abstracts presented in this period, 5 203 (50%) were reviewed. Results. Concerning methodological design, 87.5% of the abstracts referred to individual-based studies, whereas 12.5% were of aggregate studies. Concerning the general category, 41.7% referred to studies with human subjects. The remaining abstracts (58.3%) described in vitro (31.1%) or in vivo (23.6%) laboratory research and literature reviews (3.6%). Concerning the Council`s specialty categories, only five had a frequency higher than 10.0%: esthetic dentistry, periodontics, endodontics, pediatric dentistry, and population-based oral health. Conclusions. Brazil`s scientific output in the field of oral health for the period 2001-2006 was balanced, with increasing interest in the area of population-based oral health.

A Floating-point Extended Kalman Filter Implementation for Autonomous Mobile Robots

Localization and Mapping are two of the most important capabilities for autonomous mobile robots and have been receiving considerable attention from the scientific computing community over the last 10 years. One of the most efficient methods to address these problems is based on the use of the Extended Kalman Filter (EKF). The EKF simultaneously estimates a model of the environment (map) and the position of the robot based on odometric and exteroceptive sensor information. As this algorithm demands a considerable amount of computation, it is usually executed on high end PCs coupled to the robot. In this work we present an FPGA-based architecture for the EKF algorithm that is capable of processing two-dimensional maps containing up to 1.8 k features at real time (14 Hz), a three-fold improvement over a Pentium M 1.6 GHz, and a 13-fold improvement over an ARM920T 200 MHz. The proposed architecture also consumes only 1.3% of the Pentium and 12.3% of the ARM energy per feature.

Scientific program of the Radioactive Ion Beams Facility in Brasil (RIBRAS)

The RIBRAS facility (Radioactive Ion Beams in Brasil) is installed in connection with the 8MV Pelletron tandem of the University of Sao Paulo Physics Institute. It consists of two superconducting solenoids which focalize light radioactive secondary beams of low energy, produced by transfer reactions. Recent experimental results include the elastic scattering and transfer reactions of (6)He halo nucleus on (9)Be, (27)Al, (51)V and (120)Sn targets. The elastic scattering and transfer of (8)Li and (7)Be on several targets is also being studied. The transfer reaction (8)Li(p,alpha)(5)He of astrophysical interest was also Studied in the E(cm)=0.2-2.5 MeV energy range.

WRITING ELECTRONIC FERROMAGNETIC STATES IN A HIGH-TEMPERATURE PARAMAGNETIC NUCLEAR SPIN SYSTEM

In this paper, we use Nuclear Magnetic Resonance (NMR) to write electronic states of a ferromagnetic system into high-temperature paramagnetic nuclear spins. Through the control of phase and duration of radio frequency pulses, we set the NMR density matrix populations, and apply the technique of quantum state tomography to experimentally obtain the matrix elements of the system, from which we calculate the temperature dependence of magnetization for different magnetic fields. The effects of the variation of temperature and magnetic field over the populations can be mapped in the angles of spin rotations, carried out by the RF pulses. The experimental results are compared to the Brillouin functions of ferromagnetic ordered systems in the mean field approximation for two cases: the mean field is given by (i) B = B(0) + lambda M and (ii) B = B(0) + lambda M + lambda`M(3), where B(0) is the external magnetic field, and lambda, lambda` are mean field parameters. The first case exhibits second order transition, whereas the second case has first order transition with temperature hysteresis. The NMR simulations are in good agreement with the magnetic predictions.

NUCLEAR SPIN 3/2 ELECTRIC QUADRUPOLE RELAXATION AS A QUANTUM COMPUTATION PROCESS

In this work we applied a quantum circuit treatment to describe the nuclear spin relaxation. From the Redfield theory, we obtain a description of the quadrupolar relaxation as a computational process in a spin 3/2 system, through a model in which the environment is comprised by five qubits and three different quantum noise channels. The interaction between the environment and the spin 3/2 nuclei is described by a quantum circuit fully compatible with the Redfield theory of relaxation. Theoretical predictions are compared to experimental data, a short review of quantum channels and relaxation in NMR qubits is also present.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

BASS CYCLIC UNITS AS FACTORS IN A FREE GROUP IN INTEGRAL GROUP RING UNITS

Marciniak and Sehgal showed that if u is a non-trivial bicyclic unit of an integral group ring then there is a bicyclic unit v such that u and v generate a non-abelian free group. A similar result does not hold for Bass cyclic units of infinite order based on non-central elements as some of them have finite order modulo the center. We prove a theorem that suggests that this is the only limitation to obtain a non-abelian free group from a given Bass cyclic unit. More precisely, we prove that if u is a Bass cyclic unit of an integral group ring ZG of a solvable and finite group G, such that u has infinite order modulo the center of U(ZG) and it is based on an element of prime order, then there is a non-abelian free group generated by a power of u and a power of a unit in ZG which is either a Bass cyclic unit or a bicyclic unit.

REIDEMEISTER SPECTRUM FOR METABELIAN GROUPS OF THE FORM Q(n) x Z AND Z[1/p](n) x Z, p PRIME

In this paper, we study the Reidemeister spectrum for metabelian groups of the form Q(n) x Z and Z[1/p](n) x Z. Particular attention is given to the R(infinity)-property of a subfamily of these groups.