যখনএবংস্বাধীনভাবে

12

$X$ এবং স্বাধীনভাবে র্যান্ডম ভেরিয়েবলগুলি বিতরণ করা হয় যেখানে এবং । এর বিতরণ কী ? $Y$ $X\sim\chi^2_{(n-1)}$ $Y\sim\text{Beta}\left(\frac{n}{2}-1,\frac{n}{2}-1\right)$ $Z=(2Y-1)\sqrt X$

এর যৌথ ঘনত্ব দেওয়া হয়েছে $(X,Y)$

f_{X, Y} (x, y) = f_{X} (x) f_{Y} (y) = \frac{e^{- \frac{x}{2}} x^{\frac{n - 1}{2} - 1}}{2^{\frac{n - 1}{2}} Γ (\frac{n - 1}{2})} \cdot \frac{y^{\frac{n}{2} - 2} (1 - y)^{\frac{n}{2} - 2}}{B (\frac{n}{2} - 1, \frac{n}{2} - 1)} 1_{{x > 0, 0 < y < 1}}

$f_{X,Y}(x,y)=f_X(x)f_Y(y)=\frac{e^{-\frac{x}{2}}x^{\frac{n-1}{2}-1}}{2^{\frac{n-1}{2}}\Gamma\left(\frac{n-1}{2}\right)}\cdot\frac{y^{\frac{n}{2}-2}(1-y)^{\frac{n}{2}-2}}{B\left(\frac{n}{2}-1,\frac{n}{2}-1\right)}\mathbf1_{\{x>0\,,\,0<y<1\}}$

ভেরিয়েবলের পরিবর্তন যেমন এবং , $(X,Y)\mapsto(Z,W)$ $Z=(2Y-1)\sqrt X$ $W=\sqrt X$

আমি যুগ্ম ঘনত্ব পেতে যেমন $(Z,W)$

f_{Z, W} (z, w) = \frac{e^{- \frac{w^{2}}{2}} w^{n - 3} {(\frac{1}{4} - \frac{z^{2}}{4 w^{2}})}^{\frac{n}{2} - 2}}{2^{\frac{n - 1}{2}} Γ (\frac{n - 1}{2}) B (\frac{n}{2} - 1, \frac{n}{2} - 1)} 1_{{w > 0, | z | < w}}

$f_{Z,W}(z,w)=\frac{e^{-\frac{w^2}{2}}w^{n-3}\left(\frac{1}{4}-\frac{z^2}{4w^2}\right)^{\frac{n}{2}-2}}{2^{\frac{n-1}{2}}\Gamma\left(\frac{n-1}{2}\right)B\left(\frac{n}{2}-1,\frac{n}{2}-1\right)}\mathbf1_{\{w>0\,,\,|z|<w\}}$

প্রান্তিক পিডিএফ তখন , যা আমাকে কোথাও নিয়ে যায় না। $Z$ $f_Z(z)=\displaystyle\int_{|z|}^\infty f_{Z,W}(z,w)\,\mathrm{d}w$

আবার, এর বিতরণ ফাংশনটি সন্ধান করার সময় , একটি অসম্পূর্ণ বিটা / গামা ফাংশন প্রদর্শিত হবে: $Z$

$F_Z(z)=\Pr(Z\le z)$

$\quad\qquad=\Pr((2Y-1)\sqrt X\le z)=\displaystyle\iint_{(2y-1)\sqrt{x}\le z}f_{X,Y}(x,y)\,\mathrm{d}x\,\mathrm{d}y$

এখানে ভেরিয়েবলের উপযুক্ত পরিবর্তন কী? বিতরণটি খুঁজে পাওয়ার অন্য কোনও উপায় আছে কি ? $Z$

আমি চি-স্কোয়ার্ড, বিটা, 'এফ' এবং 'টি' বিতরণগুলির মধ্যে বিভিন্ন সম্পর্ক ব্যবহার করার চেষ্টা করেছি তবে কিছুই কাজ করছে বলে মনে হচ্ছে না। সম্ভবত আমি স্পষ্ট কিছু মিস করছি।

@ ফ্রেঞ্চিস দ্বারা উল্লিখিত হিসাবে, এই রূপান্তরটি বক্স-মেলার ট্রান্সফর্মের একটি সাধারণীকরণ।

— StubbornAtom
সূত্র

