Structure learning of context-specific graphical models

Abstract The ultimate problem considered in this thesis is modeling a high-dimensional joint distribution over a set of discrete variables. For this purpose, we consider classes of context-specific graphical models and the main emphasis is on learning the structure of such models from data. Traditional graphical models compactly represent a joint distribution through a factorization justi ed by statements of conditional independence which are encoded by a graph structure. Context-speci c independence is a natural generalization of conditional independence that only holds in a certain context, speci ed by the conditioning variables. We introduce context-speci c generalizations of both Bayesian networks and Markov networks by including statements of context-specific independence which can be encoded as a part of the model structures. For the purpose of learning context-speci c model structures from data, we derive score functions, based on results from Bayesian statistics, by which the plausibility of a structure is assessed. To identify high-scoring structures, we construct stochastic and deterministic search algorithms designed to exploit the structural decomposition of our score functions. Numerical experiments on synthetic and real-world data show that the increased exibility of context-specific structures can more accurately emulate the dependence structure among the variables and thereby improve the predictive accuracy of the models.

Development of Craniofacial Bone Defect Reconstruction Implant Based on Fibre-Reinforced Composite with Photopolymerisable Resin Systems. Experimental studies in vitro and in vivo.

Reconstruction of defects in the craniomaxillofacial (CMF) area has mainly been based on bone grafts or metallic fixing plates and screws. Particularly in the case of large calvarial and/or craniofacial defects caused by trauma, tumours or congenital malformations, there is a need for reliable reconstruction biomaterials, because bone grafts or metallic fixing systems do not completely fulfill the criteria for the best possible reconstruction methods in these complicated cases. In this series of studies, the usability of fibre-reinforced composite (FRC) was studied as a biostable, nonmetallic alternative material for reconstructing artificially created bone defects in frontal and calvarial areas of rabbits. The experimental part of this work describes the different stages of the product development process from the first in vitro tests with resin-impregnated fibrereinforced composites to the in vivo animal studies, in which this FRC was tested as an implant material for reconstructing different size bone defects in rabbit frontal and calvarial areas. In the first in vitro study, the FRC was polymerised in contact with bone or blood in the laboratory. The polymerised FRC samples were then incubated in water, which was analysed for residual monomer content by using high performance liquid chromatography (HPLC). It was found that this in vitro polymerisation in contact with bone and blood did not markedly increase the residual monomer leaching from the FRC. In the second in vitro study, different adhesive systems were tested in fixing the implant to bone surface. This was done to find an alternative implant fixing system to screws and pins. On the basis of this study, it was found that the surface of the calvarial bone needed both mechanical and chemical treatments before the resinimpregnated FRC could be properly fixed onto it. In three animal studies performed with rabbit frontal bone defects and critical size calvarial bone defect models, biological responses to the FRC implants were evaluated. On the basis of theseevaluations, it can be concluded that the FRC, based on E-glass (electrical glass) fibres forming a porous fibre veil enables the ingrowth of connective tissues to the inner structures of the material, as well as the bone formation and mineralization inside the fibre veil. Bone formation could be enhanced by using bioactive glass granules fixed to the FRC implants. FRC-implanted bone defects healed partly; no total healing of defects was achieved. Biological responses during the follow-up time, at a maximum of 12 weeks, to resin-impregnated composite implant seemed to depend on the polymerization time of the resin matrix of the FRC. Both of the studied resin systems used in the FRC were photopolymerised and the heat-induced postpolymerisation was used additionally.

Empirical Tests of Asset Pricing Models in Finnish Stock Market

