Catatan Kuliah MA3081 STATISTIKA MATEMATIKA Statistika Mengalahkan Matematika

disusun oleh

Khreshna I.A. Syuhada, MSc. PhD.

Kelompok Keilmuan STATISTIKA - FMIPA Institut Teknologi Bandung 2011

Daftar Isi
1 Peubah Acak dan Distribusi Kontinu 1.1 Fungsi distribusi . . . . . . . . . . . 1.2 Unsur Peluang . . . . . . . . . . . . 1.3 Ekspektasi . . . . . . . . . . . . . . . 1.4 Distribusi Bivariat . . . . . . . . . . 1.5 Distribusi Bersyarat . . . . . . . . . 1.6 Fungsi Pembangkit Momen . . . . . 1 1 4 7 9 11 14 1 1 1 5 6 6 7 7 1 1 2

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

. . . . . .

2 Distribusi Sampel, Likelihood dan Penaksir 2.1 Sampel Acak . . . . . . . . . . . . . . . . . 2.2 Likelihood . . . . . . . . . . . . . . . . . . . 2.3 Statistic Cukup . . . . . . . . . . . . . . . . 2.4 Distribusi Sampel . . . . . . . . . . . . . . . 2.5 Statistik Terurut . . . . . . . . . . . . . . . 2.6 Momen dari Mean dan Proporsi Sampel . . 2.7 Teorema Limit Pusat . . . . . . . . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

. . . . . . .

3 Penaksiran dan Selang Kepercayaan 3.1 Sifat-sifat (Kesalahan) penaksiran . . . . . . . . . . . . . . . . 3.2 Konsistensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

BAB 1 Peubah Acak dan Distribusi Kontinu

Definisi: Misalkan X peubah acak. Fungsi distribusi (kumulatif) dari X adalah FX (x) = P (X x) Contoh: 1. Misalkan X Bin(3, 0.5), maka fungsi distribusi F (x) adalah fungsi tangga berikut 0, 1/8, F (x) = 1/2, 7/8, 1, x (, 0); x [0, 1); x [1, 2); x [2, 3); x [3, ).

2. Misalkan X peubah acak dengan support S = [a, b], b > 0. Misalkan peluang X akan berada di selang S proporsional terhadap panjang selang. Dengan kata lain, P (x1 X x2 ) = (x2 x1 ),

untuk a x1 x2 b. Untuk menentukan , misalkan x1 = a dan x2 = b. Maka, P (a X b) = 1 = (b a) = 1/(b a) Fungsi distribusinya: 0, x < a; x [a, b]; x > b.

F (x) = P (X x) = P (a X x) =

x a , ba

1,

Peubah acak X dikatakan berdistribusi Uniform, X U (a, b). Sifat-sifat fungsi distribusi: F () = 0 dan F () = 1 F merupakan fungsi tidak turun; F (a) F (b) untuk a b F adalah fungsi kontinu kanan; lim0+ F (x + ) = F (x) Misalkan X peubah acak dengan fungsi distribusi F (x). Jika b a, maka P (a < X b) = F (b) F (a) Untuk setiap x, P (X = x) = lim0+ P (x < X ) = F (x) F (x) (Perhatikan notasi F (x) dan kasus apabila fungsi distribusi kontinu kiri) Definisi: Distribusi dari peubah acak X dikatakan KONTINU jika fungsi distribusi disetiap x kontinu dan fungsi distribusi tersebut dapat diturunkan. Misalkan X peubah acak kontinu dengan fungsi distribusi FX (x). Misalkan g (X ) fungsi naik satu-satu kontinu. Untuk y yang berada di daerah hasil dari g , fungsi invers x = g 1 (y ) ada. Misalkan Y = g (X ). Fungsi distribusi dari Y adalah P (Y y ) = P (g (X ) y ) = P (X g 1 (y )) = FX (g 1 (y )) Misalkan g (X ) fungsi turun satu-satu kontinu. Untuk y yang berada di daerah hasil dari g , fungsi invers x = g 1 (y ) ada. Misalkan Y = g (X ). Fungsi distribusi dari Y adalah P (Y y ) = P (g (X ) y ) = P (X > g 1 (y )) = 1 FX (g 1 (y )) MA3081 Stat.Mat. 2 K. Syuhada, PhD.

 Misalkan X U (0, 1) dan Y = g (X ) = hX + k, h < 0. Maka X = g 1 (Y ) = (Y k )/h FX (x) = FY (y ) = Y U (h + k, k ) Latihan: 1. Misalkan X peubah acak kontinu yang memiliki fungsi distribusi FX (x) yang naik murni. Misalkan Y = FX (X ). Tentukan distribusi dari Y 2. Misalkan U peubah acak berdistribusi U (0, 1). Misalkan FX (x) fungsi distribusi yang naik murni dari X . Tentukan fungsi distribusi dari peubah 1 acak FX (U ) 3. Misalkan U1 , U2 , . . . , Un sampel acak dari U (0, 1). Bangkitkan sampel acak dari FX (x) (ambil contoh misalnya untuk FX (x) = 1 e x , x > 0) Misalkan X peubah acak kontinu dengan fungsi distribusi FX (x). Misalkan Y = g (X ) fungsi kontinu tidak monoton. Kita ketahui bahwa pada fungsi yang monoton, FY (y ) = P (Y y ) = P (g (X ) y ) dimana dalam hal ini setiap solusi inverse x = g 1 (y ) digunakan untuk menentukan FY (y ) dengan menggunakan FX (g 1 (y )). Untuk X U (1, 2) dan g (X ) = Y = X 2 , kita dapatkan fungsi distribusi dari Y : FY (y ) =

