求如果n--√(X¯--√- 1 )n(X¯−1)n--√(X¯− 1 ) → N( 0 , 1 )n(X¯−1)→N(0,1)

机器算法验证 自习 收敛 渐近的 增量法
2022-04-03 19:04:34

如果你能在下面检查我的工作吗?n(X¯1)n(X¯1)N(0,1)

原则上,这里应该使用 Delta 方法。我们知道,根据该特定方法,如果 then

n(Xnθ)N(0,σ2)
n(g(Xn)g(θ))N(0,σ2g(θ)2)

我们有这意味着然后是极限分布是g(t)=tg(1)=1g(1)=12N(0,1/4)

这里一切都好吗?谢谢。

1个回答

结果是正确的(高达的一个因子,这是一个不重要的印刷遗漏)。这个答案提供了两种单独的方法来仔细检查它。σ

事实上,我们可以直接得到转换后变量的 PDF:当正好是 Normal 时(而不只是渐近如此),可以找到的 PDF通过集成作为Xnn(g(Xn)1)

1(σ/2)2πn(n+x)exp(x2(2n+x)22nσ2)

对于 (否则等于对于固定的极限值x>n0xn

1(σ/2)2πexp(8x22σ2)=1(σ/2)2πexp(x22(σ/2)2),

正态分布的 PDF。(0,σ2/4)


可以检查也可以通过模拟检查结果,例如通过以下R代码执行:

set.seed(17)
n <- 10^6
x <- sqrt(n)*(sqrt(rnorm(n, 1, 1/sqrt(n)))-1)
m <- mean(x); v <- var(x); k <- mean((x-m)^4)
se.v <- sqrt(((n-1)^2 * k - (n-1)*(n-3)*v^2) / n^3)
print(v)              # Variance
print((v - 1/4)/se.v) # Its standardized standard error

这报告了模拟方差(当时)为(单位)标准差个标准误:证明计算的方差对于如此大的并不正确。n=1060.250150.441/4n