Search

Posturography is a tool for quantitative measurement of balance in various static and dynamic test conditions. Current medical grade Posturography systems are very costly and rare outside research institutes. The Wii Balance Board is an inexpensive gaming device which works on similar principles and uses the center of pressure of an individual as an input for motion-controlled gaming. The aim of this study was to assess the validity of using a gaming platform (Wii Balance Board) as a substitute to medical grade Posturography machine. The objectives were to compare the measures of postural sway obtained from a gaming platform with those from a medical grade Posturography machine and to demonstrate a low-cost method of analysis of balance using the gaming platform. Descriptive Observational study. Individuals aged 18–65 years and having no apparent dysfunction of balance were assessed objectively using Wii Balance Board and Medical-grade Posturography machine. A program was developed in Scilab for computation of sway parameters. The results were compared in terms of correlation coefficient, Bland-Altmann plots, and ability to detect abnormal sway. A Pearson's correlation coefficient ranging from 0.53–0.91 was observed in various test conditions, with an overall value of 0.62 between the results of the two machines. The sensitivity was 69.39%, specificity was 73.16%, positive predictive value 24.29%, negative predictive value 95.07%, and accuracy 72.75%. Thus, we can conclude that the Wii Balance Board can be used as a low-cost substitute to medical grade Posturography machine for quantification of balance in situations where precise measurement of balance is not required. [ABSTRACT FROM AUTHOR]

Translating progress in neuroscience into clinical benefits for patients with psychiatric disorders is challenging because it involves the brain as the most complex organ and its interaction with a complex environment and condition. Dealing with such complexity requires powerful techniques. Computational neuroscience approach to psychiatry integrates multiple levels and types of simulation, analysis and computation according to the different types of computational models to enhance comprehending, prediction and treatment of psychiatric disorder. This approach comprises two approaches: theory-driven and data-driven. In this review, we focus on recent advances in theory-driven approaches that mathematically and mechanistically examine the relationships between disorder-related changes and behavior at different level of brain organization. We discuss recent progresses in computational neuroscience models that relate to psychiatry and show how principles of neural computational modeling can be employed to explain psychopathology. [ABSTRACT FROM AUTHOR]

Background: Circadian clock alterations were poorly reported in trauma patients, although they have a critical role in human physiology. Core body temperature is a clinical variable regulated by the circadian clock. Our objective was to identify the circadian temperature disruption in trauma patients and to determine whether these disruptions were associated with the 28-day mortality rate. Methods: A retrospective and observational single-center cohort study was conducted. All adult severe trauma patients admitted to the intensive care unit of Aix Marseille University, North Hospital, from November 2013 to February 2018, were evaluated. The variations of core body temperature for each patient were analyzed between days 2 and 3 after intensive care unit admission. Core body temperature variations were defined by three parameters: mesor, amplitude, and period. A logistic regression model was used to determine the variables influencing these three parameters. A survival analysis was performed assessing the association between core body temperature rhythm disruption and 28-day mortality rate. A post hoc subgroup analysis focused on the patients with head trauma. Results: Among the 1584 screened patients, 248 were included in this study. The period differed from 24 h in 177 (71%) patients. The mesor value (°C) was associated with body mass index and ketamine use. Amplitude (°C) was associated with ketamine use only. The 28-day mortality rate was 18%. For all trauma patients, age, body mass index, intracranial hypertension, and amplitude were independent risk factors. The patients with a mesor value < 36.9 °C (p < 0.001) and an amplitude > 0.6 °C (p < 0.001) had a higher 28-day mortality rate. Among the patients with head trauma, mesor and amplitude were identified as independent risk factors (HR = 0.40, 95% CI [0.23–0.70], p = 0.001 and HR = 4.73, 95% CI [1.38–16.22], p = 0.01). Conclusions: Our results highlight an association between core body temperature circadian alteration and 28-day mortality rate. This association was more pronounced in the head trauma patients than in the non-head trauma patients. Further studies are needed to show a causal link and consider possible interventions. [ABSTRACT FROM AUTHOR]

Mitotic apparatus, which comprises hundreds of microtubules, plays an essential role in cell division, ensuring the correct segregation of chromosomes into each daughter cell. To gain insight into its regulatory mechanisms, it is essential to detect and analyze the behavior of individual microtubule filaments. However, the discrimination of discrete microtubule filaments within the mitotic apparatus is beyond the capabilities of conventional light microscopic technologies. Recently, we detected three-dimensional (3-D) microtubule growth dynamics within the cellular cytoplasmic space using lattice light-sheet microscopy in conjunction with microtubule growth marker protein end-binding 1, a microtubule plus-end-tracking protein, which was fused to green fluorescent protein (EB1-GFP). This technique enables high-resolution 3-D imaging at subsecond intervals. We adapted mathematical computing and geometric representation techniques to analyze spatial variations in microtubule growth dynamics within the mitotic spindle apparatus. Our analytical approach enabled the different dynamic properties of individual microtubules to be determined, including the direction and speed of their growth, and their growth duration within a 3-D spatial map. Our analysis framework provides an important step toward a more comprehensive understanding of the mechanisms driving cellular machinery at the whole-cell level. [ABSTRACT FROM AUTHOR]

Objective Uncertainty evaluation of determination of vitamin B6 content in oyster soybean peptide carnitine oral liquid by high performance liquid chromatography (HPLC assay). Methods HPLC was adopted for the determination of vitamin B6 content. Mathematical model was established for the evaluation of uncertainty to quantize influential factors of uncertainty and evaluate the uncertainty factors further. Results The content of vitamin B6 in sample oral liquid was 0.498 mg/100mL, in which, the relative uncertainty of vitamin B6 content was 0.024 9, the standard uncertainty was 0.0124 mg/100 mL, the expanded uncertainty was 0.024 8 mg/100 mL(k=2). Conclusion The calculation method is practical and reliable, being access to the uncertainty evaluation of content of vitamin B6 in oyster soybean peptide carnitine oral liquid by high performance liquid chromatography. Upon the evaluation, it is found that the quantitative repetitive of the HPLC is the major impact. [ABSTRACT FROM AUTHOR]