MA3081 Stat.Mat.

K. Syuhada, PhD.

1.2

Unsur Peluang

Misalkan X peubah acak kontinu, x bilangan positif kecil. Denisikan h(a, b) =def P (a X a + b) = FX (a + b) FX (a) Untuk h(x, x) = P (x X x + x), maka deret Taylor-nya disekitar x = 0 adalah h(x, x) = F (x + x) F (x) d h(x, x) = h(x, 0) + dx = = dimana
x0

x=0

x + o(x)

lim

o(x) =0 x

Fungsi ] d F (x) x dF (x) = dx [ disebut DIFERENSIAL. Dalam statistika, diferensial dari fungsi distribusi adalah UNSUR PELUANG (yang merupakan pendekatan terhadap h(x, x)). d Unsur peluang adalah fungsi linier dari dx F (x). Contoh: Misalkan F (x) = 1 e3x untuk x 0. Apakah F (x) suatu fungsi distribusi? Hitung unsur peluang di x = 2. Cari pendekatan untuk P (2 X 2.01). Densitas rata-rata pada selang (x, x + x) didenisikan: Density rata-rata =def P (x X x + x) x

Sedangkan fungsi densitas peluang atau fungsi peluang (f.p) di x adalah limit

MA3081 Stat.Mat.

K. Syuhada, PhD.

densitas rata-rata saat x 0: f.p = f (x) =def lim = = = d F (x) dx

P (x X x + x) x

Catatan: Unsur peluang dituliskan sebagai dF (x) = f (x)x. Sifat-sifat fungsi peluang: f (x) 0 untuk semua x f (x) = 1 Hubungan antara fungsi peluang dan fungsi distribusi: d F (x) dx x F ( x) = f (u)du f (x) =

P (a < X < b) = ... = ... = ... = F (b) F (a) = Latihan:

f (x)dx

1. Misalkan bilangan riil positif. Jika F (x) = 1 ex , maka f (x) = 2. Jika X U (a, b) maka F (x) = dan f (x) = 3. *Misalkan f (x) = c/(1 + x2 ) untuk < x < dan c konstanta. Fungsi f (x) tak negatif dan (1 + x2 )1 dx = . Berapa nilai c agar f (x) menjadi fungsi peluag? Tentukan fungsi distribusinya. 4. *Pandang distribusi waktu tunggu. Misalkan T adalah waktu kedatangan kejadian ke-r dalam Proses Poisson dengan laju . Tentukan fungsi peluang dari T

MA3081 Stat.Mat.

K. Syuhada, PhD.

Misalkan X peubah acak kontinu dengan fungsi peluang f (x) dan Y = g (X ) fungsi yang terdiferensial bernilai tunggal. Maka fungsi peluang dari Y : fY (y ) = fX (g 1 (y )) d 1 g (y ) dy

untuk support Y = g (X ). Komponen J (y ) = d 1 g (y ) dy

adalah transformasi Jacobian. BUKTI: Misalkan g (x) memiliki lebih dari satu fungsi invers maka unsur peluang yang terpisah harus dihitung untuk setiap fungsi invers. Contoh, misalkan X U (1, 2) dan Y = g (X ) = X 2 . Maka untuk y [0, 1], terdapat 2 fungsi invers yaitu ?, dan satu fungsi invers untuk y (1, 4] yaitu ?. Fungsi peluang dari Y adalah: f (y ) =

MA3081 Stat.Mat.

K. Syuhada, PhD.

1.3

Ekspektasi

Misalkan X peubah acak dengan fungsi peluang f (x). Nilai harapan dari X , jika ada, adalah E (X ) = X = f (x)dx

Catatan: nilai ekspektasi dikatakan ada jika nilai integral adalah hingga. Misalkan X rv dengan pdf f (x). Maka nilai harapan dari g (X ), jika ada, adalah E [g (X )] = g (x)f (x)dx

. Operator integral bersifat linier. Jika g1 (X ) dan g2 (X ) fungsi-fungsi yang memiliki ekspektasi dan a, b, c konstanta, maka E [ag1 (X ) + bg2 (X ) + c] = aE [g1 (X )] + bE [g2 (X )] + c Contoh/Latihan: 1. Jika distribusi X simetrik di sekitar c dan nilai harapanny ada maka E (X ) = c. Bukti: E (X c) = =

(x c)f (x) dx (x c)f (x)dx +

(x c)f (x)dx

uf (c u)du + uf (c + u)du = 0 0 = u(f (c + u) f (c u)) du = 0

2. Misalkan X U (a, b). Tunjukkan bahwa distribusi tersebut simetrik disekitar (a + b)/2.

MA3081 Stat.Mat.

K. Syuhada, PhD.

Bukti:

) ( ) a+b a+b 1 f =f + = 2 2 ba [ ] a ba untuk b , 2 2 3. Misalkan X berdistribusi Cauchy dengan fungsi peluang f (x) = [ 1
(x)2 2

],

