Inferring neuronal network connectivity from spike data: A temporal data mining approach

Understanding the functioning of a neural system in terms of its underlying circuitry is an important problem in neuroscience. Recent d evelopments in electrophysiology and imaging allow one to simultaneously record activities of hundreds of neurons. Inferring the underlying neuronal connectivity patterns from such multi-neuronal spike train data streams is a challenging statistical and computational problem. This task involves finding significant temporal patterns from vast amounts of symbolic time series data. In this paper we show that the frequent episode mining methods from the field of temporal data mining can be very useful in this context. In the frequent episode discovery framework, the data is viewed as a sequence of events, each of which is characterized by an event type and its time of occurrence and episodes are certain types of temporal patterns in such data. Here we show that, using the set of discovered frequent episodes from multi-neuronal data, one can infer different types of connectivity patterns in the neural system that generated it. For this purpose, we introduce the notion of mining for frequent episodes under certain temporal constraints; the structure of these temporal constraints is motivated by the application. We present algorithms for discovering serial and parallel episodes under these temporal constraints. Through extensive simulation studies we demonstrate that these methods are useful for unearthing patterns of neuronal network connectivity.

NETGEM: Network Embedded Temporal GEnerative Model for gene expression data

Background: Temporal analysis of gene expression data has been limited to identifying genes whose expression varies with time and/or correlation between genes that have similar temporal profiles. Often, the methods do not consider the underlying network constraints that connect the genes. It is becoming increasingly evident that interactions change substantially with time. Thus far, there is no systematic method to relate the temporal changes in gene expression to the dynamics of interactions between them. Information on interaction dynamics would open up possibilities for discovering new mechanisms of regulation by providing valuable insight into identifying time-sensitive interactions as well as permit studies on the effect of a genetic perturbation. Results: We present NETGEM, a tractable model rooted in Markov dynamics, for analyzing the dynamics of the interactions between proteins based on the dynamics of the expression changes of the genes that encode them. The model treats the interaction strengths as random variables which are modulated by suitable priors. This approach is necessitated by the extremely small sample size of the datasets, relative to the number of interactions. The model is amenable to a linear time algorithm for efficient inference. Using temporal gene expression data, NETGEM was successful in identifying (i) temporal interactions and determining their strength, (ii) functional categories of the actively interacting partners and (iii) dynamics of interactions in perturbed networks. Conclusions: NETGEM represents an optimal trade-off between model complexity and data requirement. It was able to deduce actively interacting genes and functional categories from temporal gene expression data. It permits inference by incorporating the information available in perturbed networks. Given that the inputs to NETGEM are only the network and the temporal variation of the nodes, this algorithm promises to have widespread applications, beyond biological systems. The source code for NETGEM is available from https://github.com/vjethava/NETGEM

Short- and long-term temporal variations in salinity and the oxygen, carbon and hydrogen isotopic compositions of the Hooghly Estuary water, India

The Hooghly River estuary provides a unique experimental site to understand the effect of monsoonal river discharge on freshwater and seawater mixing. Water samples collected bi-weekly for a duration of 17 months were analyzed for salinity, delta O-18,delta C-13(DIC), as well as delta D to investigate the differential mixing of freshwater and seawater. The differences in salinity and delta O-18 of samples collected during low and high tides on the same day are strongly correlated suggesting a well mixed water column at our sampling site. Low salinity and depleted delta O-18 during monsoon is consistent with increased river discharge as well as high rainfall. We identified different slopes in a delta O-18 versus salinity plot for the estuary water samples collected during monsoon and non-monsoon seasons. This is driven by composition of the freshwater source which is dominated by rainwater during monsoon and rivers during non-monsoon months. Selected delta D analyses of samples indicate that groundwater contributes significantly to the Hooghly Estuary during low rainfall times of the year. delta C-13(DIC) measured in the water recorded low values towards the end of monsoon indicating low productivity (i.e. increased organic respiration) while progressively increasing delta C-13(DIC) values from October till January as well as during some of the pre-monsoon months can be explained by increasing productivity. Very low delta C-13(DIC) (similar to-20%0) suggests involvement of carbon derived from anaerobic oxidation of organics and/or methane with potential contribution from increased anthropogenic water supply. An estimate of seawater incursion into the Hooghly Estuary at different times of the year is obtained by using salinity data in a two-component mixing model. Presence of seawater was found maximum (31-37%) during February till July and lowest (less than or equal to 6%) from September till November. We notice a temporal offset between Ganges River discharge farther upstream at Farakka and salinity variation at the Hooghly Estuary. We believe that this time lag is a result of the physical distance between Farakka and Kakdweep (our sampling location) and put constraints on the travel time of river water during early monsoon. (c) 2012 Published by Elsevier B.V.

