Towards an exact reconstruction of a time-invariant model from time series data

Dynamic processes in biological systems may be profiled by measuring system properties over time. One way of representing such time series data is through weighted interaction networks, where the nodes in the network represent the measurables and the weighted edges represent interactions between any pair of nodes. Construction of these network models from time series data may involve seeking a robust data-consistent and time-invariant model to approximate and describe system dynamics. Many problems in mathematics, systems biology and physics can be recast into this form and may require finding the most consistent solution to a set of first order differential equations. This is especially challenging in cases where the number of data points is less than or equal to the number of measurables. We present a novel computational method for network reconstruction with limited time series data. To test our method, we use artificial time series data generated from known network models. We then attempt to reconstruct the original network from the time series data alone. We find good agreement between the original and predicted networks

Wave modelling - the state of the art

This paper is the product of the wave modelling community and it tries to make a picture of the present situation in this branch of science, exploring the previous and the most recent results and looking ahead towards the solution of the problems we presently face. Both theory and applications are considered. The many faces of the subject imply separate discussions. This is reflected into the single sections, seven of them, each dealing with a specific topic, the whole providing a broad and solid overview of the present state of the art. After an introduction framing the problem and the approach we followed, we deal in sequence with the following subjects: (Section) 2, generation by wind; 3, nonlinear interactions in deep water; 4, white-capping dissipation; 5, nonlinear interactions in shallow water; 6, dissipation at the sea bottom; 7, wave propagation; 8, numerics. The two final sections, 9 and 10, summarize the present situation from a general point of view and try to look at the future developments

Experimental design trade-offs for gene regulatory network inference: an in silico study of the yeast Saccharomyces cerevisiae cell cycle

Time-series of high throughput gene sequencing data intended for gene regulatory network (GRN) inference are often short due to the high costs of sampling cell systems. Moreover, experimentalists lack a set of quantitative guidelines that prescribe the minimal number of samples required to infer a reliable GRN model. We study the temporal resolution of data vs quality of GRN inference in order to ultimately overcome this deficit. The evolution of a Markovian jump process model for the Ras/cAMP/PKA pathway of proteins and metabolites in the G1 phase of the Saccharomyces cerevisiae cell cycle is sampled at a number of different rates. For each time-series we infer a linear regression model of the GRN using the LASSO method. The inferred network topology is evaluated in terms of the area under the precision-recall curve AUPR. By plotting the AUPR against the number of samples, we show that the trade-off has a, roughly speaking, sigmoid shape. An optimal number of samples corresponds to values on the ridge of the sigmoid

From Markovian to pairwise epidemic models and the performance of moment closure approximations

Many if not all models of disease transmission on networks can be linked to the exact state-based Markovian formulation. However the large number of equations for any system of realistic size limits their applicability to small populations. As a result, most modelling work relies on simulation and pairwise models. In this paper, for a simple SIS dynamics on an arbitrary network, we formalise the link between a well known pairwise model and the exact Markovian formulation. This involves the rigorous derivation of the exact ODE model at the level of pairs in terms of the expected number of pairs and triples. The exact system is then closed using two different closures, one well established and one that has been recently proposed. A new interpretation of both closures is presented, which explains several of their previously observed properties. The closed dynamical systems are solved numerically and the results are compared to output from individual-based stochastic simulations. This is done for a range of networks with the same average degree and clustering coefficient but generated using different algorithms. It is shown that the ability of the pairwise system to accurately model an epidemic is fundamentally dependent on the underlying large-scale network structure. We show that the existing pairwise models are a good fit for certain types of network but have to be used with caution as higher-order network structures may compromise their effectiveness

Blackbox Quantization of Superconducting Circuits using exact Impedance Synthesis

We propose a new quantization method for superconducting electronic circuits involving a Josephson junction device coupled to a linear microwave environment. The method is based on an exact impedance synthesis of the microwave environment considered as a blackbox with impedance function Z(s). The synthesized circuit captures dissipative dynamics of the system with resistors coupled to the reactive part of the circuit in a non-trivial way. We quantize the circuit and compute relaxation rates following previous formalisms for lumped element circuit quantization. Up to the errors in the fit our method gives an exact description of the system and its losses