The Generalized Random Variable Appears the Trace of Fractional Calculus in Statistics

Total Page:16

File Type:pdf, Size:1020Kb

The Generalized Random Variable Appears the Trace of Fractional Calculus in Statistics Appl. Math. Inf. Sci. Lett. 3, No. 2, 61-67 (2015) 61 Applied Mathematics & Information Sciences Letters An International Journal http://dx.doi.org/10.12785/amisl/030204 The Generalized Random Variable Appears the Trace of Fractional Calculus in Statistics Masoud Ganji1,∗ and Fatemeh Gharari2 1 Department of Statistics, University of Mohaghegh Ardabili, Ardabil, Iran 2 Department of Mathematics, University of Mohaghegh Ardabili, Ardabil, Iran Received: 5 Aug. 2014, Revised: 14 Mar. 2015, Accepted: 17 Mar. 2015 Published online: 1 May 2015 Abstract: In this paper, introducing a generalized random variable, we obtained a new relationship between fractional calculus and statistics. With this definition, the domain of shape parameter space for some of distributions, such as Gamma, Beta and Weibull, were expanded from (0,∞) to (−1,∞). It is shown that the expectation of such random variable ,Φα (x), coincides with the fractional integral of probability density function (PDF), at the origin, for α > 0, and also the fractional derivative of PDF, at the origin, for −1 < α < 0. We also, presented PDFs of such distributions as a product of fractional derivation of Dirac delta function of shape parameter order, α δ ( )(.). Finally, we showed that the Liouville fractional differintegral operator on the moment generating function (MGF) of positive random variable, at the origin, gives fractional moments. Keywords: Fractional calculus, Dirac delta function, Fractional moments, Weibull distribution, Gamma distribution, Beta distribution 1 Introduction Dirac delta function which is defined as: δ C The study of generalized functions is now widely used in : X → applied mathematics and engineering sciences. δ(x)ϕ(x)dx = ϕ(0) Generalized functions are defined as a linear functional Z on a space X of conveniently chosen test functions. ϕ For every locally integrable function f L 1 R , there for every test function ∈ X. The test function space X is ∈ loc( ) ∞ R exists a distribution F : X → C defined by: usually chosen as a subspace of C ( ), the space of f infinitely differentiable functions [6]. ∞ In the present work, we suppose that X is a positive Ff (ϕ)= h f ,ϕi = f (x)ϕ(x)dx (1) random variable and α is the shape parameter of ∞ Z− distribution. The new random variable is represented as α X −1 where ϕ ∈ X is test function from a suitable space X of Φα Φα + the function (x) and defined by (x) = Γ (α) . The test functions. A distribution that corresponds to functions function Φα (x) can be extended to all complex values of via equation (1) are called regular distributions. Examples α α as a pseudo function and is a distribution whose for regular distributions are the convolution kernels K± ∈ support is [0,∞) except for the case α = 0,−1,.... Since L 1 R loc( ) defined as: distributions or generalized functions in mathematical analysis are not really functions in the classical sense, we also call our proposed random variable a generalized α−1 α ±x K = H(±x) (2) random variable. The expectation of this generalized ± Γ (α) random variable coincides with Riemann-Liouville left fractional integral of the PDF, at the origin, for α > 0 and for α > 0 and H(x) is a Heaviside unit step function. Marchaud fractional derivative of the PDF, at the origin, Distributions that are not regular are sometimes called for −1 < α < 0. The generalized random variable singular. An example for a singular distribution is the appears in some distributions that belong to the ∗ Corresponding author e-mail: [email protected] c 2015 NSP Natural Sciences Publishing Cor. 62 M. Ganji, F. Gharari: The Generalized Random Variable Appears the Trace... exponential family like Weibull, Gamma, Beta, Binomial where α > −1 and β > −1 such that α + β > −1 and we and poisson distributions. By this generalized random used the star notation for the convolution operation. variable, we more expand the domain of shape parameter The proof is easy for α > 0 and β > 0, as instance, see space. For example, the domain of shape parameter which [2]. Other values of α and β can be proved using analytic has been already expanded for Gamma and Weibull continuation. distributions of α > 0 to α > −1 and for Beta distribution from (0,∞) × (0,∞) to (−1,∞) × (0,∞). In case of (e) If X1,X2,...,Xn be i.i.d. W(x;α,θ) then negative values of the shape parameter space, the Φ Γ Γ α θ relationship between fractional derivatives and statistics α+1(x) ∼ (1, ( + 1) ), theory can be obtained, this means we can write the PDF of distributions like Gamma, Weibull, Beta as a product and of fractional derivatives of Dirac delta function. The Φα+1(x) characteristics of this generalized random variable are as ∼ B(1,n − 1). Σ nΦα (x ) following: 1 +1 i Also, we demonstrate that the αth moment, α ∈ R, of the (a) By rewriting Binomial and poisson distributions in positive random variable X can be obtained directly from terms of the generalized random variable, they are the Liouville fractional derivation and integrations of the respectively transformed into Gamma and Beta MGF at the origin. distributions with a discrete parameter space, that is 2 Preliminaries [P(x;λ);Φλ (x)] = [Γ (λ;x,1);Φx(λ)], x = 0,1,... and In this section, we introduce notations, definitions and preliminary facts which are used throughout this paper. [Bin(x; p);Φp(x)] = [B(p;x,n − x + 2);Φx(p)] We need some basic definitions and properties of the fractional calculus theory and the generalized functions where x = 1,2,...,n and [.,Φ] is the rewriting form of theory which are used further in this paper. As mentioned generalized random variable. in references [1] and [3], the definitions 1 and 2 of fractional calculus are as following: (b) The expectation of this generalized random Φ variable, α (x), coincides with Riemann-Liouville left Definition 1. For a function f defined on an interval [a,b], fractional integral of the PDF at the origin for α > 0 and α α the Riemann-Liouville (R-L) integrals Ia+ f and Ib− f of Marchaud fractional derivative of the PDF at the origin order α ∈ C , (R(α) > 0) are defined, respectively, by for −1 < α < 0, that is, we have t α 1 α (I f )(t)= (t − ξ ) −1 f (ξ )dξ (7) a+ Γ (α) α α Za Φ (I− f )(0), > 0 E[ α (X)] = α α (3) ((D− f )(0), −1 < < 0 and b α 1 α where (I f )(t)= (ξ −t) −1 f (ξ )dξ . (8) b− Γ (α) ∞ Zt α 1 α−1 (I− f )(x)= t f (x +t)dt (4) α Γ (α) 0 Also, the left and right R-L fractional derivations Da+ f Z α α C R α and Db− f of order ∈ , ( ( ) > 0) are defined, is the Riemann-Liouville left fractional integral, while respectively, by ∞ α 1 α α d α (D f )(x)= t− −1{ f (x +t) − f (x)}dt (5) (D f )(t) = ( )n(In− f )(t), (9) − Γ (−α) a+ dt a+ Z0 is the Marchaud fractional derivative. and α −d n n−α (Db− f )(t) = ( ) (Ia+ f )(t). (10) (c) The integer and fractional derivatives of this dt generalized random variable are the generalized random If in above definition, respectively, a = −∞ and b = ∞, variable, too. then we get Liouville fractional differintegral. (d) Another property is α Definition 2. The Liouville fractional differintegral D± is defined by Φ Φ Φ α (t − a) ∗ β (t)= α+β (t − a), (6) c 2015 NSP Natural Sciences Publishing Cor. Appl. Math. Inf. Sci. Lett. 3, No. 2, 61-67 (2015) / www.naturalspublishing.com/Journals.asp 63 for 0 < α < 1. [6] Γ α x −α−1 , ∞ R ′ R+ α frac1 (− ) −∞(x −t) f (t)dt Now if f ∈ C0 ( ) with supp f ⊂ , then D+{ f (x)} = dn α−n ( dxn {D+ { f (Rx)} (11) α α −α 0 δ (α) D0+ f = D0+(I f ) = (K+ · K+) · f = · f (20) R α where the first expression satisfies for ( ) < 0 and for all α ∈ C. Also, the differentiation rule the second expression satisfies for R(α) > 0;n = [R(α)] + 1. Also α β β −α D0+K+ = K+ (21) β α C 1 ∞ holds for all , ∈ . It contains (x −t)−α−1 f (t)dt, α Γ α x D f x (− ) (12) β β −1 −{ ( )} = n n DK+ = K+ (−1) d R α−n n {D− { f (x)} dx for all β ∈ C as a special case [6]. with R(α) <0 for the first expression and R(α) > 0 such that n = [R(α)] + 1, for the second expression satisfies. In particular, when α = n ∈ N0, then 3 The generalized Weibull random variable 0 n (n) The Wiebull distribution was originally used for modeling D+{ f (x)} = f (x),D+{ f (x)} = f (x), (13) fatigue data and is at present used extensively in and engineering problems for modeling, the distribution of the lifetime of an object which consists of several parts and 0 n n (n) D−{ f (x)} = f (x),D−{ f (x)} = (−1) f (x). (14) that fails if any component part fails [7]. It is important to note that the fractional calculus has also ∞ R ′ Definition 3. Let f be a generalized function f ∈ C0 ( ) been applied in describing several PDFs in mathematical with supp f ⊂ R+. Then its fractional integral is the statistics in terms of fractional integral and derivative of α distribution I0+ f defined as: exponential functions. For instance, the Wiebull distribution can be writtenα in terms of the fractional α ϕ α ϕ α ϕ −θx hI0+ f , i = hI f , i = hK+ ∗ f , i (15) derivative or integral of e , as α for R(α) > 0[6].
Recommended publications
  • The Exponential Family 1 Definition
    The Exponential Family David M. Blei Columbia University November 9, 2016 The exponential family is a class of densities (Brown, 1986). It encompasses many familiar forms of likelihoods, such as the Gaussian, Poisson, multinomial, and Bernoulli. It also encompasses their conjugate priors, such as the Gamma, Dirichlet, and beta. 1 Definition A probability density in the exponential family has this form p.x / h.x/ exp >t.x/ a./ ; (1) j D f g where is the natural parameter; t.x/ are sufficient statistics; h.x/ is the “base measure;” a./ is the log normalizer. Examples of exponential family distributions include Gaussian, gamma, Poisson, Bernoulli, multinomial, Markov models. Examples of distributions that are not in this family include student-t, mixtures, and hidden Markov models. (We are considering these families as distributions of data. The latent variables are implicitly marginalized out.) The statistic t.x/ is called sufficient because the probability as a function of only depends on x through t.x/. The exponential family has fundamental connections to the world of graphical models (Wainwright and Jordan, 2008). For our purposes, we’ll use exponential 1 families as components in directed graphical models, e.g., in the mixtures of Gaussians. The log normalizer ensures that the density integrates to 1, Z a./ log h.x/ exp >t.x/ d.x/ (2) D f g This is the negative logarithm of the normalizing constant. The function h.x/ can be a source of confusion. One way to interpret h.x/ is the (unnormalized) distribution of x when 0. It might involve statistics of x that D are not in t.x/, i.e., that do not vary with the natural parameter.
    [Show full text]
  • A Skew Extension of the T-Distribution, with Applications
    J. R. Statist. Soc. B (2003) 65, Part 1, pp. 159–174 A skew extension of the t-distribution, with applications M. C. Jones The Open University, Milton Keynes, UK and M. J. Faddy University of Birmingham, UK [Received March 2000. Final revision July 2002] Summary. A tractable skew t-distribution on the real line is proposed.This includes as a special case the symmetric t-distribution, and otherwise provides skew extensions thereof.The distribu- tion is potentially useful both for modelling data and in robustness studies. Properties of the new distribution are presented. Likelihood inference for the parameters of this skew t-distribution is developed. Application is made to two data modelling examples. Keywords: Beta distribution; Likelihood inference; Robustness; Skewness; Student’s t-distribution 1. Introduction Student’s t-distribution occurs frequently in statistics. Its usual derivation and use is as the sam- pling distribution of certain test statistics under normality, but increasingly the t-distribution is being used in both frequentist and Bayesian statistics as a heavy-tailed alternative to the nor- mal distribution when robustness to possible outliers is a concern. See Lange et al. (1989) and Gelman et al. (1995) and references therein. It will often be useful to consider a further alternative to the normal or t-distribution which is both heavy tailed and skew. To this end, we propose a family of distributions which includes the symmetric t-distributions as special cases, and also includes extensions of the t-distribution, still taking values on the whole real line, with non-zero skewness. Let a>0 and b>0be parameters.
    [Show full text]
  • On a Problem Connected with Beta and Gamma Distributions by R
    ON A PROBLEM CONNECTED WITH BETA AND GAMMA DISTRIBUTIONS BY R. G. LAHA(i) 1. Introduction. The random variable X is said to have a Gamma distribution G(x;0,a)if du for x > 0, (1.1) P(X = x) = G(x;0,a) = JoT(a)" 0 for x ^ 0, where 0 > 0, a > 0. Let X and Y be two independently and identically distributed random variables each having a Gamma distribution of the form (1.1). Then it is well known [1, pp. 243-244], that the random variable W = X¡iX + Y) has a Beta distribution Biw ; a, a) given by 0 for w = 0, (1.2) PiW^w) = Biw;x,x)=\ ) u"-1il-u)'-1du for0<w<l, Ío T(a)r(a) 1 for w > 1. Now we can state the converse problem as follows : Let X and Y be two independently and identically distributed random variables having a common distribution function Fix). Suppose that W = Xj{X + Y) has a Beta distribution of the form (1.2). Then the question is whether £(x) is necessarily a Gamma distribution of the form (1.1). This problem was posed by Mauldon in [9]. He also showed that the converse problem is not true in general and constructed an example of a non-Gamma distribution with this property using the solution of an integral equation which was studied by Goodspeed in [2]. In the present paper we carry out a systematic investigation of this problem. In §2, we derive some general properties possessed by this class of distribution laws Fix).
    [Show full text]
  • Lectures on the Local Semicircle Law for Wigner Matrices
    Lectures on the local semicircle law for Wigner matrices Florent Benaych-Georges∗ Antti Knowlesy September 11, 2018 These notes provide an introduction to the local semicircle law from random matrix theory, as well as some of its applications. We focus on Wigner matrices, Hermitian random matrices with independent upper-triangular entries with zero expectation and constant variance. We state and prove the local semicircle law, which says that the eigenvalue distribution of a Wigner matrix is close to Wigner's semicircle distribution, down to spectral scales containing slightly more than one eigenvalue. This local semicircle law is formulated using the Green function, whose individual entries are controlled by large deviation bounds. We then discuss three applications of the local semicircle law: first, complete delocalization of the eigenvectors, stating that with high probability the eigenvectors are approximately flat; second, rigidity of the eigenvalues, giving large deviation bounds on the locations of the individual eigenvalues; third, a comparison argument for the local eigenvalue statistics in the bulk spectrum, showing that the local eigenvalue statistics of two Wigner matrices coincide provided the first four moments of their entries coincide. We also sketch further applications to eigenvalues near the spectral edge, and to the distribution of eigenvectors. arXiv:1601.04055v4 [math.PR] 10 Sep 2018 ∗Universit´eParis Descartes, MAP5. Email: [email protected]. yETH Z¨urich, Departement Mathematik. Email: [email protected].
    [Show full text]
  • A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess
    S S symmetry Article A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess Guillermo Martínez-Flórez 1, Víctor Leiva 2,* , Emilio Gómez-Déniz 3 and Carolina Marchant 4 1 Departamento de Matemáticas y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Montería 14014, Colombia; [email protected] 2 Escuela de Ingeniería Industrial, Pontificia Universidad Católica de Valparaíso, 2362807 Valparaíso, Chile 3 Facultad de Economía, Empresa y Turismo, Universidad de Las Palmas de Gran Canaria and TIDES Institute, 35001 Canarias, Spain; [email protected] 4 Facultad de Ciencias Básicas, Universidad Católica del Maule, 3466706 Talca, Chile; [email protected] * Correspondence: [email protected] or [email protected] Received: 30 June 2020; Accepted: 19 August 2020; Published: 1 September 2020 Abstract: In this paper, we consider skew-normal distributions for constructing new a distribution which allows us to model proportions and rates with zero/one inflation as an alternative to the inflated beta distributions. The new distribution is a mixture between a Bernoulli distribution for explaining the zero/one excess and a censored skew-normal distribution for the continuous variable. The maximum likelihood method is used for parameter estimation. Observed and expected Fisher information matrices are derived to conduct likelihood-based inference in this new type skew-normal distribution. Given the flexibility of the new distributions, we are able to show, in real data scenarios, the good performance of our proposal. Keywords: beta distribution; centered skew-normal distribution; maximum-likelihood methods; Monte Carlo simulations; proportions; R software; rates; zero/one inflated data 1.
    [Show full text]
  • Lecture 2 — September 24 2.1 Recap 2.2 Exponential Families
    STATS 300A: Theory of Statistics Fall 2015 Lecture 2 | September 24 Lecturer: Lester Mackey Scribe: Stephen Bates and Andy Tsao 2.1 Recap Last time, we set out on a quest to develop optimal inference procedures and, along the way, encountered an important pair of assertions: not all data is relevant, and irrelevant data can only increase risk and hence impair performance. This led us to introduce a notion of lossless data compression (sufficiency): T is sufficient for P with X ∼ Pθ 2 P if X j T (X) is independent of θ. How far can we take this idea? At what point does compression impair performance? These are questions of optimal data reduction. While we will develop general answers to these questions in this lecture and the next, we can often say much more in the context of specific modeling choices. With this in mind, let's consider an especially important class of models known as the exponential family models. 2.2 Exponential Families Definition 1. The model fPθ : θ 2 Ωg forms an s-dimensional exponential family if each Pθ has density of the form: s ! X p(x; θ) = exp ηi(θ)Ti(x) − B(θ) h(x) i=1 • ηi(θ) 2 R are called the natural parameters. • Ti(x) 2 R are its sufficient statistics, which follows from NFFC. • B(θ) is the log-partition function because it is the logarithm of a normalization factor: s ! ! Z X B(θ) = log exp ηi(θ)Ti(x) h(x)dµ(x) 2 R i=1 • h(x) 2 R: base measure.
    [Show full text]
  • Semi-Parametric Likelihood Functions for Bivariate Survival Data
    Old Dominion University ODU Digital Commons Mathematics & Statistics Theses & Dissertations Mathematics & Statistics Summer 2010 Semi-Parametric Likelihood Functions for Bivariate Survival Data S. H. Sathish Indika Old Dominion University Follow this and additional works at: https://digitalcommons.odu.edu/mathstat_etds Part of the Mathematics Commons, and the Statistics and Probability Commons Recommended Citation Indika, S. H. S.. "Semi-Parametric Likelihood Functions for Bivariate Survival Data" (2010). Doctor of Philosophy (PhD), Dissertation, Mathematics & Statistics, Old Dominion University, DOI: 10.25777/ jgbf-4g75 https://digitalcommons.odu.edu/mathstat_etds/30 This Dissertation is brought to you for free and open access by the Mathematics & Statistics at ODU Digital Commons. It has been accepted for inclusion in Mathematics & Statistics Theses & Dissertations by an authorized administrator of ODU Digital Commons. For more information, please contact [email protected]. SEMI-PARAMETRIC LIKELIHOOD FUNCTIONS FOR BIVARIATE SURVIVAL DATA by S. H. Sathish Indika BS Mathematics, 1997, University of Colombo MS Mathematics-Statistics, 2002, New Mexico Institute of Mining and Technology MS Mathematical Sciences-Operations Research, 2004, Clemson University MS Computer Science, 2006, College of William and Mary A Dissertation Submitted to the Faculty of Old Dominion University in Partial Fulfillment of the Requirement for the Degree of DOCTOR OF PHILOSOPHY DEPARTMENT OF MATHEMATICS AND STATISTICS OLD DOMINION UNIVERSITY August 2010 Approved by: DayanandiN; Naik Larry D. Le ia M. Jones ABSTRACT SEMI-PARAMETRIC LIKELIHOOD FUNCTIONS FOR BIVARIATE SURVIVAL DATA S. H. Sathish Indika Old Dominion University, 2010 Director: Dr. Norou Diawara Because of the numerous applications, characterization of multivariate survival dis­ tributions is still a growing area of research.
    [Show full text]
  • Delta Functions and Distributions
    When functions have no value(s): Delta functions and distributions Steven G. Johnson, MIT course 18.303 notes Created October 2010, updated March 8, 2017. Abstract x = 0. That is, one would like the function δ(x) = 0 for all x 6= 0, but with R δ(x)dx = 1 for any in- These notes give a brief introduction to the mo- tegration region that includes x = 0; this concept tivations, concepts, and properties of distributions, is called a “Dirac delta function” or simply a “delta which generalize the notion of functions f(x) to al- function.” δ(x) is usually the simplest right-hand- low derivatives of discontinuities, “delta” functions, side for which to solve differential equations, yielding and other nice things. This generalization is in- a Green’s function. It is also the simplest way to creasingly important the more you work with linear consider physical effects that are concentrated within PDEs, as we do in 18.303. For example, Green’s func- very small volumes or times, for which you don’t ac- tions are extremely cumbersome if one does not al- tually want to worry about the microscopic details low delta functions. Moreover, solving PDEs with in this volume—for example, think of the concepts of functions that are not classically differentiable is of a “point charge,” a “point mass,” a force plucking a great practical importance (e.g. a plucked string with string at “one point,” a “kick” that “suddenly” imparts a triangle shape is not twice differentiable, making some momentum to an object, and so on.
    [Show full text]
  • Distributions (3) © 2008 Winton 2 VIII
    © 2008 Winton 1 Distributions (3) © onntiW8 020 2 VIII. Lognormal Distribution • Data points t are said to be lognormally distributed if the natural logarithms, ln(t), of these points are normally distributed with mean μ ion and standard deviatσ – If the normal distribution is sampled to get points rsample, then the points ersample constitute sample values from the lognormal distribution • The pdf or the lognormal distribution is given byf - ln(x)()μ 2 - 1 2 f(x)= e ⋅ 2 σ >(x 0) σx 2 π because ∞ - ln(x)()μ 2 - ∞ -() t μ-2 1 1 2 1 2 ∫ e⋅ 2eσ dx = ∫ dt ⋅ 2σ (where= t ln(x)) 0 xσ2 π 0σ2 π is the pdf for the normal distribution © 2008 Winton 3 Mean and Variance for Lognormal • It can be shown that in terms of μ and σ σ2 μ + E(X)= 2 e 2⋅μ σ + 2 σ2 Var(X)= e( ⋅ ) e - 1 • If the mean E and variance V for the lognormal distribution are given, then the corresponding μ and σ2 for the normal distribution are given by – σ2 = log(1+V/E2) – μ = log(E) - σ2/2 • The lognormal distribution has been used in reliability models for time until failure and for stock price distributions – The shape is similar to that of the Gamma distribution and the Weibull distribution for the case α > 2, but the peak is less towards 0 © 2008 Winton 4 Graph of Lognormal pdf f(x) 0.5 0.45 μ=4, σ=1 0.4 0.35 μ=4, σ=2 0.3 μ=4, σ=3 0.25 μ=4, σ=4 0.2 0.15 0.1 0.05 0 x 0510 μ and σ are the lognormal’s mean and std dev, not those of the associated normal © onntiW8 020 5 IX.
    [Show full text]
  • MTHE/MATH 332 Introduction to Control
    1 Queen’s University Mathematics and Engineering and Mathematics and Statistics MTHE/MATH 332 Introduction to Control (Preliminary) Supplemental Lecture Notes Serdar Yuksel¨ April 3, 2021 2 Contents 1 Introduction ....................................................................................1 1.1 Introduction................................................................................1 1.2 Linearization...............................................................................4 2 Systems ........................................................................................7 2.1 System Properties...........................................................................7 2.2 Linear Systems.............................................................................8 2.2.1 Representation of Discrete-Time Signals in terms of Unit Pulses..............................8 2.2.2 Linear Systems.......................................................................8 2.3 Linear and Time-Invariant (Convolution) Systems................................................9 2.4 Bounded-Input-Bounded-Output (BIBO) Stability of Convolution Systems........................... 11 2.5 The Frequency Response (or Transfer) Function of Linear Time-Invariant Systems..................... 11 2.6 Steady-State vs. Transient Solutions............................................................ 11 2.7 Bode Plots................................................................................. 12 2.8 Interconnections of Systems..................................................................
    [Show full text]
  • Computational Methods in Nonlinear Physics 1. Review of Probability and Statistics
    Computational Methods in Nonlinear Physics 1. Review of probability and statistics P.I. Hurtado Departamento de Electromagnetismo y F´ısicade la Materia, and Instituto Carlos I de F´ısica Te´oricay Computacional. Universidad de Granada. E-18071 Granada. Spain E-mail: [email protected] Abstract. These notes correspond to the first part (20 hours) of a course on Computational Methods in Nonlinear Physics within the Master on Physics and Mathematics (FisyMat) of University of Granada. In this first chapter we review some concepts of probability theory and stochastic processes. Keywords: Computational physics, probability and statistics, Monte Carlo methods, stochastic differential equations, Langevin equation, Fokker-Planck equation, molecular dynamics. References and sources [1] R. Toral and P. Collet, Stochastic Numerical Methods, Wiley (2014). [2] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford Univ. Press (2002). [3] N.G. van Kampen, Stochastic Processes in Physics and Chemistry, Elsevier (2007). [4] Wikipedia: https://en.wikipedia.org CONTENTS 2 Contents 1 The probability space4 1.1 The σ-algebra of events....................................5 1.2 Probability measures and Kolmogorov axioms........................6 1.3 Conditional probabilities and statistical independence....................8 2 Random variables9 2.1 Definition of random variables................................. 10 3 Average values, moments and characteristic function 15 4 Some Important Probability Distributions 18 4.1 Bernuilli distribution.....................................
    [Show full text]
  • Theoretical Properties of the Weighted Feller-Pareto Distributions." Asian Journal of Mathematics and Applications, 2014: 1-12
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Georgia Southern University: Digital Commons@Georgia Southern Georgia Southern University Digital Commons@Georgia Southern Mathematical Sciences Faculty Publications Mathematical Sciences, Department of 2014 Theoretical Properties of the Weighted Feller- Pareto Distributions Oluseyi Odubote Georgia Southern University Broderick O. Oluyede Georgia Southern University, [email protected] Follow this and additional works at: https://digitalcommons.georgiasouthern.edu/math-sci-facpubs Part of the Mathematics Commons Recommended Citation Odubote, Oluseyi, Broderick O. Oluyede. 2014. "Theoretical Properties of the Weighted Feller-Pareto Distributions." Asian Journal of Mathematics and Applications, 2014: 1-12. source: http://scienceasia.asia/index.php/ama/article/view/173/ https://digitalcommons.georgiasouthern.edu/math-sci-facpubs/307 This article is brought to you for free and open access by the Mathematical Sciences, Department of at Digital Commons@Georgia Southern. It has been accepted for inclusion in Mathematical Sciences Faculty Publications by an authorized administrator of Digital Commons@Georgia Southern. For more information, please contact [email protected]. ASIAN JOURNAL OF MATHEMATICS AND APPLICATIONS Volume 2014, Article ID ama0173, 12 pages ISSN 2307-7743 http://scienceasia.asia THEORETICAL PROPERTIES OF THE WEIGHTED FELLER-PARETO AND RELATED DISTRIBUTIONS OLUSEYI ODUBOTE AND BRODERICK O. OLUYEDE Abstract. In this paper, for the first time, a new six-parameter class of distributions called weighted Feller-Pareto (WFP) and related family of distributions is proposed. This new class of distributions contains several other Pareto-type distributions such as length-biased (LB) Pareto, weighted Pareto (WP I, II, III, and IV), and Pareto (P I, II, III, and IV) distributions as special cases.
    [Show full text]