Reading [this](https://math.stackexchange.com/questions/1476568/strange-definitions-about-basic-probability-need-clarification/1580395#comment3335933_1580395) post, I realize that is possible to have *another* type of PDF (probability density function) in the special case when the sample space is an Euclidean space.


Usually, we have a probability space $(\Omega, \mathcal{F}, P)$, a random variable $X:\Omega\to\mathbb{R}$ and the PDF of $X$, $f\_X:\mathbb{R}\to[0, \infty)$. The whole point of the PDF is that we can integrate it over any borel set $A$ (contained in the codomain of $X$) and get the probability $ P[X\in A] = \int\_A f\_X\ dx$. 


What about the case when we have that $\Omega = \mathbb{R}$ and $\mathcal{F}$ is the Borel sigma algebra over $\mathbb{R}$? Looks like it's easier to define the PDF to be a function $f:\Omega\to[0,\infty)$ such that $P[A] = \int\_A f\ dx$. So we still can use this PDF to get the probability of events via integration over borel sets. The idea in the book (mentioned on the other post) is that $\Omega$ can be a manifold, so we integrate the PDF over a subset contained in this manifold to get the probability of that subset. If $X:\Omega\to\mathbb{R}$ is a random variable (not necessarily the map $x\to x$), we could talk about $f\_X$, the PDF of $X$, such as described before. But in this sample space, we can use the *other* PDF $f$. So, instead calculating $\int\_Af\_X\ dx = P[X\in A]$, we can calculate $\int\_{\{X\in A\}}f\ dx = P[X\in A]$. Note that the first integral is over the codomain of $X$, but the second is over the sample space $\Omega$. 


[This](http://www.ams.org/notices/200511/what-is.pdf) article in particular is an example of what I'm saying. Instead defining real random variables $X\_{ij}$ as the entries of a random matrix, it's possible to see the sample space as $\mathbb{R}^{n^2}$ and **define** a probability distribution on $\mathbb{R}^{n^2}$ to work with.


I want to know if there is some technical detail I should be aware of and how common is this. Thank you.



If the sample space is an Euclidean Space, we can use a different type of PDF
Your formula 
$$\int\_{\{X\in A\}}f\ dx = P[X\in A]\tag{1}$$
(I guess you wanted to say it should hold for all Borel $A\subseteq\mathbb R$) can be rewritten as $\int\_B f\ dx = P(B)$ for all $B$ in $\sigma(X)$, the smallest sigma-algebra with respect to which the random variable (r.v.) $X$ is measurable. So, condition $(1)$ means precisely that $f$ is the density (that is, the Radon--Nikodym derivative) 
of the restriction $P|\_{\sigma(X)}$ of the measure $P$ to $\sigma(X)$ with respect to the restriction $\lambda|\_{\sigma(X)}$ of the Lebesgue measure $\lambda$ to $\sigma(X)$:
$$f=\frac{dP|\_{\sigma(X)}}{d\lambda|\_{\sigma(X)}}. \tag{2}
$$
Such a density (when it exists) depends on $X$ only through $\sigma(X)$ and thus usually provides little information on the r.v. $X$, and it provides no information on its distribution. 
To illustrate this point, suppose that $(2)$ holds for a r.v. $X$ and some "alternative pdf" $f$. Let $Y\_1:=1-X$, $Y\_2:=X^3$, and $Y\_3:=\arctan X$. Then $\sigma(Y\_1)=\sigma(Y\_2)=\sigma(Y\_3)=\sigma(X)$ and hence the same "alternative pdf" $f$ will serve all these r.v.'s: $X,Y\_1,Y\_2,Y\_3$. 


Somewhat more generally, suppose in addition that the distribution function (d.f.) $F$ of $X$ is continuous and strictly increasing (c.s.i.). Let $G$ be any other c.s.i. d.f. Let $Y:=G^{-1}(F(X))$. Then the d.f. of the r.v. $Y$ is $G$. On the other hand, $Y$ is a c.s.i. transformation of $X$ and hence $\sigma(Y)=\sigma(X)$. So, $Y$ is served by the same "alternative pdf" $f$ as $X$, even though the d.f. $Y$ is an arbitrary c.s.i. d.f. $G$.


Therefore, I don't see why such a "alternative pdf" would be of use. 


As for the linked note by Diaconis, it seems to have nothing to do with formula $(1)$, even though he does let $\Omega=\mathbb R^{n^2}$ (nothing wrong or unusual with that). What Diaconis does in that note to answer the question in its title "What is a Random Matrix?" is just provide a particular construction, via Gram--Schmidt orthogonalization, of the Haar measure on the orthogonal group $O\_n$, starting with the standard Gaussian measure on $\mathbb R^{n^2}$, which latter is the product of $n^2$ copies of the standard Gaussian measure on $\mathbb R$. 