A system for real-time ultrasound (US) elastography will advance interventions for the diagnosis and treatment of cancer by advancing methods such as thermal monitoring of tissue ablation. A multi-stream graphics processing unit (GPU) based accelerated normalized cross-correlation (NCC) elastography, with a maximum frame rate of 78 frames per second, is presented in this paper. A study of NCC window size is undertaken to determine the effect on frame rate and the quality of output elastography images. This paper also presents a novel system for Online Tracked Ultrasound Elastography (O-TRuE), which extends prior work on an offline method. By tracking the US probe with an electromagnetic (EM) tracker, the system selects in-plane radio frequency (RF) data frames for generating high quality elastograms. A novel method for evaluating the quality of an elastography output stream is presented, suggesting that O-TRuE generates more stable elastograms than generated by untracked, free-hand palpation. Since EM tracking cannot be used in all systems, an integration of real-time elastography and the da Vinci Surgical System is presented and evaluated for elastography stream quality based on our metric. The da Vinci surgical robot is outfitted with a laparoscopic US probe, and palpation motions are autonomously generated by customized software. It is found that a stable output stream can be achieved, which is affected by both the frequency and amplitude of palpation. The GPU framework is validated using data from in-vivo pig liver ablation; the generated elastography images identify the ablated region, outlined more clearly than in the corresponding B-mode US images. [ABSTRACT FROM AUTHOR]

Animals may respond to habitat quality and habitat edges and these responses may affect their distribution between habitats. We studied the movement behaviour of a ground-dwelling generalist predator, the carabid beetle Pterostichus melanarius (Illiger). We performed a mark-recapture experiment in two adjacent habitats; a large plot with oilseed radish (Raphanus sativus) and a plot with rye (Secale cereale). We used model selection to identify a minimal model representing the mark-recapture data, and determine whether habitat-specific motility and boundary behaviour affected population redistribution. We determined movement characteristics of P. melanarius in laboratory arenas with the same plant species using video recording. Both the field and arena results showed preference behaviour of P. melanarius at the habitat interface. In the field, significantly more beetles moved from rye to oilseed radish than from radish to rye. In the arena, habitat entry was more frequent into oilseed radish than into rye. In the field, movement was best described by a Fokker-Planck diffusion model that contained preference behaviour at the interface and did not account for habitat specific motility. Likewise, motility calculated from movement data using the Patlak model was not different between habitats in the arena studies. Motility (m2 d−1) calculated from behavioural data resulted in estimates that were similar to those determined in the field. Thus individual behaviour explained population redistribution in the field qualitatively as well as quantitatively. The findings provide a basis for evaluating movement within and across habitats in complex agricultural landscapes with multiple habitats and habitat interfaces. [ABSTRACT FROM AUTHOR]

Coexistence of apparently similar species remains an enduring paradox in ecology. Spatial structure has been predicted to enable coexistence even when population-level models predict competitive exclusion if it causes each species to limit its own population more than that of its competitor. Nevertheless, existing hypotheses conflict with regard to whether clustering favours or precludes coexistence. The spatial segregation hypothesis predicts that in clustered populations the frequency of intra-specific interactions will be increased, causing each species to be self-limiting. Alternatively, individuals of the same species might compete over greater distances, known as heteromyopia, breaking down clusters and opening space for a second species to invade. In this study we create an individual-based model in homogeneous two-dimensional space for two putative sessile species differing only in their demographic rates and the range and strength of their competitive interactions. We fully characterise the parameter space within which coexistence occurs beyond population-level predictions, thereby revealing a region of coexistence generated by a previously-unrecognised process which we term the triadic mechanism. Here coexistence occurs due to the ability of a second generation of offspring of the rarer species to escape competition from their ancestors. We diagnose the conditions under which each of three spatial coexistence mechanisms operates and their characteristic spatial signatures. Deriving insights from a novel metric — ecological pressure — we demonstrate that coexistence is not solely determined by features of the numerically-dominant species. This results in a common framework for predicting, given any pair of species and knowledge of the relevant parameters, whether they will coexist, the mechanism by which they will do so, and the resultant spatial pattern of the community. Spatial coexistence arises from complementary combinations of traits in each species rather than solely through self-limitation. [ABSTRACT FROM AUTHOR]

It is well known that the Neolithic transition spread across Europe at a speed of about 1 km/yr. This result has been previously interpreted as a range expansion of the Neolithic driven mainly by demic diffusion (whereas cultural diffusion played a secondary role). However, a long-standing problem is whether this value (1 km/yr) and its interpretation (mainly demic diffusion) are characteristic only of Europe or universal (i.e. intrinsic features of Neolithic transitions all over the world). So far Neolithic spread rates outside Europe have been barely measured, and Neolithic spread rates substantially faster than 1 km/yr have not been previously reported. Here we show that the transition from hunting and gathering into herding in southern Africa spread at a rate of about 2.4 km/yr, i.e. about twice faster than the European Neolithic transition. Thus the value 1 km/yr is not a universal feature of Neolithic transitions in the world. Resorting to a recent demic-cultural wave-of-advance model, we also find that the main mechanism at work in the southern African Neolithic spread was cultural diffusion (whereas demic diffusion played a secondary role). This is in sharp contrast to the European Neolithic. Our results further suggest that Neolithic spread rates could be mainly driven by cultural diffusion in cases where the final state of this transition is herding/pastoralism (such as in southern Africa) rather than farming and stockbreeding (as in Europe). [ABSTRACT FROM AUTHOR]

In this paper, a fractional complex transform (FCT) is used to convert the given fractional partial differential equations (FPDEs) into corresponding partial differential equations (PDEs) and subsequently Reduced Differential Transform Method (RDTM) is applied on the transformed system of linear and nonlinear time-fractional PDEs. The results so obtained are re-stated by making use of inverse transformation which yields it in terms of original variables. It is observed that the proposed algorithm is highly efficient and appropriate for fractional PDEs and hence can be extended to other complex problems of diversified nonlinear nature. [ABSTRACT FROM AUTHOR]

Reconstructing regulatory and signaling response networks is one of the major goals of systems biology. While several successful methods have been suggested for this task, some integrating large and diverse datasets, these methods have so far been applied to reconstruct a single response network at a time, even when studying and modeling related conditions. To improve network reconstruction we developed MT-SDREM, a multi-task learning method which jointly models networks for several related conditions. In MT-SDREM, parameters are jointly constrained across the networks while still allowing for condition-specific pathways and regulation. We formulate the multi-task learning problem and discuss methods for optimizing the joint target function. We applied MT-SDREM to reconstruct dynamic human response networks for three flu strains: H1N1, H5N1 and H3N2. Our multi-task learning method was able to identify known and novel factors and genes, improving upon prior methods that model each condition independently. The MT-SDREM networks were also better at identifying proteins whose removal affects viral load indicating that joint learning can still lead to accurate, condition-specific, networks. Supporting website with MT-SDREM implementation: [ABSTRACT FROM AUTHOR]

