A boundary preserving numerical algorithm for the Wright-Fisher model with mutation

The Wright-Fisher model is an Itô stochastic differential equation that was originally introduced to model genetic drift within finite populations and has recently been used as an approximation to ion channel dynamics within cardiac and neuronal cells. While analytic solutions to this equation remain within the interval [0,1], current numerical methods are unable to preserve such boundaries in the approximation. We present a new numerical method that guarantees approximations to a form of Wright-Fisher model, which includes mutation, remain within [0,1] for all time with probability one. Strong convergence of the method is proved and numerical experiments suggest that this new scheme converges with strong order 1/2. Extending this method to a multidimensional case, numerical tests suggest that the algorithm still converges strongly with order 1/2. Finally, numerical solutions obtained using this new method are compared to those obtained using the Euler-Maruyama method where the Wiener increment is resampled to ensure solutions remain within [0,1].

Stochastic models and simulation of ion channel dynamics

The behaviour of ion channels within cardiac and neuronal cells is intrinsically stochastic in nature. When the number of channels is small this stochastic noise is large and can have an impact on the dynamics of the system which is potentially an issue when modelling small neurons and drug block in cardiac cells. While exact methods correctly capture the stochastic dynamics of a system they are computationally expensive, restricting their inclusion into tissue level models and so approximations to exact methods are often used instead. The other issue in modelling ion channel dynamics is that the transition rates are voltage dependent, adding a level of complexity as the channel dynamics are coupled to the membrane potential. By assuming that such transition rates are constant over each time step, it is possible to derive a stochastic differential equation (SDE), in the same manner as for biochemical reaction networks, that describes the stochastic dynamics of ion channels. While such a model is more computationally efficient than exact methods we show that there are analytical problems with the resulting SDE as well as issues in using current numerical schemes to solve such an equation. We therefore make two contributions: develop a different model to describe the stochastic ion channel dynamics that analytically behaves in the correct manner and also discuss numerical methods that preserve the analytical properties of the model.

Probabilité de fixation dans des modèles génétiques de populations à plusieurs allèles

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Processus de coalescence dans une population subdivisée avec possibilité de coalescences multiples

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Power law behavior of the isotope yield distributions in the multifragmentation regime of heavy ion reactions

Isotope yield distributions in the multifragmentation regime were studied with high-quality isotope identification, focusing on the intermediate mass fragments (IMFs) produced in semiviolent collisions. The yields were analyzed within the framework of a modified Fisher model. Using the ratio of the mass-dependent symmetry energy coefficient relative to the temperature, a(sym)/T, extracted in previous work and that of the pairing term, a(p)/T, extracted from this work, and assuming that both reflect secondary decay processes, the experimentally observed isotope yields were corrected for these effects. For a given I = N - Z value, the corrected yields of isotopes relative to the yield of C-12 show a power law distribution Y (N, Z)/Y(C-12) similar to A(-tau) in the mass range 1 <= A <= 30, and the distributions are almost identical for the different reactions studied. The observed power law distributions change systematically when I of the isotopes changes and the extracted tau value decreases from 3.9 to 1.0 as I increases from -1 to 3. These observations are well reproduced by a simple deexcitation model, with which the power law distribution of the primary isotopes is determined to be tau(prim) = 2.4 +/- 0.2, suggesting that the disassembling system at the time of the fragment formation is indeed at, or very near, the critical point.

Isobaric yield ratios and the symmetry energy in heavy-ion reactions near the Fermi energy

The relative isobaric yields of fragments produced in a series of heavy-ion-induced multifragmentation reactions have been analyzed in the framework of a modified Fisher model, primarily to determine the ratio of the symmetry energy coefficient to the temperature, a(sym)/T, as a function of fragment mass A. The extracted values increase from 5 to similar to 16 as A increases from 9 to 37. These values have been compared to the results of calculations using the antisymmetrized molecular dynamics (AMD) model together with the statistical decay code GEMINI. The calculated ratios are in good agreement with those extracted from the experiment. In contrast, the values extracted from the ratios of the primary isobars from the AMD model calculation are similar to 4 to 5 and show little variation with A. This observation indicates that the value of the symmetry energy coefficient derived from final fragment observables may be significantly different than the actual value at the time of fragment formation. The experimentally observed pairing effect is also studied within the same simulations. The Coulomb coefficient is also discussed.

A novel approach to isoscaling: The role of the order parameter m = N-f-Z(f)/A(f)

Isoscaling is derived within a recently proposed modified Fisher model where the free energy near the critical point is described by the Landau O(m(6)) theory. In this model m = N-f-Z(f)/A(f) is the order parameter, a consequence of (one of) the symmetries of the nuclear Hamiltonian. Within this framework we show that isoscaling depends mainly on this order parameter through the 'external (conjugate) field' H. The external field is just given by the difference in chemical potentials of the neutrons and protons of the two sources. To distinguish from previously employed isoscaling relationships, this approach is dubbed: m-scaling. We discuss the relationship between this framework and the standard isoscaling formalism and point out some substantial differences in interpretation of experimental results which might result. These should be investigated further both theoretically and experimentally. (C) 2010 Elsevier B.V. All rights reserved.

Comparative evaluation of a new effective population size estimator based on approximate Bayesian computation

We describe and evaluate a new estimator of the effective population size (N-e), a critical parameter in evolutionary and conservation biology. This new "SummStat" N-e. estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer N-e. Simulations of a Wright-Fisher population with known N-e show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and N-e values. We also address the paucity of information about the relative performance of N-e estimators by comparing the SUMMStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated rising initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and Ne less than or equal to 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true N-e. The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any, potentially informative summary statistic from Population genetic data.

