Gibbs Voronoi tessellations: Modeling, Simulation, Estimation.
|
|
- Eric Green
- 5 years ago
- Views:
Transcription
1 Gibbs Voronoi tessellations: Modeling, Simulation, Estimation. Frédéric Lavancier, Laboratoire Jean Leray, Nantes (France) Work with D. Dereudre (LAMAV, Valenciennes, France). SPA, Berlin, July 27-31, 2009.
2 Introduction Voronoï tessellations : applications in Astronomy, Biology, Physics, etc.
3 Introduction Voronoï tessellations : applications in Astronomy, Biology, Physics, etc. Studied as a random object : based on Poisson point Processes.
4 Introduction Voronoï tessellations : applications in Astronomy, Biology, Physics, etc. Studied as a random object : based on Poisson point Processes.
5 Introduction Voronoï tessellations : applications in Astronomy, Biology, Physics, etc. Studied as a random object : based on Poisson point Processes. Drawback : Strong independent structures coming from the Poisson process
6 Introduction Voronoï tessellations : applications in Astronomy, Biology, Physics, etc. Studied as a random object : based on Poisson point Processes. Drawback : Strong independent structures coming from the Poisson process Interactions between the cells?
7 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations.
8 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions?
9 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction.
10 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction. Hardcore interaction (some Voronoi tessellations are forbidden)
11 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction. Hardcore interaction (some Voronoi tessellations are forbidden) Existence of models.
12 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction. Hardcore interaction (some Voronoi tessellations are forbidden) Existence of models. Unicity of Gibbs measures.
13 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction. Hardcore interaction (some Voronoi tessellations are forbidden) Existence of models. Unicity of Gibbs measures. Simulations.
14 Introduction One solution : Gibbs modifications of Poisson Voronoï tessellations. Questions : What kind of interactions? Smooth interaction. Hardcore interaction (some Voronoi tessellations are forbidden) Existence of models. Unicity of Gibbs measures. Simulations. Parametric estimations.
15 2 Definitions
16 Notations B(R 2 ) denotes the space of bounded sets in R 2.
17 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <.
18 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <. γ Λ is the restriction of γ on Λ : γ Λ = γ Λ.
19 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <. γ Λ is the restriction of γ on Λ : γ Λ = γ Λ. Vor(γ) : Voronoï tessellation coming from γ
20 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <. γ Λ is the restriction of γ on Λ : γ Λ = γ Λ. Vor(γ) : Voronoï tessellation coming from γ λ is the Lebesgue measure on R 2.
21 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <. γ Λ is the restriction of γ on Λ : γ Λ = γ Λ. Vor(γ) : Voronoï tessellation coming from γ λ is the Lebesgue measure on R 2. For z > 0, π z : Poisson point process with intensity zλ
22 Notations B(R 2 ) denotes the space of bounded sets in R 2. M(R 2 ) is the space of locally finite point configurations γ in R 2 : γ R 2, such that for all Λ B(R 2 ), card(γ Λ) <. γ Λ is the restriction of γ on Λ : γ Λ = γ Λ. Vor(γ) : Voronoï tessellation coming from γ λ is the Lebesgue measure on R 2. For z > 0, π z : Poisson point process with intensity zλ πλ z : πz restricted on Λ.
23 Gibbs measures Let (H Λ ) Λ B(R 2 ) be a family of energies H Λ : M(Λ) M(Λ c ) R {+ } (γ Λ, γ Λ c) H Λ (γ Λ γ Λ c) We suppose that it is compatible. For every Λ Λ H Λ (γ Λ γ Λ c) = H Λ (γ Λ γ Λ c) + ϕ Λ,Λ (γ Λ c).
24 Gibbs measures Let (H Λ ) Λ B(R 2 ) be a family of energies H Λ : M(Λ) M(Λ c ) R {+ } (γ Λ, γ Λ c) H Λ (γ Λ γ Λ c) We suppose that it is compatible. For every Λ Λ H Λ (γ Λ γ Λ c) = H Λ (γ Λ γ Λ c) + ϕ Λ,Λ (γ Λ c). Definition A probability measure P on M(R 2 ) is a Gibbs measure for z > 0 and (H Λ ) if for every Λ B(R 2 ) and P -almost every γ Λ c P (dγ Λ γ Λ c) = 1 Z Λ (γ Λ c) e H Λ(γ Λ γ Λc ) π z Λ(dγ Λ ), where Z Λ (γ Λ c) = e H Λ(γ Λ γ Λ c ) π Λ (dγ Λ ).
25 A typical energy of a Voronoï tessellation : H Λ (γ Λ γ Λ c) = C Vor(γ) C Λ V 1 (C) + C,C Vor(γ) C and C are neighbors (C C ) Λ V 2 (C, C ).
26 A typical energy of a Voronoï tessellation : H Λ (γ Λ γ Λ c) = C Vor(γ) C Λ V 1 (C) + C,C Vor(γ) C and C are neighbors (C C ) Λ V 2 (C, C ). Our guiding example : + if h min (C) ε + if h V 1 (C) = max (C) α + if h 2 max(c)/vol(c) B 0 otherwise x h min C h max 0 < ε < α, B > 1/2 3 ;
27 A typical energy of a Voronoï tessellation : H Λ (γ Λ γ Λ c) = C Vor(γ) C Λ V 1 (C) + C,C Vor(γ) C and C are neighbors (C C ) Λ V 2 (C, C ). Our guiding example : + if h min (C) ε + if h V 1 (C) = max (C) α + if h 2 max(c)/vol(c) B 0 otherwise x h min C h max 0 < ε < α, B > 1/2 3 ; ( max(vol(c), Vol(C V 2 (C, C ) 1 )) ) = θ min(vol(c), Vol(C )) 1 2, θ R
28 Existence results First existence results (bounded interactions) : Bertin, Billiot and Drouilhet, : 1) Existence of nearest-neighbors spatial Gibbs models, Adv. Appl. Prob. (SGSA) (1999) 31, ) Existence of Delaunay pairwise Gibbs point process with superstable component, J. Stat. phys. (1999) Vol 95 3/
29 Existence results First existence results (bounded interactions) : Bertin, Billiot and Drouilhet, : 1) Existence of nearest-neighbors spatial Gibbs models, Adv. Appl. Prob. (SGSA) (1999) 31, ) Existence of Delaunay pairwise Gibbs point process with superstable component, J. Stat. phys. (1999) Vol 95 3/ Existence results with hardcore interactions (B = + ) : Dereudre, Gibbs Delaunay tessellations with geometric hardcore conditions, J. Stat. Phys.,(2008) 131,
30 Existence results First existence results (bounded interactions) : Bertin, Billiot and Drouilhet, : 1) Existence of nearest-neighbors spatial Gibbs models, Adv. Appl. Prob. (SGSA) (1999) 31, ) Existence of Delaunay pairwise Gibbs point process with superstable component, J. Stat. phys. (1999) Vol 95 3/ Existence results with hardcore interactions (B = + ) : Dereudre, Gibbs Delaunay tessellations with geometric hardcore conditions, J. Stat. Phys.,(2008) 131, In the general case : Dereudre, Drouilhet and Georgii, Existence of Gibbs process with stable cluster interactions, preprint.
31 Unicity For the interaction given before : A Gibbs measure exists but we don t know if it is unique or not (phase transition problem!)
32 3 Simulation
33 Simulations Strong hardcore interaction
34 Simulations Strong hardcore interaction Rigidity of the tessellation
35 Simulations Strong hardcore interaction Rigidity of the tessellation Several difficulties for the simulations.
36 Simulations Strong hardcore interaction Rigidity of the tessellation Several difficulties for the simulations. Birth-death-move MCMC algorithm on [0, 1] 2 : 1 Draw independently a and b uniformly on [0, 1]. 2 If a < 1/3 then generate x uniformly on [0, 1] 2 and if b < f(γ + x)z, then γ+x γ otherwise "do nothing". (n + 1)f(γ) 3 If 1/3 < a < 2/3 then generate x on γ and if b < nf(γ x), then γ x γ otherwise "do nothing". f(γ)z 4 If a > 2/3 then generate x on γ, y N (x, σ 2 ) and if b < f(γ x + y), then γ x+y γ otherwise "do nothing". f(γ)
37 Introduction Definitions Simulation Estimation Theory Examples of simulations We fix z = 100, ε = 0, α = 0.05 : B = +, θ = 0.5 B = 1, θ = 0.5 B = 0.625, θ = 0.5 B = +, θ = 0.5 B = 1, θ = 0.5 B = 0.625, θ = 0.5
38 Introduction Definitions Simulation Estimation Theory Monitoring control B = +, θ = 0.5 B = 1, θ = 0.5 B = 0.625, θ = 0.5 B = +, θ = 0.5 B = 1, θ = 0.5 B = 0.625, θ = 0.5
39 4 Estimation
40 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure.
41 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure. Hardcore parameters : ε, α and B. Empirical extremum hardcore parameters.
42 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure. Hardcore parameters : ε, α and B. Empirical extremum hardcore parameters. Smooth parameters : z and θ. Pseudolikelihood procedure.
43 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure. Hardcore parameters : ε, α and B. Empirical extremum hardcore parameters. Smooth parameters : z and θ. Pseudolikelihood procedure. Why the pseudo and not the MLE?
44 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure. Hardcore parameters : ε, α and B. Empirical extremum hardcore parameters. Smooth parameters : z and θ. Pseudolikelihood procedure. Why the pseudo and not the MLE? MLE is too time consuming (because of the estimation by simulations of the normalizing constant).
45 Pseudo-likelihood Estimation The aim : Estimate the parameters of the interaction from one realization γ of the Gibbs measure. Hardcore parameters : ε, α and B. Empirical extremum hardcore parameters. Smooth parameters : z and θ. Pseudolikelihood procedure. Why the pseudo and not the MLE? MLE is too time consuming (because of the estimation by simulations of the normalizing constant). Pseudo is proved to be asymptotically consistent and normal in most cases. Bibliography : Besag (1975), Jensen and Moller (1991), Jensen and Kunsch (1994), Mase (1995), Billiot, Coeurjolly and Drouilhet (2008), Dereudre and L. (2009).
46 Practical estimation procedures Let Λ n = [ n, n] 2 be the observation window and γ a realization of the Gibbs measure P.
47 Practical estimation procedures Let Λ n = [ n, n] 2 be the observation window and γ a realization of the Gibbs measure P. -Hardcore parameter estimators : ˆε = min{h min (C), C V or(γ) and C Λ n }, ˆα = max{h max (C), C V or(γ) and C Λ n }, ˆB = max{h 2 max(c)/vol(c), C V or(γ) and C Λ n }.
48 Practical estimation procedures Let Λ n = [ n, n] 2 be the observation window and γ a realization of the Gibbs measure P. -Hardcore parameter estimators : ˆε = min{h min (C), C V or(γ) and C Λ n }, ˆα = max{h max (C), C V or(γ) and C Λ n }, ˆB = max{h 2 max(c)/vol(c), C V or(γ) and C Λ n }. -Smooth parameter estimators : (ẑ, ˆθ) = argmin z,θ P LL Λn (γ, z, θ, ˆε, ˆα, ˆB),
49 Practical estimation procedures Let Λ n = [ n, n] 2 be the observation window and γ a realization of the Gibbs measure P. -Hardcore parameter estimators : ˆε = min{h min (C), C V or(γ) and C Λ n }, ˆα = max{h max (C), C V or(γ) and C Λ n }, ˆB = max{h 2 max(c)/vol(c), C V or(γ) and C Λ n }. -Smooth parameter estimators : (ẑ, ˆθ) = argmin z,θ P LL Λn (γ, z, θ, ˆε, ˆα, ˆB), with P LL Λn () = z exp ( h(x, γ)) dx+ Λ n where h(x, γ) = H Λn (γ + x) H Λn (γ). x γ Λn H Λn (γ x)< ( h(x, γ x) ln(z) ),
50 Theoretical results For the hardcore parameters : Theorem (Dereudre-L. (2009)) For P -almost all γ lim (ˆε, ˆα, ˆB) = (ε, α, B). n
51 Theoretical results For the hardcore parameters : Theorem (Dereudre-L. (2009)) For P -almost all γ For the smooth parameters : Theorem (Dereudre-L. (2009)) For P -almost all γ lim (ˆε, ˆα, ˆB) = (ε, α, B). n lim (ẑ, ˆθ) = (z, θ). n (ẑ, ˆθ) are asymptotic normal if ε, α and B are supposed to be known.
52 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5.
53 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5. Typical tessellation : Hardcore parameter estimators : ˆα ˆB
54 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5. Smooth parameter estimators : ˆθ when z is known ˆθ when z is estimated ẑ
55 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5.
56 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5. Typical tessellation : Hardcore parameter estimators : ˆα ˆB
57 Estimation results The true parameters : ε = 0, α = 0.05, B = 0.625, z = 100 and θ = 0.5. Smooth parameter estimators : ˆθ when z is known ˆθ when z is estimated ẑ
58 Conclusion Our Gibbs Voronoi model : forces the shape and the maximal size of the cells provides some repulsive or attractive interaction between two neighbour cells.
59 Conclusion Our Gibbs Voronoi model : forces the shape and the maximal size of the cells provides some repulsive or attractive interaction between two neighbour cells. The simulation can be achieved by a Birth-Death-Move MCMC algorithm very time consuming because of the hardcore interactions.
60 Conclusion Our Gibbs Voronoi model : forces the shape and the maximal size of the cells provides some repulsive or attractive interaction between two neighbour cells. The simulation can be achieved by a Birth-Death-Move MCMC algorithm very time consuming because of the hardcore interactions. A two-step estimation procedure can be applied 1 the hardcore parameters are estimated in a natural way, 2 the smooth parameters are estimated by pseudo-likelihood where the the hardcore parameters are plugged in.
61 Conclusion Our Gibbs Voronoi model : forces the shape and the maximal size of the cells provides some repulsive or attractive interaction between two neighbour cells. The simulation can be achieved by a Birth-Death-Move MCMC algorithm very time consuming because of the hardcore interactions. A two-step estimation procedure can be applied 1 the hardcore parameters are estimated in a natural way, 2 the smooth parameters are estimated by pseudo-likelihood where the the hardcore parameters are plugged in. This is consistent and allows to distinguish between the repulsive and the attractive case in a non-trivial situation.
62 A. Baddeley, R. Turner, J. Moller and M. Hazelton, Residual analysis for spatial point processes, J. R. Statist. Soc. B, 65, (2005). E. Bertin, J.M. Billiot, R. Drouilhet, (1999)Existence of nearest neighbors spatial Gibbs models, Adv. Appl. Prob. (SGSA) 31, J.-M. Billiot, J.-F. Coeurjolly and R. Drouilhet, (2008) Maximum pseudolikelihood estimator for exponential family models of marked Gibbs point processes, Electronic J. Statistics, 2, D. Dereudre, (2008) Gibbs Delaunay tessellations with geometric hardcore conditions, J. Stat. Phys., 131, D. Dereudre, R. Drouilhet and H.-O. Georgii, (2009) Existence of Gibbs process with stable cluster interactions, preprint. D. Dereudre, F. Lavancier, (2009) Campbell equilibrium equation and pseudo-likelihood estimation for non-hereditary Gibbs point processes, to appear in Bernoulli. D. Dereudre, F. Lavancier, (2009) Practical simulation and estimation for Gibbs Delaunay-Voronoi tessellations with geometric hardcore interaction, preprint. J.L. Jensen and J. Moller (1991) Pseudolikelihood for exponential family models of spatial point processes, Ann. Appl. Probab. 1, C. Preston, Random Fields, LNM Vol. 534, Springer.
63 5 Some theoretical points
64 The problem of heredity Definition The family of energies (H Λ ) Λ is said hereditary if for every Λ, every γ M(R 2 ) and every x Λ H Λ (γ) = + H Λ (γ + δ x ) = +.
65 The problem of heredity Definition The family of energies (H Λ ) Λ is said hereditary if for every Λ, every γ M(R 2 ) and every x Λ H Λ (γ) = + H Λ (γ + δ x ) = +. γ is forbidden γ + δ x is forbidden
66 The problem of heredity Definition The family of energies (H Λ ) Λ is said hereditary if for every Λ, every γ M(R 2 ) and every x Λ H Λ (γ) = + H Λ (γ + δ x ) = +. γ is forbidden γ + δ x is forbidden γ + δ x is allowed γ is allowed
67 The problem of heredity Definition The family of energies (H Λ ) Λ is said hereditary if for every Λ, every γ M(R 2 ) and every x Λ H Λ (γ) = + H Λ (γ + δ x ) = +. γ is forbidden γ + δ x is forbidden γ + δ x is allowed γ is allowed It is a standard assumption in classical statistical mechanics. (Example : The classical hard ball model is hereditary.)
68 The problem of heredity Definition The family of energies (H Λ ) Λ is said hereditary if for every Λ, every γ M(R 2 ) and every x Λ H Λ (γ) = + H Λ (γ + δ x ) = +. γ is forbidden γ + δ x is forbidden γ + δ x is allowed γ is allowed It is a standard assumption in classical statistical mechanics. (Example : The classical hard ball model is hereditary.) The Gibbs Voronoi Tessellations are not hereditary. When one adds a point in a too large cell, the new tessellation may be allowed.
69 Nguyen-Zessin equilibrium equation Theorem (Nguyen-Zessin (1979), hereditary case) Suppose that the energy (H Λ ) Λ is hereditary. P is Gibbs measure with intensity measure ν if and only if, for every bounded non negative measurable function ψ from R 2 M(R 2 ) to R, ( ) ( ) E P ψ(x, γ x) = E P ψ(x, γ)e h(x,γ) ν(dx), x γ R 2 where h(x, γ) = H Λn (γ + x) H Λn (γ). Proposition (Dereudre, L. (2009), general case) Let P be a Gibbs measure with intensity measure ν, then ( ) E P ψ(x, γ x) = E P ψ(x, γ)e h(x,γ) ν(dx). R 2 x γ Λn H Λn (γ x)<
70 Validation : residuals process We can extend the concept of residuals (see Baddeley et al., 2005) to the non-hereditary setting. The residuals process on a set is defined for any function ψ by ) R (, ψ, ĥ, ˆν = x γ H (γ x)< ψ(x, γ x) ψ(x, γ)e ĥ(x,γ)ˆν(dx), From the equilibrium equation given before, under the true model, ) R (, ψ, ĥ, ˆν 0 ) R (, ψ, ĥ, ˆν is approximatively gaussian. Several diagnostic tools can then be applied when fitting a Gibbs Voronoi model
71 Validation : residuals process We can extend the concept of residuals (see Baddeley et al., 2005) to the non-hereditary setting. The residuals process on a set is defined for any function ψ by ) R (, ψ, ĥ, ˆν = x γ H (γ x)< ψ(x, γ x) ψ(x, γ)e ĥ(x,γ)ˆν(dx), From the equilibrium equation given before, under the true model, ) R (, ψ, ĥ, ˆν 0 ) R (, ψ, ĥ, ˆν is approximatively gaussian. Several diagnostic tools can then be applied when fitting a Gibbs Voronoi model For further asymptotic results on the residuals process R : See the talk of J.-F. Coeurjolly on Friday morning.
Gibbs point processes : modelling and inference
Gibbs point processes : modelling and inference J.-F. Coeurjolly (Grenoble University) et J.-M Billiot, D. Dereudre, R. Drouilhet, F. Lavancier 03/09/2010 J.-F. Coeurjolly () Gibbs models 03/09/2010 1
More informationJean-Michel Billiot, Jean-François Coeurjolly and Rémy Drouilhet
Colloque International de Statistique Appliquée pour le Développement en Afrique International Conference on Applied Statistics for Development in Africa Sada 07 nn, 1?? 007) MAXIMUM PSEUDO-LIKELIHOOD
More informationStatistical study of spatial dependences in long memory random fields on a lattice, point processes and random geometry.
Statistical study of spatial dependences in long memory random fields on a lattice, point processes and random geometry. Frédéric Lavancier, Laboratoire de Mathématiques Jean Leray, Nantes 9 décembre 2011.
More informationGibbs point processes with geometrical interactions. Models for strongly structured patterns
Gibbs point processes with geometrical interactions. Models for strongly structured patterns David Dereudre, Laboratoire Paul Painlevé, Lille (joint work with R. Drouilhet and H.-O. Georgii) 8th September
More informationResiduals and Goodness-of-fit tests for marked Gibbs point processes
Residuals and Goodness-of-fit tests for marked Gibbs point processes Frédéric Lavancier (Laboratoire Jean Leray, Nantes, France) Joint work with J.-F. Coeurjolly (Grenoble, France) 09/06/2010 F. Lavancier
More informationTakacs Fiksel Method for Stationary Marked Gibbs Point Processes
Scandinavian Journal of Statistics doi: 10.1111/j.1467-9469.2011.00738.x Published by Blackwell Publishing Ltd. Takacs Fiksel Method for Stationary Marked Gibbs Point Processes JEAN-FRANÇOIS COEURJOLLY
More informationarxiv:math/ v1 [math.st] 4 Jan 2006
arxiv:math/060065v [math.st] 4 Jan 2006 Maximum pseudo-likelihood estimator for nearest-neighbours Gibbs point processes By Jean-Michel Billiot, Jean-François Coeurjolly,2 and Rémy Drouilhet Labsad, University
More informationarxiv: v2 [math.pr] 6 Apr 2018
Introduction to the theory of Gibbs point processes DEREUDRE David arxiv:70.0805v2 [math.pr] 6 Apr 208 Abstract The Gibbs point processes GPP constitute a large class of point processes with interaction
More informationPerfect simulation for repulsive point processes
Perfect simulation for repulsive point processes Why swapping at birth is a good thing Mark Huber Department of Mathematics Claremont-McKenna College 20 May, 2009 Mark Huber (Claremont-McKenna College)
More informationAn introduction to spatial point processes
An introduction to spatial point processes Jean-François Coeurjolly 1 Examples of spatial data 2 Intensities and Poisson p.p. 3 Summary statistics 4 Models for point processes A very very brief introduction...
More informationRandom graphs: Random geometric graphs
Random graphs: Random geometric graphs Mathew Penrose (University of Bath) Oberwolfach workshop Stochastic analysis for Poisson point processes February 2013 athew Penrose (Bath), Oberwolfach February
More informationBootstrap Approximation of Gibbs Measure for Finite-Range Potential in Image Analysis
Bootstrap Approximation of Gibbs Measure for Finite-Range Potential in Image Analysis Abdeslam EL MOUDDEN Business and Management School Ibn Tofaïl University Kenitra, Morocco Abstract This paper presents
More informationLikelihood Inference for Lattice Spatial Processes
Likelihood Inference for Lattice Spatial Processes Donghoh Kim November 30, 2004 Donghoh Kim 1/24 Go to 1234567891011121314151617 FULL Lattice Processes Model : The Ising Model (1925), The Potts Model
More informationEconometrics I, Estimation
Econometrics I, Estimation Department of Economics Stanford University September, 2008 Part I Parameter, Estimator, Estimate A parametric is a feature of the population. An estimator is a function of the
More informationChapter 2. Poisson point processes
Chapter 2. Poisson point processes Jean-François Coeurjolly http://www-ljk.imag.fr/membres/jean-francois.coeurjolly/ Laboratoire Jean Kuntzmann (LJK), Grenoble University Setting for this chapter To ease
More informationSpatial point processes: Theory and practice illustrated with R
Spatial point processes: Theory and practice illustrated with R Department of Mathematical Sciences Aalborg University Lecture IV, February 24, 2011 1/28 Contents of lecture IV Papangelou conditional intensity.
More informationPOWER DIAGRAMS AND INTERACTION PROCESSES FOR UNIONS OF DISCS
30 July 2007 POWER DIAGRAMS AND INTERACTION PROCESSES FOR UNIONS OF DISCS JESPER MØLLER, Aalborg University KATEŘINA HELISOVÁ, Charles University in Prague Abstract We study a flexible class of finite
More informationEstimating Unnormalised Models by Score Matching
Estimating Unnormalised Models by Score Matching Michael Gutmann Probabilistic Modelling and Reasoning (INFR11134) School of Informatics, University of Edinburgh Spring semester 2018 Program 1. Basics
More informationRESEARCH REPORT. A logistic regression estimating function for spatial Gibbs point processes.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 2013 Adrian Baddeley, Jean-François Coeurjolly, Ege Rubak and Rasmus Waagepetersen A logistic regression estimating function
More informationLikelihood inference in the presence of nuisance parameters
Likelihood inference in the presence of nuisance parameters Nancy Reid, University of Toronto www.utstat.utoronto.ca/reid/research 1. Notation, Fisher information, orthogonal parameters 2. Likelihood inference
More informationKernel families of probability measures. Saskatoon, October 21, 2011
Kernel families of probability measures Saskatoon, October 21, 2011 Abstract The talk will compare two families of probability measures: exponential, and Cauchy-Stjelties families. The exponential families
More informationi=1 h n (ˆθ n ) = 0. (2)
Stat 8112 Lecture Notes Unbiased Estimating Equations Charles J. Geyer April 29, 2012 1 Introduction In this handout we generalize the notion of maximum likelihood estimation to solution of unbiased estimating
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationHigh dimensional Ising model selection
High dimensional Ising model selection Pradeep Ravikumar UT Austin (based on work with John Lafferty, Martin Wainwright) Sparse Ising model US Senate 109th Congress Banerjee et al, 2008 Estimate a sparse
More informationR. Lachieze-Rey Recent Berry-Esseen bounds obtained with Stein s method andgeorgia PoincareTech. inequalities, 1 / with 29 G.
Recent Berry-Esseen bounds obtained with Stein s method and Poincare inequalities, with Geometric applications Raphaël Lachièze-Rey, Univ. Paris 5 René Descartes, Georgia Tech. R. Lachieze-Rey Recent Berry-Esseen
More informationMinimax estimators of the coverage probability of the impermissible error for a location family
Minimax estimators of the coverage probability of the impermissible error for a location family by Miguel A. Arcones Binghamton University arcones@math.binghamton.edu Talk based on: Arcones, M. A. (2008).
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationThe extremal elliptical model: Theoretical properties and statistical inference
1/25 The extremal elliptical model: Theoretical properties and statistical inference Thomas OPITZ Supervisors: Jean-Noel Bacro, Pierre Ribereau Institute of Mathematics and Modeling in Montpellier (I3M)
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationInvariant HPD credible sets and MAP estimators
Bayesian Analysis (007), Number 4, pp. 681 69 Invariant HPD credible sets and MAP estimators Pierre Druilhet and Jean-Michel Marin Abstract. MAP estimators and HPD credible sets are often criticized in
More informationLIMIT THEORY FOR PLANAR GILBERT TESSELLATIONS
PROBABILITY AND MATHEMATICAL STATISTICS Vol. 31, Fasc. 1 (2011), pp. 149 160 LIMIT THEORY FOR PLANAR GILBERT TESSELLATIONS BY TOMASZ S C H R E I B E R (TORUŃ) AND NATALIA S O JA (TORUŃ) Abstract. A Gilbert
More informationLecture 2: Poisson point processes: properties and statistical inference
Lecture 2: Poisson point processes: properties and statistical inference Jean-François Coeurjolly http://www-ljk.imag.fr/membres/jean-francois.coeurjolly/ 1 / 20 Definition, properties and simulation Statistical
More informationSemi-Parametric Importance Sampling for Rare-event probability Estimation
Semi-Parametric Importance Sampling for Rare-event probability Estimation Z. I. Botev and P. L Ecuyer IMACS Seminar 2011 Borovets, Bulgaria Semi-Parametric Importance Sampling for Rare-event probability
More informationA = {(x, u) : 0 u f(x)},
Draw x uniformly from the region {x : f(x) u }. Markov Chain Monte Carlo Lecture 5 Slice sampler: Suppose that one is interested in sampling from a density f(x), x X. Recall that sampling x f(x) is equivalent
More informationAdvances and Applications in Perfect Sampling
and Applications in Perfect Sampling Ph.D. Dissertation Defense Ulrike Schneider advisor: Jem Corcoran May 8, 2003 Department of Applied Mathematics University of Colorado Outline Introduction (1) MCMC
More informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
More informationOn Ergodic Impulse Control with Constraint
On Ergodic Impulse Control with Constraint Maurice Robin Based on joint papers with J.L. Menaldi University Paris-Sanclay 9119 Saint-Aubin, France (e-mail: maurice.robin@polytechnique.edu) IMA, Minneapolis,
More informationRandom Geometric Graphs
Random Geometric Graphs Mathew D. Penrose University of Bath, UK Networks: stochastic models for populations and epidemics ICMS, Edinburgh September 2011 1 MODELS of RANDOM GRAPHS Erdos-Renyi G(n, p):
More informationA PARAMETRIC MODEL FOR DISCRETE-VALUED TIME SERIES. 1. Introduction
tm Tatra Mt. Math. Publ. 00 (XXXX), 1 10 A PARAMETRIC MODEL FOR DISCRETE-VALUED TIME SERIES Martin Janžura and Lucie Fialová ABSTRACT. A parametric model for statistical analysis of Markov chains type
More informationLearning MN Parameters with Approximation. Sargur Srihari
Learning MN Parameters with Approximation Sargur srihari@cedar.buffalo.edu 1 Topics Iterative exact learning of MN parameters Difficulty with exact methods Approximate methods Approximate Inference Belief
More informationBayesian model selection in graphs by using BDgraph package
Bayesian model selection in graphs by using BDgraph package A. Mohammadi and E. Wit March 26, 2013 MOTIVATION Flow cytometry data with 11 proteins from Sachs et al. (2005) RESULT FOR CELL SIGNALING DATA
More informationStatistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation.
Statistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation. 1. Questions and code from last time. The difference between ETAS
More informationRandom Geometric Graphs
Random Geometric Graphs Mathew D. Penrose University of Bath, UK SAMBa Symposium (Talk 2) University of Bath November 2014 1 MODELS of RANDOM GRAPHS Erdos-Renyi G(n, p): start with the complete n-graph,
More informationarxiv: v1 [stat.me] 27 May 2015
Modelling aggregation on the large scale and regularity on the small scale in spatial point pattern datasets arxiv:1505.07215v1 [stat.me] 27 May 2015 Frédéric Lavancier 1,2 and Jesper Møller 3 1 Laboratoire
More informationTwo step estimation for Neyman-Scott point process with inhomogeneous cluster centers. May 2012
Two step estimation for Neyman-Scott point process with inhomogeneous cluster centers Tomáš Mrkvička, Milan Muška, Jan Kubečka May 2012 Motivation Study of the influence of covariates on the occurrence
More informationPoisson line processes. C. Lantuéjoul MinesParisTech
Poisson line processes C. Lantuéjoul MinesParisTech christian.lantuejoul@mines-paristech.fr Bertrand paradox A problem of geometrical probability A line is thrown at random on a circle. What is the probability
More informationWhere Probability Meets Combinatorics and Statistical Mechanics
Where Probability Meets Combinatorics and Statistical Mechanics Daniel Ueltschi Department of Mathematics, University of Warwick MASDOC, 15 March 2011 Collaboration with V. Betz, N.M. Ercolani, C. Goldschmidt,
More informationNon-parametric inference for stationary pairwise interaction point process
Non-parametric inference for stationary pairwise interaction point process Nadia Morsli To cite this version: Nadia Morsli. Non-parametric inference for stationary pairwise interaction point process. 2016.
More informationOn the entropy flows to disorder
CHAOS 2009 Charnia, Crete 1-5 June 2009 On the entropy flows to disorder C.T.J. Dodson School of Mathematics University of Manchester, UK Abstract Gamma distributions, which contain the exponential as
More informationRandom hyperplane tessellations and dimension reduction
Random hyperplane tessellations and dimension reduction Roman Vershynin University of Michigan, Department of Mathematics Phenomena in high dimensions in geometric analysis, random matrices and computational
More informationLecture 4 Lebesgue spaces and inequalities
Lecture 4: Lebesgue spaces and inequalities 1 of 10 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 4 Lebesgue spaces and inequalities Lebesgue spaces We have seen how
More informationGenerative Models and Stochastic Algorithms for Population Average Estimation and Image Analysis
Generative Models and Stochastic Algorithms for Population Average Estimation and Image Analysis Stéphanie Allassonnière CIS, JHU July, 15th 28 Context : Computational Anatomy Context and motivations :
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationLog-Density Estimation with Application to Approximate Likelihood Inference
Log-Density Estimation with Application to Approximate Likelihood Inference Martin Hazelton 1 Institute of Fundamental Sciences Massey University 19 November 2015 1 Email: m.hazelton@massey.ac.nz WWPMS,
More informationAsymptotic standard errors of MLE
Asymptotic standard errors of MLE Suppose, in the previous example of Carbon and Nitrogen in soil data, that we get the parameter estimates For maximum likelihood estimation, we can use Hessian matrix
More information12 - Nonparametric Density Estimation
ST 697 Fall 2017 1/49 12 - Nonparametric Density Estimation ST 697 Fall 2017 University of Alabama Density Review ST 697 Fall 2017 2/49 Continuous Random Variables ST 697 Fall 2017 3/49 1.0 0.8 F(x) 0.6
More informationOn probability... Mark Huber (Duke University) Swap moves for spatial point processes Graduate/Faculty Seminar 1 / 42
On probability... I think you re begging the question, said Haydock, and I can see looming ahead one of those terrible exercises in probability where six men have white hats and six men have black hats
More informationAN EFFICIENT ESTIMATOR FOR GIBBS RANDOM FIELDS
K Y B E R N E T I K A V O L U M E 5 0 ( 2 0 1 4, N U M B E R 6, P A G E S 8 8 3 8 9 5 AN EFFICIENT ESTIMATOR FOR GIBBS RANDOM FIELDS Martin Janžura An efficient estimator for the expectation R f dp is
More informationEstimation theory and information geometry based on denoising
Estimation theory and information geometry based on denoising Aapo Hyvärinen Dept of Computer Science & HIIT Dept of Mathematics and Statistics University of Helsinki Finland 1 Abstract What is the best
More informationStatistical Estimation: Data & Non-data Information
Statistical Estimation: Data & Non-data Information Roger J-B Wets University of California, Davis & M.Casey @ Raytheon G.Pflug @ U. Vienna, X. Dong @ EpiRisk, G-M You @ EpiRisk. a little background Decision
More informationFrom independent component analysis to score matching
From independent component analysis to score matching Aapo Hyvärinen Dept of Computer Science & HIIT Dept of Mathematics and Statistics University of Helsinki Finland 1 Abstract First, short introduction
More informationLecture 28: Asymptotic confidence sets
Lecture 28: Asymptotic confidence sets 1 α asymptotic confidence sets Similar to testing hypotheses, in many situations it is difficult to find a confidence set with a given confidence coefficient or level
More informationSAMPLE PATH PROPERIES OF RANDOM TRANSFORMATIONS.
SAMPLE PATH PROPERIES OF RANDOM TRANSFORMATIONS. DMITRY DOLGOPYAT 1. Models. Let M be a smooth compact manifold of dimension N and X 0, X 1... X d, d 2, be smooth vectorfields on M. (I) Let {w j } + j=
More informationLarge Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials
Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials Maxim L. Yattselev joint work with Christopher D. Sinclair International Conference on Approximation
More informationAsymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ Lawrence D. Brown University
More informationarxiv:math/ v2 [math.pr] 1 Oct 2008
The nnals of Probability 2008, Vol. 36, No. 4, 209 220 DOI: 0.24/07-OP352 c Institute of Mathematical Statistics, 2008 arxiv:math/06369v2 [math.pr] Oct 2008 SUBCRITICL REGIMES IN THE POISSON BOOLEN MODEL
More informationOn thinning a spatial point process into a Poisson process using the Papangelou intensity
On thinning a spatial point process into a Poisson process using the Papangelou intensity Frederic Paik choenberg Department of tatistics, University of California, Los Angeles, CA 90095-1554, UA and Jiancang
More informationarxiv: v2 [math.pr] 26 Jun 2017
Existence of an unbounded vacant set for subcritical continuum percolation arxiv:1706.03053v2 math.pr 26 Jun 2017 Daniel Ahlberg, Vincent Tassion and Augusto Teixeira Abstract We consider the Poisson Boolean
More informationOikos. Appendix 1 and 2. o20751
Oikos o20751 Rosindell, J. and Cornell, S. J. 2013. Universal scaling of species-abundance distributions across multiple scales. Oikos 122: 1101 1111. Appendix 1 and 2 Universal scaling of species-abundance
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationis a Borel subset of S Θ for each c R (Bertsekas and Shreve, 1978, Proposition 7.36) This always holds in practical applications.
Stat 811 Lecture Notes The Wald Consistency Theorem Charles J. Geyer April 9, 01 1 Analyticity Assumptions Let { f θ : θ Θ } be a family of subprobability densities 1 with respect to a measure µ on a measurable
More informationA Bayesian perspective on GMM and IV
A Bayesian perspective on GMM and IV Christopher A. Sims Princeton University sims@princeton.edu November 26, 2013 What is a Bayesian perspective? A Bayesian perspective on scientific reporting views all
More informationIntroduction to graphical models: Lecture III
Introduction to graphical models: Lecture III Martin Wainwright UC Berkeley Departments of Statistics, and EECS Martin Wainwright (UC Berkeley) Some introductory lectures January 2013 1 / 25 Introduction
More informationON COMPOUND POISSON POPULATION MODELS
ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre
More informationLogistic regression for spatial Gibbs point processes
Biometrika Advance Access published March 6, 2014 Biometrika (2014), pp. 1 16 Printed in Great Britain doi: 10.1093/biomet/ast060 Logistic regression for spatial Gibbs point processes BY ADRIAN BADDELEY
More information1 Mixed effect models and longitudinal data analysis
1 Mixed effect models and longitudinal data analysis Mixed effects models provide a flexible approach to any situation where data have a grouping structure which introduces some kind of correlation between
More informationModelling Non-linear and Non-stationary Time Series
Modelling Non-linear and Non-stationary Time Series Chapter 2: Non-parametric methods Henrik Madsen Advanced Time Series Analysis September 206 Henrik Madsen (02427 Adv. TS Analysis) Lecture Notes September
More informationConvergence rates in weighted L 1 spaces of kernel density estimators for linear processes
Alea 4, 117 129 (2008) Convergence rates in weighted L 1 spaces of kernel density estimators for linear processes Anton Schick and Wolfgang Wefelmeyer Anton Schick, Department of Mathematical Sciences,
More informationStochastic Proximal Gradient Algorithm
Stochastic Institut Mines-Télécom / Telecom ParisTech / Laboratoire Traitement et Communication de l Information Joint work with: Y. Atchade, Ann Arbor, USA, G. Fort LTCI/Télécom Paristech and the kind
More information3 : Representation of Undirected GM
10-708: Probabilistic Graphical Models 10-708, Spring 2016 3 : Representation of Undirected GM Lecturer: Eric P. Xing Scribes: Longqi Cai, Man-Chia Chang 1 MRF vs BN There are two types of graphical models:
More informationGeneralized Linear Models
Generalized Linear Models David Rosenberg New York University April 12, 2015 David Rosenberg (New York University) DS-GA 1003 April 12, 2015 1 / 20 Conditional Gaussian Regression Gaussian Regression Input
More informationApproximate Bayesian Computation and Model Assessment for Repulsive Spatial Point Processes
Approximate Bayesian Computation and Model Assessment for Repulsive Spatial Point Processes arxiv:1604.07027v2 [stat.co] 26 Aug 2016 Shinichiro Shirota Department of Statistical Science, Duke University,
More informationBayesian Modeling of Inhomogeneous Marked Spatial Point Patterns with Location Dependent Mark Distributions
Inhomogeneous with Location Dependent Mark Distributions Matt Bognar Department of Statistics and Actuarial Science University of Iowa December 4, 2014 A spatial point pattern describes the spatial location
More informationMachine Learning Basics Lecture 2: Linear Classification. Princeton University COS 495 Instructor: Yingyu Liang
Machine Learning Basics Lecture 2: Linear Classification Princeton University COS 495 Instructor: Yingyu Liang Review: machine learning basics Math formulation Given training data x i, y i : 1 i n i.i.d.
More informationGeneralised AR and MA Models and Applications
Chapter 3 Generalised AR and MA Models and Applications 3.1 Generalised Autoregressive Processes Consider an AR1) process given by 1 αb)x t = Z t ; α < 1. In this case, the acf is, ρ k = α k for k 0 and
More informationCauchy-Stieltjes kernel families
Cauchy-Stieltjes kernel families Wlodek Bryc Fields Institute, July 23, 2013 NEF versus CSK families The talk will switch between two examples of kernel families K(µ) = {P θ (dx) : θ Θ} NEF versus CSK
More informationChapters 9. Properties of Point Estimators
Chapters 9. Properties of Point Estimators Recap Target parameter, or population parameter θ. Population distribution f(x; θ). { probability function, discrete case f(x; θ) = density, continuous case The
More informationThis paper is not to be removed from the Examination Halls
~~ST104B ZA d0 This paper is not to be removed from the Examination Halls UNIVERSITY OF LONDON ST104B ZB BSc degrees and Diplomas for Graduates in Economics, Management, Finance and the Social Sciences,
More informationarxiv: v1 [math.pr] 17 Jan 2019
Perfect Sampling for Gibbs Point Processes Using Partial Rejection Sampling Sarat B. Moka and Dirk P. Kroese arxiv:90.05624v [math.pr] 7 Jan 209 School of Mathematics and Physics The University of Queensland,
More informationA tailor made nonparametric density estimate
A tailor made nonparametric density estimate Daniel Carando 1, Ricardo Fraiman 2 and Pablo Groisman 1 1 Universidad de Buenos Aires 2 Universidad de San Andrés School and Workshop on Probability Theory
More informationOn the entropy flows to disorder
On the entropy flows to disorder C.T.J. Dodson School of Mathematics, University of Manchester, Manchester M 9PL, UK ctdodson@manchester.ac.uk May, 9 Abstract Gamma distributions, which contain the exponential
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationChapter 20. Deep Generative Models
Peng et al.: Deep Learning and Practice 1 Chapter 20 Deep Generative Models Peng et al.: Deep Learning and Practice 2 Generative Models Models that are able to Provide an estimate of the probability distribution
More informationNeighbourhoods of Randomness and Independence
Neighbourhoods of Randomness and Independence C.T.J. Dodson School of Mathematics, Manchester University Augment information geometric measures in spaces of distributions, via explicit geometric representations
More informationBayesian Inference by Density Ratio Estimation
Bayesian Inference by Density Ratio Estimation Michael Gutmann https://sites.google.com/site/michaelgutmann Institute for Adaptive and Neural Computation School of Informatics, University of Edinburgh
More informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
More informationSCHWARTZ SPACES ASSOCIATED WITH SOME NON-DIFFERENTIAL CONVOLUTION OPERATORS ON HOMOGENEOUS GROUPS
C O L L O Q U I U M M A T H E M A T I C U M VOL. LXIII 1992 FASC. 2 SCHWARTZ SPACES ASSOCIATED WITH SOME NON-DIFFERENTIAL CONVOLUTION OPERATORS ON HOMOGENEOUS GROUPS BY JACEK D Z I U B A Ń S K I (WROC
More information