Refine
Has Fulltext
- yes (520) (remove)
Year of publication
Document Type
- Preprint (374)
- Doctoral Thesis (64)
- Postprint (32)
- Article (31)
- Monograph/Edited Volume (7)
- Master's Thesis (7)
- Bachelor Thesis (2)
- Conference Proceeding (2)
- Course Material (1)
Language
- English (476)
- German (40)
- French (3)
- Multiple languages (1)
Keywords
- random point processes (19)
- statistical mechanics (19)
- stochastic analysis (19)
- index (11)
- boundary value problems (10)
- elliptic operators (9)
- Fredholm property (8)
- cluster expansion (8)
- K-theory (7)
- manifolds with singularities (7)
Institute
- Institut für Mathematik (520) (remove)
The aim of this paper is to build and compare estimators of the infection parameter in the different phases of an epidemic (growth and extinction phases). The epidemic is modeled by a Markovian process of order d > 1 (allowing non-Markovian life spans), and can be written as a multitype branching process. We propose three estimators suitable for the different classes of criticality of the process, in particular for the subcritical case corresponding to the extinction phase. We prove their consistency and asymptotic normality for two asymptotics, when the number of ancestors (resp. number of generations) tends to infinity. We illustrate the asymptotic properties with simulated examples, and finally use our estimators to study the infection intensity in the extinction phase of the BSE epidemic in Great-Britain.
Estimation and testing of distributions in metric spaces are well known. R.A. Fisher, J. Neyman, W. Cochran and M. Bartlett achieved essential results on the statistical analysis of categorical data. In the last 40 years many other statisticians found important results in this field. Often data sets contain categorical data, e.g. levels of factors or names. There does not exist any ordering or any distance between these categories. At each level there are measured some metric or categorical values. We introduce a new method of scaling based on statistical decisions. For this we define empirical probabilities for the original observations and find a class of distributions in a metric space where these empirical probabilities can be found as approximations for equivalently defined probabilities. With this method we identify probabilities connected with the categorical data and probabilities in metric spaces. Here we get a mapping from the levels of factors or names into points of a metric space. This mapping yields the scale for the categorical data. From the statistical point of view we use multivariate statistical methods, we calculate maximum likelihood estimations and compare different approaches for scaling.
Estimation and testing of distributions in metric spaces are well known. R.A. Fisher, J. Neyman, W. Cochran and M. Bartlett achieved essential results on the statistical analysis of categorical data. In the last 40 years many other statisticians found important results in this field. Often data sets contain categorical data, e.g. levels of factors or names. There does not exist any ordering or any distance between these categories. At each level there are measured some metric or categorical values. We introduce a new method of scaling based on statistical decisions. For this we define empirical probabilities for the original observations and find a class of distributions in a metric space where these empirical probabilities can be found as approximations for equivalently defined probabilities. With this method we identify probabilities connected with the categorical data and probabilities in metric spaces. Here we get a mapping from the levels of factors or names into points of a metric space. This mapping yields the scale for the categorical data. From the statistical point of view we use multivariate statistical methods, we calculate maximum likelihood estimations and compare different approaches for scaling.
We consider an infinite system of non overlaping globules undergoing Brownian motions in R3. The term globules means that the objects we are dealing with are spherical, but with a radius which is random and time-dependent. The dynamics is modelized by an infinitedimensional Stochastic Differential Equation with local time. Existence and uniqueness of a strong solution is proven for such an equation with fixed deterministic initial condition. We also find a class of reversible measures.
Multitype branching processes and Feller diffusion processes are conditioned on very late extinction. The conditioned laws are expressed as Doob h-transforms of the unconditioned laws, and an interpretation of the conditioned paths for the branching process is given, via the immortal particle. We study different limits for the conditioned process (increasing delay of extinction, long-time behavior, scaling limit) and provide an exhaustive list of exchangeability results.
Zufällige Punktprozesse beschreiben eine (zufällige) zeitliche Abfolge von Ereignissen oder eine (zufällige) räumliche Anordnung von Objekten. Deren wichtigster Vertreter ist der Poissonprozess. Der Poissonprozess zum Intensitätsmaß, das Lebesgue-Maß ordnet jedem Gebiet sein Volumen zu, erzeugt lokal, d.h in einem beschränkten Gebiet B, gerade eine mit dem Volumen von B poissonverteilte Anzahl von Punkten, die identisch und unabhängig voneinander in B plaziert werden; im Mittel ist diese Anzahl (B). Ersetzt man durch ein Vielfaches a, so wird diese Anzahl mit dem a-fachen Mittelwert erzeugt. Poissonprozesse, die im gesamten Raum unendlich viele Punkte realisieren, enthalten bereits in einer einzigen Stichprobe genügend Informationen, um Statistik betreiben zu können: Bedingt man lokal bzgl. der Anzahl der Teilchen einer Stichprobe, so fragt man nach allen Punktprozessen, die eine solche Beobachtung hätten liefern können. Diese sind Limespunktprozesse zu dieser Beobachtung. Kommt mehr als einer in Frage, spricht man von einem Phasenübergang. Da die Menge dieser Limespunktprozesse konvex ist, fragt man nach deren Extremalpunkten, dem Rand. Im ersten Teil wird ein Poissonprozess für ein physikalisches Teilchenmodell für Bosonen konstruiert. Dieses erzeugt sogenannte Loops, das sind geschlossene Polygonzüge, die dadurch charakterisiert sind, dass man an einem Ort mit einem Punkt startet, den mit einem normalverteilten Schritt läuft und dabei nach einer gegebenen, aber zufälligen Anzahl von Schritten zum Ausgangspunkt zurückkehrt. Für verschiedene Beobachtungen von Stichproben werden zugehörige Limespunktprozesse diskutiert. Diese Beobachtungen umfassen etwa das Zählen der Loops gemäaß ihrer Länge, das Zählen der Loops insgesamt, oder das Zählen der von den Loops gemachten Schritte. Jede Wahl zieht eine charakteristische Struktur der invarianten Punktprozesse nach sich. In allen hiesigen Fällen wird ein charakteristischer Phasenübergang gezeigt und Extremalpunkte werden als spezielle Poissonprozesse identifiziert. Insbesondere wird gezeigt, wie die Wahl der Beobachtung die Länge der Loops beeinflusst. Geometrische Eigenschaften dieser Poissonprozesse sind der Gegenstand des zweiten Teils der Arbeit. Die Technik der Palmschen Verteilungen eines Punktprozesses erlaubt es, unter den unendlich vielen Loops einer Realisierung den typischen Loop herauszupicken, dessen Geometrie dann untersucht wird. Eigenschaften sind unter anderem die euklidische Länge eines Schrittes oder, nimmt man mehrere aufeinander folgende Schritte, das Volumen des von ihnen definierten Simplex. Weiterhin wird gezeigt, dass der Schwerpunkt eines typischen Loops normalverteilt ist mit einer festen Varianz. Der dritte und letzte Teil befasst sich mit der Konstruktion, den Eigenschaften und der Statistik eines neuartigen Punktprozesses, der Polyascher Summenprozess genannt wird. Seine Konstruktion verallgemeinert das Prinzip der Polyaschen Urne: Im Gegensatz zum Poissonprozess, der alle Punkte unabhängig und vor allem identisch verteilt, werden hier die Punkte nacheinander derart verteilt, dass der Ort, an dem ein Punkt plaziert wird, eine Belohnung auf die Wahrscheinlichkeit bekommt, nach der nachfolgende Punkte verteilt werden. Auf diese Weise baut der Polyasche Summenprozess "Türmchen", indem sich verschiedene Punkte am selben Ort stapeln. Es wird gezeigt, dass dennoch grundlegende Eigenschaften mit denjenigen des Poissonprozesses übereinstimmen, dazu gehören unendliche Teilbarkeit sowie Unabhängigkeit der Zuwächse. Zudem werden sein Laplace-Funktional sowie seine Palmsche Verteilung bestimmt. Letztere zeigt, dass die Höhe der Türmchen gerade geometrisch verteilt ist. Abschließend werden wiederum Statistiken, nun für den Summenprozess, diskutiert. Je nach Art der Beobachtung von der Stichprobe, etwa Anzahl, Gesamthöhe der Türmchen oder beides, gibt es in jedem der drei Fälle charakteristische Limespunktprozesse und es stellt sich heraus, dass die zugehörigen Extremalverteilungen wiederum Polyasche Summenprozesse sind.
We consider a class of infinite-dimensional diffusions where the interaction between the components is both spatial and temporal. We start the system from a Gibbs measure with finiterange uniformly bounded interaction. Under suitable conditions on the drift, we prove that there exists t0 > 0 such that the distribution at time t = t0 is a Gibbs measure with absolutely summable interaction. The main tool is a cluster expansion of both the initial interaction and certain time-reversed Girsanov factors coming from the dynamics.
Studying the influence of the updating scheme for MCMC algorithm on spatially extended models is a well known problem. For discrete-time interacting particle systems we study through simulations the effectiveness of a synchronous updating scheme versus the usual sequential one. We compare the speed of convergence of the associated Markov chains from the point of view of the time-to-coalescence arising in the coupling-from-the-past algorithm. Unlike the intuition, the synchronous updating scheme is not always the best one. The distribution of the time-to-coalescence for these spatially extended models is studied too.
Aus dem Inhalt: Einleitung Kapitel 1. Starke Gesetze der Grossen Zahlen 1. SGGZ unter Wachstumsbedingungen an die p-ten Momente 2. SGGZ für identisch verteilte Zufallsvariablen 3. SGGZ für Prozesse mit *-mixing-Eigenschaft Kapitel 2. Einführung zu diskreten (Sub-,Super-)Martingalen 1. Vorhersagbarkeit 2. gestoppte (Sub-,Super-)Martingale 3. Upcrossings 4. Konvergenzsätze 5. Doob-Zerlegung 6. Eine äquivalente Definition eines (Sub-)Martingals Kapitel 3. Martingale und gleichgradige Integrierbarkeit 1. Gleichmäßige(-f¨ormige,-gradige) Integrierbarkeit 2. gleichgradig integrierbare Martingale Kapitel 4. Martingale und das SGGZ Kapitel 5.”reversed“ (Sub-,Super-)Martingale 1. Konvergenzsätze Kapitel 6. (Sub-,Super-)Martingale mit gerichteter Indexmenge 1. Äquivalente Formulierung eines (Sub-)Martingals 2. Konvergenzsätze Kapitel 7. Quasimartingale,Amarts und Semiamarts 1. Konvergenzsätze 2. Riesz-Zerlegung 3. Doob-Zerlegung Kapitel 8. Amarts und das SGGZ Kapitel 9.”reversed“ Amarts und Semiamarts 1. Konvergenzsätze 2.”Aufwärts“- gegen ”Abwärts“-Adaptiertheit 3. Riesz-Zerlegung 4. Stabilitätsanalyse Kapitel 10. Amarts mit gerichteter Indexmenge 1. Konvergenzsätze 2. Riesz-Zerlegung Anhang A. zur Existenz einer Folge unabhängiger Zufallsvariablen B. Konvergenz
We give the explicit solution for the minimax linear estimate. For scale dependent models an empirical minimax linear estimates is de¯ned and we prove that these estimates are Stein's estimates.
Aus dem Inhalt: 0.1 Danksagung 0.2 Einleitung 1 Allgemeines und Grundlagen 1.1 Die Brownsche Bewegung 2 Die Dualitätsformel des Wienermaßes 2.1 Wienermaß erfüllt Dualitätsformel 2.2 Dualitätsformel charakterisiert Wienermaß 3 Die diskrete Dualitätsformel der Irrfahrt 3.1 Verallgemeinerte symmetrische Irrfahrt erfüllt diskrete Dualitätsformel 3.2 Diskrete Dualitätsformel charakterisiert verallgemeinerte symmetrische Irrfahrt 4 Donskers Theorem und die Dualitätsformeln 4.1 Straffheit der renormierten stetigen Irrfahrt 4.2 Konvergenz der Irrfahrt 5 Anhang
Aus dem Inhalt: 1 Abraham Wald (1902-1950) 2 Einführung der Grundbegriffe. Einige technische bekannte Ergebnisse 2.1 Martingal und Doob-Ungleichung 2.2 Brownsche Bewegung und spezielle Martingale 2.3 Gleichgradige Integrierbarkeit von Prozessen 2.4 Gestopptes Martingal 2.5 Optionaler Stoppsatz von Doob 2.6 Lokales Martingal 2.7 Quadratische Variation 2.8 Die Dichte der ersten einseitigen Überschreitungszeit der Brown- schen Bewegung 2.9 Waldidentitäten für die Überschreitungszeiten der Brownschen Bewegung 3 Erste Waldidentität 3.1 Burkholder, Gundy und Davis Ungleichungen der gestoppten Brown- schen Bewegung 3.2 Erste Waldidentität für die Brownsche Bewegung 3.3 Verfeinerungen der ersten Waldidentität 3.4 Stärkere Verfeinerung der ersten Waldidentität für die Brown- schen Bewegung 3.5 Verfeinerung der ersten Waldidentität für spezielle Stoppzeiten der Brownschen Bewegung 3.6 Beispiele für lokale Martingale für die Verfeinerung der ersten Waldidentität 3.7 Überschreitungszeiten der Brownschen Bewegung für nichtlineare Schranken 4 Zweite Waldidentität 4.1 Zweite Waldidentität für die Brownsche Bewegung 4.2 Anwendungen der ersten und zweitenWaldidentität für die Brown- schen Bewegung 5 Dritte Waldidentität 5.1 Dritte Waldidentität für die Brownsche Bewegung 5.2 Verfeinerung der dritten Waldidentität 5.3 Eine wichtige Voraussetzung für die Verfeinerung der drittenWal- didentität 5.4 Verfeinerung der dritten Waldidentität für spezielle Stoppzeiten der Brownschen Bewegung 6 Waldidentitäten im Mehrdimensionalen 6.1 Erste Waldidentität im Mehrdimensionalen 6.2 Zweite Waldidentität im Mehrdimensionalen 6.3 Dritte Waldidentität im Mehrdimensionalen 7 Appendix
We study resonances for the generator of a diffusion with small noise in R(d) : L = -∈∆ + ∇F * ∇, when the potential F grows slowly at infinity (typically as a square root of the norm). The case when F grows fast is well known, and under suitable conditions one can show that there exists a family of exponentially small eigenvalues, related to the wells of F. We show that, for an F with a slow growth, the spectrum is R+, but we can find a family of resonances whose real parts behave as the eigenvalues of the "quick growth" case, and whose imaginary parts are small.
We consider a class of ergodic Hamilton-Jacobi-Bellman (HJB) equations, related to large time asymptotics of non-smooth multiplicative functional of difusion processes. Under suitable ergodicity assumptions on the underlying difusion, we show existence of these asymptotics, and that they solve the related HJB equation in the viscosity sense.
A multitype Dawson-Watanabe process is conditioned, in subcritical and critical cases, on non-extinction in the remote future. On every nite time interval, its distribution law is absolutely continuous with respect to the law of the unconditioned process. A martingale problem characterization is also given. The explicit form of the Laplace functional of the conditioned process is used to obtain several results on the long time behaviour of the mass of the conditioned and unconditioned processes. The general case is considered first, where the mutation matrix which modelizes the interaction between the types, is irreducible. Several two-type models with decomposable mutation matrices are also analysed.
The accelerated life time model is considered. First, test procedures for testing the parameter of a parametric acceleration function is investigated; this is done under the assumption of parametric and nonparametric baseline distribution. Further, based on nonparametric estimators for regression functions tests are proposed for checking whether a parametric acceleration function is appropriate to model the influence of the covariates. Resampling procedures are discussed for the realization of these methods. Simulations complete the considerations.
We consider the problem of testing whether the density of a mul- tivariate random variable can be expressed by a prespecified copula function and the marginal densities. The proposed test procedure is based on the asymptotic normality of the properly standardized integrated squared distance between a multivariate kernel density estimator and an estimator of its expectation under the hypothesis. The test of independence is a special case of this approach.
We consider an infinite system of hard balls in Rd undergoing Brownian motions and submitted to a pair potential with infinite range and quasi polynomial decay. It is modelized by an infinite-dimensional Stochastic Differential Equation with an infinite-dimensional local time term. Existence and uniqueness of a strong solution is proven for such an equation with deterministic initial condition. We also show that the set of all equilibrium measures, solution of a Detailed Balance Equation, coincides with the set of canonical Gibbs measures associated to the hard core potential.
We consider a system of infinitely many hard balls in Rd undergoing Brownian motions and submitted to a smooth pair potential. It is modelized by an infinite-dimensional Stochastic Differential Equation with a local time term. We prove that the set of all equilibrium measures, solution of a Detailed Balance Equation, coincides with the set of canonical Gibbs measures associated to the hard core potential added to the smooth interaction potential.
The aim of this paper is to demonstrate that nonparametric smoothing methods for estimating functions can be an useful tool in the analysis of life time data. After stating some basic notations we will present a data example. Applying standard parametric methods to these data we will see that this approach fails - basic features of the underlying functions are not reflected by their estimates. Our proposal is to use nonparametric estimation methods. These methods are explained in section 2. Nonparametric approaches are better in the sense that they are more flexible, and misspecifications of the model are avoided. But, parametric models have the advantage that the parameters can be interpreted. So, finally, we will formulate a test procedure to check whether a parametric or a nonparametric model is appropriate.
We give a necessary and sufficient condition for the existence of an increasing coupling of N (N >= 2) synchronous dynamics on S-Zd (PCA). Increasing means the coupling preserves stochastic ordering. We first present our main construction theorem in the case where S is totally ordered; applications to attractive PCAs are given. When S is only partially ordered, we show on two examples that a coupling of more than two synchronous dynamics may not exist. We also prove an extension of our main result for a particular class of partially ordered spaces.
We first introduce some coupling of a finite number of Probabilistic Cellular Automata dynamics (PCA), preserving the stochastic ordering. Using this tool, for a general attractive probabilistic cellular automata on SZd, where S is finite, we prove that a condition (A) is equivalent to the (time-) convergence towards equilibrium of this Markovian parallel dynamics, in the uniform norm, exponentially fast. This condition (A) means the exponential decay of the influence from the boundary for the invariant measures of the system restricted to finite ‘box’-volume. For a class of reversible PCA dynamics on {−1, +1}Zd , with a naturally associated Gibbsian potential ϕ, we prove that a Weak Mixing condition for ϕ implies the validity of the assumption (A); thus the ‘exponential ergodicity’ of the dynamics towards the unique Gibbs measure associated to ϕ holds. On some particular examples of this PCA class, we verify that our assumption (A) is weaker than the Dobrushin-Vasershtein ergodicity condition. For some special PCA, the ‘exponential ergodicity’ holds as soon as there is no phase transition.
The two and k-sample tests of equality of the survival distributions against the alternatives including cross-effects of survival functions, proportional and monotone hazard ratios, are given for the right censored data. The asymptotic power against approaching alternatives is investigated. The tests are applied to the well known chemio and radio therapy data of the Gastrointestinal Tumor Study Group. The P-values for both proposed tests are much smaller then in the case of other known tests. Differently from the test of Stablein and Koutrouvelis the new tests can be applied not only for singly but also to randomly censored data.
We consider a nonparametric survival model with random censoring. To test whether the hazard rate has a parametric form the unknown hazard rate is estimated by a kernel estimator. Based on a limit theorem stating the asymptotic normality of the quadratic distance of this estimator from the smoothed hypothesis an asymptotic ®-test is proposed. Since the test statistic depends on the maximum likelihood estimator for the unknown parameter in the hypothetical model properties of this parameter estimator are investigated. Power considerations complete the approach.
The dependence between survival times and covariates is described e.g. by proportional hazard models. We consider partly parametric Cox models and discuss here the estimation of interesting parameters. We represent the ma- ximum likelihood approach and extend the results of Huang (1999) from linear to nonlinear parameters. Then we investigate the least squares esti- mation and formulate conditions for the a.s. boundedness and consistency of these estimators.
We give a survey on procedures for testing functions which are based on quadratic deviation measures. The following problems are considered: Testing whether a density function lies in a parametric class of functions, whether continuous random variables are independent; testing cell probabilities and independence in sparse data sets; testing the parametric fit of a regression homoscedasticity in a regression model and testing the hazard rate in survival models with censoring and with and without covariates.
In diesem Beitrag wird der Zusammenhang zwischen Algebrodifferentialgleichungen (ADGL) und Vektorfeldern auf Mannigfaltigkeiten untersucht. Dazu wird zunächst der Begriff der regulären ADGL eingeführt, wobei unter eirter regulären ADGL eine ADGL verstanden wird, deren Lösungsmenge identisch mit der Lösungsmenge eines Vektorfeldes ist. Ausgehend von bekannten Aussagen über die Lösungsmenge eines Vektorfeldes werden analoge Aussagen für die Lösungsmenge einer regulären ADGL abgeleitet. Es wird eine Reduktionsmethode angegeben, die auf ein Kriterium für die Begularität einer ADGL und auf die Definition des Index einer nichtlinearen ADGL führt. Außerdem wird gezeigt, daß beliebige Vektorfelder durch reguläre ADGL so realisiert werden können, daß die Lösungsmenge des Vektorfeldes mit der der realisierenden ADGL identisch ist. Abschließend werden die für autonome ADGL gewonnenen Aussagen auf den Fall der nichtautonomen ADGL übertragen.
A theoretical famework for the investigation of the qualitative behavior of differential-algebraic equations (DAEs) near an equilibrium point is established. The key notion of our approach is the notion of regularity. A DAE is called regular locally around an equilibrium point if there is a unique vector field such that the solutions of the DAE and the vector field are in one-to-one correspondence in a neighborhood of this equili Drium point. Sufficient conditions for the regularity of an equilibrium point are stated. This in turn allows us to translate several local results, as formulated for vector fields, to DAEs that are regular locally around a g: ven equilibrium point (e.g. Local Stable and Unstable Manifold Theorem, Hopf theorem). It is important that ihese theorems are stated in terms of the given problem and not in terms of the corresponding vector field.
An existence and uniqueness theory is developed for general nonlinear and nonautonomous differential-algebraic equations (DAEs) by exploiting their underlying differential-geometric structure. A DAE is called regular if there is a unique nonautonomous vector field such that the solutions of the DAE and the solutions of the vector field are in one-to-one correspondence. Sufficient conditions for regularity of a DAE are derived in terms of constrained manifolds. Based on this differential-geometric characterization, existence and uniqueness results are stated for regular DAEs. Furthermore, our not ons are compared with techniques frequently used in the literature such as index and solvability. The results are illustrated in detail by means of a simple circuit example.
The subject of this paper is the relation of differential-algebraic equations (DAEs) to vector fields on manifolds. For that reason, we introduce the notion of a regular DAE as a DAE to which a vector field uniquely corresponds. Furthermore, a technique is described which yields a family of manifolds for a given DAE. This socalled family of constraint manifolds allows in turn the formulation of sufficient conditions for the regularity of a DAE. and the definition of the index of a regular DAE. We also state a method for the reduction of higher-index DAEs to lowsr-index ones that can be solved without introducing additional constants of integration. Finally, the notion of realizability of a given vector field by a regular DAE is introduced, and it is shown that any vector field can be realized by a regular DAE. Throughout this paper the problem of path-tracing is discussed as an illustration of the mathematical phenomena.
Technical and physical systems, especially electronic circuits, are frequently modeled as a system of differential and nonlinear implicit equations. In the literature such systems of equations are called differentialalgebraic equations (DAEs). It turns out that the numerical and analytical properties of a DAE depend on an integer called the index of the problem. For example, the well-known BDF method of Gear can be applied, in general, to a DAE only if the index does not exceed one. In this paper we give a geometric interpretation of higherindex DAEs and indicate problems arising in connection with such DAEs by means of several examples.
This thesis is concerned with the issue of extinction of populations composed of different types of individuals, and their behavior before extinction and in case of a very late extinction. We approach this question firstly from a strictly probabilistic viewpoint, and secondly from the standpoint of risk analysis related to the extinction of a particular model of population dynamics. In this context we propose several statistical tools. The population size is modeled by a branching process, which is either a continuous-time multitype Bienaymé-Galton-Watson process (BGWc), or its continuous-state counterpart, the multitype Feller diffusion process. We are interested in different kinds of conditioning on non-extinction, and in the associated equilibrium states. These ways of conditioning have been widely studied in the monotype case. However the literature on multitype processes is much less extensive, and there is no systematic work establishing connections between the results for BGWc processes and those for Feller diffusion processes. In the first part of this thesis, we investigate the behavior of the population before its extinction by conditioning the associated branching process X_t on non-extinction (X_t≠0), or more generally on non-extinction in a near future 0≤θ<∞ (X_{t+θ}≠0), and by letting t tend to infinity. We prove the result, new in the multitype framework and for θ>0, that this limit exists and is non-degenerate. This reflects a stationary behavior for the dynamics of the population conditioned on non-extinction, and provides a generalization of the so-called Yaglom limit, corresponding to the case θ=0. In a second step we study the behavior of the population in case of a very late extinction, obtained as the limit when θ tends to infinity of the process conditioned by X_{t+θ}≠0. The resulting conditioned process is a known object in the monotype case (sometimes referred to as Q-process), and has also been studied when X_t is a multitype Feller diffusion process. We investigate the not yet considered case where X_t is a multitype BGWc process and prove the existence of the associated Q-process. In addition, we examine its properties, including the asymptotic ones, and propose several interpretations of the process. Finally, we are interested in interchanging the limits in t and θ, as well as in the not yet studied commutativity of these limits with respect to the high-density-type relationship between BGWc processes and Feller processes. We prove an original and exhaustive list of all possible exchanges of limit (long-time limit in t, increasing delay of extinction θ, diffusion limit). The second part of this work is devoted to the risk analysis related both to the extinction of a population and to its very late extinction. We consider a branching population model (arising notably in the epidemiological context) for which a parameter related to the first moments of the offspring distribution is unknown. We build several estimators adapted to different stages of evolution of the population (phase growth, decay phase, and decay phase when extinction is expected very late), and prove moreover their asymptotic properties (consistency, normality). In particular, we build a least squares estimator adapted to the Q-process, allowing a prediction of the population development in the case of a very late extinction. This would correspond to the best or to the worst-case scenario, depending on whether the population is threatened or invasive. These tools enable us to study the extinction phase of the Bovine Spongiform Encephalopathy epidemic in Great Britain, for which we estimate the infection parameter corresponding to a possible source of horizontal infection persisting after the removal in 1988 of the major route of infection (meat and bone meal). This allows us to predict the evolution of the spread of the disease, including the year of extinction, the number of future cases and the number of infected animals. In particular, we produce a very fine analysis of the evolution of the epidemic in the unlikely event of a very late extinction.
On a method for solution of the ordinary differential equations connected with Huygens' equations
(2010)
We establish elements of a new approach to ellipticity and parametrices within operator algebras on manifolds with higher singularities, only based on some general axiomatic requirements on parameter-dependent operators in suitable scales of spaes. The idea is to model an iterative process with new generations of parameter-dependent operator theories, together with new scales of spaces that satisfy analogous requirements as the original ones, now on a corresponding higher level. The "full" calculus involves two separate theories, one near the tip of the corner and another one at the conical exit to infinity. However, concerning the conical exit to infinity, we establish here a new concrete calculus of edge-degenerate operators which can be iterated to higher singularities.
This thesis considers on the one hand the construction of point processes via conditional intensities, motivated by the partial Integration of the Campbell measure of a point process. Under certain assumptions on the intensity the existence of such a point process is shown. A fundamental example turns out to be the Pólya sum process, whose conditional intensity is a generalisation of the Pólya urn dynamics. A Cox process representation for that point process is shown. A further process considered is a Poisson process of Gaussian loops, which represents a noninteracting particle system derived from the discussion of indistinguishable particles. Both processes are used to define particle systems locally, for which thermodynamic limits are determined.
We construct elliptic elements in the algebra of (classical pseudo-differential) operators on a manifold M with conical singularities. The ellipticity of any such operator A refers to a pair of principal symbols (σ0, σ1) where σ0 is the standard (degenerate) homogeneous principal symbol, and σ1 is the so-called conormal symbol, depending on the complex Mellin covariable z. The conormal symbol, responsible for the conical singularity, is operator-valued and acts in Sobolev spaces on the base X of the cone. The σ1-ellipticity is a bijectivity condition for all z of real part (n + 1)/2 − γ, n = dimX, for some weight γ. In general, we have to rule out a discrete set of exceptional weights that depends on A. We show that for every operator A which is elliptic with respect to σ0, and for any real weight γ there is a smoothing Mellin operator F in the cone algebra such that A + F is elliptic including σ1. Moreover, we apply the results to ellipticity and index of (operator-valued) edge symbols from the calculus on manifolds with edges.
The Cauchy problem of the vacuum Einstein's equations aims to find a semimetric g(αβ) of a spacetime with vanishing Ricci curvature Rα,β and prescribed initial data. Under the harmonic gauge condition, the equations Rα,β = 0 are transferred into a system of quasi-linear wave equations which are called the reduced Einstein equations. The initial data for Einstein's equations are a proper Riemannian metric h(αβ) and a second fundamental form K(αβ). A necessary condition for the reduced Einstein equation to satisfy the vacuum equations is that the initial data satisfy Einstein constraint equations. Hence the data (h(αβ),K(αβ)) cannot serve as initial data for the reduced Einstein equations. Previous results in the case of asymptotically flat spacetimes provide a solution to the constraint equations in one type of Sobolev spaces, while initial data for the evolution equations belong to a different type of Sobolev spaces. The goal of our work is to resolve this incompatibility and to show that under the harmonic gauge the vacuum Einstein equations are well-posed in one type of Sobolev spaces.
We give a survey on the calculus of (pseudo-differential) boundary value problems with the transmision property at the boundary, and ellipticity in the Shapiro-Lopatinskij sense. Apart from the original results of the work of Boutet de Monvel we present an approach based on the ideas of the edge calculus. In a final section we introduce symbols with the anti-transmission property.
We give a brief survey on some new developments on elliptic operators on manifolds with polyhedral singularities. The material essentially corresponds to a talk given by the author during the Conference “Elliptic and Hyperbolic Equations on Singular Spaces”, October 27 - 31, 2008, at the MSRI, University of Berkeley.
We prove a local in time existence and uniqueness theorem of classical solutions of the coupled Einstein{Euler system, and therefore establish the well posedness of this system. We use the condition that the energy density might vanish or tends to zero at infinity and that the pressure is a certain function of the energy density, conditions which are used to describe simplified stellar models. In order to achieve our goals we are enforced, by the complexity of the problem, to deal with these equations in a new type of weighted Sobolev spaces of fractional order. Beside their construction, we develop tools for PDEs and techniques for hyperbolic and elliptic equations in these spaces. The well posedness is obtained in these spaces.
We consider a mixed problem for a degenerate differentialoperator equation of higher order. We establish some embedding theorems in weighted Sobolev spaces and show existence and uniqueness of the generalized solution of this problem. We also give a description of the spectrum for the corresponding operator.
The ellipticity of boundary value problems on a smooth manifold with boundary relies on a two-component principal symbolic structure (σψ; σ∂), consisting of interior and boundary symbols. In the case of a smooth edge on manifolds with boundary we have a third symbolic component, namely the edge symbol σ∧, referring to extra conditions on the edge, analogously as boundary conditions. Apart from such conditions in integral form' there may exist singular trace conditions, investigated in [6] on closed' manifolds with edge. Here we concentrate on the phenomena in combination with boundary conditions and edge problem.
We establish elements of a new approch to ellipticity and parametrices within operator algebras on a manifold with higher singularities, only based on some general axiomatic requirements on parameter-dependent operators in suitable scales of spaces. The idea is to model an iterative process with new generations of parameter-dependent operator theories, together with new scales of spaces that satisfy analogous requirements as the original ones, now on a corresponding higher level. The “full” calculus is voluminous; so we content ourselves here with some typical aspects such as symbols in terms of order reducing families, classes of relevant examples, and operators near the conical exit to infinity.