Background: Predicting species’ potential geographical range by species distribution models (SDMs) is central to understand their ecological requirements. However, the effects of using different modeling techniques need further investigation. In order to improve the prediction effect, we need to assess the predictive performance and stability of different SDMs. Methodology: We collected the distribution data of five common tree species (Pinus massoniana, Betula platyphylla, Quercus wutaishanica, Quercus mongolica and Quercus variabilis) and simulated their potential distribution area using 13 environmental variables and six widely used SDMs: BIOCLIM, DOMAIN, MAHAL, RF, MAXENT, and SVM. Each model run was repeated 100 times (trials). We compared the predictive performance by testing the consistency between observations and simulated distributions and assessed the stability by the standard deviation, coefficient of variation, and the 99% confidence interval of Kappa and AUC values. Results: The mean values of AUC and Kappa from MAHAL, RF, MAXENT, and SVM trials were similar and significantly higher than those from BIOCLIM and DOMAIN trials (p<0.05), while the associated standard deviations and coefficients of variation were larger for BIOCLIM and DOMAIN trials (p<0.05), and the 99% confidence intervals for AUC and Kappa values were narrower for MAHAL, RF, MAXENT, and SVM. Compared to BIOCLIM and DOMAIN, other SDMs (MAHAL, RF, MAXENT, and SVM) had higher prediction accuracy, smaller confidence intervals, and were more stable and less affected by the random variable (randomly selected pseudo-absence points). Conclusions: According to the prediction performance and stability of SDMs, we can divide these six SDMs into two categories: a high performance and stability group including MAHAL, RF, MAXENT, and SVM, and a low performance and stability group consisting of BIOCLIM, and DOMAIN. We highlight that choosing appropriate SDMs to address a specific problem is an important part of the modeling process. [ABSTRACT FROM AUTHOR]

SLAM is a popular task used by robots and autonomous vehicles to build a map of an unknown environment and, at the same time, to determine their location within the map. This paper describes a SLAM-based, probabilistic robotic system able to learn the essential features of different parts of its environment. Some previous SLAM implementations had computational complexities ranging from O(Nlog(N)) to O(N2), where N is the number of map features. Unlike these methods, our approach reduces the computational complexity to O(N) by using a model to fuse the information from the sensors after applying the Bayesian paradigm. Once the training process is completed, the robot identifies and locates those areas that potentially match the sections that have been previously learned. After the training, the robot navigates and extracts a three-dimensional map of the environment using a single laser sensor. Thus, it perceives different sections of its world. In addition, in order to make our system able to be used in a low-cost robot, low-complexity algorithms that can be easily implemented on embedded processors or microcontrollers are used. [ABSTRACT FROM AUTHOR]

The steady two-dimensional flow and heat transfer over a stretching/shrinking sheet in a nanofluid is investigated using Buongiorno’s nanofluid model. Different from the previously published papers, in the present study we consider the case when the nanofluid particle fraction on the boundary is passively rather than actively controlled, which make the model more physically realistic. The governing partial differential equations are transformed into nonlinear ordinary differential equations by a similarity transformation, before being solved numerically by a shooting method. The effects of some governing parameters on the fluid flow and heat transfer characteristics are graphically presented and discussed. Dual solutions are found to exist in a certain range of the suction and stretching/shrinking parameters. Results also indicate that both the skin friction coefficient and the local Nusselt number increase with increasing values of the suction parameter. [ABSTRACT FROM AUTHOR]

A method of calculating inductances based on first principles is presented, which has the advantage over the more popular simulators in that fundamental formulas are explicitly used so that a deeper understanding of the inductance calculation is obtained with no need for explicit discretization of the inductor. It also has the advantage over the traditional method of formulas or table lookups in that it can be used for a wider range of configurations. It relies on the use of fast computers with a sophisticated mathematical computing language such as Mathematica to perform the required integration numerically so that the researcher can focus on the physics of the inductance calculation and not on the numerical integration. [ABSTRACT FROM AUTHOR]

Protein carbonylation is one of the most pervasive oxidative stress-induced post-translational modifications (PTMs), which plays a significant role in the etiology and progression of several human diseases. It has been regarded as a biomarker of oxidative stress due to its relatively early formation and stability compared with other oxidative PTMs. Only a subset of proteins is prone to carbonylation and most carbonyl groups are formed from lysine (K), arginine (R), threonine (T) and proline (P) residues. Recent advancements in analysis of the PTM by mass spectrometry provided new insights into the mechanisms of protein carbonylation, such as protein susceptibility and exact modification sites. However, the experimental approaches to identifying carbonylation sites are costly, time-consuming and capable of processing a limited number of proteins, and there is no bioinformatics method or tool devoted to predicting carbonylation sites of human proteins so far. In the paper, a computational method is proposed to identify carbonylation sites of human proteins. The method extracted four kinds of features and combined the minimum Redundancy Maximum Relevance (mRMR) feature selection criterion with weighted support vector machine (WSVM) to achieve total accuracies of 85.72%, 85.95%, 83.92% and 85.72% for K, R, T and P carbonylation site predictions respectively using 10-fold cross-validation. The final optimal feature sets were analysed, the position-specific composition and hydrophobicity environment of flanking residues of modification sites were discussed. In addition, a software tool named CarSPred has been developed to facilitate the application of the method. Datasets and the software involved in the paper are available at . [ABSTRACT FROM AUTHOR]

In this paper we analyse the word frequency profiles of a set of works from the Shakespearean era to uncover patterns of relationship between them, highlighting the connections within authorial canons. We used a text corpus comprising 256 plays and poems from the 16th and 17th centuries, with 17 works of uncertain authorship. Our clustering approach is based on the Jensen-Shannon divergence and a graph partitioning algorithm, and our results show that authors' characteristic styles are very powerful factors in explaining the variation of word use, frequently transcending cross-cutting factors like the differences between tragedy and comedy, early and late works, and plays and poems. Our method also provides an empirical guide to the authorship of plays and poems where this is unknown or disputed. [ABSTRACT FROM AUTHOR]

Purpose: Compressed sensing (CS) has been widely applied to prospective cardiac cine MRI. The aim of this work is to study the benefits obtained by including motion estimation in the CS framework for small-animal retrospective cardiac cine. Methods: We propose a novel B-spline-based compressed sensing method (SPLICS) that includes motion estimation and generalizes previous spatiotemporal total variation (ST-TV) methods by taking into account motion between frames. In addition, we assess the effect of an optimum weighting between spatial and temporal sparsity to further improve results. Both methods were implemented using the efficient Split Bregman methodology and were evaluated on rat data comparing animals with myocardial infarction with controls for several acceleration factors. Results: ST-TV with optimum selection of the weighting sparsity parameter led to results similar to those of SPLICS; ST-TV with large relative temporal sparsity led to temporal blurring effects. However, SPLICS always properly corrected temporal blurring, independently of the weighting parameter. At acceleration factors of 15, SPLICS did not distort temporal intensity information but led to some artefacts and slight over-smoothing. At an acceleration factor of 7, images were reconstructed without significant loss of quality. Conclusion: We have validated SPLICS for retrospective cardiac cine in small animal, achieving high acceleration factors. In addition, we have shown that motion modelling may not be essential for retrospective cine and that similar results can be obtained by using ST-TV provided that an optimum selection of the spatiotemporal sparsity weighting parameter is performed. [ABSTRACT FROM AUTHOR]