Constraints in relativistic Hamiltonian mechanics

We propose a new scheme for the use of constraints in setting up classical, Hamiltonian, relativistic, interacting particle theories. We show that it possesses both Poincaré invariance and invariance of world lines. We discuss the transition to the physical phase space and the nonrelativistic limit.

An Algorithm to Decide Feasibility of Linear Integer Constraints Occurrng in Decision Tables

Abstract-To detect errors in decision tables one needs to decide whether a given set of constraints is feasible or not. This paper describes an algorithm to do so when the constraints are linear in variables that take only integer values. Decision tables with such constraints occur frequently in business data processing and in nonnumeric applications. The aim of the algorithm is to exploit. the abundance of very simple constraints that occur in typical decision table contexts. Essentially, the algorithm is a backtrack procedure where the the solution space is pruned by using the set of simple constrains. After some simplications, the simple constraints are captured in an acyclic directed graph with weighted edges. Further, only those partial vectors are considered from extension which can be extended to assignments that will at least satisfy the simple constraints. This is how pruning of the solution space is achieved. For every partial assignment considered, the graph representation of the simple constraints provides a lower bound for each variable which is not yet assigned a value. These lower bounds play a vital role in the algorithm and they are obtained in an efficient manner by updating older lower bounds. Our present algorithm also incorporates an idea by which it can be checked whether or not an (m - 2)-ary vector can be extended to a solution vector of m components, thereby backtracking is reduced by one component.

Monoblock cylinder internally pressurised under asymmetric end constraints Some aspects of design

The finite element method is used to analyse stresses and displacements in a monoblock cylinder open at one end only. The cylinder is internally pressurised. The analysis shows that the minimum pressure required to cause yield in the cylinder decreases rapidly with increasing cylinder height until the height is about the same as the outer radius of the cylinder, beyond which the decrease is marginal. Introduction of a fillet at the internal corner enhances the design pressure substantially while a fillet at the outer corner affects this pressure only marginally.

Ionized Gas towards Galactic Centre-Constraints from Low-Frequency Recombination Lines

Observations of the H272α recombination line towards the galactic centre show features near VLSR= 0, -50 and + 36 kms-1 . We have combined the parameters of these features with the available -166 measurements to obtain the properties of the ionized gas present along the line of sight and also in the -3 kpc arm-. For the line-of-sight ionized gas we get an electron density around 7 cm-3 and a pathlength through it 10-60 pc. The emission measure and the electron temperature are in the range 500-2900 pc cm-6 and 2000-6000 - respectively. The ionized gas in the 3 kpc arm has an electron density of 30 cm-3 and extends over 9 pc along the line of sight if we assume an electron temperature of 104 K. Using the available upper limit to the intensity of the H351α recombination line, we show that the distributed ionized gas responsible for the dispersion of pulsar signals should have a temperature > 4500 - and a minimum filling factor of 20 per cent. We also show that recombination lines from the -warm ionized- gas proposed by McKee & Ostriker (1977) should be detectable in the frequency range 100-150 MHz towards the galactic centre with the sensitivity available at present.

Reassessing the earthquake hazard in Kerala based on the historical and current seismicity

Given the lack of proper constraints in understanding earthquake mechanisms in the cratonic interiors and the general absence of good quality database, here we reassess the seismic hazard in the province of Kerala, a part of the aEuro cent stable continental interioraEuro cent, based on an improved historical and instrumental database. The temporal pattern of the current seismicity suggests that > 60% of the microtremors in Kerala occurs with a time lag after the peak rainfall, indicating that hydroseismicity may be a plausible model to explain the low-level seismicity in this region. Further, an increment in overall seismicity rate in the region in the recent years is explained as due to increased anthropogenic activities, which includes changes in hydrological pathways as a consequence of rapid landscape changes. Our analyses of the historical database eliminate a few events that are ascribed to this region; this exercise has also led to identification of a few events, not previously noted. The improved historical database essentially suggests that the central midland region is more prone to seismic activity compared to other parts of Kerala. This region appears to have generated larger number of significant earthquakes; the most prominent being the multiple events (doublets) of 1856 and 1953, whose magnitudes are comparable to that of the 2000/2001 (central Kerala) events. Occurrences of these historical events and the recent earthquakes, and the local geology indicative of pervasive faulting as shown by widely distributed pseudotachylite veins suggest that the NNW-SSE trending faults in central midland Kerala may host discrete potentially active sources that may be capable of generating light to moderate size earthquakes. The frequency of earthquakes in central Kerala evident from the historical database requires that the seismic codes stipulated for this region are made mandatory.