Tämän tutkielman tavoitteena on selvittää mitkä riskitekijät vaikuttavat osakkeiden tuottoihin. Arvopapereina käytetään kuutta portfoliota, jotka ovat jaoteltu markkina-arvon mukaan. Aikaperiodi on vuoden 1987 alusta vuoden 2004 loppuun. Malleina käytetään pääomamarkkinoiden hinnoittelumallia, arbitraasihinnoitteluteoriaa sekä kulutuspohjaista pääomamarkkinoiden hinnoittelumallia. Riskifaktoreina kahteen ensimmäiseen malliin käytetään markkinariskiä sekä makrotaloudellisia riskitekijöitä. Kulutuspohjaiseen pääomamarkkinoiden hinnoinoittelumallissa keskitytään estimoimaan kuluttajien riskitottumuksia sekä diskonttaustekijää, jolla kuluttaja arvostavat tulevaisuuden kulutusta. Tämä työ esittelee momenttiteorian, jolla pystymme estimoimaan lineaarisia sekä epälineaarisia yhtälöitä. Käytämme tätä menetelmää testaamissamme malleissa. Yhteenvetona tuloksista voidaan sanoa, että markkinabeeta onedelleen tärkein riskitekijä, mutta löydämme myös tukea makrotaloudellisille riskitekijöille. Kulutuspohjainen mallimme toimii melko hyvin antaen teoreettisesti hyväksyttäviä arvoja.

Modeling of non-isothermalvapor membrane separation with thermodynamic models and generalized mass transfer equations

A rigorous unit operation model is developed for vapor membrane separation. The new model is able to describe temperature, pressure, and concentration dependent permeation as wellreal fluid effects in vapor and gas separation with hydrocarbon selective rubbery polymeric membranes. The permeation through the membrane is described by a separate treatment of sorption and diffusion within the membrane. The chemical engineering thermodynamics is used to describe the equilibrium sorption of vapors and gases in rubbery membranes with equation of state models for polymeric systems. Also a new modification of the UNIFAC model is proposed for this purpose. Various thermodynamic models are extensively compared in order to verify the models' ability to predict and correlate experimental vapor-liquid equilibrium data. The penetrant transport through the selective layer of the membrane is described with the generalized Maxwell-Stefan equations, which are able to account for thebulk flux contribution as well as the diffusive coupling effect. A method is described to compute and correlate binary penetrant¿membrane diffusion coefficients from the experimental permeability coefficients at different temperatures and pressures. A fluid flow model for spiral-wound modules is derived from the conservation equation of mass, momentum, and energy. The conservation equations are presented in a discretized form by using the control volume approach. A combination of the permeation model and the fluid flow model yields the desired rigorous model for vapor membrane separation. The model is implemented into an inhouse process simulator and so vapor membrane separation may be evaluated as an integralpart of a process flowsheet.

The Value System of Consumer Board Packaging

Olemassa olevien teollisuuden rakenteiden ja arvoketjujen muutos on usein onnistuneiden liiketoimintainnovaatioiden takana. Yritysten tulisikin olla tietoisia rakenteista, jotka määrittelevät arvon jakaantumista niiden arvoketjuissa. Tämän työn tavoitteena oli tutkia ja arvioida kartongin pakkausketjujen rakenteita kartongintuottajalta aina kauppiaalle asti. Työn kuluessa yritysketjujen arviointia varten oli kehitettävä uudenlainen analysointityökalu. Työssä käsitellyt kartonkipakkausketjut olivat nestekartonkipakkaaminen, elintarvikepakkaaminen, lääkepakkaaminen sekä graafisten kartonkien ketju. Jokaista edellä mainituista ketjuista arvioitiin käyttämällä kehitettyä analysointityökalua. Työn merkittävin tulos on kehitetty työkalu, jolla pakkausketjuja sekä myös muita yritysketjuja on mahdollista analysoida. Kehitettyä työkalua käyttämällä oli mahdollista tehdä johtopäätöksiä vertikaalin integraation mahdollisuudesta sekä pakkausketjujen suotuisuudesta uusille bisnesmalleille. Sekä analysointityökalua että saatuja johtopäätöksiä voidaan käyttää hyväksi, kun kartongintuottaja suunnittelee tulevaisuuden toimenpiteitä.

Development of Customer Satisfaction Models for Electronics Engineering Company