Parameter identifiability problems can plague biomodelers when they reach the quantification stage of development, even for relatively simple models. Structural identifiability (SI) is the primary question, usually understood as knowing which of P unknown biomodel parameters p1,…, pi,…, pP are-and which are not-quantifiable in principle from particular input-output (I-O) biodata. It is not widely appreciated that the same database also can provide quantitative information about the structurally identifiable ( quantifiable) subset, in the form of explicit algebraic relationships among unidentifiable pi. Importantly, this is a first step toward finding what else is needed to quantify particular identifiable parameters of interest from new I–O experiments. We further develop, implement and exemplify novel algorithms that address and solve the SI problem for a practical class of ordinary differential equation (ODE) systems biology models, as a user-friendly and universally-accessible web application (app)–COMBOS. Users provide the structural ODE and output measurement models in one of two standard forms to a remote server via their web browser. COMBOS provides a list of uniquely and non-uniquely SI model parameters, and–importantly-the combinations of parameters not individually SI. If non-uniquely SI, it also provides the maximum number of different solutions, with important practical implications. The behind-the-scenes symbolic differential algebra algorithms are based on computing Gröbner bases of model attributes established after some algebraic transformations, using the computer-algebra system Maxima. COMBOS was developed for facile instructional and research use as well as modeling. We use it in the classroom to illustrate SI analysis; and have simplified complex models of tumor suppressor p53 and hormone regulation, based on explicit computation of parameter combinations. It’s illustrated and validated here for models of moderate complexity, with and without initial conditions. Built-in examples include unidentifiable 2 to 4-compartment and HIV dynamics models. [ABSTRACT FROM AUTHOR]

Finding solutions to the classical transportation problem is of great importance, since this optimization problem arises in many engineering and computer science applications. Especially the Earth Mover's Distance is used in a plethora of applications ranging from content-based image retrieval, shape matching, fingerprint recognition, object tracking and phishing web page detection to computing color differences in linguistics and biology. Our starting point is the well-known revised simplex algorithm, which iteratively improves a feasible solution to optimality. that we propose substantially reduces the number of candidates inspected for improving the solution, while at the same time balancing the number of pivots required. Tests on simulated benchmarks demonstrate a considerable reduction in computation time for the new method as compared to the usual revised simplex algorithm implemented with state-of-the-art initialization and pivot strategies. As a consequence, the Shortlist Method facilitates the computation of large scale transportation problems in viable time. In addition we describe a novel method for finding an initial feasible solution which we coin Modified Russell's Method. [ABSTRACT FROM AUTHOR]

An aggregate signature scheme allows anyone to compress multiple individual signatures from various users into a single compact signature. The main objective of such a scheme is to reduce the costs on storage, communication and computation. However, among existing aggregate signature schemes in the identity-based setting, some of them fail to achieve constant-length aggregate signature or require a large amount of pairing operations which grows linearly with the number of signers, while others have some limitations on the aggregated signatures. The main challenge in building efficient aggregate signature scheme is to compress signatures into a compact, constant-length signature without any restriction. To address the above drawbacks, by using the bilinear pairings, we propose an efficient unrestricted identity-based aggregate signature. Our scheme achieves both full aggregation and constant pairing computation. We prove that our scheme has existential unforgeability under the computational Diffie-Hellman assumption. [ABSTRACT FROM AUTHOR]

In this study combined heat and mass transfer by mixed convective flow along a moving vertical flat plate with hydrodynamic slip and thermal convective boundary condition is investigated. Using similarity variables, the governing nonlinear partial differential equations are converted into a system of coupled nonlinear ordinary differential equations. The transformed equations are then solved using a semi-numerical/analytical method called the differential transform method and results are compared with numerical results. Close agreement is found between the present method and the numerical method. Effects of the controlling parameters, including convective heat transfer, magnetic field, buoyancy ratio, hydrodynamic slip, mixed convective, Prandtl number and Schmidt number are investigated on the dimensionless velocity, temperature and concentration profiles. In addition effects of different parameters on the skin friction factor, , local Nusselt number, , and local Sherwood number are shown and explained through tables. [ABSTRACT FROM AUTHOR]

In this paper, the heat transfer effect on the unsteady boundary layer flow of a Casson fluid past an infinite oscillating vertical plate with Newtonian heating is investigated. The governing equations are transformed to a systems of linear partial differential equations using appropriate non-dimensional variables. The resulting equations are solved analytically by using the Laplace transform method and the expressions for velocity and temperature are obtained. They satisfy all imposed initial and boundary conditions and reduce to some well-known solutions for Newtonian fluids. Numerical results for velocity, temperature, skin friction and Nusselt number are shown in various graphs and discussed for embedded flow parameters. It is found that velocity decreases as Casson parameters increases and thermal boundary layer thickness increases with increasing Newtonian heating parameter. [ABSTRACT FROM AUTHOR]

Modelled as finite homogeneous Markov chains, probabilistic cellular automata with local transition probabilities in (0, 1) always posses a stationary distribution. This result alone is not very helpful when it comes to predicting the final configuration; one needs also a formula connecting the probabilities in the stationary distribution to some intrinsic feature of the lattice configuration. Previous results on the asynchronous cellular automata have showed that such feature really exists. It is the number of zero-one borders within the automaton's binary configuration. An exponential formula in the number of zero-one borders has been proved for the 1-D, 2-D and 3-D asynchronous automata with neighborhood three, five and seven, respectively. We perform computer experiments on a synchronous cellular automaton to check whether the empirical distribution obeys also that theoretical formula. The numerical results indicate a perfect fit for neighbourhood three and five, which opens the way for a rigorous proof of the formula in this new, synchronous case. [ABSTRACT FROM AUTHOR]

In this article, we have studied the flow and heat transfer in Sisko fluid with convective boundary condition over a non-isothermal stretching sheet. The flow is influenced by non-linearly stretching sheet in the presence of a uniform transverse magnetic field. The partial differential equations governing the problem have been reduced by similarity transformations into the ordinary differential equations. The transformed coupled ordinary differential equations are then solved analytically by using the homotopy analysis method (HAM) and numerically by the shooting method. Effects of different parameters like power-law index , magnetic parameter , stretching parameter , generalized Prandtl number Pr and generalized Biot number are presented graphically. It is found that temperature profile increases with the increasing value of and whereas it decreases for . Numerical values of the skin-friction coefficient and local Nusselt number are tabulated at various physical situations. In addition, a comparison between the HAM and exact solutions is also made as a special case and excellent agreement between results enhance a confidence in the HAM results. [ABSTRACT FROM AUTHOR]

