ไฟล์ PDF ของ

สมมติว่าเป็น iid จากโดยไม่ทราบและ $X_1, X_2,...,X_n$ $N(\mu,\sigma^2)$ $\mu \in \mathcal R$ $\sigma^2>0$

ให้ S คือค่าเบี่ยงเบนมาตรฐานที่นี่ $Z=\frac{X_1-\bar{X}}{S},$

มันสามารถแสดงให้เห็นว่า มีไฟล์ Lebesgue $Z$

f (z) = \frac{\sqrt{n} Γ (\frac{n - 1}{2})}{\sqrt{π} (n - 1) Γ (\frac{n - 2}{2})} {[1 - \frac{n z^{2}}{(n - 1)^{2}}]}^{n / 2 - 2} I_{(0, (n - 1) / \sqrt{n})} (| Z |)

$f(z)=\frac{\sqrt{n} \Gamma\left(\frac{n-1}{2}\right)}{\sqrt{\pi}(n-1)\Gamma\left(\frac{n-2}{2}\right)}\left[1-\frac{nz^2}{(n-1)^2}\right]^{n/2-2}I_{(0,(n-1)/\sqrt{n})}(|Z|)$

คำถามของฉันคือวิธีการรับ pdf นี้

คำถามคือจากที่นี่ในตัวอย่าง 3.3.4 เพื่อหา UMVUE ของ $P(X_1 \le c)$ ค) ฉันเข้าใจตรรกะและขั้นตอนเพื่อค้นหา UMVUE แต่ไม่รู้วิธีรับ PDF

ผมคิดว่าคำถามนี้ยังเกี่ยวข้องกับเรื่องนี้อย่างใดอย่างหนึ่ง

ขอบคุณมากสำหรับความช่วยเหลือหรือชี้ไปที่การอ้างอิงใด ๆ ที่เกี่ยวข้องจะได้รับการจัดสรร

self-study umvue

— ภาคเหนือตอนล่าง
แหล่งที่มา

คำตอบ:

สิ่งที่น่าสนใจมากเกี่ยวกับผลลัพธ์นี้คือดูเหมือนว่าการกระจายของสัมประสิทธิ์สหสัมพันธ์ มีเหตุผล

สมมติว่าเป็นค่าปกติ bivariate ที่ไม่มีสหสัมพันธ์และความแปรปรวนทั่วไปสำหรับทั้งสองตัวแปร วาดตัวอย่าง IIDy_n) มันเป็นที่รู้จักกันดีและสร้างขึ้นอย่างง่ายดายในเชิงเรขาคณิต (เช่นเดียวกับฟิชเชอร์เมื่อหนึ่งศตวรรษก่อน) ว่าการกระจายตัวของสัมประสิทธิ์สหสัมพันธ์ตัวอย่าง $(X,Y)$ $\sigma^2$ $(x_1,y_1), \ldots, (x_n,y_n)$

r = \frac{\sum_{i = 1}^{n} (x_{i} - \bar{x}) (y_{i} - \bar{y})}{(n - 1) S_{x} S_{y}}

$r = \frac{\sum_{i=1}^n(x_i - \bar x)(y_i - \bar y)}{(n-1) S_x S_y}$

คือ

f (r) = \frac{1}{B (\frac{1}{2}, \frac{n}{2} - 1)} {(1 - r^{2})}^{n / 2 - 2}, - 1 \leq r \leq 1.

$f(r) = \frac{1}{B\left(\frac{1}{2}, \frac{n}{2}-1\right)}\left(1-r^2\right)^{n/2-2},\ -1 \le r \le 1.$

(ที่นี่ตามปกติและเป็นค่าเฉลี่ยตัวอย่างและและเป็นรากที่สองของการประมาณค่าความแปรปรวนแบบไม่เอนเอียง) คือฟังก์ชันเบต้าซึ่ง $\bar x$ $\bar y$ $S_x$ $S_y$ $B$

\begin{matrix} (1) & \frac{1}{B (\frac{1}{2}, \frac{n}{2} - 1)} = \frac{Γ (\frac{n - 1}{2})}{Γ (\frac{1}{2}) Γ (\frac{n}{2} - 1)} = \frac{Γ (\frac{n - 1}{2})}{\sqrt{π} Γ (\frac{n}{2} - 1)} . \end{matrix}

$\frac{1}{B\left(\frac{1}{2}, \frac{n}{2}-1\right)} = \frac{\Gamma\left(\frac{n-1}{2}\right)}{\Gamma\left(\frac{1}{2}\right)\Gamma\left(\frac{n}{2}-1\right)} = \frac{\Gamma\left(\frac{n-1}{2}\right)}{\sqrt{\pi}\Gamma\left(\frac{n}{2}-1\right)} . \tag{1}$