Diplomityön tarkoituksena oli luoda ja kehittää kaksi asiakastyytyväisyysmallia asiakastyytyväisyyden mittaamisen aloittamiseksi ja toteuttamiseksi kohdeyrityksessä. Työ pohjautuu nykyisten tyytyväisyysprosessien analysointiin sekä työn teoriaosaan, joka käsittelee yksityiskohtaisesti niitä asioita, joita asiakastyytyväisyyden mittaamisessa ja prosessissa tulisi huomioida. Työssä tehdyn mallien tarkoituksen on auttaa kohdeyritystä hyödyntämään asiakastyytyväisyysmittauksen tuloksia paremmin liiketoiminnassa, sekä asiakkaiden keskuudessa. Työn yhtenä tavoitteena oli myös sopivan mittaustyökalun löytäminen ja suositteleminen kohdeyritykselle.Teorian ja analysoinnin pohjalta luotiin molemmat asiakastyytyväisyysmallit vastamaan kohdeyksiköiden tarpeita. Kun ulkoiset seikat, kuten mittaustavat, mittausinstrumentit, kyselylomakkeet ja vastaajaryhmät oli määritelty, keskityttiin tulosten analysointiin ja hyödyntämiseen, mikä korostui asiakassuuntautuneessa organisaatiossa. Työssä pohdittiin myös yhtenäisen asiakastyytyväisyysprosessin merkitystä ja etuja kohdeyrityksessä.

Business models enhancing the diffusion of mobile Internet in emerging markets

Tämän tutkimuksen päätavoitteena oli selvittää, millaiset liiketoimintamallit soveltuvat mobiilin internet-liiketoiminnan harjoittamiseen kehittyvillä markkinoilla. Tavoitteena oli myös selvittää tekijöitä, jotka vaikuttavat mobiilin internetin diffuusioon. Tutkimus tehtiin käyttäen sekä kvantitatiivista että kvalitatiivista tutkimusmenetelmää. Klusterianalyysin avulla 40 Euroopan maasta muodostettiin sisäisesti homogeenisiä maaklustereita. Näiden klustereiden avulla oli mahdollista suunnitella erityyppisille markkinoille soveltuvat liiketoimintamallit. Haastatteluissa selvitettiin asiantuntijoiden näkemyksiä tekijöistä, jotka vaikuttavat mobiilin internetin diffuusioon kehittyvillä markkinoilla. Tutkimuksessa saatiin selville, että tärkeimmät liiketoimintamallin elementit kehittyvillä markkinoilla ovat hinnoittelu, arvotarjooma ja arvoverkko. Puutteellisen kiinteän verkon todettiin olevan yksi tärkeimmistä mobiilin internetin diffuusiota edistävistä tekijöistä kehittyvillä markkinoilla.

A comparison of electricity spot prices simulation using ARMA-GARCH and mean-reverting models

The aim of this work is to compare two families of mathematical models for their respective capability to capture the statistical properties of real electricity spot market time series. The first model family is ARMA-GARCH models and the second model family is mean-reverting Ornstein-Uhlenbeck models. These two models have been applied to two price series of Nordic Nord Pool spot market for electricity namely to the System prices and to the DenmarkW prices. The parameters of both models were calibrated from the real time series. After carrying out simulation with optimal models from both families we conclude that neither ARMA-GARCH models, nor conventional mean-reverting Ornstein-Uhlenbeck models, even when calibrated optimally with real electricity spot market price or return series, capture the statistical characteristics of the real series. But in the case of less spiky behavior (System prices), the mean-reverting Ornstein-Uhlenbeck model could be seen to partially succeeded in this task.

Innovations of ageing and societal transition : Dynamics of change of the socio-technical regime of ageing