Multi-compartmental models of neurons provide insight into the complex, integrative properties of dendrites. Because it is not feasible to experimentally determine the exact density and kinetics of each channel type in every neuronal compartment, an essential goal in developing models is to help characterize these properties. To address biological variability inherent in a given neuronal type, there has been a shift away from using hand-tuned models towards using ensembles or populations of models. In collectively capturing a neuron's output, ensemble modeling approaches uncover important conductance balances that control neuronal dynamics. However, conductances are never entirely known for a given neuron class in terms of its types, densities, kinetics and distributions. Thus, any multi-compartment model will always be incomplete. In this work, our main goal is to use ensemble modeling as an investigative tool of a neuron's biophysical balances, where the cycling between experiment and model is a design criterion from the start. We consider oriens-lacunosum/moleculare (O-LM) interneurons, a prominent interneuron subtype that plays an essential gating role of information flow in hippocampus. O-LM cells express the hyperpolarization-activated current (Ih). Although dendritic Ih could have a major influence on the integrative properties of O-LM cells, the compartmental distribution of Ih on O-LM dendrites is not known. Using a high-performance computing cluster, we generated a database of models that included those with or without dendritic Ih. A range of conductance values for nine different conductance types were used, and different morphologies explored. Models were quantified and ranked based on minimal error compared to a dataset of O-LM cell electrophysiological properties. Co-regulatory balances between conductances were revealed, two of which were dependent on the presence of dendritic Ih. These findings inform future experiments that differentiate between somatic and dendritic Ih, thereby continuing a cycle between model and experiment. [ABSTRACT FROM AUTHOR]

The 2D Wavelet-Transform Modulus Maxima (WTMM) method was used to detect microcalcifications (MC) in human breast tissue seen in mammograms and to characterize the fractal geometry of benign and malignant MC clusters. This was done in the context of a preliminary analysis of a small dataset, via a novel way to partition the wavelet-transform space-scale skeleton. For the first time, the estimated 3D fractal structure of a breast lesion was inferred by pairing the information from two separate 2D projected mammographic views of the same breast, i.e. the cranial-caudal (CC) and mediolateral-oblique (MLO) views. As a novelty, we define the “CC-MLO fractal dimension plot”, where a “fractal zone” and “Euclidean zones” (non-fractal) are defined. 118 images (59 cases, 25 malignant and 34 benign) obtained from a digital databank of mammograms with known radiologist diagnostics were analyzed to determine which cases would be plotted in the fractal zone and which cases would fall in the Euclidean zones. 92% of malignant breast lesions studied (23 out of 25 cases) were in the fractal zone while 88% of the benign lesions were in the Euclidean zones (30 out of 34 cases). Furthermore, a Bayesian statistical analysis shows that, with 95% credibility, the probability that fractal breast lesions are malignant is between 74% and 98%. Alternatively, with 95% credibility, the probability that Euclidean breast lesions are benign is between 76% and 96%. These results support the notion that the fractal structure of malignant tumors is more likely to be associated with an invasive behavior into the surrounding tissue compared to the less invasive, Euclidean structure of benign tumors. Finally, based on indirect 3D reconstructions from the 2D views, we conjecture that all breast tumors considered in this study, benign and malignant, fractal or Euclidean, restrict their growth to 2-dimensional manifolds within the breast tissue. [ABSTRACT FROM AUTHOR]

Chronic diseases impose a tremendous global health problem of the 21st century. Epidemiological and public health models help to gain insight into the distribution and burden of chronic diseases. Moreover, the models may help to plan appropriate interventions against risk factors. To provide accurate results, models often need to take into account three different time-scales: calendar time, age, and duration since the onset of the disease. Incidence and mortality often change with age and calendar time. In many diseases such as, for example, diabetes and dementia, the mortality of the diseased persons additionally depends on the duration of the disease. The aim of this work is to describe an algorithm and a flexible software framework for the simulation of populations moving in an illness-death model that describes the epidemiology of a chronic disease in the face of the different times-scales. We set up a discrete event simulation in continuous time involving competing risks using the freely available statistical software R. Relevant events are birth, the onset (or diagnosis) of the disease and death with or without the disease. The Lexis diagram keeps track of the different time-scales. Input data are birth rates, incidence and mortality rates, which can be given as numerical values on a grid. The algorithm manages the complex interplay between the rates and the different time-scales. As a result, for each subject in the simulated population, the algorithm provides the calendar time of birth, the age of onset of the disease (if the subject contracts the disease) and the age at death. By this means, the impact of interventions may be estimated and compared. [ABSTRACT FROM AUTHOR]

In computerized adaptive testing (CAT), examinees are presented with various sets of items chosen from a precalibrated item pool. Consequently, the attrition speed of the items is extremely fast, and replenishing the item pool is essential. Therefore, item calibration has become a crucial concern in maintaining item banks. In this study, a two-parameter logistic model is used. We applied optimal designs and adaptive sequential analysis to solve this item calibration problem. The results indicated that the proposed optimal designs are cost effective and time efficient. [ABSTRACT FROM AUTHOR]

We present a novel approach for analyzing biological time-series data using a context-free language (CFL) representation that allows the extraction and quantification of important features from the time-series. This representation results in Hierarchically AdaPtive (HAP) analysis, a suite of multiple complementary techniques that enable rapid analysis of data and does not require the user to set parameters. HAP analysis generates hierarchically organized parameter distributions that allow multi-scale components of the time-series to be quantified and includes a data analysis pipeline that applies recursive analyses to generate hierarchically organized results that extend traditional outcome measures such as pharmacokinetics and inter-pulse interval. Pulsicons, a novel text-based time-series representation also derived from the CFL approach, are introduced as an objective qualitative comparison nomenclature. We apply HAP to the analysis of 24 hours of frequently sampled pulsatile cortisol hormone data, which has known analysis challenges, from 14 healthy women. HAP analysis generated results in seconds and produced dozens of figures for each participant. The results quantify the observed qualitative features of cortisol data as a series of pulse clusters, each consisting of one or more embedded pulses, and identify two ultradian phenotypes in this dataset. HAP analysis is designed to be robust to individual differences and to missing data and may be applied to other pulsatile hormones. Future work can extend HAP analysis to other time-series data types, including oscillatory and other periodic physiological signals. [ABSTRACT FROM AUTHOR]