1 +

dengan , konstanta yang memenuhi || < dan (0, ). Tunjukkan bahwa fungsi peluang simetrik di sekitar namun ekspektasinya bukanlah . 4. Misalkan X Exp(). Nilai harapan dari X adalah...

MA3081 Stat.Mat.

K. Syuhada, PhD.

1.4

Distribusi Bivariat

Suatu fungsi fX,Y (x, y ) dikatakan fungsi peluang bivariat jika fX,Y (x, y ) 0, untuk semua x, y fX,Y (x, y ) dxdy = 1 Jika fX,Y (x, y ) fungsi peluang bivariat maka x FX,Y (x, y ) = P (X x, Y y ) =

fX,Y (u, v ) dvdu

Sifat-sifat fungsi distribusi bivariat: 1. FX,Y (x, ) = FX (x) 2. FX,Y (, y ) = FY (y ) 3. FX,Y (, ) = 1 4. FX,Y (, y ) = FX,Y (x, ) = FX,Y (, ) = 0 5. fX,Y (x, y ) =
2 xy

FX,Y (x, y )

fX,Y (x, y )xy adalah unsur peluang bersama, P (x X x + x, y Y y + y ) = fX,Y (x, y )xy + o(xy ) Contoh/Latihan: 1. Jika (X, Y ) U (a, b, c, d) maka fX,Y (x, y ) = 1 , x (a, b), y (c, d) (b a)(d c)

2. Untuk soal no 1 di atas, misalkan a = c = 0, b = 4, d = 6 maka P (2.5 X 3.5, 1 Y 4) = 3/24 P (X 2 + Y 2 > 16) = 1 P (X 2 + Y 2 16) = 1 /6

MA3081 Stat.Mat.

K. Syuhada, PhD.

3. Jika fX,Y (x, y ) = (6/5) (x + y 2 ) untuk x (0, 1) dan y (0, 1). Tentukan P (X + Y < 1). 1 1y P (X + y < 1) = P (X < 1 Y ) = fX,Y (x, y ) dx dy
0 0

= = 3/10

Untuk menentukan fungsi peluang marginal, integralkan peubah yang tidak diinginkan: fX (x) = fX,Y (x, y ) dy fY (y ) =