ในการคำนวณเราอาจใช้ความไม่แปรเปลี่ยนของมันภายใต้การหมุนในรอบเส้นที่สร้างโดยพร้อมกับความแปรปรวนของการกระจายตัวของตัวอย่างภายใต้การหมุนเดียวกัน และเลือกเป็นเวกเตอร์หน่วยใดก็ได้ที่มีผลรวมองค์ประกอบเป็นศูนย์ หนึ่งเวกเตอร์ดังกล่าวเป็นสัดส่วนกับ-1) ค่าเบี่ยงเบนมาตรฐานคือ $r$ $\mathbb{R}^n$ $(1,1,\ldots, 1)$ $y_i/S_y$ $v = (n-1, -1, \ldots, -1)$

S_{v} = \sqrt{\frac{1}{n - 1} ((n - 1)^{2} + (- 1)^{2} + \dots + (- 1)^{2})} = \sqrt{n} .

$S_v = \sqrt{\frac{1}{n-1}\left((n-1)^2 + (-1)^2 + \cdots + (-1)^2\right)} = \sqrt{n}.$

ดังนั้นต้องมีการแจกแจงแบบเดียวกับ $r$

\frac{\sum_{i = 1}^{n} (x_{i} - \bar{x}) (v_{i} - \bar{v})}{(n - 1) S_{x} S_{v}} = \frac{(n - 1) x_{1} - x_{2} - \dots - x_{n}}{(n - 1) S_{x} \sqrt{n}} = \frac{n (x_{1} - \bar{x})}{(n - 1) S_{x} \sqrt{n}} = \frac{\sqrt{n}}{n - 1} Z .

$\frac{\sum_{i=1}^n(x_i - \bar x)(v_i - \bar v)}{(n-1) S_x S_v} = \frac{(n-1)x_1 - x_2-\cdots-x_n}{(n-1) S_x \sqrt{n}} = \frac{n(x_1 - \bar x)}{(n-1) S_x \sqrt{n}} = \frac{\sqrt{n}}{n-1}Z.$

ดังนั้นสิ่งที่เราต้องทำคือ rescaleเพื่อค้นหาการกระจายตัวของ : $r$ $Z$

f_{Z} (z) = | \frac{\sqrt{n}}{n - 1} | f (\frac{\sqrt{n}}{n - 1} z) = \frac{1}{B (\frac{1}{2}, \frac{n}{2} - 1)} \frac{\sqrt{n}}{n - 1} {(1 - \frac{n}{(n - 1)^{2}} z^{2})}^{n / 2 - 2}

$f_Z(z) = \big|\frac{\sqrt{n}}{n-1}\big| f\left(\frac{\sqrt{n}}{n-1}z\right) = \frac{1}{B\left(\frac{1}{2}, \frac{n}{2}-1\right)} \frac{\sqrt{n}}{n-1}\left(1- \frac{n}{(n-1)^2}z^2\right)^{n/2-2}$

สำหรับ{n}} สูตร (1) แสดงสิ่งนี้เหมือนกับคำถาม $|z| \le \frac{n-1}{\sqrt{n}}$

ไม่มั่นใจอย่างสมบูรณ์? นี่คือผลลัพธ์ของการจำลองสถานการณ์นี้ 100,000 ครั้ง (ด้วยที่การแจกแจงสม่ำเสมอ) $n=4$

ครั้งแรกที่แปลงกราฟค่าสัมประสิทธิ์สหสัมพันธ์ของในขณะที่แปลง histogram สองค่าสัมประสิทธิ์สหสัมพันธ์ของสำหรับการสุ่มเลือก vectorที่คงที่สำหรับการวนซ้ำทั้งหมด พวกเขาทั้งสองเหมือนกัน พล็อต QQ- ทางด้านขวาเป็นการยืนยันว่าการแจกแจงเหล่านี้เหมือนกัน $(x_i,y_i),i=1,\ldots,4$ $(x_i,v_i),i=1,\ldots,4)$ $v_i$

นี่คือRรหัสที่สร้างพล็อต

n <- 4
n.sim <- 1e5
set.seed(17)
par(mfrow=c(1,3))
#
# Simulate spherical bivariate normal samples of size n each.
#
x <- matrix(rnorm(n.sim*n), n)
y <- matrix(rnorm(n.sim*n), n)
#
# Look at the distribution of the correlation of `x` and `y`.
#
sim <- sapply(1:n.sim, function(i) cor(x[,i], y[,i]))
hist(sim)
#
# Specify *any* fixed vector in place of `y`.
#
v <- c(n-1, rep(-1, n-1)) # The case in question
v <- rnorm(n)             # Can use anything you want
#
# Look at the distribution of the correlation of `x` with `v`.
#
sim2 <- sapply(1:n.sim, function(i) cor(x[,i], v))
hist(sim2)
#
# Compare the two distributions.
#
qqplot(sim, sim2, main="QQ Plot")

