Misalkan x1, x2, …xn adalah variable random dimana kita akan mencari ekspektasi mean dengan menggunakan metode bootstrap maka hal yang perlu kita lakukan adalah:
Pertama kita harus mengambil sampling dengan pengembalian seukuran n yang dilakukan sebanyak b kali (biasanya b=1000). Setiap pengambilan sampling yang berukuran n, kita harus menghitung mean sehingga kita akan mendapatkan mean sebanyak b. Kemudian mean sebanyak b tersebut kita cari rata-ratanya (mean) dan rata-rata itulah yang dinamakan ekspektasi mean dengan metode bootstrap.
Berikut adalah macro minitab ekspektasi mean dengan metode bootstrap
Macro
bootmean x
mconstant i n b lower upper stat_b se_bmcolumn x y stat
let n=count(x)
let b=1000
do i=1:b
sample n x y;
replacement.
let stat(i)=mean(y)
enddo
let stat_b=mean(stat)
let se_b=stde(stat)
histo stat
sort stat stat
let lower=stat(25)
let upper=stat(975)
print stat_b se_b lower upper
endmacro
Berikut adalah pemrograman matlab ekspektasi mean dengan menggunakan metode bootstrap
function f=bmean(x)
n=length(x);
b=1000;f
or i=1:b
y=randsample(x,n,true);
stat(i)=mean(y);
end
stat_b=mean(stat)
se_b=std(stat)^2
hist(stat)stat=sort(stat);
lower=stat(25)
upper=stat(975)
0 komentar:
Posting Komentar