A recently proposed methodology called the Horizontal Visibility Graph (HVG) [Luque et al., Phys. Rev. E., 80, 046103 (2009)] that constitutes a geometrical simplification of the well known Visibility Graph algorithm [Lacasa et al., Proc. Natl. Sci. U.S.A. 105, 4972 (2008)], has been used to study the distinction between deterministic and stochastic components in time series [L. Lacasa and R. Toral, Phys. Rev. E., 82, 036120 (2010)]. Specifically, the authors propose that the node degree distribution of these processes follows an exponential functional of the form , in which is the node degree and is a positive parameter able to distinguish between deterministic (chaotic) and stochastic (uncorrelated and correlated) dynamics. In this work, we investigate the characteristics of the node degree distributions constructed by using HVG, for time series corresponding to chaotic maps, 2 chaotic flows and different stochastic processes. We thoroughly study the methodology proposed by Lacasa and Toral finding several cases for which their hypothesis is not valid. We propose a methodology that uses the HVG together with Information Theory quantifiers. An extensive and careful analysis of the node degree distributions obtained by applying HVG allow us to conclude that the Fisher-Shannon information plane is a remarkable tool able to graphically represent the different nature, deterministic or stochastic, of the systems under study. [ABSTRACT FROM AUTHOR]

This article studies the viscous flow and heat transfer over a plane horizontal surface stretched non-linearly in two lateral directions. Appropriate wall conditions characterizing the non-linear variation in the velocity and temperature of the sheet are employed for the first time. A new set of similarity variables is introduced to reduce the boundary layer equations into self-similar forms. The velocity and temperature distributions are determined by two methods, namely (i) optimal homotopy analysis method (OHAM) and (ii) fourth-fifth-order Runge-Kutta integration based shooting technique. The analytic and numerical solutions are compared and these are found in excellent agreement. Influences of embedded parameters on momentum and thermal boundary layers are sketched and discussed. [ABSTRACT FROM AUTHOR]

We quantify the notion of pattern and formalize the process of pattern discovery under the framework of binary bipartite networks. Patterns of particular focus are interrelated global interactions between clusters on its row and column axes. A binary bipartite network is built into a thermodynamic system embracing all up-and-down spin configurations defined by product-permutations on rows and columns. This system is equipped with its ferromagnetic energy ground state under Ising model potential. Such a ground state, also called a macrostate, is postulated to congregate all patterns of interest embedded within the network data in a multiscale fashion. A new computing paradigm for indirect searching for such a macrostate, called Data Mechanics, is devised by iteratively building a surrogate geometric system with a pair of nearly optimal marginal ultrametrics on row and column spaces. The coupling measure minimizing the Gromov-Wasserstein distance of these two marginal geometries is also seen to be in the vicinity of the macrostate. This resultant coupling geometry reveals multiscale block pattern information that characterizes multiple layers of interacting relationships between clusters on row and on column axes. It is the nonparametric information content of a binary bipartite network. This coupling geometry is then demonstrated to shed new light and bring resolution to interaction issues in community ecology and in gene-content-based phylogenetics. Its implied global inferences are expected to have high potential in many scientific areas. [ABSTRACT FROM AUTHOR]

This article addresses the steady three-dimensional flow of an Oldroyd-B nanofluid over a bidirectional stretching surface with heat generation/absorption effects. Suitable similarity transformations are employed to reduce the governing partial differential equations into coupled nonlinear ordinary differential equations. These nonlinear ordinary differential equations are then solved analytically by using the homotpy analysis method (HAM). Graphically results are presented and discussed for various parameters, namely, Deborah numbers and , heat generation/absorption parameter Prandtl parameter , Brownian motion parameters, thermophoresis parameter and Lewis number . We have seen that the increasing values of the Brownian motion parameter and thermophoresis parameter leads to an increase in the temperature field and thermal boundary layer thickness while the opposite behavior is observed for concentration field and concentration boundary layer thickness. To see the validity of the present work, the numerical results are compared with the analytical solutions obtained by Homotopy analysis method and noted an excellent agreement for the limiting cases. [ABSTRACT FROM AUTHOR]

In this paper, a Semantic Web service matchmaker called UltiMatch-NL is presented. UltiMatch-NL applies two filters namely Signature-based and Description-based on different abstraction levels of a service profile to achieve more accurate results. More specifically, the proposed filters rely on semantic knowledge to extract the similarity between a given pair of service descriptions. Thus it is a further step towards fully automated Web service discovery via making this process more semantic-aware. In addition, a new technique is proposed to weight and combine the results of different filters of UltiMatch-NL, automatically. Moreover, an innovative approach is introduced to predict the relevance of requests and Web services and eliminate the need for setting a threshold value of similarity. In order to evaluate UltiMatch-NL, the repository of OWLS-TC is used. The performance evaluation based on standard measures from the information retrieval field shows that semantic matching of OWL-S services can be significantly improved by incorporating designed matching filters. [ABSTRACT FROM AUTHOR]

We develop a detailed theoretical framework for various types of transcription factor gene oscillators. We further demonstrate that one can build genetic-oscillators which are tunable and robust against perturbations in the critical control parameters by coupling two or more independent Goodwin-Griffith oscillators through either -OR- or -AND- type logic. Most of the coupled oscillators constructed in the literature so far seem to be of -OR- type. When there are transient perturbations in one of the -OR- type coupled-oscillators, then the overall period of the system remains constant (period-buffering) whereas in case of -AND- type coupling the overall period of the system moves towards the perturbed oscillator. Though there is a period-buffering, the amplitudes of oscillators coupled through -OR- type logic are more sensitive to perturbations in the parameters associated with the promoter state dynamics than -AND- type. Further analysis shows that the period of -AND- type coupled dual-feedback oscillators can be tuned without conceding on the amplitudes. Using these results we derive the basic design principles governing the robust and tunable synthetic gene oscillators without compromising on their amplitudes. [ABSTRACT FROM AUTHOR]

An individual-tree diameter growth model was developed for Cunninghamia lanceolata in Fujian province, southeast China. Data were obtained from 72 plantation-grown China-fir trees in 24 single-species plots. Ordinary non-linear least squares regression was used to choose the best base model from among 5 theoretical growth equations; selection criteria were the smallest absolute mean residual and root mean square error and the largest adjusted coefficient of determination. To account for autocorrelation in the repeated-measures data, we developed one-level and nested two-level nonlinear mixed-effects (NLME) models, constructed on the selected base model; the NLME models incorporated random effects of the tree and plot. The best random-effects combinations for the NLME models were identified by Akaike's information criterion, Bayesian information criterion and −2 logarithm likelihood. Heteroscedasticity was reduced with two residual variance functions, a power function and an exponential function. The autocorrelation was addressed with three residual autocorrelation structures: a first-order autoregressive structure [AR(1)], a combination of first-order autoregressive and moving average structures [ARMA(1,1)] and a compound symmetry structure (CS). The one-level (tree) NLME model performed best. Independent validation data were used to test the performance of the models and to demonstrate the advantage of calibrating the NLME models. [ABSTRACT FROM AUTHOR]