การอ้างอิง

ชาวประมง RA, ความถี่การกระจายของค่าสัมประสิทธิ์สหสัมพันธ์ในตัวอย่างจากประชากรที่มีขนาดใหญ่ไปเรื่อย ๆ Biometrika , 10 , 507 ดูหัวข้อ 3 (อ้างถึงในทฤษฎีขั้นสูงทางสถิติของเคนดัลล์ , 5 เอ็ด, มาตรา 16.24)

— whuber
แหล่งที่มา

ลิงก์ไปยังข้อมูลอ้างอิงใช้งานไม่ได้

— Sextus Empiricus

@Martijn ขอบคุณสำหรับการตรวจสอบ ฉันเห็นสิ่งที่คุณหมายถึง - การเชื่อมโยงใช้งานได้ แต่มันไม่ได้เกี่ยวข้องอะไรเลย! ฉันได้แก้ไขแล้ว

— whuber

ฉันอยากจะแนะนำวิธีนี้ในการรับไฟล์ PDF ของ Z โดยการคำนวณ MVUE ของโดยตรงโดยใช้ทฤษฎีบทของเบย์แม้ว่ามันจะมีจำนวนไม่มากและซับซ้อน $P(X\leq c)$

เนื่องจากและ ,เป็นสถิติที่สมบูรณ์ที่สมบูรณ์ร่วม MVUE ของ จะเป็นเช่นนี้: $E[I_{(-\infty,c)}(X_1)]=P(X_1\leq c)$ $Z_1=\bar X$ $Z_2=S^2$ $P(X\leq c)$

ψ (z_{1}, z_{2}) = E [I_{(- \infty, c)} (X_{1}) | z_{1}, z_{2}] = \int_{- \infty}^{\infty} I_{(- \infty, c)} f_{X | Z_{1}, Z_{2}} (x_{1} | z_{1}, z_{2}) d x_{1}

$\psi(z_1,z_2)=E[I_{(-\infty,c)}(X_1)|z_1,z_2]=\int_{-\infty}^{\infty}I_{(-\infty,c)}f_{X|Z_1,Z_2}(x_1|z_1,z_2)dx_1$

ตอนนี้ใช้ทฤษฎีบทของเบย์เราจะได้

f_{X | Z_{1}, Z_{2}} (x_{1} | z_{1}, z_{2}) = \frac{f_{Z_{1}, Z_{2} | X_{1}} (z_{1}, z_{2} | x_{1}) f_{X_{1}} (x_{1})}{f_{Z_{1}, Z_{2}} (z_{1}, z_{2})}

$f_{X|Z_1,Z_2}(x_1|z_1,z_2)={{f_{Z_1,Z_2|X_1}(z_1,z_2|x_1)f_{X_1}(x_1)}\over{f_{Z_1,Z_2}(z_1,z_2)}}$

ตัวส่วนสามารถเขียนในรูปแบบปิดเนื่องจาก ,เป็นอิสระต่อกัน $f_{Z_1,Z_2}(z_1,z_2)=f_{Z_1}(z_1)f_{Z_2}(z_2)$ $Z_1 \sim N(\mu,\frac{\sigma^2}{n})$ $Z_2 \sim \Gamma({n-1\over 2},{2 \sigma^2\over n-1})$

เพื่อให้ได้รูปแบบตัวเลขที่ปิดเราสามารถนำสถิติเหล่านี้มาใช้:

W_{1} = \frac{\sum_{i = 2}^{n} X_{i}}{n - 1}

$W_1 = {\sum_{i=2}^n X_i \over n-1}$

W_{2} = \frac{\sum_{i = 2}^{n} X_{i}^{2} - (n - 1) W_{1}^{2}}{(n - 1) - 1}

$W_2 = {\sum_{i=2}^n X_i^2 -(n-1) W_1^2 \over (n-1)-1}$

ซึ่งเป็นค่าเฉลี่ยและความแปรปรวนของตัวอย่างและพวกเขามีความเป็นอิสระของแต่ละอื่น ๆ และยังเป็นอิสระจากX_1เราสามารถแสดงเหล่านี้ในแง่ของการZ_2 $X_2, X_3, ..., X_n$ $X_1$ $Z_1, Z_2$