The purpose of this dissertation is to examine the dynamics of the socio-technical system in the field of ageing. The study stems from the notion that the ageing of the population as a powerful megatrend has wide societal effects, and is not just a matter for the social and health sector. The central topic in the study is change: not only the age structures and structures of society are changing, but also at the same time there is constant development, for instance, in technologies, infrastructures and cultural perceptions. The changing concept of innovation has widened the understanding of innovations related to ageing from medical and assistive technological innovations to service and social innovations, as well as systemic innovations at different levels, which means the intertwined and co-evolutionary change in technologies, structures, services and thinking models. By the same token, the perceptions of older people and old age are becoming more multi-faceted: old age is no longer equated to illnesses and decline, but visions of active ageing and a third age have emerged, which are framed by choices, opportunities, resources and consumption in later life. The research task in this study is to open up the processes and mechanisms of change in the field of ageing, which are studied as a complex, multi-level and interrelated socio-technical system. The question is about co-effective elements consisting of macro-level landscape changes, the existing socio-technical regime (the rule system, practices and structures) and bottom-up niche-innovations. Societal transitions do not account for the things inside the regime alone, or for the long-term changes in the landscape, nor for the radical innovations, but for the interplay between all these levels. The research problem is studied through five research articles, which offer micro-level case studies to macro-level phenomenon. Each of the articles focus on different aspects related to ageing and change, and utilise various datasets. The framework of this study leans on the studies of socio-technical systems and multi-level perspective on transitions mainly developed by Frank Geels. Essential factors in transition from one socio-technological regime to another are the co-evolutionary processes between landscape changes, regime level and experimental niches. Landscape level changes, like the ageing of the population, destabilise the regime in the forms of coming pressures. This destabilization offers windows for opportunity to niche-innovations outside or at fringe of the regime, which, through their breakthrough, accelerate the transition process. However, the change is not easy because of various kinds of lock-ins and inertia, which tend to maintain the stability of the regime. In this dissertation, a constructionist approach of society is applied leaning mainly to the ideas of Anthony Giddens’ theory of structuration, with the dual nature of structures. The change is taking place in the interplay between actors and structures: structures shape people’s practices, but at the same time these practices constitute and reproduce social systems. Technology and other material aspects, as part of socio-technical systems, and the use of them, also take part in the structuration process. The findings of the study point out that co-evolutionary and co-effective relationships between economic, cultural, technological and institutional fields, as well as relationships between landscape changes, changes in the local and regime-level practices and rule systems, are a very complex and multi-level dynamic socio-technical phenomenon. At the landscape level of ageing, which creates the pressures and triggers to the regime change, there are three remarkable megatrends: demographic change, changes in the global economy and the development of technologies. These exert pressures to the socio-technical regime, which as a rule system is experiencing changes in the form of new markets and consumer habits, new ways of perceiving ageing, new models of organising the health care and other services and as new ways of considering innovation and innovativeness. There are also inner dynamics in the relationships between these aspects within the regime. These are interrelated and coconstructed: the prevailing perceptions of ageing and innovation, for instance, reflect the ageing policies, innovation policies, societal structures, organising models, technology and scientific discussion, and vice versa. Technology is part of the inner dynamics of the sociotechnological regime. Physical properties of the artefacts set limitations and opportunities with regard to their functions and uses. The use of and discussion about technology, contributes producing and reproducing the perceptions of old age. For societal transition, micro-level changes are also needed, in form of niche-innovations, for instance new services, organisational models or new technologies, Regimes, as stabilitystriven systems, tend to generate incremental innovations, but radically new innovations are generated in experimental niches protected from ‘normal’ market selection. The windows of opportunity for radical novelties may be opened if the circumstances are favourable for instance by tensions in the socio-technical regime affected by landscape level changes. This dissertation indicates that a change is taking place, firstly, in the dynamic interactionbetween levels, as a result of purposive action and governance to some extent. Breaking the inertia and using the window of opportunity for change and innovation offered by dynamics between levels, presupposes the actors’ special capabilities and actions such as dynamic capabilities and distance management. Secondly, the change is taking place the socio-technological negotiations inside the regime: interaction between technological and social, which is embodied in the use of technology. The use of technology includes small-level contextual scripts that also participate in forming broader societal scripts (for instance defining old age at the society level), which in their turn affect the formation of policies for innovation and ageing. Thirdly, the change is taking place by the means of active formation of the multi-actor innovation networks, where the role of distance management is crucial to facilitate the communication between actors coming from different backgrounds as well as to help the niches born outside the regime to utilise the window of opportunity offered by regime destabilisation. This dissertation has both theoretical and practical contributions. This study participates in the discussion of action-oriented view on transition by opening up of the socio-technological, coevolutionary processes of the multi-faceted phenomenon of ageing, which has lacked systematic analyses. The focus of this study, however, is not on the large-scale coordination and governance, but rather on opening up the incremental elements and structuration processes, which contribute to the transition little by little, and which can be affected to. This increases the practical importance of this dissertation, by highlighting the importance of very tiny, everyday elements in the change processes in the long run.

