拓端tecdat|R語言可視化漸近正態(tài)性、收斂性:大數(shù)定律、中心極限定理、經(jīng)驗累積分布函
原文鏈接: http://tecdat.cn/?p=23777
原文出處:拓端數(shù)據(jù)部落公眾號
在我們的數(shù)理統(tǒng)計課程中,已經(jīng)看到了大數(shù)定律(這在概率課程中已經(jīng)被證明),證明
?

給出一組i.i.d.隨機變量?

?,其中有

為了直觀地看到這種收斂性,我們可以使用
> for(i in 1:20)B[,i]=mean_samples(i*10)
> boxplot(B)
也可以直觀地看到邊界?

?(用于中心極限定理,獲得極限的非退化分布)。

我們一直在討論經(jīng)驗累積分布函數(shù)的特點。

我們已經(jīng)看到了格利文科-坎特利定理,該定理指出

為了直觀地看到這種收斂。這里我使用了一個技巧可視化

獲得兩個矩陣之間的最大值(分量)。
+ Df=(D1+D2)/2+abs(D2-D1)/2
> boxplot(B)

我們還討論了經(jīng)驗累積分布函數(shù)的逐點漸近正態(tài)性

在這里,又可以把它形象化。第一步是計算經(jīng)驗累積分布函數(shù)的幾條軌跡
> plot(u,u)
請注意,我們可以計算(逐點)置信帶
> lines(u,apply(M,1,function(x) quantile(x,.05)
> lines(u,apply(M,1,function(x) quantile(x,.95)

現(xiàn)在,如果我們專注于一個特定的點,我們可以直觀地看到漸近正態(tài)性(即當我們有一個大小為100的樣本時,幾乎是正態(tài)的)。
> hist(y)
> lines(vu,dnorm(vu,pnorm(x0)
+ sqrt((pnorm(x0)*(1-pnorm(x0)))/100)


最受歡迎的見解
1.Matlab馬爾可夫鏈蒙特卡羅法(MCMC)估計隨機波動率(SV,Stochastic Volatility) 模型
2.基于R語言的疾病制圖中自適應(yīng)核密度估計的閾值選擇方法
3.WinBUGS對多元隨機波動率模型:貝葉斯估計與模型比較
4.R語言回歸中的hosmer-lemeshow擬合優(yōu)度檢驗
5.matlab實現(xiàn)MCMC的馬爾可夫切換ARMA – GARCH模型估計
6.R語言區(qū)間數(shù)據(jù)回歸分析
7.R語言WALD檢驗 VS 似然比檢驗
8.python用線性回歸預(yù)測股票價格
9.
R語言如何在生存分析與Cox回歸中計算IDI,NRI指標