4

বক্স-মুলার রূপান্তরকে সাধারণীকরণের মতো দেখায়

— ফ্রান্সিস

10

এখানে একটি বীজগণিতিক প্রমাণ রয়েছে। আমি এর পরিবর্তে (স্কোয়ার নয়) করতে যাচ্ছি যাতে আমাদের । এগুলি বৈধ ঘনত্বের নিশ্চয়তাযুক্ত তাই আমি স্বাভাবিককরণের ধ্রুবকগুলি ট্র্যাক করব না track আমাদের কাছে যাক এবং তাই বিপরীত রূপান্তরগুলির হয় এবং । এটি আমাদের দেয় । এটি আমাদের দিকে নিয়ে যায় $X \sim \chi_{n-1}$ $Z := (2Y - 1)X$

f_{X, Y} (x, y) \propto x^{n - 2} e^{- x^{2} / 2} {[y (1 - y)]}^{n / 2 - 2} 1_{{0 < x, 0 < y < 1}} .

$f_{X,Y}(x,y) \propto x^{n-2}e^{-x^2/2}\left[y(1-y)\right]^{n/2-2} \mathbf 1_{\{0 < x, \, 0 < y < 1\}}.$

Z = (2 y - 1) X

$Z = (2y-1)X$

W = X

$W=X$

x (z, w) = w

$x(z,w) = w$

y (z, w) = \frac{z + w}{2 w} = \frac{z}{2 w} + \frac{1}{2}

$y(z,w) = \frac{z+w}{2w} = \frac{z}{2w} + \frac 12$

| J | = \frac{1}{2 w}

$|J| = \frac 1{2w}$

f_{Z, W} (z, w) \propto w^{n - 1} e^{- w^{2} / 2} {[\frac{z + w}{2 w} (1 - \frac{z + w}{2 w})]}^{n / 2 - 2} 1_{{0 < w, - 1 < \frac{z}{w} < 1}}

$f_{Z,W}(z,w) \propto w^{n-1}e^{-w^2/2} \left[\frac{z+w}{2w} \left( 1 - \frac{z+w}{2w}\right)\right]^{n/2-2} \mathbf 1_{\{0 < w, \, -1 < \frac zw < 1\}}$

\propto w e^{- w^{2} / 2} {(w^{2} - z^{2})}^{n / 2 - 2} 1_{{| z | < w}} .

$\propto w e^{-w^2/2} \left(w^2-z^2\right)^{n/2-2} \mathbf 1_{\{|z| < w\}}.$ সুতরাং

f_{Z} (z) \propto \int_{w > | z |} w e^{- w^{2} / 2} {(w^{2} - z^{2})}^{n / 2 - 2} d w .

$f_Z(z) \propto \int_{w > |z|} w e^{-w^2/2} \left(w^2-z^2\right)^{n/2-2} \,\text dw.$

সুবিধার জন্য । পেতে উভয় পক্ষকে এখন তাই । এটি আমাদেরকে এই চূড়ান্ত অবিচ্ছেদ্য উপর নির্ভর করে না , তাই আমরা দেখিয়েছি যে , সুতরাং therefore $m = n/2-2$ $e^{z^2/2}$

e^{z^{2} / 2} f_{Z} (z) \propto \int_{| z |}^{\infty} w e^{- (w^{2} - z^{2}) / 2} (w^{2} - z^{2})^{m} d w .

$e^{z^2/2} f_Z(z) \propto \int_{|z|}^\infty w e^{-(w^2-z^2)/2}(w^2-z^2)^m \,\text dw.$

2 u = w^{2} - z^{2}

$2u = w^2 - z^2$

d u = w d w

$\text du = w\,\text dw$

e^{z^{2} / 2} f_{Z} (z) \propto 2^{m} \int_{0}^{\infty} u^{m} e^{- u} d u = 2^{m} Γ (m + 1) .

$e^{z^2/2} f_Z(z) \propto 2^m \int_0^\infty u^m e^{-u} \,\text du = 2^m \Gamma(m+1).$

z

$z$

e^{z^{2} / 2} f_{Z} (z) \propto 1

$e^{z^2/2} f_Z(z) \propto 1$

Z \sim N (0, 1) .

$Z \sim \mathcal N(0, 1).$

— jld
সূত্র

1

