Bu yazı, soruya yapılan yorumların cevaplarını detaylandırmaktadır.
Let . Herhangi bir saptamak E 1 ∈ R, n birim uzunluğu. Böyle bir vektör, her zaman ortonormal esasına tamamlanabilir ( e 1 , e 2 , ... , e n ) (vasıtasıyla Gram-Schmidt işlemi , örneğin). Bu temel değişim (olağan olandan) diktir: uzunlukları değiştirmez. Böylece dağılımıX= ( X1, X2, … , Xn)e1∈ Rn( e1, e2, … , En)
( e1⋅ X)2| | X| |2= ( e1⋅ X)2X21+ X22+ ⋯ + X2n
bağlı değildir . Alarak e 1 = ( 1 , 0 , 0 , ... , 0 ) gösterir, bu aynı dağılımı olarak vare1e1= ( 1 , 0 , 0 , … , 0 )
X21X21+ X22+ ⋯ + X2n.(1)
Yana olarak Normal IID edilir, bunlar yazılmış olabilir σ kez standart normal değişkenler IID Y 1 , ... , Y n ve onların kareler vardır σ 2 kez Γ ( 1 / 2 ) dağılımları. Yana toplamı n - 1 bağımsız Γ ( 1 / 2 ) dağılımları olan Γ ( ( N - 1 ) / 2 )XbenσY1, … , Ynσ2Γ ( 1 / 2 )n - 1Γ ( 1 / 2 )Γ ( ( n - 1 ) / 2 ), in dağılımının( 1 )
σ2Uσ2U+ σ2V= UU+ V
burada ve V = ( x 2 2 + ⋯ + x 2 , n ) / σ 2 ~ Γ ( ( n - 1 ) / 2 ) bağımsızdır. Olduğu iyi bilinen bu oran bir Beta sahip olduğu ( 1 / 2 , ( n - 1U= X21/ σ2~ Γ ( 1 / 2 )V= ( X22+ ⋯ + X2n) / σ2∼ Γ ( ( n - 1 ) / 2 ) dağıtım. (Aynı zamanda yakından ilgili iplik bakınızDağıtım X , Y ise X ~ Beta ( 1 , K - 1 ) ve Y ~ ki-kare ile 2 K derece).( 1 / 2 , ( n - 1 ) / 2 )XYX~( 1 , K- 1 )Y~2 K
Yana
X1+ ⋯ + Xn= ( 1 , 1 , … , 1 ) ⋅ ( X1, X2, ⋯ , Xn) = n--√e1⋅ X
birim vektör , biz sonucunaZ'ninise( √e1= ( 1 , 1 , … , 1 ) / n--√Zkez beta(1/2,(n-1)/2)değişken. ( n--√)2= n( 1 / 2 , ( n - 1 ) / 2 ) İçinbu nedenle yoğunluğu fonksiyonu vardırn ≥ 2
fZ( z) = n1 - n / 2B ( 12, n - 12)( n - z)n−3z−−−−−−−−−√
aralığında ( ve aksi halde sıfır).( 0 , n )
Bir kontrol olarak, σ = 1 ve n = 2 , 3 , 10 için bağımsız Z gerçekleştirmesi simüle ettim , histogramlarını çizdim ve karşılık gelen Beta yoğunluğunun grafiğini (kırmızı olarak) üst üste koydum. Anlaşmalar mükemmel.100 , 000Zσ= 1n = 2 , 3 , 10
İşte R
kod. Bu formül aracılığıyla simülasyonu gerçekleştirir sum(x)^2 / sum(x^2)
için , uzunluğunun bir vektör tarafından oluşturulan . Gerisi sadece döngü ( , ) ve çizim ( , ).Zx
n
rnorm
for
apply
hist
curve
for (n in c(2, 3, 10)) {
z <- apply(matrix(rnorm(n*1e5), nrow=n), 2, function(x) sum(x)^2 / sum(x^2))
hist(z, freq=FALSE, breaks=seq(0, n, length.out=50), main=paste("n =", n), xlab="Z")
curve(dbeta(x/n, 1/2, (n-1)/2)/n, add=TRUE, col="Red", lwd=2)
}