Pemrograman Matlab dan Macro Minitab Ekspektasi Mean dengan Menggunakan Metode Bootstrap

Kamis, 29 April 2010

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

Share

| More