$W_1={n Z_1 - X_1\over n-1}$ , $W_2={(n-1)Z_2+nZ_1^2-X_1^2-(n-1)W_1^2 \over n-2}$

เราสามารถใช้การแปลงในขณะที่ , $X_1=x_1$

f_{Z_{1}, Z_{2} | X_{1}} (z_{1}, z_{2} | x_{1}) = \frac{n}{n - 2} f_{W_{1}, W_{2}} (w_{1}, w_{2}) = \frac{n}{n - 2} f_{W_{1}} (w_{1}) f_{W_{2}} (w_{2})

$f_{Z_1,Z_2|X_1}(z_1,z_2|x_1)={n \over n-2}f_{W_1,W_2}(w_1,w_2)={n \over n-2}f_{W_1}(w_1)f_{W_2}(w_2)$

ตั้งแต่ ,เราสามารถรับรูปแบบปิดนี้ โปรดทราบว่าสิ่งนี้มีไว้สำหรับซึ่ง จำกัดถึงz_2} $W_1 \sim N(\mu,\frac{\sigma^2}{n-1})$ $W_2 \sim \Gamma({n-2\over 2},{2 \sigma^2\over n-2})$ $w_2 \geq 0$ $x_1$ $z_1-{n-1 \over \sqrt n}\sqrt{z_2} \leq x_1 \leq z_1+{n-1 \over \sqrt n}\sqrt{z_2}$

ดังนั้นให้รวมเข้าด้วยกันเทอมชี้แจงจะหายไปและคุณจะได้

f_{X | Z_{1}, Z_{2}} (x_{1} | z_{1}, z_{2}) = \frac{Γ (\frac{n - 1}{2})}{\sqrt{π} Γ (\frac{n - 2}{2})} \frac{\sqrt{n}}{\sqrt{z_{2}} (n - 1)} (1 - {(\frac{\sqrt{n} (x_{1} - z_{1})}{\sqrt{z_{2}} (n - 1)})}^{2})

$f_{X|Z_1,Z_2}(x_1|z_1,z_2)={\Gamma({n-1 \over 2}) \over \sqrt{\pi} \Gamma({n-2 \over 2})} {\sqrt{n} \over \sqrt{z_2} (n-1)} (1-{({\sqrt{n} (x_1 -z_1) \over \sqrt{z_2} (n-1) })}^2)$ โดยที่และศูนย์อื่น ๆ

z_{1} - \frac{n - 1}{\sqrt{n}} \sqrt{z_{2}} \leq x_{1} \leq z_{1} + \frac{n - 1}{\sqrt{n}} \sqrt{z_{2}}

$z_1-{n-1 \over \sqrt n}\sqrt{z_2} \leq x_1 \leq z_1+{n-1 \over \sqrt n}\sqrt{z_2}$

จากจุดนี้ ณ จุดนี้เราสามารถรับไฟล์ PDF ของโดยใช้การแปลง $Z={X_1- z_1 \over \sqrt{z_2}}$

โดยวิธีการ MVUE จะเป็นเช่นนี้: ในขณะที่และจะเป็น 1 ถ้า

ψ (z_{1}, z_{2}) = \frac{Γ (\frac{n - 1}{2})}{\sqrt{π} Γ (\frac{n - 2}{2})} \int_{- \frac{π}{2}}^{θ_{c}} c o s^{n - 3} θ d θ

$\psi(z_1,z_2)={\Gamma({n-1 \over 2}) \over \sqrt{\pi} \Gamma({n-2 \over 2})} \int ^{\theta_c} _{-{\pi \over2}} cos^{n-3} \theta d\theta$

θ_{c} = s i n^{- 1} (\frac{\sqrt{n} (c - z_{1})}{(n - 1) \sqrt{z_{1}}})

$\theta_c = sin^{-1} ({\sqrt{n}(c-z_1)\over(n-1)\sqrt{z_1}})$

c \geq z_{1} + \frac{n - 1}{\sqrt{n} \sqrt{z_{2}}}

$c \geq z_1+{n-1 \over \sqrt{n} \sqrt{z_2} }$

ฉันไม่ใช่คนพูดภาษาอังกฤษและอาจมีประโยคที่น่าอึดอัดใจ ฉันกำลังศึกษาสถิติด้วยตัวเองด้วยการแนะนำหนังสือตำราสถิติทางคณิตศาสตร์โดย Hogg ดังนั้นอาจมีข้อผิดพลาดทางแนวคิดทางไวยากรณ์หรือคณิตศาสตร์ มันจะได้รับการชื่นชมถ้ามีคนแก้ไขพวกเขา

ขอบคุณสำหรับการอ่าน.

— KDG
แหล่งที่มา