Trajectory Sensitivity Method and Master-Slave Synchronization to Estimate Parameters of Nonlinear Systems

A combination of trajectory sensitivity method and master-slave synchronization was proposed to parameter estimation of nonlinear systems. It was shown that master-slave coupling increases the robustness of the trajectory sensitivity algorithm with respect to the initial guess of parameters. Since synchronization is not a guarantee that the estimation process converges to the correct parameters, a conditional test that guarantees that the new combined methodology estimates the true values of parameters was proposed. This conditional test was successfully applied to Lorenz's and Chua's systems, and the proposed parameter estimation algorithm has shown to be very robust with respect to parameter initial guesses and measurement noise for these examples. Copyright (C) 2009 Elmer P. T. Cari et al.

Synchronization of Discrete-Time Chaotic Systems in Bandlimited Channels

Over the last couple of decades, many methods for synchronizing chaotic systems have been proposed with communications applications in view. Yet their performance has proved disappointing in face of the nonideal character of usual channels linking transmitter and receiver, that is, due to both noise and signal propagation distortion. Here we consider a discrete-time master-slave system that synchronizes despite channel bandwidth limitations and an allied communication system. Synchronization is achieved introducing a digital filter that limits the spectral content of the feedback loop responsible for producing the transmitted signal. Copyright (C) 2009 Marcio Eisencraft et al.

Geographic Information Systems (GIS) in Assessing Dental Health

The present study investigated the distribution profile of dental caries and its association with areas of social deprivation at the individual and contextual level. The cluster sample consisted of 1,002 12-year-old schoolchildren from Piracicaba, SP, Brazil. The DMFT Index was used for dental caries and the Care Index was used to determine access to dental services. On the individual level, variables were associated with a better oral status. On the contextual level, areas were not associated with oral status. However, maps enabled determining that the central districts have better social and oral conditions than the deprived outlying districts.

Socioeconomic Differences in Preferences and Willingness-to-Pay for Insulin Delivery Systems in Type 1 and Type 2 Diabetes

Background: An evaluation of patients' preferences is necessary to understand the demand for different insulin delivery systems. The aim of this study was to investigate the association between socioeconomic status (SES) and patients' preferences and willingness to pay (WTP) for various attributes of insulin administration for diabetes management. Methods: We conducted a discrete choice experiment (DCE) to determine patients' preferences and their WTP for hypothetical insulin treatments. Both self-reported annual household income and education completed were used to explore differences in treatment preferences and WTP for different attributes of treatment across different levels of SES. Results: The DCE questionnaire was successfully completed by 274 patients. Overall, glucose control was the most valued attribute by all socioeconomic groups, while route of insulin delivery was not as important. Patients with higher incomes were willing to pay significantly more for better glucose control and to avoid adverse events compared to lower income groups. In addition, they were willing to pay more for an oral short-acting insulin ($Can 71.65 [95% confidence interval, $40.68, $102.62]) compared to the low income group ($Can 9.85 [95% confidence interval, 14.86, 34.56; P < 0.01]). Conversely, there were no differences in preferences when the sample was stratified by level of education. Conclusions: This study revealed that preferences and WTP for insulin therapy are influenced by income but not by level of education. Specifically, the higher the income, the greater desire for an oral insulin delivery system, whereas an inhaled route becomes less important for patients.

Crossover between the extended and localized regimes in stochastic partially self-avoiding walks in one-dimensional disordered systems

Consider N sites randomly and uniformly distributed in a d-dimensional hypercube. A walker explores this disordered medium going to the nearest site, which has not been visited in the last mu (memory) steps. The walker trajectory is composed of a transient part and a periodic part (cycle). For one-dimensional systems, travelers can or cannot explore all available space, giving rise to a crossover between localized and extended regimes at the critical memory mu(1) = log(2) N. The deterministic rule can be softened to consider more realistic situations with the inclusion of a stochastic parameter T (temperature). In this case, the walker movement is driven by a probability density function parameterized by T and a cost function. The cost function increases as the distance between two sites and favors hops to closer sites. As the temperature increases, the walker can escape from cycles that are reminiscent of the deterministic nature and extend the exploration. Here, we report an analytical model and numerical studies of the influence of the temperature and the critical memory in the exploration of one-dimensional disordered systems.

Sensitization Prevalence, Antibody Cross-Reactivity and Immunogenic Peptide Profile of Api g 2, the Non-Specific Lipid Transfer Protein 1 of Celery