A Farewell to Flat Biology. Three-dimensional Cell Culture Models in Cancer Drug Target Identification and Validation

Cells of epithelial origin, e.g. from breast and prostate cancers, effectively differentiate into complex multicellular structures when cultured in three-dimensions (3D) instead of conventional two-dimensional (2D) adherent surfaces. The spectrum of different organotypic morphologies is highly dependent on the culture environment that can be either non-adherent or scaffold-based. When embedded in physiological extracellular matrices (ECMs), such as laminin-rich basement membrane extracts, normal epithelial cells differentiate into acinar spheroids reminiscent of glandular ductal structures. Transformed cancer cells, in contrast, typically fail to undergo acinar morphogenic patterns, forming poorly differentiated or invasive multicellular structures. The 3D cancer spheroids are widely accepted to better recapitulate various tumorigenic processes and drug responses. So far, however, 3D models have been employed predominantly in the Academia, whereas the pharmaceutical industry has yet to adopt a more widely and routine use. This is mainly due to poor characterisation of cell models, lack of standardised workflows and high throughput cell culture platforms, and the availability of proper readout and quantification tools. In this thesis, a complete workflow has been established entailing well-characterised 3D cell culture models for prostate cancer, a standardised 3D cell culture routine based on high-throughput-ready platform, automated image acquisition with concomitant morphometric image analysis, and data visualisation, in order to enable large-scale high-content screens. Our integrated suite of software and statistical analysis tools were optimised and validated using a comprehensive panel of prostate cancer cell lines and 3D models. The tools quantify multiple key cancer-relevant morphological features, ranging from cancer cell invasion through multicellular differentiation to growth, and detect dynamic changes both in morphology and function, such as cell death and apoptosis, in response to experimental perturbations including RNA interference and small molecule inhibitors. Our panel of cell lines included many non-transformed and most currently available classic prostate cancer cell lines, which were characterised for their morphogenetic properties in 3D laminin-rich ECM. The phenotypes and gene expression profiles were evaluated concerning their relevance for pre-clinical drug discovery, disease modelling and basic research. In addition, a spontaneous model for invasive transformation was discovered, displaying a highdegree of epithelial plasticity. This plasticity is mediated by an abundant bioactive serum lipid, lysophosphatidic acid (LPA), and its receptor LPAR1. The invasive transformation was caused by abrupt cytoskeletal rearrangement through impaired G protein alpha 12/13 and RhoA/ROCK, and mediated by upregulated adenylyl cyclase/cyclic AMP (cAMP)/protein kinase A, and Rac/ PAK pathways. The spontaneous invasion model tangibly exemplifies the biological relevance of organotypic cell culture models. Overall, this thesis work underlines the power of novel morphometric screening tools in drug discovery.

Analysis and validation of space averaged drag model for numerical simulations of gas-solid flows in fluidized beds

