Scholarly Works - Isaac Newton Institute for Mathematical Sciences
Permanent URI for this collection
Browse
Recent Submissions
Item Open Access Published version Peer-reviewed An integral equation method for the homogenization of unidirectional fibre-reinforced media; antiplane elasticity and other potential problems(The Royal Society Publishing, 2017-05) Joyce, D; Parnell, WJ; Assier, RC; Abrahams, ID; Abrahams, Ian [0000-0002-4134-4438]In Parnell & Abrahams (2008 Proc. R. Soc. A 464, 1461–1482. (doi:10.1098/rspa.2007.0254)), a homogenization scheme was developed that gave rise to explicit forms for the effective antiplane shear moduli of a periodic unidirectional fibre-reinforced medium where fibres have non-circular cross section. The explicit expressions are rational functions in the volume fraction. In that scheme, a (non-dilute) approximation was invoked to determine leading-order expressions. Agreement with existing methods was shown to be good except at very high volume fractions. Here, the theory is extended in order to determine higher-order terms in the expansion. Explicit expressions for effective properties can be derived for fibres with non-circular cross section, without recourse to numerical methods. Terms appearing in the expressions are identified as being associated with the lattice geometry of the periodic fibre distribution, fibre cross-sectional shape and host/fibre material properties. Results are derived in the context of antiplane elasticity but the analogy with the potential problem illustrates the broad applicability of the method to, e.g. thermal, electrostatic and magnetostatic problems. The efficacy of the scheme is illustrated by comparison with the well-established method of asymptotic homogenization where for fibres of general cross section, the associated cell problem must be solved by some computational scheme.Item Open Access Accepted version Peer-reviewed Delay-independent asymptotic stability in monotone systems(IEEE, 2015) Devane, E; Lestas, IMonotone systems comprise an important class of dynamical systems that are of interest both for their wide applicability and because of their interesting mathematical properties. It is known that under the property of quasimono-tonicity time-delayed systems become monotone, and some remarkable properties have been reported for such systems. These include, for example, the fact that for linear systems global asymptotic stability of the undelayed system implies global asymptotic stability for the delayed system under arbitrary bounded delays. Nevertheless, extensions to nonlinear systems have thus far relied on various restrictive conditions, such as homogeneity and subhomogeneity, and it has been conjectured that these can be relaxed. Our aim in this paper is to show that this is feasible for a general class of nonlinear monotone systems, by deriving asymptotic stability results in which simple properties of the undelayed system lead to delay-independent stability. In particular, one of our results is to show that if the undelayed system has a convergent trajectory that is unbounded in all components as t → -∞ then the system is globally asymptotically stable for arbitrary time-varying delays. This follows from a more general result derived in the paper where delay-independent regions of attraction are quantified from the asymptotic behavior of individual trajectories of the undelayed system. This result recovers various known delay-independent stability results, and several examples are included in the paper to illustrate the significance of the proposed stability conditions.Item Open Access BayesPeak: Bayesian analysis of ChIP-seq data.(Springer Science and Business Media LLC, 2009-09-21) Spyrou, Christiana; Stark, Rory; Lynch, Andy G; Tavaré, Simon; Stark, Rory [0000-0002-1790-5469]; Lynch, Andy [0000-0002-7876-7338]BACKGROUND: High-throughput sequencing technology has become popular and widely used to study protein and DNA interactions. Chromatin immunoprecipitation, followed by sequencing of the resulting samples, produces large amounts of data that can be used to map genomic features such as transcription factor binding sites and histone modifications. METHODS: Our proposed statistical algorithm, BayesPeak, uses a fully Bayesian hidden Markov model to detect enriched locations in the genome. The structure accommodates the natural features of the Solexa/Illumina sequencing data and allows for overdispersion in the abundance of reads in different regions. Moreover, a control sample can be incorporated in the analysis to account for experimental and sequence biases. Markov chain Monte Carlo algorithms are applied to estimate the posterior distributions of the model parameters, and posterior probabilities are used to detect the sites of interest. CONCLUSION: We have presented a flexible approach for identifying peaks from ChIP-seq reads, suitable for use on both transcription factor binding and histone modification data. Our method estimates probabilities of enrichment that can be used in downstream analysis. The method is assessed using experimentally verified data and is shown to provide high-confidence calls with low false positive rates.