fX,Y (x, y ) dx

fX,Y (x, y ) =

fW,X,Y,Z (w, x, y, z ) dwdz

Pada fungsi peluang fX,Y (x, y ) = 6/5(x + y 2 ) diperoleh fX (x) = fY (y ) = 6x + 2 , x (0, 1) 5

6 y2 + 3 , y (0, 1) 5 dan nilai harapan E (g (X, Y )) = E (X ) =

g (x, y ) fX,Y (x, y ) dxdy = = 3/5

MA3081 Stat.Mat.

10

K. Syuhada, PhD.

1.5

Distribusi Bersyarat

Misalkan fX,Y (x, y ) adalah fungsi peluang bersama, maka fungsi peluang Y , diberikan X = x, adalah fY |X (y |x) =def asalkan fX (x) > 0. Contoh: Misalkan X dan Y memiliki distribusi bersama fX,Y (x, y ) = 8xy, 0 < x < y < 1, maka fX (x) = 4 x 4 x3 , 0 < x < 1 E (X r ) = 8 (r + 2)(r + 4) fX,Y (x, y ) , fX (x)

fY (y ) = 4 y 3 , 0 < y < 1 4 r+4 2x fX |Y (x|y ) = 2 , 0 < x < y y 2y fY |X (y |x) = ,x<y<1 1 x2 2 yr E (X r |Y = y ) = r+2 2 (1 xr+2 ) r E (Y |X = x) = (r + 2)(1 x2 ) E (Y r ) = Misalkan (X, Y ) adalah peubah acak berpasangan dengan fungsi peluang bersama fX,Y (x, y ). Pandang persoalan memprediksi Y setelah X = x terobservasi. Prediktor dinotasikan sebagai y (x). Prediktor terbaik didenisikan sebagai (X ) yang meminimumkan fungsi Y [ ]2 (X ) = (y y (x))2 fX,Y (x, y ) dydx E Y Y

Prediktor terbaik adalah y (x) = E (Y |X = x). BUKTI: MA3081 Stat.Mat. 11 K. Syuhada, PhD.

Contoh/Latihan: 1. Misalkan X dan Y memiliki distribusi bersama fX,Y (x, y ) = 8xy, 0 < x < y < 1, maka fY |X (y |x) = 2y ,x<y<1 1 x2 2 (1 x3 ) 3 (1 x2 )

y (x) = E (Y |X = x) =

2. Misalkan (Y, X ) berdistribusi normal bivariat dengan E (Y ) = Y , E (X ) = 2 2 , V ar(X ) = X , Cov (X, Y ) = X,Y X Y . Distribusi X , V ar(Y ) = Y bersyarat Y , diberikan X , adalah (Y |X = x) 3. Tunjukkan bahwa [ ] EX fY |X (y |X ) = fY (y ) 4. Buktikan { [ ]} [ ] EX E h(Y )|X = E h(Y ) 5. Buktikan

[ ] [ ] V ar(Y ) = EX V ar(Y |X ) + V ar E (Y |X )

6. Misalkan X dan Y memiliki distribusi bersama 3y 2 fX,Y (x, y ) = 3 , 0 < y < x < 1 x Maka fY (y ) = 1.5 (1 y 2 ), 0 < y < 1 E (Y r ) = 3 , E (Y ) = 3/8, V ar(Y ) = 19/320 (r + 1)(r + 3)

fX (x) = 1, 0 < x < 1

MA3081 Stat.Mat.

12

K. Syuhada, PhD.

3 y2 , 0<y<x<1 x3 3 xr 3x 3 x2 E (Y r |X = x) = , E (Y |X = x) = , V ar(Y |X = x) = 3+r 4 80 V ar(E (Y |X )) = 3/64 fY |X (y |x) = E (V ar(Y |X )) = 1/80

MA3081 Stat.Mat.

13

K. Syuhada, PhD.

1.6

Fungsi Pembangkit Momen

Misalkan X peubah acak kontinu, fungsi pembangkit momen dari X adalah tX MX (t) = E (e ) = etx f (x)dx,

