\Only if part" Suppose U is a normal random vector. Assuming that we can compute \(F^{-1}\), the previous exercise shows how we can simulate a distribution with distribution function \(F\). Work on the task that is enjoyable to you. If x_mean is the mean of my first normal distribution, then can the new mean be calculated as : k_mean = x . 2. Similarly, \(V\) is the lifetime of the parallel system which operates if and only if at least one component is operating. The result in the previous exercise is very important in the theory of continuous-time Markov chains. Let \(Y = a + b \, X\) where \(a \in \R\) and \(b \in \R \setminus\{0\}\). Subsection 3.3.3 The Matrix of a Linear Transformation permalink. Random variable \( V = X Y \) has probability density function \[ v \mapsto \int_{-\infty}^\infty f(x, v / x) \frac{1}{|x|} dx \], Random variable \( W = Y / X \) has probability density function \[ w \mapsto \int_{-\infty}^\infty f(x, w x) |x| dx \], We have the transformation \( u = x \), \( v = x y\) and so the inverse transformation is \( x = u \), \( y = v / u\). Find the probability density function of each of the follow: Suppose that \(X\), \(Y\), and \(Z\) are independent, and that each has the standard uniform distribution. Both results follows from the previous result above since \( f(x, y) = g(x) h(y) \) is the probability density function of \( (X, Y) \). Suppose that \(T\) has the exponential distribution with rate parameter \(r \in (0, \infty)\). Moreover, this type of transformation leads to simple applications of the change of variable theorems. \( g(y) = \frac{3}{25} \left(\frac{y}{100}\right)\left(1 - \frac{y}{100}\right)^2 \) for \( 0 \le y \le 100 \). The change of temperature measurement from Fahrenheit to Celsius is a location and scale transformation. In this section, we consider the bivariate normal distribution first, because explicit results can be given and because graphical interpretations are possible. Then \( X + Y \) is the number of points in \( A \cup B \). Linear Transformation of Gaussian Random Variable Theorem Let , and be real numbers . Set \(k = 1\) (this gives the minimum \(U\)). The Jacobian is the infinitesimal scale factor that describes how \(n\)-dimensional volume changes under the transformation. Let \(U = X + Y\), \(V = X - Y\), \( W = X Y \), \( Z = Y / X \). The Erlang distribution is studied in more detail in the chapter on the Poisson Process, and in greater generality, the gamma distribution is studied in the chapter on Special Distributions. Suppose that \(\bs X\) is a random variable taking values in \(S \subseteq \R^n\), and that \(\bs X\) has a continuous distribution with probability density function \(f\). Transform a normal distribution to linear. (z - x)!} Then \( (R, \Theta, Z) \) has probability density function \( g \) given by \[ g(r, \theta, z) = f(r \cos \theta , r \sin \theta , z) r, \quad (r, \theta, z) \in [0, \infty) \times [0, 2 \pi) \times \R \], Finally, for \( (x, y, z) \in \R^3 \), let \( (r, \theta, \phi) \) denote the standard spherical coordinates corresponding to the Cartesian coordinates \((x, y, z)\), so that \( r \in [0, \infty) \) is the radial distance, \( \theta \in [0, 2 \pi) \) is the azimuth angle, and \( \phi \in [0, \pi] \) is the polar angle. Then the lifetime of the system is also exponentially distributed, and the failure rate of the system is the sum of the component failure rates. When V and W are finite dimensional, a general linear transformation can Algebra Examples. It is also interesting when a parametric family is closed or invariant under some transformation on the variables in the family. Bryan 3 years ago The Cauchy distribution is studied in detail in the chapter on Special Distributions. Hence \[ \frac{\partial(x, y)}{\partial(u, v)} = \left[\begin{matrix} 1 & 0 \\ -v/u^2 & 1/u\end{matrix} \right] \] and so the Jacobian is \( 1/u \). If \( a, \, b \in (0, \infty) \) then \(f_a * f_b = f_{a+b}\). Suppose that \(X\) and \(Y\) are random variables on a probability space, taking values in \( R \subseteq \R\) and \( S \subseteq \R \), respectively, so that \( (X, Y) \) takes values in a subset of \( R \times S \). The independence of \( X \) and \( Y \) corresponds to the regions \( A \) and \( B \) being disjoint. \( f \) is concave upward, then downward, then upward again, with inflection points at \( x = \mu \pm \sigma \). In the reliability setting, where the random variables are nonnegative, the last statement means that the product of \(n\) reliability functions is another reliability function. The following result gives some simple properties of convolution. Let X N ( , 2) where N ( , 2) is the Gaussian distribution with parameters and 2 . The formulas in last theorem are particularly nice when the random variables are identically distributed, in addition to being independent. These can be combined succinctly with the formula \( f(x) = p^x (1 - p)^{1 - x} \) for \( x \in \{0, 1\} \). Our team is available 24/7 to help you with whatever you need. = f_{a+b}(z) \end{align}. \(Y\) has probability density function \( g \) given by \[ g(y) = \frac{1}{\left|b\right|} f\left(\frac{y - a}{b}\right), \quad y \in T \]. (1) (1) x N ( , ). Recall that for \( n \in \N_+ \), the standard measure of the size of a set \( A \subseteq \R^n \) is \[ \lambda_n(A) = \int_A 1 \, dx \] In particular, \( \lambda_1(A) \) is the length of \(A\) for \( A \subseteq \R \), \( \lambda_2(A) \) is the area of \(A\) for \( A \subseteq \R^2 \), and \( \lambda_3(A) \) is the volume of \(A\) for \( A \subseteq \R^3 \). A fair die is one in which the faces are equally likely. Thus we can simulate the polar radius \( R \) with a random number \( U \) by \( R = \sqrt{-2 \ln(1 - U)} \), or a bit more simply by \(R = \sqrt{-2 \ln U}\), since \(1 - U\) is also a random number. Hence the PDF of \( V \) is \[ v \mapsto \int_{-\infty}^\infty f(u, v / u) \frac{1}{|u|} du \], We have the transformation \( u = x \), \( w = y / x \) and so the inverse transformation is \( x = u \), \( y = u w \). Show how to simulate, with a random number, the exponential distribution with rate parameter \(r\). Suppose that \(X\) has a continuous distribution on an interval \(S \subseteq \R\) Then \(U = F(X)\) has the standard uniform distribution. In many cases, the probability density function of \(Y\) can be found by first finding the distribution function of \(Y\) (using basic rules of probability) and then computing the appropriate derivatives of the distribution function. (2) (2) y = A x + b N ( A + b, A A T). We also acknowledge previous National Science Foundation support under grant numbers 1246120, 1525057, and 1413739. Then \(Y\) has a discrete distribution with probability density function \(g\) given by \[ g(y) = \int_{r^{-1}\{y\}} f(x) \, dx, \quad y \in T \]. Note that \( \P\left[\sgn(X) = 1\right] = \P(X \gt 0) = \frac{1}{2} \) and so \( \P\left[\sgn(X) = -1\right] = \frac{1}{2} \) also. Suppose that \(X\) and \(Y\) are independent random variables, each with the standard normal distribution. In the order statistic experiment, select the exponential distribution. For example, recall that in the standard model of structural reliability, a system consists of \(n\) components that operate independently. \(f^{*2}(z) = \begin{cases} z, & 0 \lt z \lt 1 \\ 2 - z, & 1 \lt z \lt 2 \end{cases}\), \(f^{*3}(z) = \begin{cases} \frac{1}{2} z^2, & 0 \lt z \lt 1 \\ 1 - \frac{1}{2}(z - 1)^2 - \frac{1}{2}(2 - z)^2, & 1 \lt z \lt 2 \\ \frac{1}{2} (3 - z)^2, & 2 \lt z \lt 3 \end{cases}\), \( g(u) = \frac{3}{2} u^{1/2} \), for \(0 \lt u \le 1\), \( h(v) = 6 v^5 \) for \( 0 \le v \le 1 \), \( k(w) = \frac{3}{w^4} \) for \( 1 \le w \lt \infty \), \(g(c) = \frac{3}{4 \pi^4} c^2 (2 \pi - c)\) for \( 0 \le c \le 2 \pi\), \(h(a) = \frac{3}{8 \pi^2} \sqrt{a}\left(2 \sqrt{\pi} - \sqrt{a}\right)\) for \( 0 \le a \le 4 \pi\), \(k(v) = \frac{3}{\pi} \left[1 - \left(\frac{3}{4 \pi}\right)^{1/3} v^{1/3} \right]\) for \( 0 \le v \le \frac{4}{3} \pi\). . In particular, suppose that a series system has independent components, each with an exponentially distributed lifetime. \sum_{x=0}^z \frac{z!}{x! In the continuous case, \( R \) and \( S \) are typically intervals, so \( T \) is also an interval as is \( D_z \) for \( z \in T \). For the following three exercises, recall that the standard uniform distribution is the uniform distribution on the interval \( [0, 1] \). \(g(u, v, w) = \frac{1}{2}\) for \((u, v, w)\) in the rectangular region \(T \subset \R^3\) with vertices \(\{(0,0,0), (1,0,1), (1,1,0), (0,1,1), (2,1,1), (1,1,2), (1,2,1), (2,2,2)\}\). With \(n = 4\), run the simulation 1000 times and note the agreement between the empirical density function and the probability density function. Note that the joint PDF of \( (X, Y) \) is \[ f(x, y) = \phi(x) \phi(y) = \frac{1}{2 \pi} e^{-\frac{1}{2}\left(x^2 + y^2\right)}, \quad (x, y) \in \R^2 \] From the result above polar coordinates, the PDF of \( (R, \Theta) \) is \[ g(r, \theta) = f(r \cos \theta , r \sin \theta) r = \frac{1}{2 \pi} r e^{-\frac{1}{2} r^2}, \quad (r, \theta) \in [0, \infty) \times [0, 2 \pi) \] From the factorization theorem for joint PDFs, it follows that \( R \) has probability density function \( h(r) = r e^{-\frac{1}{2} r^2} \) for \( 0 \le r \lt \infty \), \( \Theta \) is uniformly distributed on \( [0, 2 \pi) \), and that \( R \) and \( \Theta \) are independent. Show how to simulate the uniform distribution on the interval \([a, b]\) with a random number. Suppose that \( X \) and \( Y \) are independent random variables, each with the standard normal distribution, and let \( (R, \Theta) \) be the standard polar coordinates \( (X, Y) \). The expectation of a random vector is just the vector of expectations. Let \(\bs Y = \bs a + \bs B \bs X\) where \(\bs a \in \R^n\) and \(\bs B\) is an invertible \(n \times n\) matrix. Suppose that \(Z\) has the standard normal distribution. Sketch the graph of \( f \), noting the important qualitative features. As usual, the most important special case of this result is when \( X \) and \( Y \) are independent. Recall that a standard die is an ordinary 6-sided die, with faces labeled from 1 to 6 (usually in the form of dots). In the discrete case, \( R \) and \( S \) are countable, so \( T \) is also countable as is \( D_z \) for each \( z \in T \). The next result is a simple corollary of the convolution theorem, but is important enough to be highligted. The critical property satisfied by the quantile function (regardless of the type of distribution) is \( F^{-1}(p) \le x \) if and only if \( p \le F(x) \) for \( p \in (0, 1) \) and \( x \in \R \). Suppose that \(\bs X = (X_1, X_2, \ldots)\) is a sequence of independent and identically distributed real-valued random variables, with common probability density function \(f\). In many respects, the geometric distribution is a discrete version of the exponential distribution. Using the random quantile method, \(X = \frac{1}{(1 - U)^{1/a}}\) where \(U\) is a random number. As before, determining this set \( D_z \) is often the most challenging step in finding the probability density function of \(Z\). Suppose that \( X \) and \( Y \) are independent random variables with continuous distributions on \( \R \) having probability density functions \( g \) and \( h \), respectively. Probability, Mathematical Statistics, and Stochastic Processes (Siegrist), { "3.01:_Discrete_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.02:_Continuous_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.03:_Mixed_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.04:_Joint_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.05:_Conditional_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.06:_Distribution_and_Quantile_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.07:_Transformations_of_Random_Variables" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.08:_Convergence_in_Distribution" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.09:_General_Distribution_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.10:_The_Integral_With_Respect_to_a_Measure" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.11:_Properties_of_the_Integral" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.12:_General_Measures" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.13:_Absolute_Continuity_and_Density_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "3.14:_Function_Spaces" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, { "00:_Front_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "01:_Foundations" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "02:_Probability_Spaces" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "03:_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "04:_Expected_Value" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "05:_Special_Distributions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "06:_Random_Samples" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "07:_Point_Estimation" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "08:_Set_Estimation" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "09:_Hypothesis_Testing" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "10:_Geometric_Models" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "11:_Bernoulli_Trials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "12:_Finite_Sampling_Models" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "13:_Games_of_Chance" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "14:_The_Poisson_Process" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "15:_Renewal_Processes" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "16:_Markov_Processes" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "17:_Martingales" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "18:_Brownian_Motion" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "zz:_Back_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, [ "article:topic", "license:ccby", "authorname:ksiegrist", "licenseversion:20", "source@http://www.randomservices.org/random" ], https://stats.libretexts.org/@app/auth/3/login?returnto=https%3A%2F%2Fstats.libretexts.org%2FBookshelves%2FProbability_Theory%2FProbability_Mathematical_Statistics_and_Stochastic_Processes_(Siegrist)%2F03%253A_Distributions%2F3.07%253A_Transformations_of_Random_Variables, \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}}}\) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\), \(\renewcommand{\P}{\mathbb{P}}\) \(\newcommand{\R}{\mathbb{R}}\) \(\newcommand{\N}{\mathbb{N}}\) \(\newcommand{\Z}{\mathbb{Z}}\) \(\newcommand{\bs}{\boldsymbol}\) \( \newcommand{\sgn}{\text{sgn}} \), Transformed Variables with Discrete Distributions, Transformed Variables with Continuous Distributions, http://mathworld.wolfram.com/PolarCoordinates.html, source@http://www.randomservices.org/random, status page at https://status.libretexts.org, \(g(y) = f\left[r^{-1}(y)\right] \frac{d}{dy} r^{-1}(y)\).