학술/수학

TRANSFORMING DENSITY FUNCTIONS

ksyoon 2019. 4. 22. 17:58

어떤 변수의 probability density function을 알고 있다면 transformation function을 사용하여 다른 것도 쉽게 구할 수있다.

일반적인 적분에서 변수 치환의 방법은 흔히 쓰는 방법인데

$$ \int^b _a f(x)dx=  \int^{y(b)} _ {y(a)} f(x(y)) {dx \over dy} dy $$

가 된다.

이 공식은 probability density function를 변환할 때 사용 된다.

\( X \)

라 는 랜덤 변수가 있고 이것의 density function \( f(x) \) 라고 하면

정의에 의하여

$$ P(a \leq X  < b) = \int^b _a f(x)dx $$

가  된다.

어떤 랜덤 변수의 함수는 또한 랜덤 변수인데, 만일 \( y \) 가 어떤 transformation function이면 \( y(X) \) 랜덤 변수가 된다.

가령 \( Y=y(X)\) 라고 하자.

만일 \( X=a\) 이면 \(Y=y(a), X=b\) 이면 \(Y=y(b)\)가 된다. 또한 \(a \leq X < b\) 라면 \(y(a) \leq Y < y(b)\) 가 되고 \(P(y(a) \leq Y < y(b)) = P(a\leq X < b) \)이다.

$$ P(y(a) \leq Y < y(b)) = P(a \leq X < b) =  \int^b _a f(x)dx=  \int^{y(b)} _ {y(a)} f(x(y)) {dx \over dy} dy$$

여기서 \(f(x(y)) {dx \over dy} dy \)는 \( y\)로 모두 표시되어 있다 이 적분을 \(  g(y) \) 라고 하면

$$ P(y(a) \leq Y < y(b)) =   \int^{y(b)} _ {y(a)} g(y)dy$$

 

Observation and Constraints

여기에는 중요한 조건이 있는데 다음과 같다.

  • \( f(x) \)는 모든 \( x\)에 대하여 single value 를 가져야한다.
  • \( f(x) \geq 0 \) for all \( x\)
  • \( \int ^{+\infty} _{-\infty} f(x)dx = 1 \)

많은 경우에 \( f(x) \)는 정해진 영역에서만 정의 되고 그 외의 영역에서는 0인 경우로 정의 된다. 예를 들어 \( 2xcos \,x^2 \) 는 다음과 같이 정의 된다.

$$ f(x) =  \begin{cases} 2xcos \,x^2 , \quad 0 \leq x < 2 , 0 \quad otherwise \end{cases} $$

이 함수는 single value를 가지고 non-negative이고 \( -\infty \) 에서 \( +\infty \) 까지의 적분 값이 1이다.

 

모든 transformation function \( y(x) \)는 암묵적으로 다음과 같은 제한 조건을 갖는다.

  • \( x\)의 유효 범위에서 \( y(x) \)와 역함수 \( x(y) \)는 정의 되어야 하고 single-value를 가져야 한다.
  • 이 범위에서 \( {dx \over dy} \) 가 정의 되어야하고 \( {dx \over dy} \ge 0\) 또는 \( {dx \over dy} \leq 0\)

만약 \( {dx \over dy} \)의 \( \pm \)이 바뀌면 \( y(x) \) 에 대한 \( x\)의값이 여러개 존재하게 된다. 즉 이러한 조건은 \( y(x) \)  이 단조 증가 또는 단조 감소 함수이어야 한다. 그러므로

$$ g(y) = f(x(y)) |{dx \over dy}| $$