asalkan ekspektasi ada untuk t disekitar 0. Jika semua momen dari X tidak ada, maka fungsi pembangkit momen juga tidak ada. Fungsi pembangkit momen berkaitan dengan fungsi pembangkit peluang MX (t) = GX (et ) asalkan GX (t) ada untuk t disekitar 1. Jika MX (t) adalah fungsi pembangkit peluang maka MX (0) = 1. Contoh/Latihan: 1. Jika fX (x) = ex I0, (x), maka MX (t) = 2. Jika MX (t) ada maka Ma+bX (t) = 3. Jika Xi , i = 1, . . . , n saling bebas, MXi (t) ada untuk setiap i, dan S = Xi , maka MS (t) = 4. Fungsi pembangkit momen bersifat unik. Setiap distribusi memiliki fungsi pembangkit momen yang unik, dan setiap fungsi pembangkit momen berkorespondensi dengan tepat satu distribusi. Akibatnya, jika fungsi pembangkit momen ada maka fungsi pembangkit momen tersebut secara unik menentukan distribusinya. Beri contoh. 5. Pandang turunan dari MX (t) yang kemudian dievaluasi di t = 0. Apa yang dapat anda katakan? Dapatkah kita mendapatkan momen orde tinggi? 6. Dapatkah hasil diatas digunakan untuk distribusi diskrit? Ambil contoh distribusi Geometrik dengan parameter p.

MA3081 Stat.Mat.

14

K. Syuhada, PhD.

7. Misalkan Y U (a, b). Gunakan fungsi pembangkit momen untuk mendapatkan momen pusat (( )r ) a+b 2 E ((Y Y ) ) = E Y 2

MA3081 Stat.Mat.

15

K. Syuhada, PhD.

BAB 2 Distribusi Sampel, Likelihood dan Penaksir

Misalkan X1 , X2 , . . . , Xn sampel acak berukuran n (random sample of size n). Fungsi peluang n-variat nya adalah fX1 ,X2 , ,Xn (x1 , x2 , . . . , xn ) =
n i=1

fXi (xi )

Contoh/Latihan: 1. Misalkan X1 , X2 , . . . , Xn sampel acak dari distribusi Eksponensial dengan parameter . Fungsi peluang n-variatnya adalah... 2. Misalkan X1 , X2 , . . . , Xn sampel acak dari distribusi Uniform pada selang (a, b). Fungsi peluang n-variatnya adalah...

2.2

Likelihood

Misalkan fungsi peluang n-variat bergantung pada parameter yang tidak diketahui . Fungsi peluang tersebut ditulis sebagai fX1 ,X2 ,...,Xn (x1 , . . . , xn |1 , . . . , k ) atau fX (x|) 1

Contoh/Latihan: 1. Misalkan X1 , X2 , . . . , Xn sampel acak dari distribusi N (, 2 ). Fungsi peluang n-variat yang bergantung pada parameternya ditulis sebagai... Definisi Fungsi likelihood adalah ukuran yang menyatakan sebarapa sering nilai , diberikan bahwa x telah terobservasi. Fungsi likelihood BUKAN suatu peluang. Fungsi likelihood diperoleh dengan (i) menukar peran dan x dalam fungsi peluang n-variat, dan (ii) membuang suku yang tidak bergantung pada . Notasi: L() = L(|x) fX (x|) Contoh/Latihan: 1. Misalkan X1 , X2 , . . . , Xn sampel acak dari distribusi Eksponensial dengan parameter . Fungsi likelihoodnya adalah... function likefunction; % this function calculates the likelihood function of distribution % % created by K Syuhada, 14/3/2011 clear clc n = input(n = ); % size of random sample % data x = exprnd(0.5,n,1); sumx = sum(x); % parameter of exponential distribution lambda = 0.5:0.05:5; for i = 1:length(lambda) L(i) = (lambda(i)^n)*exp(-lambda(i)*sumx); end plot(lambda,L) MA3081 Stat.Mat. 2 K. Syuhada, PhD.