Background: Celery (Apium graveolens) represents a relevant allergen source that can elicit severe reactions in the adult population. To investigate the sensitization prevalence and cross-reactivity of Api g 2 from celery stalks in a Mediterranean population and in a mouse model. Methodology: 786 non-randomized subjects from Italy were screened for IgE reactivity to rApi g 2, rArt v 3 (mugwort pollen LTP) and nPru p 3 (peach LTP) using an allergen microarray. Clinical data of 32 selected patients with reactivity to LTP under investigation were evaluated. Specific IgE titers and cross-inhibitions were performed in ELISA and allergen microarray. Balb/c mice were immunized with purified LTPs; IgG titers were determined in ELISA and mediator release was examined using RBL-2H3 cells. Simulated endolysosomal digestion was performed using microsomes obtained from human DCs. Results: IgE testing showed a sensitization prevalence of 25.6% to Api g 2, 18.6% to Art v 3, and 28.6% to Pru p 3 and frequent co-sensitization and correlating IgE-reactivity was observed. 10/32 patients suffering from LTP-related allergy reported symptoms upon consumption of celery stalks which mainly presented as OAS. Considerable IgE cross-reactivity was observed between Api g 2, Art v 3, and Pru p 3 with varying inhibition degrees of individual patients' sera. Simulating LTP mono-sensitization in a mouse model showed development of more congruent antibody specificities between Api g 2 and Art v 3. Notably, biologically relevant murine IgE cross-reactivity was restricted to the latter and diverse from Pru p 3 epitopes. Endolysosomal processing of LTP showed generation of similar clusters, which presumably represent T-cell peptides. Conclusions: Api g 2 represents a relevant celery stalk allergen in the LTP-sensitized population. The molecule displays common B cell epitopes and endolysosomal peptides that encompass T cell epitopes with pollen and plant-food derived LTP.

Integration and visualization of systems biology data in context of the genome

Background: High-density tiling arrays and new sequencing technologies are generating rapidly increasing volumes of transcriptome and protein-DNA interaction data. Visualization and exploration of this data is critical to understanding the regulatory logic encoded in the genome by which the cell dynamically affects its physiology and interacts with its environment. Results: The Gaggle Genome Browser is a cross-platform desktop program for interactively visualizing high-throughput data in the context of the genome. Important features include dynamic panning and zooming, keyword search and open interoperability through the Gaggle framework. Users may bookmark locations on the genome with descriptive annotations and share these bookmarks with other users. The program handles large sets of user-generated data using an in-process database and leverages the facilities of SQL and the R environment for importing and manipulating data. A key aspect of the Gaggle Genome Browser is interoperability. By connecting to the Gaggle framework, the genome browser joins a suite of interconnected bioinformatics tools for analysis and visualization with connectivity to major public repositories of sequences, interactions and pathways. To this flexible environment for exploring and combining data, the Gaggle Genome Browser adds the ability to visualize diverse types of data in relation to its coordinates on the genome. Conclusions: Genomic coordinates function as a common key by which disparate biological data types can be related to one another. In the Gaggle Genome Browser, heterogeneous data are joined by their location on the genome to create information-rich visualizations yielding insight into genome organization, transcription and its regulation and, ultimately, a better understanding of the mechanisms that enable the cell to dynamically respond to its environment.

Continental passive microwave-based rainfall estimation algorithm: Application to the Amazon Basin

This paper presents a new statistical algorithm to estimate rainfall over the Amazon Basin region using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm relies on empirical relationships derived for different raining-type systems between coincident measurements of surface rainfall rate and 85-GHz polarization-corrected brightness temperature as observed by the precipitation radar (PR) and TMI on board the TRMM satellite. The scheme includes rain/no-rain area delineation (screening) and system-type classification routines for rain retrieval. The algorithm is validated against independent measurements of the TRMM-PR and S-band dual-polarization Doppler radar (S-Pol) surface rainfall data for two different periods. Moreover, the performance of this rainfall estimation technique is evaluated against well-known methods, namely, the TRMM-2A12 [ the Goddard profiling algorithm (GPROF)], the Goddard scattering algorithm (GSCAT), and the National Environmental Satellite, Data, and Information Service (NESDIS) algorithms. The proposed algorithm shows a normalized bias of approximately 23% for both PR and S-Pol ground truth datasets and a mean error of 0.244 mm h(-1) ( PR) and -0.157 mm h(-1)(S-Pol). For rain volume estimates using PR as reference, a correlation coefficient of 0.939 and a normalized bias of 0.039 were found. With respect to rainfall distributions and rain area comparisons, the results showed that the formulation proposed is efficient and compatible with the physics and dynamics of the observed systems over the area of interest. The performance of the other algorithms showed that GSCAT presented low normalized bias for rain areas and rain volume [0.346 ( PR) and 0.361 (S-Pol)], and GPROF showed rainfall distribution similar to that of the PR and S-Pol but with a bimodal distribution. Last, the five algorithms were evaluated during the TRMM-Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) 1999 field campaign to verify the precipitation characteristics observed during the easterly and westerly Amazon wind flow regimes. The proposed algorithm presented a cumulative rainfall distribution similar to the observations during the easterly regime, but it underestimated for the westerly period for rainfall rates above 5 mm h(-1). NESDIS(1) overestimated for both wind regimes but presented the best westerly representation. NESDIS(2), GSCAT, and GPROF underestimated in both regimes, but GPROF was closer to the observations during the easterly flow.