This thesis presents an approach for formulating and validating a space averaged drag model for coarse mesh simulations of gas-solid flows in fluidized beds using the two-fluid model. Proper modeling for fluid dynamics is central in understanding any industrial multiphase flow. The gas-solid flows in fluidized beds are heterogeneous and usually simulated with the Eulerian description of phases. Such a description requires the usage of fine meshes and small time steps for the proper prediction of its hydrodynamics. Such constraint on the mesh and time step size results in a large number of control volumes and long computational times which are unaffordable for simulations of large scale fluidized beds. If proper closure models are not included, coarse mesh simulations for fluidized beds do not give reasonable results. The coarse mesh simulation fails to resolve the mesoscale structures and results in uniform solids concentration profiles. For a circulating fluidized bed riser, such predicted profiles result in a higher drag force between the gas and solid phase and also overestimated solids mass flux at the outlet. Thus, there is a need to formulate the closure correlations which can accurately predict the hydrodynamics using coarse meshes. This thesis uses the space averaging modeling approach in the formulation of closure models for coarse mesh simulations of the gas-solid flow in fluidized beds using Geldart group B particles. In the analysis of formulating the closure correlation for space averaged drag model, the main parameters for the modeling were found to be the averaging size, solid volume fraction, and distance from the wall. The closure model for the gas-solid drag force was formulated and validated for coarse mesh simulations of the riser, which showed the verification of this modeling approach. Coarse mesh simulations using the corrected drag model resulted in lowered values of solids mass flux. Such an approach is a promising tool in the formulation of appropriate closure models which can be used in coarse mesh simulations of large scale fluidized beds.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Future of Newspaper Business Models in Finland

The purpose of this Master’s thesis was to study the business model development in Finnish newspaper industry during the next then years through scenario planning. The objective was to see how will the business models develop amidst the many changes in the industry, what factors are affecting the change, what are the implications of these changes for the players in the industry and how should the Finnish newspaper companies evolve in order to succeed in the future. In this thesis the business model change is studied based on all the elements of business models, as it was discovered that the industry is too often focusing on changes in only few of those elements and a more broader view can provide valuable information for the companies. The results revealed that the industry is affected by many changes during the next ten years. Scenario planning provides a good tool for analyzing this change and for developing valuable options for businesses. After conducting series of interviews and discovering forces affecting the change, four different scenarios were developed centered on the role that newspaper will take and the level at which they are providing the content in the future. These scenarios indicated that there are varieties of options in the way the business models may develop and that companies should start making decisions proactively in order to succeed. As the business model elements are interdepended, changes made in the other elements will affect the whole model, making these decisions about the role and level of content important for the companies. In the future, it is likely that the Finnish newspaper industry will include many different kinds of business models, some of which can be drastically different from the current ones and some of which can still be similar, but take better into account the new kind of media environment.

Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

Bayesian analysis of SEIR epidemic models

This thesis concerns the analysis of epidemic models. We adopt the Bayesian paradigm and develop suitable Markov Chain Monte Carlo (MCMC) algorithms. This is done by considering an Ebola outbreak in the Democratic Republic of Congo, former Zaïre, 1995 as a case of SEIR epidemic models. We model the Ebola epidemic deterministically using ODEs and stochastically through SDEs to take into account a possible bias in each compartment. Since the model has unknown parameters, we use different methods to estimate them such as least squares, maximum likelihood and MCMC. The motivation behind choosing MCMC over other existing methods in this thesis is that it has the ability to tackle complicated nonlinear problems with large number of parameters. First, in a deterministic Ebola model, we compute the likelihood function by sum of square of residuals method and estimate parameters using the LSQ and MCMC methods. We sample parameters and then use them to calculate the basic reproduction number and to study the disease-free equilibrium. From the sampled chain from the posterior, we test the convergence diagnostic and confirm the viability of the model. The results show that the Ebola model fits the observed onset data with high precision, and all the unknown model parameters are well identified. Second, we convert the ODE model into a SDE Ebola model. We compute the likelihood function using extended Kalman filter (EKF) and estimate parameters again. The motivation of using the SDE formulation here is to consider the impact of modelling errors. Moreover, the EKF approach allows us to formulate a filtered likelihood for the parameters of such a stochastic model. We use the MCMC procedure to attain the posterior distributions of the parameters of the SDE Ebola model drift and diffusion parts. In this thesis, we analyse two cases: (1) the model error covariance matrix of the dynamic noise is close to zero , i.e. only small stochasticity added into the model. The results are then similar to the ones got from deterministic Ebola model, even if methods of computing the likelihood function are different (2) the model error covariance matrix is different from zero, i.e. a considerable stochasticity is introduced into the Ebola model. This accounts for the situation where we would know that the model is not exact. As a results, we obtain parameter posteriors with larger variances. Consequently, the model predictions then show larger uncertainties, in accordance with the assumption of an incomplete model.