2. Misalkan X1 , X2 , . . . , Xn sampel acak dari distribusi Uniform pada selang (, b). Fungsi likelihoodnya adalah... Prinsip Likelihood Jika dua percobaan, yang melibatkan model dengan parameter , memberikan likelihood yang sama, maka inferensi terhadap haruslah sama. Ilustrasi Pandang percobaan 1 dimana sebuah koin dilantunkan sebanyak n kali secara bebas. Misalkan p adalah peluang muncul MUKA dan X peubah acak yang menyatakan banyaknya MUKA yang muncul. Fungsi peluang dari X dan fungsi likelihoodnya adalah... Pandang percobaan 2 dimana sebuah koin dilantunkan hingga diperoleh MUKA sebanyak 6 kali secara bebas. Misalkan Y peubah acak yang menyatakan banyaknya lantunan yang dibutuhkan agar diperoleh enam MUKA. Fungsi peluang dari X dan fungsi likelihoodnya adalah... Dari 2 percobaan diatas, misalkan kita ingin melakukan uji hipotesis: H0 : p = 0.5 versus H0 : p < 0.5 Nilai signkansinya atau p-value adalah... Penaksir Likelihood Maksimum Misalkan L() adalah fungsi likelihood (fungsi dari parameter ). Kita dapat menentukan nilai yang memaksimumkan L(). Penaksir untuk , yaitu disebut Penaksir Likelihood Maksimum (maximum likelihood estimator, MLE). Penaksir suatu parameter adalah fungsi dari peubah acak. Contoh/Latihan: 1. Diketahui sampel acak berukuran n dari distribusi Bernoulli (p). Fungsi likelihoodnya: L() =
xi

(1 )n

xi

, 0 < < 1.

Untuk menentukan nilai yang memaksimumkan L(), transformasikan L() menjadi log L(): () = log L() = xi log() + (n xi ) log(1 ), kemudian hitung turunan pertama () terhadap : d() xi xi n = . d 1 MA3081 Stat.Mat. 3 K. Syuhada, PhD.

Normalisasi dari turunan pertama tersebut memberikan nilai xi = , n yang mana sebagai penaksir ditulis sebagai berikut: Xi = = X. n (Pr: Tunjukkan bahwa ini memaksimumkan L() dengan menghitung turunan kedua). 2. Misalkan X1 , . . . , Xn sampel acak berdistribusi U (0, ). Tentukan yang untuk memaksimumkan L(0, ). Dengan kata lain, tentukan penaksir .