MULTILINE DOPPLER IMAGING OF MR SER IN HIGH STATE

Doppler images in Balmer, He I, He II and C II lines, and simultaneous I-band photometry of the polar MR Ser are presented and analyzed. The Balmer and Helium Doppler tomograms, of this bright polar at high mass transfer state show the emission from the accretion flow and the heated surface of the companion star. As a result of a comparison between the Doppler tomograms, the ionization structure of the flow could be constrained. The highest ionization region was found in the vicinity of the magnetospheric radius. Photoionization modeling of the accretion column indicates that the Balmer and Helium emission line production in this system can be explained only by the central soft X-ray illumination. The orbital ephemeris of MR Ser has been revised.

A NEW DIAGNOSTIC OF THE RADIAL DENSITY STRUCTURE OF Be DISKS

We analyze the intrinsic polarization of two classical Be stars in the process of losing their circumstellar disks via a Be to normal B star transition originally reported by Wisniewski et al. During each of five polarimetric outbursts which interrupt these disk-loss events, we find that the ratio of the polarization across the Balmer jump (BJ+/BJ-) versus the V-band polarization traces a distinct loop structure as a function of time. Since the polarization change across the Balmer jump is a tracer of the innermost disk density whereas the V-band polarization is a tracer of the total scattering mass of the disk, we suggest that such correlated loop structures in Balmer jump-V-band polarization diagrams (BJV diagrams) provide a unique diagnostic of the radial distribution of mass within Be disks. We use the three-dimensional Monte Carlo radiation transfer code HDUST to reproduce the observed clockwise loops simply by turning ""on/off"" the mass decretion from the disk. We speculate that counterclockwise loop structures we observe in BJV diagrams might be caused by the mass decretion rate changing between subsequent ""on/off"" sequences. Applying this new diagnostic to a larger sample of Be disk systems will provide insight into the time-dependent nature of each system's stellar decretion rate.

The role of magnetic reconnection on jet/accretion disk systems

Context. It was proposed earlier that the relativistic ejections observed in microquasars could be produced by violent magnetic reconnection episodes at the inner disk coronal region (de Gouveia Dal Pino & Lazarian 2005). Aims. Here we revisit this model, which employs a standard accretion disk description and fast magnetic reconnection theory, and discuss the role of magnetic reconnection and associated heating and particle acceleration in different jet/disk accretion systems, namely young stellar objects (YSOs), microquasars, and active galactic nuclei (AGNs). Methods. In microquasars and AGNs, violent reconnection episodes between the magnetic field lines of the inner disk region and those that are anchored in the black hole are able to heat the coronal/disk gas and accelerate the plasma to relativistic velocities through a diffusive first-order Fermi-like process within the reconnection site that will produce intermittent relativistic ejections or plasmons. Results. The resulting power-law electron distribution is compatible with the synchrotron radio spectrum observed during the outbursts of these sources. A diagram of the magnetic energy rate released by violent reconnection as a function of the black hole (BH) mass spanning 10(9) orders of magnitude shows that the magnetic reconnection power is more than sufficient to explain the observed radio luminosities of the outbursts from microquasars to low luminous AGNs. In addition, the magnetic reconnection events cause the heating of the coronal gas, which can be conducted back to the disk to enhance its thermal soft X-ray emission as observed during outbursts in microquasars. The decay of the hard X-ray emission right after a radio flare could also be explained in this model due to the escape of relativistic electrons with the evolving jet outburst. In the case of YSOs a similar magnetic configuration can be reached that could possibly produce observed X-ray flares in some sources and provide the heating at the jet launching base, but only if violent magnetic reconnection events occur with episodic, very short-duration accretion rates which are similar to 100-1000 times larger than the typical average accretion rates expected for more evolved (T Tauri) YSOs.

The BRAzilian Seismographic Integrated Systems (BRASIS): infrastructure and data management