It is a classic topic of social network analysis to evaluate the importance of nodes and identify the node that takes on the role of core or bridge in a network. Because a single indicator is not sufficient to analyze multiple characteristics of a node, it is a natural solution to apply multiple indicators that should be selected carefully. An intuitive idea is to select some indicators with weak correlations to efficiently assess different characteristics of a node. However, this paper shows that it is much better to select the indicators with strong correlations. Because indicator correlation is based on the statistical analysis of a large number of nodes, the particularity of an important node will be outlined if its indicator relationship doesn't comply with the statistical correlation. Therefore, the paper selects the multiple indicators including degree, ego-betweenness centrality and eigenvector centrality to evaluate the importance and the role of a node. The importance of a node is equal to the normalized sum of its three indicators. A candidate for core or bridge is selected from the great degree nodes or the nodes with great ego-betweenness centrality respectively. Then, the role of a candidate is determined according to the difference between its indicators' relationship with the statistical correlation of the overall network. Based on 18 real networks and 3 kinds of model networks, the experimental results show that the proposed methods perform quite well in evaluating the importance of nodes and in identifying the node role. [ABSTRACT FROM AUTHOR]

In this paper a novel application of a particular type of spiking neural network, a Polychronous Spiking Network, was used for financial time series prediction. It is argued that the inherent temporal capabilities of this type of network are suited to non-stationary data such as this. The performance of the spiking neural network was benchmarked against three systems: two “traditional”, rate-encoded, neural networks; a Multi-Layer Perceptron neural network and a Dynamic Ridge Polynomial neural network, and a standard Linear Predictor Coefficients model. For this comparison three non-stationary and noisy time series were used: IBM stock data; US/Euro exchange rate data, and the price of Brent crude oil. The experiments demonstrated favourable prediction results for the Spiking Neural Network in terms of Annualised Return and prediction error for 5-Step ahead predictions. These results were also supported by other relevant metrics such as Maximum Drawdown and Signal-To-Noise ratio. This work demonstrated the applicability of the Polychronous Spiking Network to financial data forecasting and this in turn indicates the potential of using such networks over traditional systems in difficult to manage non-stationary environments. [ABSTRACT FROM AUTHOR]

The automatic detection of bilateral symmetry is a challenging task in computer vision and pattern recognition. This paper presents an approach for the detection of bilateral symmetry in digital single object images. Our method relies on the extraction of Scale Invariant Feature Transform (SIFT) based feature points, which serves as the basis for the ascertainment of the centroid of the object; the latter being the origin under the Cartesian coordinate system to be converted to the polar coordinate system in order to facilitate the selection symmetric coordinate pairs. This is followed by comparing the gradient magnitude and orientation of the corresponding points to evaluate the amount of symmetry exhibited by each pair of points. The experimental results show that our approach draw the symmetry line accurately, provided that the observed centroid point is true. [ABSTRACT FROM AUTHOR]

Electroporation-based therapies are powerful biotechnological tools for enhancing the delivery of exogeneous agents or killing tissue with pulsed electric fields (PEFs). Electrochemotherapy (ECT) and gene therapy based on gene electrotransfer (EGT) both use reversible electroporation to deliver chemotherapeutics or plasmid DNA into cells, respectively. In both ECT and EGT, the goal is to permeabilize the cell membrane while maintaining high cell viability in order to facilitate drug or gene transport into the cell cytoplasm and induce a therapeutic response. Irreversible electroporation (IRE) results in cell kill due to exposure to PEFs without drugs and is under clinical evaluation for treating otherwise unresectable tumors. These PEF therapies rely mainly on the electric field distributions and do not require changes in tissue temperature for their effectiveness. However, in immediate vicinity of the electrodes the treatment may results in cell kill due to thermal damage because of the inhomogeneous electric field distribution and high current density during the electroporation-based therapies. Therefore, the main objective of this numerical study is to evaluate the influence of pulse number and electrical conductivity in the predicted cell kill zone due to irreversible electroporation and thermal damage. Specifically, we simulated a typical IRE protocol that employs ninety 100-µs PEFs. Our results confirm that it is possible to achieve predominant cell kill due to electroporation if the PEF parameters are chosen carefully. However, if either the pulse number and/or the tissue conductivity are too high, there is also potential to achieve cell kill due to thermal damage in the immediate vicinity of the electrodes. Therefore, it is critical for physicians to be mindful of placement of electrodes with respect to critical tissue structures and treatment parameters in order to maintain the non-thermal benefits of electroporation and prevent unnecessary damage to surrounding healthy tissue, critical vascular structures, and/or adjacent organs. [ABSTRACT FROM AUTHOR]

Considering that many natural stimuli are sparse, can a sensory system evolve to take advantage of this sparsity? We explore this question and show that significant downstream reductions in the numbers of neurons transmitting stimuli observed in early sensory pathways might be a consequence of this sparsity. First, we model an early sensory pathway using an idealized neuronal network comprised of receptors and downstream sensory neurons. Then, by revealing a linear structure intrinsic to neuronal network dynamics, our work points to a potential mechanism for transmitting sparse stimuli, related to compressed-sensing (CS) type data acquisition. Through simulation, we examine the characteristics of networks that are optimal in sparsity encoding, and the impact of localized receptive fields beyond conventional CS theory. The results of this work suggest a new network framework of signal sparsity, freeing the notion from any dependence on specific component-space representations. We expect our CS network mechanism to provide guidance for studying sparse stimulus transmission along realistic sensory pathways as well as engineering network designs that utilize sparsity encoding. [ABSTRACT FROM AUTHOR]

In a previous work by Alvarez-Martinez et al. (2011), the authors pointed out some fallacies in the mainstream interpretation of the prion amyloid formation. It appeared necessary to propose an original hypothesis able to reconcile the in vitro data with the predictions of a mathematical model describing the problem. Here, a model is developed accordingly with the hypothesis that an intermediate on-pathway leads to the conformation of the prion protein into an amyloid competent isoform thanks to a structure, called micelles, formed from hydrodynamic interaction. The authors also compare data to the prediction of their model and propose a new hypothesis for the formation of infectious prion amyloids. [ABSTRACT FROM AUTHOR]

: A small molecule inhibitor of the malarial protease Plasmepsin V impairs protein export and cellular remodeling, reducing parasite survival in human erythrocytes. [ABSTRACT FROM AUTHOR]

