This is not homework. Case I is mostly for background. The real questions are in Case II.

Case I (one dimension):

a. Suppose X is a continuous r.v. with pdf f

[tex]f_Y(y) = f_X(g^{-1}(y) | (g^{-1})'(y) |[/tex].

All of the above is well-known (http://en.wikipedia.org/wiki/Density_function" [Broken]).

b. Now suppose the transformation g is not one-to-one. Denote the real roots of y = g(x) by x

[tex]f_Y(y) = \sum_k = \frac{f_X(x_k)}{|g'(x_k)|}[/tex].

Again, this is (relatively) well-known.

Case II (n dimensions):

a. The case for n functions of n random variables where g is one-to-one is well-known (see link above):

[tex]f_Y(y) = f_X(g^{-1}(y))| J_{g^{-1}}(g(x)) |[/tex].

b. What about the case for n functions of n random variables where g is not one-to-one? Is there an analogous result to that in Case Ib? Can someone provide a reference?

c. What about the case for Y = g(X

Thanks for any insights you can provide.

Case I (one dimension):

a. Suppose X is a continuous r.v. with pdf f

_{X}(x), y = g(x) is one-to-one, and the inverse x = g^{-1}(y) exists. Then the pdf of Y = g(X) is found by[tex]f_Y(y) = f_X(g^{-1}(y) | (g^{-1})'(y) |[/tex].

All of the above is well-known (http://en.wikipedia.org/wiki/Density_function" [Broken]).

b. Now suppose the transformation g is not one-to-one. Denote the real roots of y = g(x) by x

_{k}, i.e., y = g(x_{1}) = ... = g(x_{k}). Then the pdf of Y = g(X) is[tex]f_Y(y) = \sum_k = \frac{f_X(x_k)}{|g'(x_k)|}[/tex].

Again, this is (relatively) well-known.

Case II (n dimensions):

a. The case for n functions of n random variables where g is one-to-one is well-known (see link above):

[tex]f_Y(y) = f_X(g^{-1}(y))| J_{g^{-1}}(g(x)) |[/tex].

b. What about the case for n functions of n random variables where g is not one-to-one? Is there an analogous result to that in Case Ib? Can someone provide a reference?

c. What about the case for Y = g(X

_{1}, x_{2}, ..., x_{n}) where g : R^{n}-> R? [Note: Though g might still be one-to-one, we suppose that it is not.] The approach to Case IIc with which I am familiar is to construct an additional n-1 functions (often just using the identity function, but sometimes more effort is needed to choose functions so the Jacobian is nonzero), and then use the known method of Case IIa. Is there an approach analogous to that in Case Ib or Case IIb? Can someone provide a reference?Thanks for any insights you can provide.

Last edited by a moderator: