how can i calculate z^2 distibution?

조회 수: 1 (최근 30일)
jinseok shin
jinseok shin 2021년 2월 26일
댓글: jinseok shin 2021년 3월 3일
Hi, I am currently trying to calculate z^2 distribution when z is folloing (0,a) gaussian distribution
how can i calculate z^2 distribution?
  댓글 수: 1
Hernia Baby
Hernia Baby 2021년 2월 26일
Do you want to calculate Chi-square distribution?
The distribution depends on the degree of freedom.

댓글을 달려면 로그인하십시오.

채택된 답변

David Goodmanson
David Goodmanson 2021년 3월 2일
편집: David Goodmanson 2021년 3월 2일
Hello JS
By gaussian I assume you mean the normal distribution, mean 0, std deviation 'a'
f(z) = N*exp(-(z/a)^2/2)
with N = 1/(sqrt(2)*sqrt(pi)*a)
y = z^2 has one sign while z can have either sign. Since the normal distribution is symmetric, to make life easier assume the initial distribution is one-sided, with z>=0 only, and double the height of the distribution to keep the normalization correct:
f(z) = 2*N*exp(-(z/a)^2/2) z>=0
Now that z >= 0 we can use z = sqrt(y). The the idea with pdfs is to find a function g(y) such that
f(z) dz = g(y) dy
with z = y^(1/2)
which means that dz = (1/2) y^(-1/2) dy
plug that in
2*N*exp(-(y/(2*a^2))) (1/2) y(-1/2) dy = g(y) dy
g(y) = N*exp(-(y/(2*a^2)))/sqrt(y)
which is a chi-squared probability distribution function (with 1 degree of freeedom) as noted by Jeff. This function has the disadvantage of being unbounded as y --> 0, although its area is still 1 as required.
If you meant y = (z/a)^2 instead of y = z^2, then the expression would be the standard chi-squared
g(y) = N*exp(-(y/2))/sqrt(y) N = 1/(sqrt(2)*sqrt(pi))

추가 답변 (1개)

Jeff Miller
Jeff Miller 2021년 2월 27일
The distribution of z^2/a^2 is a chi-square with 1 df.
  댓글 수: 2
Jeff Miller
Jeff Miller 2021년 3월 2일
To normalize. If y is normal(0,1), then the distribution of y^ is chi-square with 1 df.
For your z ~ N(0,a), y = z/a is N(0,1), so y^2 = z^2/a^2 is chi-square(1).
jinseok shin
jinseok shin 2021년 3월 2일
Then how can i make pdf? Can I divide the Gaussian distribution by the variance?

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Hypothesis Tests에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by