+1 টি। আমি আনন্দিত আপনি এই উত্তরটি পুনরুদ্ধার করেছেন, কারণ এটি সমস্ত মানকে অন্তর্ভুক্ত করে , কেবলমাত্র অবিচ্ছেদ্য নয়।

n

$n$

— হোবার

@whuber ধন্যবাদ, আমি একরকম করা এর পরিবর্তে এবং এটা আমার কিছুটা সময় লেগেছিল জিনিসটা কেন আমি অদ্ভুত আচরণ যখন পেয়ে ছিল বিজোড়

z^{2} - w^{2}

$z^2-w^2$

w^{2} - z^{2}

$w^2-z^2$

n

$n$

— jld

9

$2Y-1$ গোলকটিতে অভিন্ন বিতরণের এক স্থানাঙ্কের $n-1$ মতো বিতরণ করা হয় ; এর আইআইডি স্ট্যান্ডার্ড নরমাল ভেরিয়েটের বর্গের যোগফলের বিতরণ রয়েছে ; এবং এই দুটি পরিমাণ স্বতন্ত্র। জ্যামিতিকভাবে one এর একটি স্থানাঙ্কের বিতরণ রয়েছে: এটির একটি অবশ্যই সাধারণ বন্টন থাকতে হবে। $X$ $n-1$ $(2Y-1)\sqrt{X}$

(এই যুক্তিটি অবিচ্ছেদ্য ক্ষেত্রে প্রযোজ্য )) $n=2,3,4,\ldots$

আপনার যদি কিছু সংখ্যক বোঝার প্রয়োজন হয় (যা সর্বদা জ্ঞানী, কারণ এটি যুক্তি এবং গণনার ক্ষেত্রে ত্রুটিগুলি উদঘাটন করতে পারে), অনুকরণ করুন:

সিমুলেটেড ফলাফল এবং দাবি করা মান সাধারণ বিতরণের মধ্যে চুক্তিটি এর মানগুলির এই পরিসীমা জুড়ে দুর্দান্ত । $n$

Rআপনি যদি চান তবে এই প্লটগুলি তৈরি করে এমন কোডটি নিয়ে আরও পরীক্ষা করুন ।

n.sim <- 1e5
n <- 2:5
X <- data.frame(Z = c(sapply(n, function(n){
  y <- rbeta(n.sim, n/2-1, n/2-1)  # Generate values of Y
  x <- rchisq(n.sim, n-1)          # Generate values of X
  (2*y - 1) * sqrt(x)              # Return the values of Z
})), n=factor(rep(n, each=n.sim)))

library(ggplot2)
#--Create points along the graph of a standard Normal density
i <- seq(min(z), max(z), length.out=501)
U <- data.frame(X=i, Z=dnorm(i))

#--Plot histograms on top of the density graphs
ggplot(X, aes(Z, ..density..)) + 
  geom_path(aes(X,Z), data=U, size=1) +
  geom_histogram(aes(fill=n), bins=50, alpha=0.5) + 
  facet_wrap(~ n) + 
  ggtitle("Histograms of Simulated Values of Z",
          paste0("Sample size ", n.sim))

— whuber
সূত্র

1

ধন্যবাদ, স্টাবর্ন প্যারামিটারগুলি সামঞ্জস্যপূর্ণ তা বিবেচনা করে না, অন্যথায় উপসংহারটি ভুল। আমি এটা ঠিক করব.

— হোবার

3

ব্যবহারকারী @ চ্যাঙ্কন ইতিমধ্যে কাজটি সম্পন্ন করার পরে, আমি এই নির্দিষ্ট রূপান্তরের সাথে একটি বীজগণিতীয় প্রমাণ সরবরাহ করতে সক্ষম হয়েছি। আমি কোন বিবরণ এড়ানো হয়নি।

( এর ঘনত্বটি বৈধ হওয়ার জন্য আমাদের ইতিমধ্যে )। $n>2$ $Y$

আসুন রূপান্তরটি যেমন এবং । $(X,Y)\mapsto (U,V)$ $U=(2Y-1)\sqrt{X}$ $V=X$

এটি এবং বোঝায় । $x=v$ $y=\frac{1}{2}\left(\frac{u}{\sqrt{v}}+1\right)$

এখন, এবং , $x>0\implies v>0$ $0<y<1\implies-\sqrt{v}<u<\sqrt{v}$