In geophysics and seismology, raw data need to be processed to generate useful information that can be turned into knowledge by researchers. The number of sensors that are acquiring raw data is increasing rapidly. Without good data management systems, more time can be spent in querying and preparing datasets for analyses than in acquiring raw data. Also, a lot of good quality data acquired at great effort can be lost forever if they are not correctly stored. Local and international cooperation will probably be reduced, and a lot of data will never become scientific knowledge. For this reason, the Seismological Laboratory of the Institute of Astronomy, Geophysics and Atmospheric Sciences at the University of Sao Paulo (IAG-USP) has concentrated fully on its data management system. This report describes the efforts of the IAG-USP to set up a seismology data management system to facilitate local and international cooperation.

Modeling Two-Oscillator Circadian Systems Entrained by Two Environmental Cycles

Several experimental studies have altered the phase relationship between photic and non-photic environmental, 24 h cycles (zeitgebers) in order to assess their role in the synchronization of circadian rhythms. To assist in the interpretation of the complex activity patterns that emerge from these ""conflicting zeitgeber'' protocols, we present computer simulations of coupled circadian oscillators forced by two independent zeitgebers. This circadian system configuration was first employed by Pittendrigh and Bruce (1959), to model their studies of the light and temperature entrainment of the eclosion oscillator in Drosophila. Whereas most of the recent experiments have restricted conflicting zeitgeber experiments to two experimental conditions, by comparing circadian oscillator phases under two distinct phase relationships between zeitgebers (usually 0 and 12 h), Pittendrigh and Bruce compared eclosion phase under 12 distinct phase relationships, spanning the 24 h interval. Our simulations using non-linear differential equations replicated complex non-linear phenomena, such as ""phase jumps'' and sudden switches in zeitgeber preferences, which had previously been difficult to interpret. Our simulations reveal that these phenomena generally arise when inter-oscillator coupling is high in relation to the zeitgeber strength. Manipulations in the structural symmetry of the model indicated that these results can be expected to apply to a wide range of system configurations. Finally, our studies recommend the use of the complete protocol employed by Pittendrigh and Bruce, because different system configurations can generate similar results when a ""conflicting zeitgeber experiment'' incorporates only two phase relationships between zeitgebers.

Activation of endogenous p53 by combined p19Arf gene transfer and nutlin-3 drug treatment modalities in the murine cell lines B16 and C6

Background: Reactivation of p53 by either gene transfer or pharmacologic approaches may compensate for loss of p19Arf or excess mdm2 expression, common events in melanoma and glioma. In our previous work, we constructed the pCLPG retroviral vector where transgene expression is controlled by p53 through a p53-responsive promoter. The use of this vector to introduce p19Arf into tumor cells that harbor p53wt should yield viral expression of p19Arf which, in turn, would activate the endogenous p53 and result in enhanced vector expression and tumor suppression. Since nutlin-3 can activate p53 by blocking its interaction with mdm2, we explored the possibility that the combination of p19Arf gene transfer and nutlin-3 drug treatment may provide an additive benefit in stimulating p53 function. Methods: B16 (mouse melanoma) and C6 (rat glioma) cell lines, which harbor p53wt, were transduced with pCLPGp19 and these were additionally treated with nutlin-3 or the DNA damaging agent, doxorubicin. Viral expression was confirmed by Western, Northern and immunofluorescence assays. p53 function was assessed by reporter gene activity provided by a p53-responsive construct. Alterations in proliferation and viability were measured by colony formation, growth curve, cell cycle and MTT assays. In an animal model, B16 cells were treated with the pCLPGp19 virus and/or drugs before subcutaneous injection in C57BL/6 mice, observation of tumor progression and histopathologic analyses. Results: Here we show that the functional activation of endogenous p53wt in B16 was particularly challenging, but accomplished when combined gene transfer and drug treatments were applied, resulting in increased transactivation by p53, marked cell cycle alteration and reduced viability in culture. In an animal model, B16 cells treated with both p19Arf and nutlin-3 yielded increased necrosis and decreased BrdU marking. In comparison, C6 cells were quite susceptible to either treatment, yet p53 was further activated by the combination of p19Arf and nutlin-3. Conclusions: To the best of our knowledge, this is the first study to apply both p19Arf and nutlin-3 for the stimulation of p53 activity. These results support the notion that a p53 responsive vector may prove to be an interesting gene transfer tool, especially when combined with p53- activating agents, for the treatment of tumors that retain wild-type p53.

SOLITON SOLUTIONS TO SYSTEMS OF COUPLED SCHRODINGER EQUATIONS OF HAMILTONIAN TYPE

We study the existence of positive solutions of Hamiltonian-type systems of second-order elliptic PDE in the whole space. The systems depend on a small parameter and involve a potential having a global well structure. We use dual variational methods, a mountain-pass type approach and Fourier analysis to prove positive solutions exist for sufficiently small values of the parameter.