Sifat Penaksir , kita dapat menentukan sifat baik peSetelah kita mendapatkan penaksir dikatakan tak bias naksir. Salah satunya adalah sifat TAK BIAS. Penaksir apabila ) = . E ( Untuk contoh sampel acak Bernoulli, ( ) X1 + + Xn E () = E n 1 = E (X1 + + Xn ) n ) 1( E (X1 ) + + E (Xn ) = n 1 = ( + + ) n = = X adalah penaksir tak bias untuk . Jadi, penaksir bersifat bias maka selisih nilai ekspektasi dan Catatan: Jika suatu penaksir nilai tidak nol, atau ) = 0. E (

MA3081 Stat.Mat.

K. Syuhada, PhD.

2.3

Statistic Cukup

Definisi -1 Suatu statistik T = t(X) adalah CUKUP atau sucient untuk suatu keluarga distribusi fX (x|) JIKA dan HANYA JIKA fungsi likelihoodnya bergantung terhadap X hanya melalui T: L() = h(t(X), ) Definisi -2 Suatu statistik T = t(X) adalah CUKUP untuk suatu keluarga distribusi fX (x|) JIKA dan HANYA JIKA distribusi bersyarat dari X TIDAK BERGANTUNG pada : fX|T (x|t, ) = h(x) Definisi -3 Suatu statistik T = t(X) adalah CUKUP untuk suatu keluarga distribusi fX (x|) JIKA dan HANYA JIKA fungsi peluangnya dapat difaktorkan sebagai: fX (x|) = g (t(x)|) h(x) Contoh/Latihan: 1. Misalkan Xi untuk i = 1, . . . , n saling bebas dan berdistribusi identik Bernoulli(p). Tunjukkan bahwa Y = n i=1 Xi adalah statistik cukup. 2. Misalkan X1 , . . . , Xn sampel acak n berdistribusi Poisson dengan parameter . Tunjukkan bahwa T = i=1 Xi adalah statistik cukup. 3. Misalkan Xi untuk i = 1, . . . , n saling bebas dan berdistribusi identik adalah statistik cukup. N(, 1). Tunjukkan bahwa Y = X 4. Misalkan X1 , . . . , Xn sampel acak berdistribusi Gamma dengan paramn eter (, ). Tunjukkan bahwa T = i=1 ln(Xi ) adalah statistik cukup. 5. Pandang sampel acak berukuran n dari U (a, b), dengan a diketahui. Tunjukkan bahwa T = X(n) adalah statistik cukup. 6. Pandang sampel acak berukuran n dari N (, 2 ), dengan , 2 tidak diketahui. Tunjukkan bahwa statistik T berikut adalah cukup: ( 2 ) SX T = X MA3081 Stat.Mat. 5 K. Syuhada, PhD.

2.4

Distribusi Sampel

Misalkan X1 , X2 , . . . , Xn sampel acak berukuran n dari distribusi Poisson dengan parameter . Peubah acak Xi , i = 1, . . . , n saling bebas dan berdistribusi identik dengan fungsi peluang n-variat: P (X = x) =
n e xi en y = n , x ! x ! i i i =1 i=1

dengan y = xi . Dapat ditunjukkan juga Y = Xi cukup. Distribusi sampel dari Y adalah fY (y|) = en (n)y . y!

Misalkan Xi U (0, ). Peubah acak-peubah acak Xi tersebut saling bebas dan berdistribusi identik, dengan fungsi peluang: fX (x|) = Statistik T = X(n) cukup dan memiliki fungsi distribusi: P (X(n) x) = dan fungsi peluang: f (x) =

2.5

Statistik Terurut

Misalkan X1 , . . . , Xn sampel acak berukuran n dari suatu populasi yang berdistribusi tertentu, dengan fungsi peluang fX dan fungsi distribusi FX . Pandang X(k) , statistik terurut ke-k . Untuk menentukan fX(k) (x), pertama partisikan I1 = (, x]; I2 = (x, x + dx]; I3 = (x + dx, ). Fungsi peluang fX(k) (x) adalah peluang mengamati sejumlah k 1 dari X di I1 , tepat sebuah X di I2 , dan sejumlah n k dari X di I3 : ( ) ( )k 1 ( )1 ( )nk n FX ( x ) fX (x)dx 1 FX (x) fX(k) (x) k 1, 1, n k

MA3081 Stat.Mat.

K. Syuhada, PhD.

yang dengan metode diferensial maka kita peroleh ( ) ( )k1 ( )nk n fX(k) (x) = FX (x) 1 FX (x) fX (x) k 1, 1, n k Contoh/Latihan: 1. Fungsi peluang dari statistik terurut terkecil/terbesar adalah... 2. Statistik terurut ke-k pada distribusi U (0, 1) memiliki fungsi peluang...

2.6 2.7

Momen dari Mean dan Proporsi Sampel Teorema Limit Pusat

Teorema Misalkan X1 , . . . , Xn sampel acak berukuran n dari populasi dengan mean X 2 dan variansi X . Distribusi dari Zn = X X X / n

konvergen ke N (0, 1) untuk n . Catatan: Hal penting dari Teorema Limit Pusat (Central Limit Theorem) adalah bahwa kekonvergenan dari Zn ke distribusi normal akan terjadi apapun bentuk distribusi dari X . Kita dapat memanipulasi sedemikian hingga N (X , 2 /n), X X asalkan n besar. Ekspresi lain dari TLP adalah ( ) n (X X ) lim P c = (c) n X

MA3081 Stat.Mat.

K. Syuhada, PhD.

 Pandang: X1 + + Xn , ( n ) E Xi = n X ,
i=1

V ar

( n
i=1

) Xi
i=1 2 = n X ,

(n

lim P

Xi n X c n X

) = (c)

berdistribusi normal? n = 1? Seberapa besar n harus kita pilih agar X Bergantung pada distribusi dari data (parent distribution)! Misalkan X Exp(). Distribusi ini memiliki kemencengan (skewness) dan kelancipan (kurtosis): E (X X )3 3 = = 2, 3 X dan 4 = E (X X )4 3 = 6, 4 X

2 berdistribusi Ga(n, n). dengan X = 1/ dan X = 1/2 . Mean sampel X adalah Kemencengan (skewness) dan kelancipan (kurtosis) dari X 3 X 2 E (X ) = , 3 = 3 n X

dan 4 =
4 X E (X ) 3 = 6/n, 4 X

Perhatikan plot berikut:

MA3081 Stat.Mat.

K. Syuhada, PhD.

Misalkan X B (n, p) (ingat bahwa distribusi X tersebut sama dengan distribudi dari sejumlah n peubah acak Bernoulli(p)). Untuk n besar, ( ) p(1 p) p N p, n ( ) n(c p) P ( p c) p(1 p) X N (np, np(1 p)) ( ) 1 1 P (X = x) = P x X x + , x = 0, 1 , . . . , n 2 2 ) ( ) ( x 0.5 np x + 0.5 np , np(1 p) np(1 p) dimana menambah dan mengurangi dengan 0.5 disebut continuity correction. Koreksi kekontinuan untuk pendekatan normal terhadap fungsi distribusi dari X dan p adalah ( ) 1 P (X c) = P X x + , x = 0, 1, . . . , n 2 ) ( x + 0.5 np np(1 p) dan ) 1 P ( p c) = P p c+ , c = 0/n, 1/n, . . . , n/n 2n ( ) n(c + 0.5/n p) p(1 p) (

MA3081 Stat.Mat.

K. Syuhada, PhD.

BAB 3 Penaksiran dan Selang Kepercayaan

adalah fungsi peubah acak. Pada penaksiran parameter , misalnya, penaksir Nilai taksirannya TIDAK akan pernah sama dengan nilai parameternya. Misalkan T = T (X) adalah penaksir untuk . Didenisikan: bT = E (T ) = E (T ) , dan
2 V ar(T ) = T = E (T T )2 = E (T ) ; T = E (T ),

adalah bias dan variansi dari penaksir T . Selain itu, didenisikan pula, MSE atau Mean Square Error, MSET () = E (T )2 = V ar(T ) + b2 T, Misalkan X1 , . . . , Xn sampel acak dari N (, 2 ). Penaksir untuk 2 adalah S2 = dan/atau
n 1 )2 , V = (Xi X n i=1 n 1 )2 , (Xi X n 1 i=1

Bias and MSE dari kedua penaksir adalah bS 2 = bV = dan MSES 2 = MSEV = Catatan: Penaksir dari deviasi standar dari suatu penaksir disebut standard error atau SE. Apakah SE dari jenis pengambilan sampel (sampling): Apapun asalkan tanpa pengembalian? Bernoulli tanpa pengembalian?

3.2

Konsistensi

Salah satu sifat dari penaksir yang baik adalah sifat tak bias. Kita akan mempelajari sifat baik yang lain yaitu konsisten. Namun sebelumnya, perhatikan Ketaksamaan Chebyshev. Misalkan X peubah acak dengan fungsi peluang fX (x). Misalkan h(X ) fungsi non-negatif dari X dan ekpektasinya ada serta k adalah konstanta positif. Maka P (h(X ) k ) Bukti: Aplikasi 1 Ketaksamaan Chebyshev. Misalkan E (X ) = X dan V ar(X ) = 2 X < . Maka [ ] |X X |2 1 2 P k 2. 2 X k Bukti: Aplikasi 2 Ketaksamaan Chebyshev. Misalkan T peubah acak (penaksir dari 2 parameter ) dengan E (T ) = T dan V ar(T ) = T < . Maka P [|X | < ] 1 Bukti: MSEX () . 2 E (h(X )) . k

MA3081 Stat.Mat.

K. Syuhada, PhD.

Konsistensi Denisi: Barisan dari penaksir-penaksir, {Tn }, disebut KONSISTEN untuk jika
n

lim P (|Tn | < ) = 1,

untuk setiap > 0. Konvergen dalam Peluang Denisi: Barisan dari penaksir-penaksir, {Tn }, KONVERGEN dalam PELUANG untuk jika barisan tersebut konsisten untuk . Notasi: Tn prob . Contoh/Latihan: adalah mean sampel dari suatu s.a 1. (Hukum Bilangan Besar) Jika X berukuran n dengan mean X , maka prob X . X Bukti: 2. Sebuah penaksir untuk dikatakan Mean Square Consistent jika
n

lim MSETn () = 0.

Buktikan bahwa jika sebuah penaksir memiliki sifat MSC maka penaksir tersebut konsiten.

MA3081 Stat.Mat.

K. Syuhada, PhD.