যাতে এর bivariate সমর্থন সহজভাবে হয় । $(U,V)$ $S=\{(u,v):0<u^2<v<\infty,\,u\in\mathbb{R}\}$

রূপান্তরের Jacobian এর পরম মান। $|J|=\frac{1}{2\sqrt{v}}$

এর যৌথ ঘনত্ব এইভাবে $(U,V)$

f_{U, V} (u, v) = \frac{e^{- \frac{v}{2}} v^{\frac{n - 1}{2} - 1} {(\frac{u}{\sqrt{v}} + 1)}^{\frac{n}{2} - 2} {(\frac{1}{2} - \frac{u}{2 \sqrt{v}})}^{\frac{n}{2} - 2} Γ (n - 2)}{(2 \sqrt{v}) 2^{\frac{n - 1}{2} + \frac{n}{2} - 2} Γ (\frac{n - 1}{2}) {(Γ (\frac{n}{2} - 1))}^{2}} 1_{S}

$f_{U,V}(u,v)=\frac{e^{-\frac{v}{2}}v^{\frac{n-1}{2}-1}\left(\frac{u}{\sqrt{v}}+1\right)^{\frac{n}{2}-2}\left(\frac{1}{2}-\frac{u}{2\sqrt{v}}\right)^{\frac{n}{2}-2}\Gamma(n-2)}{(2\sqrt{v})\,2^{\frac{n-1}{2}+\frac{n}{2}-2}\,\Gamma\left(\frac{n-1}{2}\right)\left(\Gamma\left(\frac{n}{2}-1\right)\right)^2}\mathbf1_{S}$

= \frac{e^{- \frac{v}{2}} v^{\frac{n - 4}{2}} (\sqrt{v} + u)^{\frac{n}{2} - 2} (\sqrt{v} - u)^{\frac{n}{2} - 2} Γ (n - 2)}{2^{\frac{2 n - 3}{2} + \frac{n}{2} - 2} (\sqrt{v})^{n - 4} Γ (\frac{n - 1}{2}) {(Γ (\frac{n - 2}{2}))}^{2}} 1_{S}

$=\frac{e^{-\frac{v}{2}}v^{\frac{n-4}{2}}(\sqrt{v}+u)^{\frac{n}{2}-2}(\sqrt{v}-u)^{\frac{n}{2}-2}\,\Gamma(n-2)}{2^{\frac{2n-3}{2}+\frac{n}{2}-2}\,(\sqrt{v})^{n-4}\,\Gamma\left(\frac{n-1}{2}\right)\left(\Gamma\left(\frac{n-2}{2}\right)\right)^2}\mathbf1_{S}$

এখন, লেজেন্ডারের নকল সূত্র ব্যবহার করে,

$\Gamma(n-2)=\frac{2^{n-3}}{\sqrt \pi}\Gamma\left(\frac{n-2}{2}\right)\Gamma\left(\frac{n-2}{2}+\frac{1}{2}\right)=\frac{2^{n-3}}{\sqrt \pi}\Gamma\left(\frac{n-2}{2}\right)\Gamma\left(\frac{n-1}{2}\right)$ যেখানে । $n>2$

সুতরাং , $n>2$

f_{U, V} (u, v) = \frac{2^{n - 3} e^{- \frac{v}{2}} (v - u^{2})^{\frac{n}{2} - 2}}{\sqrt{π} 2^{\frac{3 n - 7}{2}} Γ (\frac{n}{2} - 1)} 1_{S}

$f_{U,V}(u,v)=\frac{2^{n-3}\,e^{-\frac{v}{2}}(v-u^2)^{\frac{n}{2}-2}}{\sqrt \pi\,2^{\frac{3n-7}{2}}\,\Gamma\left(\frac{n}{2}-1\right)}\mathbf1_{S}$

প্রান্তিক পিডিএফ তারপর দেওয়া হয় $U$

f_{U} (u) = \frac{1}{2^{\frac{n - 1}{2}} \sqrt{π} Γ (\frac{n}{2} - 1)} \int_{u^{2}}^{\infty} e^{- \frac{v}{2}} (v - u^{2})^{\frac{n}{2} - 2} d v