Essays on monetary and banking theory

Esta tese de Doutorado é dedicada ao estudo de instabilidade financeira e dinâmica em Teoria Monet ária. E demonstrado que corridas banc árias são eliminadas sem custos no modelo padrão de teoria banc ária quando a popula ção não é pequena. É proposta uma extensão em que incerteza agregada é mais severa e o custo da estabilidade financeira é relevante. Finalmente, estabelece-se otimalidade de transições na distribui ção de moeda em economias em que oportunidades de trocas são escassas e heterogêneas. Em particular, otimalidade da inflação depende dos incentivos dinâmicos proporcionados por tais transi ções. O capí tulo 1 estabelece o resultado de estabilidade sem custos para economias grandes ao estudar os efeitos do tamanho populacional na an álise de corridas banc árias de Peck & Shell. No capí tulo 2, otimalidade de dinâmica é estudada no modelo de monet ário de Kiyotaki & Wright quando a sociedade é capaz de implementar uma polí tica inflacion ária. Apesar de adotar a abordagem de desenho de mecanismos, este capí tulo faz um paralelo com a an álise de Sargent & Wallace (1981) ao destacar efeitos de incentivos dinâmicos sobre a interação entre as polí ticas monet ária e fiscal. O cap ítulo 3 retoma o tema de estabilidade fi nanceira ao quanti car os custos envolvidos no desenho ótimo de um setor bancário à prova de corridas e ao propor uma estrutura informacional alternativa que possibilita bancos insolventes. A primeira an álise mostra que o esquema de estabilidade ótima exibe altas taxas de juros de longo prazo e a segunda que monitoramento imperfeito pode levar a corridas bancárias com insolvência.

Seawater carbonate chemistry and shell weights of Laternula elliptica, Yoldia eightsi, Nacella concinna and Liothyrella uva during experiments, 2009

Antarctic calcified macroorganisms are particularly vulnerable to ocean acidification because many are weakly calcified, the dissolution rates of calcium carbonate are inversely related to temperature, and high latitude seas are predicted to become undersaturated in aragonite by the year 2100. We examined the post-mortem dissolution rates of aragonitic and calcitic shells from four species of Antarctic benthic marine invertebrates (two bivalves, one limpet, one brachiopod) and the thallus of a limpet shell-encrusting coralline alga exposed to acidified pH (7.4) or non-acidified pH (8.2) seawater at a constant temperature of 4 C. Within a period of only 14-35 days, shells of all four species held in pH 7.4 seawater had suffered significant dissolution. Despite calcite being 35% less soluble in seawater than aragonite, there was surprisingly, no consistent pattern of calcitic shells having slower dissolution rates than aragonitic shells. Outer surfaces of shells held in pH 7.4 seawater exhibited deterioration by day 35, and by day 56 there was exposure of aragonitic or calcitic prisms within the shell architecture of three of the macroinvertebrate species. Dissolution of coralline algae was confirmed by differences in weight loss in limpet shells with and without coralline algae. By day 56, thalli of the coralline alga held in pH 7.4 displayed a loss of definition of the conceptacle pores and cracking was evident at the zone of interface with limpet shells. Experimental studies are needed to evaluate whether there are adequate compensatory mechanisms in these and other calcified Antarctic benthic macroorganisms to cope with anticipated ocean acidification. In their absence, these organisms, and the communities they comprise, are likely to be among the first to experience the cascading impacts of ocean acidification.

Estimating cell diffusivity and cell proliferation rate by interpreting IncuCyte ZOOM™ assay data using the Fisher-Kolmogorov model

Background: Standard methods for quantifying IncuCyte ZOOM™ assays involve measurements that quantify how rapidly the initially-vacant area becomes re-colonised with cells as a function of time. Unfortunately, these measurements give no insight into the details of the cellular-level mechanisms acting to close the initially-vacant area. We provide an alternative method enabling us to quantify the role of cell motility and cell proliferation separately. To achieve this we calibrate standard data available from IncuCyte ZOOM™ images to the solution of the Fisher-Kolmogorov model. Results: The Fisher-Kolmogorov model is a reaction-diffusion equation that has been used to describe collective cell spreading driven by cell migration, characterised by a cell diffusivity, D, and carrying capacity limited proliferation with proliferation rate, λ, and carrying capacity density, K. By analysing temporal changes in cell density in several subregions located well-behind the initial position of the leading edge we estimate λ and K. Given these estimates, we then apply automatic leading edge detection algorithms to the images produced by the IncuCyte ZOOM™ assay and match this data with a numerical solution of the Fisher-Kolmogorov equation to provide an estimate of D. We demonstrate this method by applying it to interpret a suite of IncuCyte ZOOM™ assays using PC-3 prostate cancer cells and obtain estimates of D, λ and K. Comparing estimates of D, λ and K for a control assay with estimates of D, λ and K for assays where epidermal growth factor (EGF) is applied in varying concentrations confirms that EGF enhances the rate of scratch closure and that this stimulation is driven by an increase in D and λ, whereas K is relatively unaffected by EGF. Conclusions: Our approach for estimating D, λ and K from an IncuCyte ZOOM™ assay provides more detail about cellular-level behaviour than standard methods for analysing these assays. In particular, our approach can be used to quantify the balance of cell migration and cell proliferation and, as we demonstrate, allow us to quantify how the addition of growth factors affects these processes individually.

Search for the standard model Higgs boson decaying to bottom quarks in pp collisions at root s=7 TeV

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Search for physics beyond the standard model using multilepton signatures in pp collisions at root s=7 TeV

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Search for physics beyond the standard model in events with a Z boson, jets, and missing transverse energy in pp collisions at root s=7 TeV

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)