The ancient Greek astronomical calculating machine, known as the Antikythera Mechanism, predicted eclipses, based on the 223-lunar month Saros cycle. Eclipses are indicated on a four-turn spiral Saros Dial by glyphs, which describe type and time of eclipse and include alphabetical index letters, referring to solar eclipse inscriptions. These include Index Letter Groups, describing shared eclipse characteristics. The grouping and ordering of the index letters, the organization of the inscriptions and the eclipse times have previously been unsolved. A new reading and interpretation of data from the back plate of the Antikythera Mechanism, including the glyphs, the index letters and the eclipse inscriptions, has resulted in substantial changes to previously published work. Based on these new readings, two arithmetical models are presented here that explain the complete eclipse prediction scheme. The first model solves the glyph distribution, the grouping and anomalous ordering of the index letters and the structure of the inscriptions. It also implies the existence of lost lunar eclipse inscriptions. The second model closely matches the glyph times and explains the four-turn spiral of the Saros Dial. Together, these models imply a surprisingly early epoch for the Antikythera Mechanism. The ancient Greeks built a machine that can predict, for many years ahead, not only eclipses but also a remarkable array of their characteristics, such as directions of obscuration, magnitude, colour, angular diameter of the Moon, relationship with the Moon’s node and eclipse time. It was not entirely accurate, but it was an astonishing achievement for its era. [ABSTRACT FROM AUTHOR]

This study investigates the unsteady flow of Powell-Eyring fluid past an inclined stretching sheet. Unsteadiness in the flow is due to the time-dependence of the stretching velocity and wall temperature. Mathematical analysis is performed in the presence of thermal radiation and non-uniform heat source/sink. The relevant boundary layer equations are reduced into self-similar forms by suitable transformations. The analytic solutions are constructed in a series form by homotopy analysis method (HAM). The convergence interval of the auxiliary parameter is obtained. Graphical results displaying the influence of interesting parameters are given. Numerical values of skin friction coefficient and local Nusselt number are computed and analyzed. [ABSTRACT FROM AUTHOR]

Background: Confident identification of microRNA-target interactions is significant for studying the function of microRNA (miRNA). Although some computational miRNA target prediction methods have been proposed for plants, results of various methods tend to be inconsistent and usually lead to more false positive. To address these issues, we developed an integrated model for identifying plant miRNA–target interactions. Results: Three online miRNA target prediction toolkits and machine learning algorithms were integrated to identify and analyze Arabidopsis thaliana miRNA-target interactions. Principle component analysis (PCA) feature extraction and self-training technology were introduced to improve the performance. Results showed that the proposed model outperformed the previously existing methods. The results were validated by using degradome sequencing supported Arabidopsis thaliana miRNA-target interactions. The proposed model constructed on Arabidopsis thaliana was run over Oryza sativa and Vitis vinifera to demonstrate that our model is effective for other plant species. Conclusions: The integrated model of online predictors and local PCA-SVM classifier gained credible and high quality miRNA-target interactions. The supervised learning algorithm of PCA-SVM classifier was employed in plant miRNA target identification for the first time. Its performance can be substantially improved if more experimentally proved training samples are provided. [ABSTRACT FROM AUTHOR]

In this paper, a randomized numerical approach is used to obtain approximate solutions for a class of nonlinear Fredholm integral equations of the second kind. The proposed approach contains two steps: at first, we define a discretized form of the integral equation by quadrature formula methods and solution of this discretized form converges to the exact solution of the integral equation by considering some conditions on the kernel of the integral equation. And then we convert the problem to an optimal control problem by introducing an artificial control function. Following that, in the next step, solution of the discretized form is approximated by a kind of Monte Carlo (MC) random search algorithm. Finally, some examples are given to show the efficiency of the proposed approach. [ABSTRACT FROM AUTHOR]

Information theory allows us to investigate information processing in neural systems in terms of information transfer, storage and modification. Especially the measure of information transfer, transfer entropy, has seen a dramatic surge of interest in neuroscience. Estimating transfer entropy from two processes requires the observation of multiple realizations of these processes to estimate associated probability density functions. To obtain these necessary observations, available estimators typically assume stationarity of processes to allow pooling of observations over time. This assumption however, is a major obstacle to the application of these estimators in neuroscience as observed processes are often non-stationary. As a solution, Gomez-Herrero and colleagues theoretically showed that the stationarity assumption may be avoided by estimating transfer entropy from an ensemble of realizations. Such an ensemble of realizations is often readily available in neuroscience experiments in the form of experimental trials. Thus, in this work we combine the ensemble method with a recently proposed transfer entropy estimator to make transfer entropy estimation applicable to non-stationary time series. We present an efficient implementation of the approach that is suitable for the increased computational demand of the ensemble method's practical application. In particular, we use a massively parallel implementation for a graphics processing unit to handle the computationally most heavy aspects of the ensemble method for transfer entropy estimation. We test the performance and robustness of our implementation on data from numerical simulations of stochastic processes. We also demonstrate the applicability of the ensemble method to magnetoencephalographic data. While we mainly evaluate the proposed method for neuroscience data, we expect it to be applicable in a variety of fields that are concerned with the analysis of information transfer in complex biological, social, and artificial systems. [ABSTRACT FROM AUTHOR]

The Library of Integrated Network-Based Cellular Signatures (LINCS) project aims to create a network-based understanding of biology by cataloging changes in gene expression and signal transduction that occur when cells are exposed to a variety of perturbations. It is helpful for understanding cell pathways and facilitating drug discovery. Here, we developed a novel approach to infer cell-specific pathways and identify a compound's effects using gene expression and phosphoproteomics data under treatments with different compounds. Gene expression data were employed to infer potential targets of compounds and create a generic pathway map. Binary linear programming (BLP) was then developed to optimize the generic pathway topology based on the mid-stage signaling response of phosphorylation. To demonstrate effectiveness of this approach, we built a generic pathway map for the MCF7 breast cancer cell line and inferred the cell-specific pathways by BLP. The first group of 11 compounds was utilized to optimize the generic pathways, and then 4 compounds were used to identify effects based on the inferred cell-specific pathways. Cross-validation indicated that the cell-specific pathways reliably predicted a compound's effects. Finally, we applied BLP to re-optimize the cell-specific pathways to predict the effects of 4 compounds (trichostatin A, MS-275, staurosporine, and digoxigenin) according to compound-induced topological alterations. Trichostatin A and MS-275 (both HDAC inhibitors) inhibited the downstream pathway of HDAC1 and caused cell growth arrest via activation of p53 and p21; the effects of digoxigenin were totally opposite. Staurosporine blocked the cell cycle via p53 and p21, but also promoted cell growth via activated HDAC1 and its downstream pathway. Our approach was also applied to the PC3 prostate cancer cell line, and the cross-validation analysis showed very good accuracy in predicting effects of 4 compounds. In summary, our computational model can be used to elucidate potential mechanisms of a compound's efficacy. [ABSTRACT FROM AUTHOR]