New constraints on the pion EM form factor using Pi `(-Q(2))

We study the constraints arising on the expansion parameters c and d of the pion electromagnetic form factor from the inclusion of pure spacelike data and the phase of timelike data along with one spacelike datum, using as input the first derivative of the QCD polarization amplitude Pi'(-Q(2)). These constraints when combined with other analyses, provide a valuable check on a determination of c due to Guo et al. and on our previous work where pionic contribution to the (g - 2) of the muon was used as the input. This work further illustrates the power of analyticity techniques in form factor analysis.

Derivation of the decoupling relations for composite models from anomaly constraints and n-independence

It is shown, in the composite fermion models studied by 't Hooft and others, that the requirements of Adler-Bell-Jackiw anomaly matching and n-independence are sufficient to fix the indices of composite representations. The third requirement, namely that of decoupling relations, follows from these two constraints in such models and hence is inessential.

Aspects of tautomerism. 13. Alkaline hydrolysis of .gamma.-, .delta.-, and .epsilon.-keto esters and their desoxy analogs. Geometrical constraints on keto participation

The rates of alkaline hydrolysis of methyl &benzoylpropionate (I), methyl y-benzoylbutyrate (11) and methyll6-benzoylvalerate (In) decrease in the order I > I1 > III. Keto participation is the predominant pathway in the case of y-keto esters. Evidence has also been obtained for keto participation in the case of 6-keto esters, whereas no such evidence is available in the case of r-keto esters studied.

Structures of isopropylidene nucleoside derivatives - implications for ribose ring flexibility under external cyclic constraints

Crystal structures of six isopropylidene nucleoside derivatives are described. The results show that, under external cyclic constraints, the ribose assumes a variety of unusual conformations. In those compounds which possess a base-to-sugar cyclization through the C(4′) atom, the furanose pucker is predominantly C(4′)-endo, O(4′)-exo. The possible relevance of the sulphur geometry in two of the compounds to certain structural aspects of the action of the enzyme thymidylate synthetase is also pointed out.

Symmetries and constraints in generalized Hamiltonian dynamics

A general analysis of symmetries and constraints for singular Lagrangian systems is given. It is shown that symmetry transformations can be expressed as canonical transformations in phase space, even for such systems. The relation of symmetries to generators, constraints, commutators, and Dirac brackets is clarified.

Temporal variation of tonal spectrum of tambura

Tambura is an essential drone accompaniment used in Indian music concerts. It acts as an immediate reference of pitch for both the artists and listeners. The four strings of Tambura are tuned to the frequency ratio :1:1: . Careful listening to Tambura sound reveals that the tonal spectrum is not stationary but is time varying. The object of this study is to make a detailed spectrum analysis to find out the nature of temporal variation of the tonal spectrum of Tambura sound. Results of the analysis are correlated with perceptual evaluation conducted in a controlled acoustic environment. A significant result of this study is to demonstrate the presence of several notes which are normally not noticed even by a professional artist. The effect of bridge in Tambura in producing the so called “live tone” is explained through time and frequency parameters of Tambura sounds.

Determination of optimal trajectory under design constraints for a satellite launch vehicle

An optimal pitch steering programme of a solid-fuel satellite launch vehicle to maximize either (1) the injection velocity at a given altitude, or (2) the size of circular orbit, for a given payload is presented. The two-dimensional model includes the rotation of atmosphere with the Earth, the vehicle's lift and drag, variation of thrust with time and altitude, inverse-square gravitational field, and the specified initial vertical take-off. The inequality constraints on the aerodynamic load, control force, and turning rates are also imposed. Using the properties of the central force motion the terminal constraint conditions at coast apogee are transferred to the penultimate stage burnout. Such a transformation converts a time-free problem into a time-fixed one, reduces the number of terminal constraints, improves accuracy, besides demanding less computer memory and time. The adjoint equations are developed in a compact matrix form. The problem is solved on an IBM 360/44 computer using a steepest ascent algorithm. An illustrative analysis of a typical launch vehicle establishes the speed of convergence, and accuracy and applicability of the algorithm.