linear transformation of normal distribution

In a normal distribution, data is symmetrically distributed with no skew. On the other hand, the uniform distribution is preserved under a linear transformation of the random variable. Thus, $ X $ also has the standard Cauchy distribution. Then, any linear transformation of x x is also multivariate normally distributed: y = Ax+ b N (A+ b,AAT). If x_mean is the mean of my first normal distribution, then can the new mean be calculated as : k_mean = x . $g(y) = \frac{1}{8 \sqrt{y}}, \quad 0 \lt y \lt 16$, $g(y) = \frac{1}{4 \sqrt{y}}, \quad 0 \lt y \lt 4$, $g(y) = \begin{cases} \frac{1}{4 \sqrt{y}}, & 0 \lt y \lt 1 \\ \frac{1}{8 \sqrt{y}}, & 1 \lt y \lt 9 \end{cases}$. . Since $1 - U$ is also a random number, a simpler solution is $X = -\frac{1}{r} \ln U$. Note that since $r$ is one-to-one, it has an inverse function $r^{-1}$. Hence the inverse transformation is $ x = (y - a) / b $ and $ dx / dy = 1 / b $. The result now follows from the change of variables theorem. \, ds = e^{-t} \frac{t^n}{n!} It is widely used to model physical measurements of all types that are subject to small, random errors. $ g(y) = \frac{3}{25} \left(\frac{y}{100}\right)\left(1 - \frac{y}{100}\right)^2 $ for $ 0 \le y \le 100 $. Suppose that $ (X, Y, Z) $ has a continuous distribution on $ \R^3 $ with probability density function $ f $, and that $ (R, \Theta, Z) $ are the cylindrical coordinates of $ (X, Y, Z) $. Now we can prove that every linear transformation is a matrix transformation, and we will show how to compute the matrix. Set $k = 1$ (this gives the minimum $U$). 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. For $y \in T$. Linear transformations (or more technically affine transformations) are among the most common and important transformations. I have a pdf which is a linear transformation of the normal distribution: T = 0.5A + 0.5B Mean_A = 276 Standard Deviation_A = 6.5 Mean_B = 293 Standard Deviation_A = 6 How do I calculate the probability that T is between 281 and 291 in Python? probability - Linear transformations in normal distributions Note that $Y$ takes values in $T = \{y = a + b x: x \in S\}$, which is also an interval. It's best to give the inverse transformation: $ x = r \cos \theta $, $ y = r \sin \theta $. The associative property of convolution follows from the associate property of addition: $ (X + Y) + Z = X + (Y + Z) $. The number of bit strings of length $ n $ with 1 occurring exactly $ y $ times is $ \binom{n}{y} $ for $y \in \{0, 1, \ldots, n\}$. The formulas for the probability density functions in the increasing case and the decreasing case can be combined: If $r$ is strictly increasing or strictly decreasing on $S$ then the probability density function $g$ of $Y$ is given by \[ g(y) = f\left[ r^{-1}(y) \right] \left| \frac{d}{dy} r^{-1}(y) \right| \]. (z - x)!} In this case, $ D_z = \{0, 1, \ldots, z\} $ for $ z \in \N $. Using your calculator, simulate 5 values from the Pareto distribution with shape parameter $a = 2$. The linear transformation of a normally distributed random variable is still a normally distributed random variable: . Suppose that $(X_1, X_2, \ldots, X_n)$ is a sequence of indendent real-valued random variables and that $X_i$ has distribution function $F_i$ for $i \in \{1, 2, \ldots, n\}$. 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 $. For $i \in \N_+$, the probability density function $f$ of the trial variable $X_i$ is $f(x) = p^x (1 - p)^{1 - x}$ for $x \in \{0, 1\}$. Related. Linear Algebra - Linear transformation question A-Z related to countries Lots of pick movement . Suppose that $(X_1, X_2, \ldots, X_n)$ is a sequence of independent real-valued random variables, with a common continuous distribution that has probability density function $f$. Suppose that two six-sided dice are rolled and the sequence of scores $(X_1, X_2)$ is recorded. I want to compute the KL divergence between a Gaussian mixture distribution and a normal distribution using sampling method. Note that since $ V $ is the maximum of the variables, $\{V \le x\} = \{X_1 \le x, X_2 \le x, \ldots, X_n \le x\}$. In particular, the $ n $th arrival times in the Poisson model of random points in time has the gamma distribution with parameter $ n $. Suppose that $X$ has a continuous distribution on a subset $S \subseteq \R^n$ and that $Y = r(X)$ has a continuous distributions on a subset $T \subseteq \R^m$. $\left|X\right|$ has probability density function $g$ given by $g(y) = f(y) + f(-y)$ for $y \in [0, \infty)$. $ f $ increases and then decreases, with mode $ x = \mu $. Then $\bs Y$ is uniformly distributed on $T = \{\bs a + \bs B \bs x: \bs x \in S\}$. Uniform distributions are studied in more detail in the chapter on Special Distributions. $\left|X\right|$ has probability density function $g$ given by $g(y) = 2 f(y)$ for $y \in [0, \infty)$. Transform Data to Normal Distribution in R: Easy Guide - Datanovia Normal Distribution | Examples, Formulas, & Uses - Scribbr With $n = 5$ run the simulation 1000 times and compare the empirical density function and the probability density function. In the order statistic experiment, select the uniform distribution. An ace-six flat die is a standard die in which faces 1 and 6 occur with probability $\frac{1}{4}$ each and the other faces with probability $\frac{1}{8}$ each. Note that he minimum on the right is independent of $T_i$ and by the result above, has an exponential distribution with parameter $\sum_{j \ne i} r_j$. Recall that the (standard) gamma distribution with shape parameter $n \in \N_+$ has probability density function \[ g_n(t) = e^{-t} \frac{t^{n-1}}{(n - 1)! 24/7 Customer Support. Location transformations arise naturally when the physical reference point is changed (measuring time relative to 9:00 AM as opposed to 8:00 AM, for example). Tour Start here for a quick overview of the site Help Center Detailed answers to any questions you might have Meta Discuss the workings and policies of this site Simple addition of random variables is perhaps the most important of all transformations. Also, for $ t \in [0, \infty) $, \[ g_n * g(t) = \int_0^t g_n(s) g(t - s) \, ds = \int_0^t e^{-s} \frac{s^{n-1}}{(n - 1)!} Thus suppose that $\bs X$ is a random variable taking values in $S \subseteq \R^n$ and that $\bs X$ has a continuous distribution on $S$ with probability density function $f$. Random variable $T$ has the (standard) Cauchy distribution, named after Augustin Cauchy. Assuming that we can compute $F^{-1}$, the previous exercise shows how we can simulate a distribution with distribution function $F$. Let X be a random variable with a normal distribution f ( x) with mean X and standard deviation X : PDF Basic Multivariate Normal Theory - Duke University The independence of $ X $ and $ Y $ corresponds to the regions $ A $ and $ B $ being disjoint. Recall that the Poisson distribution with parameter $t \in (0, \infty)$ has probability density function $f$ given by \[ f_t(n) = e^{-t} \frac{t^n}{n! A possible way to fix this is to apply a transformation. If $ (X, Y) $ takes values in a subset $ D \subseteq \R^2 $, then for a given $ v \in \R $, the integral in (a) is over $ \{x \in \R: (x, v / x) \in D\} $, and for a given $ w \in \R $, the integral in (b) is over $ \{x \in \R: (x, w x) \in D\} $. The first image below shows the graph of the distribution function of a rather complicated mixed distribution, represented in blue on the horizontal axis. The Irwin-Hall distributions are studied in more detail in the chapter on Special Distributions. Clearly convolution power satisfies the law of exponents: $ f^{*n} * f^{*m} = f^{*(n + m)} $ for $ m, \; n \in \N $. Then, with the aid of matrix notation, we discuss the general multivariate distribution. Conversely, any continuous distribution supported on an interval of $\R$ can be transformed into the standard uniform distribution. Types Of Transformations For Better Normal Distribution SummaryThe problem of characterizing the normal law associated with linear forms and processes, as well as with quadratic forms, is considered. Check if transformation is linear calculator - Math Practice On the other hand, $W$ has a Pareto distribution, named for Vilfredo Pareto. Suppose again that $(T_1, T_2, \ldots, T_n)$ is a sequence of independent random variables, and that $T_i$ has the exponential distribution with rate parameter $r_i \gt 0$ for each $i \in \{1, 2, \ldots, n\}$. Thus, suppose that random variable $X$ has a continuous distribution on an interval $S \subseteq \R$, with distribution function $F$ and probability density function $f$. Let M Z be the moment generating function of Z . As usual, we will let $G$ denote the distribution function of $Y$ and $g$ the probability density function of $Y$. Suppose that $Z$ has the standard normal distribution, and that $\mu \in (-\infty, \infty)$ and $\sigma \in (0, \infty)$. The distribution arises naturally from linear transformations of independent normal variables. We introduce the auxiliary variable $ U = X $ so that we have bivariate transformations and can use our change of variables formula. The result in the previous exercise is very important in the theory of continuous-time Markov chains. How to Transform Data to Better Fit The Normal Distribution Here we show how to transform the normal distribution into the form of Eq 1.1: Eq 3.1 Normal distribution belongs to the exponential family. Suppose that $r$ is strictly decreasing on $S$. Work on the task that is enjoyable to you. $\left|X\right|$ has distribution function $G$ given by$G(y) = 2 F(y) - 1$ for $y \in [0, \infty)$. Scale transformations arise naturally when physical units are changed (from feet to meters, for example). Vary $n$ with the scroll bar and note the shape of the probability density function. probability - Normal Distribution with Linear Transformation Recall that $ F^\prime = f $. 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. The distribution of $ Y_n $ is the binomial distribution with parameters $n$ and $p$. Show how to simulate, with a random number, the exponential distribution with rate parameter $r$. The change of temperature measurement from Fahrenheit to Celsius is a location and scale transformation. Suppose that $(T_1, T_2, \ldots, T_n)$ is a sequence of independent random variables, and that $T_i$ has the exponential distribution with rate parameter $r_i \gt 0$ for each $i \in \{1, 2, \ldots, n\}$. Hence the PDF of W is \[ w \mapsto \int_{-\infty}^\infty f(u, u w) |u| du \], Random variable $ V = X Y $ has probability density function \[ v \mapsto \int_{-\infty}^\infty g(x) h(v / x) \frac{1}{|x|} dx \], Random variable $ W = Y / X $ has probability density function \[ w \mapsto \int_{-\infty}^\infty g(x) h(w x) |x| dx \]. a^{x} b^{z - x} \\ & = e^{-(a+b)} \frac{1}{z!} The Jacobian of the inverse transformation is the constant function $\det (\bs B^{-1}) = 1 / \det(\bs B)$. Systematic component - $x$ is the explanatory variable (can be continuous or discrete) and is linear in the parameters. \exp\left(-e^x\right) e^{n x}\) for $x \in \R$. Suppose that $X$ has the probability density function $f$ given by $f(x) = 3 x^2$ for $0 \le x \le 1$. Next, for $ (x, y, z) \in \R^3 $, let $ (r, \theta, z) $ denote the standard cylindrical coordinates, so that $ (r, \theta) $ are the standard polar coordinates of $ (x, y) $ as above, and coordinate $ z $ is left unchanged. But a linear combination of independent (one dimensional) normal variables is another normal, so aTU is a normal variable. Clearly we can simulate a value of the Cauchy distribution by $ X = \tan\left(-\frac{\pi}{2} + \pi U\right) $ where $ U $ is a random number. calculus - Linear transformation of normal distribution - Mathematics Using the change of variables theorem, the joint PDF of $ (U, V) $ is $ (u, v) \mapsto f(u, v / u)|1 /|u| $. This is known as the change of variables formula. For $ u \in (0, 1) $ recall that $ F^{-1}(u) $ is a quantile of order $ u $. Both distributions in the last exercise are beta distributions. Note that $ Z $ takes values in $ T = \{z \in \R: z = x + y \text{ for some } x \in R, y \in S\} $. As we all know from calculus, the Jacobian of the transformation is $ r $. We will solve the problem in various special cases. $g(u) = \frac{a / 2}{u^{a / 2 + 1}}$ for $ 1 \le u \lt \infty$, $h(v) = a v^{a-1}$ for $ 0 \lt v \lt 1$, $k(y) = a e^{-a y}$ for $ 0 \le y \lt \infty$, Find the probability density function $ f $ of $X = \mu + \sigma Z$. Using the theorem on quotient above, the PDF $ f $ of $ T $ is given by \[f(t) = \int_{-\infty}^\infty \phi(x) \phi(t x) |x| dx = \frac{1}{2 \pi} \int_{-\infty}^\infty e^{-(1 + t^2) x^2/2} |x| dx, \quad t \in \R\] Using symmetry and a simple substitution, \[ f(t) = \frac{1}{\pi} \int_0^\infty x e^{-(1 + t^2) x^2/2} dx = \frac{1}{\pi (1 + t^2)}, \quad t \in \R \]. Then $U$ is the lifetime of the series system which operates if and only if each component is operating. (These are the density functions in the previous exercise). These can be combined succinctly with the formula $ f(x) = p^x (1 - p)^{1 - x} $ for $ x \in \{0, 1\} $. Distributions with Hierarchical models. Convolution can be generalized to sums of independent variables that are not of the same type, but this generalization is usually done in terms of distribution functions rather than probability density functions. Linear transformation of multivariate normal random variable is still multivariate normal. $\P(Y \in B) = \P\left[X \in r^{-1}(B)\right]$ for $B \subseteq T$. Initialy, I was thinking of applying "exponential twisting" change of measure to y (which in this case amounts to changing the mean from $\mathbf{0}$ to $\mathbf{c}$) but this requires taking . If $ A \subseteq (0, \infty) $ then \[ \P\left[\left|X\right| \in A, \sgn(X) = 1\right] = \P(X \in A) = \int_A f(x) \, dx = \frac{1}{2} \int_A 2 \, f(x) \, dx = \P[\sgn(X) = 1] \P\left(\left|X\right| \in A\right) \], The first die is standard and fair, and the second is ace-six flat. Suppose that $r$ is strictly increasing on $S$. From part (a), note that the product of $n$ distribution functions is another distribution function. Suppose that $X$ has the exponential distribution with rate parameter $a \gt 0$, $Y$ has the exponential distribution with rate parameter $b \gt 0$, and that $X$ and $Y$ are independent. from scipy.stats import yeojohnson yf_target, lam = yeojohnson (df ["TARGET"]) Yeo-Johnson Transformation To check if the data is normally distributed I've used qqplot and qqline . $ \P\left(\left|X\right| \le y\right) = \P(-y \le X \le y) = F(y) - F(-y) $ for $ y \in [0, \infty) $. \Only if part" Suppose U is a normal random vector. This follows from part (a) by taking derivatives with respect to $ y $ and using the chain rule. $U = \min\{X_1, X_2, \ldots, X_n\}$ has distribution function $G$ given by $G(x) = 1 - \left[1 - F_1(x)\right] \left[1 - F_2(x)\right] \cdots \left[1 - F_n(x)\right]$ for $x \in \R$. 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)$. Find linear transformation associated with matrix | Math Methods Linear Transformations - gatech.edu Recall that $ \frac{d\theta}{dx} = \frac{1}{1 + x^2} $, so by the change of variables formula, $ X $ has PDF $g$ given by \[ g(x) = \frac{1}{\pi \left(1 + x^2\right)}, \quad x \in \R \]. Note that $\bs Y$ takes values in $T = \{\bs a + \bs B \bs x: \bs x \in S\} \subseteq \R^n$. The Jacobian is the infinitesimal scale factor that describes how $n$-dimensional volume changes under the transformation. Bryan 3 years ago The inverse transformation is $\bs x = \bs B^{-1}(\bs y - \bs a)$. 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 \]. Random component - The distribution of $Y$ is Poisson with mean $\lambda$. Using the change of variables theorem, If $ X $ and $ Y $ have discrete distributions then $ Z = X + Y $ has a discrete distribution with probability density function $ g * h $ given by \[ (g * h)(z) = \sum_{x \in D_z} g(x) h(z - x), \quad z \in T \], If $ X $ and $ Y $ have continuous distributions then $ Z = X + Y $ has a continuous distribution with probability density function $ g * h $ given by \[ (g * h)(z) = \int_{D_z} g(x) h(z - x) \, dx, \quad z \in T \], In the discrete case, suppose $ X $ and $ Y $ take values in $ \N $. Vary $n$ with the scroll bar and set $k = n$ each time (this gives the maximum $V$). In the continuous case, $ R $ and $ S $ are typically intervals, so $ T $ is also an interval as is $ D_z $ for $ z \in T $.