$f_U(u)=\frac{1}{2^{\frac{n-1}{2}}\sqrt \pi\,\Gamma\left(\frac{n}{2}-1\right)}\int_{u^2}^\infty e^{-\frac{v}{2}}(v-u^2)^{\frac{n}{2}-2}\,\mathrm{d}v$

= \frac{e^{- \frac{u^{2}}{2}}}{2^{\frac{n - 1}{2}} \sqrt{π} Γ (\frac{n}{2} - 1)} \int_{0}^{\infty} e^{- \frac{t}{2}} t^{(\frac{n}{2} - 1 - 1)} d t

$=\frac{e^{-\frac{u^2}{2}}}{2^{\frac{n-1}{2}}\sqrt \pi\,\Gamma\left(\frac{n}{2}-1\right)}\int_0^\infty e^{-\frac{t}{2}}\,t^{(\frac{n}{2}-1-1)}\,\mathrm{d}t$

= \frac{1}{2^{\frac{n - 1}{2}} \sqrt{π} {(\frac{1}{2})}^{\frac{n}{2} - 1}} e^{- \frac{u^{2}}{2}}

$=\frac{1}{2^{\frac{n-1}{2}}\sqrt \pi\,\left(\frac{1}{2}\right)^{\frac{n}{2}-1}}e^{-\frac{u^2}{2}}$

= \frac{1}{\sqrt{2 π}} e^{- u^{2} / 2}, u \in R

$=\frac{1}{\sqrt{2\pi}}e^{-u^2/2}\,,u\in\mathbb{R}$

— StubbornAtom
সূত্র

2

এই কালো বাক্স উত্তর (অর্থাত, বীজগাণিতিক বিবরণ অনুপস্থিত) ব্যবহার বেশি হয় ম্যাথামেটিকাল । সংক্ষেপে @ শুভর জবাব অনুসারে উত্তরটি হ'ল বিতরণ একটি স্ট্যান্ডার্ড স্বাভাবিক বিতরণ। $Z$

(* Transformation *)
f = {(2 y - 1) Sqrt[x], Sqrt[x]};
sol = Solve[{z == (2 y - 1) Sqrt[x], w == Sqrt[x]}, {x, y}][[1]]
(*{x -> w^2,y -> (w+z)/(2 w)} *)
(* Jacobian *)
J = D[f, {{x, y}}]

(* Joint pdf of Z and W *)
{jointpdf, conditions} = FullSimplify[PDF[BetaDistribution[n/2 - 1, n/2 - 1], y] 
  PDF[ChiSquareDistribution[n - 1], x] Abs[Det[J]] /. sol,
  Assumptions -> {w >= 0, 0 <= y <= 1}][[1, 1]]

(* Integrate over W *)
Integrate[jointpdf Boole[conditions], {w, 0, \[Infinity]}, Assumptions -> {n > 2, z == 0}]
(* 1/Sqrt[2 \[Pi]] *)

Integrate[jointpdf Boole[conditions], {w, 0, \[Infinity]}, Assumptions -> {n > 2, z > 0}]
(* Exp[-(z^2/2)]/Sqrt[2 \[Pi]] *)

Integrate[jointpdf Boole[conditions], {w, 0, \[Infinity]}, Assumptions -> {n > 2, z < 0}]
(* Exp[-(z^2/2)]/Sqrt[2 \[Pi]] *)

— JimB
সূত্র

1

প্রতি সেবার কোনও উত্তর নয় , তবে বক্স-মুলার রূপান্তরের সংযোগটি উল্লেখ করা সার্থক হতে পারে।

বাক্স-মুলার রূপান্তরটি , যেখানে । আমরা এটি দেখাতে পারি , অর্থাৎ । অন্যদিকে, আমরা দেখাতে পারি যে এর লোকেশন-স্কেল আরকসিন বিতরণ রয়েছে , যা এর বিতরণে সম্মত হয় । এর অর্থ বক্স-মুলার রূপান্তর একটি বিশেষ ক্ষেত্রে যখন । $Z=\sqrt{-2\ln U}\sin(2\pi V)$ $U,V\sim U(0,1)$ $-\ln U\sim\operatorname{Exp}(1)$ $-2\ln U\sim\chi^2_2$ $\sin(2\pi V)$ $2\operatorname{B}(1/2,1/2)-1$ $(2Y-1)\sqrt{X}$ $n=3$

সম্পর্কিত :

— ফ্রান্সিস
সূত্র