II. TINJAUAN PUSTAKA
Dalam bab ini akan dibahas beberapa konsep dasar, definisi-definisi serta teorema yang berkaitan dalam hal pendugaan parameter pada model linier campuran ini, yaitu sebagai berikut :
2.1
Model Linear Umum
Model Linear merupakan pemodelan khusus dengan ciri linear dalam parameter.
Definisi 2.8 ๐ = ๐ฟ๐ท + ๐บ
(Persamaan 1)
Dimana : Y
: vektor peubah acak ๐ ๐ฅ1 yang teramati;
X
: mariks ๐ ๐ฅ ๐ (๐ > ๐) dengan unsur-unsurnya adalah bilangan tertentu yang diketahui;
๐ท
: vektor parameter ๐ ๐ฅ1 yang tidak diketahui nilainya;
๐บ
: vektor peubah acak ๐ ๐ฅ1 yang tidak teramati,
Dengan, E(๐) = 0 dan cov (๐บ) = โ
Sehingga persamaan linear pada persamaan 1 dinamakan model linear umum. Model tersebut mempunyai pengertian-pengertian khusus tergantung pada : 1.
Distribusi dari ๐บ
2.
Struktur matriks kovarians โ
Peringkat struktur dari matriks X, jika peringkat atau rank dari matriks X sama dengan jumlah kolomnya, maka matriks tersebut dinamakan matriks berpangkat penuh dan jika peringkat matriksnya tidak penuh maka modelnya dinamakan model tidak penuh (non-full rank model) (Usman dan Warsono, 2009).
2.2
Model Linear Campuran
Model linear campuran merupakan model yang terdiri dari tetap (fixed effect) dan faktor acak (random effect). Suatu faktor dikatakan tetap (fixed) jika dalam suatu penelitian faktor tersebut telah mewakili semua tingkat kemungkinan faktor yang ada yang telah ditetapkan sebelumnya. Suatu faktor dikatakan acak (random) jika faktor tersebut merupakan sampel yang diambil secara acak pada suatu populasi yang kiranya dapat mewakili karakteristik dari populasi tersebut, (Little dkk,1996). Model linear campuran dapat ditulis sebagai berikut : ๐ = ๐ฟ๐ท + ๐๐ธ + ๐บ
(Persamaan 2)
Dimana, X
: Matriks pada pengaruh tetap pada model
Z
: Matriks pada pengaruh acak pada model
ฮฒ
: Vektor pada pengaruh tetap
๐ธ
: Vektor pada pengaruh acak
5
dan dengan asumsi
: ๐ธ ~ N (๐, ๐ฎ) ๐บ ~ N (๐, ๐น)
serta cov (Y) =
= ๐น + ๐ ๐ฎ๐โฒ
Model linear campuran ini biasa digunakan pada data yang tak saling bebas serta data yang tak homogen (Henderson, 1984). Dalam teori model linear campuran, terdapat persamaan model linear campuran (Mixed Model Equation).
2.2.1
Persamaan Model linear Campuran (Mixed Model Equation)
Persamaan model linear campuran merupakan persamaan yang terbentuk dari penurunan fungsi gabungan dari model linear campuran yang diturunkan masingmasing terhadap kedua parameter yang terdapat pada model linear campuran. Persamaan tersebut sama halnya seperti persamaan normal seperti biasanya, sejauh ini persamaan model linear campuran digunakan secara serempak untuk menentukan BLUE dan BLUP. Kita misalkan R adalah matriks nonsingular. Jika ๐ธ acak, persamaan normal untuk model linear campuran dapat ditulis sebagai berikut : ๐ฟโฒ ๐นโ๐ ๐ฟ, ๐ ๐โฒ
๐ท = ๐ฟโฒ ๐นโ๐ ๐ ๐ธ ๐โฒ
Atau, ๐ฟโฒ๐นโ๐ ๐ฟ ๐ฟโฒ๐นโ๐ ๐ ๐โฒ๐นโ๐ ๐ฟ ๐โฒ๐นโ๐ ๐
โ๐ ๐ท = ๐ฟโฒ๐นโ๐ ๐ ๐ธ ๐โฒ๐น ๐
6
Untuk G nonsingular, persamaan dalam model linear campuran dapat didefinisikan sebagai berikut : ๐ฟโฒ๐นโ๐ ๐ฟ ๐ฟโฒ๐นโ๐ ๐ ๐โฒ๐นโ๐ ๐ฟ ๐ฎโ๐ + ๐โฒ๐นโ๐ ๐
โ๐ ๐ท = ๐ฟโฒ๐นโ๐ ๐ ๐ธ ๐โฒ๐น ๐
(Milliken dan Littel, 1998) 2.3
Dalam
Pendugaan Parameter
statistika
inferensial,
sangat
dibutuhkan
pemahaman-pemahaman
mengenai kaidah-kaidah pengambilan kesimpulan tentang suatu parameter populasi berdasarkan karakteristik sampel. Hal ini membangun apa yang disebut dengan pendugaan titik dari suatu fungsi kepekatan peluang parameter yang tidak diketahui. Dalam model linear, terdapat beberapa metode yang digunakan dalam pendugaan parameter, seperti Metode Kuadrat Terkecil (Leaqst Square), MLE (Maximum likelihood Estimation) dan RMLE (Restricted Maximum likelihood Estimation) (Hoog dan Craig, 1995).
2.3.1
Maximum Likelihood Estimation (MLE)
MLE (Maximum Likelihood Estimation) merupakan metode pendugaan parameter yang menggunakan pendekatan distribusi dari data yang dimiliki serta asumsi distribusi yang diberlakukan oleh data tersebut, sehingga dapat diperoleh fungsi likelihood dari suatu data tersebut. MLE (Maximum Likelihood Estimation) ini menyediakan metode umum, dimana suatu kondisi dalam suatu sampel acak
7
terdapat penduga, dimana penduga tersebut mempunyai sifat yang konsisten dan sifat-sifat lain yang sangat diperlukan sebagai suatu penduga. Dalam menggunakan metode MLE, pertama kita misalkan bahwa peubah acak dari suatu populasi adalah X, dimana X mempunyai fungsi peluang yang mewakili beberapa parameter ๐ โถ ๐ท๐ ๐ = ๐ฅ = ๐(๐ฅ; ๐). Lalu kita misalkan bahwa fungsi f diketahui tetapi nilai ๐ tidak diketahui. Fungsi peluang bersama dari peubah acak (๐ฅ1, ๐ฅ2, โฆ โฆ ๐ฅ๐ ) dapat ditulis menjadi :
๐
๐ฟ ๐ = ๐ (๐ฅ1, ๐ฅ2, โฆ โฆ , ๐ฅ๐ ; ๐) =
๐ (๐ฅ๐ ; ๐) ๐=1
Fungsi diatas tersebut lebih dikenal dengan sebutan likelihood function dari suatu sampel. Sifat dari MLE ini diperlukan untuk memilih penduga dari parameter yang tidak diketahui. Jika suatu kelompok distribusi ingin menentukkan dua atau lebih dari parameter yang tidak diketahui, yaitu ๐1, ๐2, โฆ โฆ ๐๐ maka fungsi likelihood dapat ditulis dalam bentuk berikut :
๐ฟ ๐1, ๐2, โฆ โฆ ๐๐
= ๐ ๐ฅ1, โฆ โฆ , ๐ฅ๐ ; ๐1, , โฆ โฆ , ๐๐ =
๐ ๐=1 ๐
(๐ฅ๐ ; ๐1, , โฆ โฆ , ๐๐ ) (Brunk, 1975)
8
2.4
Distribusi Normal Multivariate
Fungsi kepekatan normal ganda (multivariate normal) adalah bentuk generalisasi dari fungsi kepekatan univariate normal dengan p โฅ 2 dimensi. Fungsi kepekatan dari peubah acak Y yang menyebar normal dengan nilai tengah ๐ dan ragam ๐ 2 , atau fungsi kepekatan normal tersebut dapat dinotasikan ๐ ~ ๐(๐; ๐ 2 ). Sehingga fungsi kepekatan normal ganda adalah: 1
๐ ๐ฆ =
2๐๐ 2
exp โ
(๐ฆ โ ๐)2 2๐ 2
Misal, diberikan peubah acak ๐โฒ = {Y1. Y2,........Yp} dimana Yi menyebar normal dan saling bebas (independent) dengan mean ๐ dan ragam ๐ 2 . Sehingga fungsi kepekatan bersama dari peubah acak Y yaitu : ๐๐ (๐)
=
๐ ๐=1 ๐๐ (๐)
=
1 ๐ ๐=1 2๐๐ 2
= (2๐)โ๐/2 = (2๐)โ๐/2
exp โ 1 ๐๐
(๐ฆ โ๐ )2 2๐ 2
exp โ
๐ 2 ๐๐
โ1/2
๐ ๐=1(๐ฆ
โ ๐)2 /2๐ 2
exp โ ๐ฆ โ ๐
โฒ
๐ 2 ๐๐
โ1
(๐ฆ โ ๐)/2
Fungsi kepekatan bersama yang terbentuk dari distribusi normal multivariate yang saling bebas dapat ditulis ๐ ~ ๐(๐; ๐ 2 ๐), dimana vektor nilai tengah dan matriks peragam adalah sebagai berikut : ๐1 ๐2 E (Y) = ๐ = โฎ ๐๐
๐2 ,dan cov (Y) = 0 โฎ 0
0 ๐2 โฎ 0
โฆ 0 โฆ 0 = ๐2๐ ๐ โฎ โฑ 2 โฆ ๐ ๐
9
Bentuk [(
๐ฆโ๐ ๐
)]2 dari eksponen fungsi sebaran normal mengukur jarak kuadrat
dari ๐ฆ๐ ke ๐ dalam unit simpangan baku. Bentuk ini dapat digeneralisasikan untuk px1 vektor Y dari pengamatan beberapa peubah sebagai : ๐ฆ โ ๐ โฒ ๐โ1 (๐ฆ โ ๐) Dimana, ๐= matriks ragam peragam (covariance) ๐2 ๐= 0 โฎ 0
0 ๐2 โฎ 0
โฆ 0 โฆ 0 โฎ โฑ 2 ๐ ๐ โฆ
Sehingga, fungsi kepekatan peluang bersama dari distribusi normal multivariate dapat ditulis sebagai berikut :
๐ ๐ฆ =
1 (2๐)๐/2
๐
1/2
exp โ
1 2
๐ฆ โ ๐ โฒ ๐โ1 (๐ฆ โ ๐) (Timm, 2002)
2.5
Aplikasi Model linear campuran pada rancangan Faktorial dalam RAL (Rancangan Acak Lengkap)
Banyak Penelitian dilakukan untuk mengetahui pengaruh dua atau lebih factor. Rancangan faktorial merupakan rancangan yang sangat efisien digunakan dalam penelitian demikian. Pada faktorial, semua kombinasi perlakuan dari setiap level/taraf suatu faktor dikombinasikan secara penuh dengan dengan level/taraf
10
faktor lainnya. Dengan demikian setiap pengamatan menyediakan informasi mengenai semua pengaruh dan dapat mengetahui respon level pada level lainnya.
Struktur rancangan perlakuan faktorial dapat diterapkan pada rancangan acak lengkap. Dalam model campuran, misal terdapat dua faktor dalam suatu percobaan yaitu A dan B, bila salah satu dari A atau B yang digunakan dalam penelitian diambil secara acak dari populasinya masing-masing sedangkan lainnya ditetapkan, maka penelitian tersebut haruslah dianalisis dengan model linear campuran. Model linear campuran untuk factorial dalam rancangan acak lengkap:
๐ = ๐๐ + ๐ฝ๐ + ๐พ๐ + (๐ฝ๐พ)๐๐ + ๐๐๐๐
dimana
: i = 1,2,โฆ.b , j = 1,2,โฆ.t
dengan
: ๐ฝ ~ ๐ 0, ๐ 2
, k = 1,2,โฆ.r
๐๐๐๐ ~ ๐ 0, ๐ 2
2.6
Karaktersistik Estimabilitas
Untuk melakukan uji hipotesis pada suatu model linear dalam bentuk ๐ณ๐ท, haruslah estimable. Matriks ๐ณ๐ท dikatakan estimable, jika dan hanya jika salah satu dari kondisi berikut terpenuhi : 1.
L = BX, untuk sebarang matriks B
2.
๐
3.
๐ ๐ฟ ๐ฐ โ ๐ณโ ๐ณ
๐ฟ = ๐(๐ฟ) ๐ณ = ๐ ๐ฟ โ ๐ ๐ณ , untuk suatu g-invers ๐จโ
11
4.
๐ณ๐ฟโ๐ฟ = ๐ณ untuk suatu g-invers ๐ฟโ
5.
โ ๐ณ๐ฟโ ๐ adalah invariant untuk setiap kuadrat terkecil g-inverse ๐ฟ๐
6.
โ ๐(๐ณ๐ฟโ ๐ ) adalah invariant untuk setiap kuadrat terkecil g-inverse ๐ฟ๐
7.
โ ๐ ๐ณ๐ฟโ ๐ = ๐(๐ณ) untuk setiap kuadrat terkecil g-inverse ๐ฟ๐
dimana g-inverse itu sendiri didefinisikan : Jika matriks X berukuran mxn dan jika ๐ฟโ ada dan memenuhi empat syarat berikut,yaitu : 1.
๐ฟ๐ฟโ simetrik;
2.
๐ฟโ๐ฟ simetrik;
3.
๐ฟ๐ฟโ๐ฟ = ๐ฟ; dan
4.
๐ฟโ๐ฟ๐ฟโ = ๐ฟโ.
maka ๐ฟโ dikatakan g-inverse X. (Usman dan Warsono, 2009).
2.7
Uji Hipotesis
Hipotesis pada dasarnya merupakan suatu proporsi atau anggapan yang mungkin benar, dan sering digunakan sebagai dasar pembuatan keputusan atau pemecahan persoalan ataupun untuk dasar penelitian lebih lanjut. Hipotesis yang akan diuji biasanya disimbolkan H0 (Hipotesis nol) dan disertai dengan H1 (Hipotesis alternatif). H1 akan secara otomatis diterima jika H0 ditolak. Cara untuk merumuskan H0 dan H1 tergantung pada jenis parameter yang akan diuji dan jenis data (informasi yang dimiliki oleh peneliti atau menurut rencananya akan dipilih). Pada pengujian hipotesis terdapat dua macam kesalahan yang secara ringkas akan disajikan dalam Tabel 1 :
12
Tabel 1. Kesalahan Pada Uji Hipotesis
H0 Kesimpulan
Hipotesis Benar
Hipotesis Salah
Tidak tolak H0
1-๐ผ
Kesalahan Tipe II (Beta)
Tolak H0
Kesalahan Tipe I Taraf nyata (ฮฑ)
1-๐ฝ
Dari Tabel 1 dapat dijelaskan bahwa didalam uji hipotesis ini terdapat dua jenis kesalahan yang dapat terjadi di dalam pengujian hipotesis. Kesalahan itu bisa terjadi karena kita menolak hipotesis padahal hipotesis itu benar atau kita menerima hipotesis padahal hipotesis itu salah. Kesalahan yang disebabkan karena kita menolak hipotesis padahal hipotesis tersebut benar, disebut kesalahan jenis I, sebaliknya kesalahan yang disebabkan karena kita menerima hipotesis padahal hipotesis itu salah maka disebut kesalahan jenis II. Sehingga power / kuasa uji adalah peluang menolak H0 dimana H0 tidak benar atau sama saja dengan peluang (statistik uji akan jatuh dalam penolakan wilayah ketika hipotesis nol salah). Kuasa uji = 1 - beta. Fungsi dari uji hipotesis ini adalah untuk menguji suatu teori, menyarankan teori baru apabila hasil pengujian hipotesis dapat membentuk proposisi asumsi atau penjelasan tentang suatu peristiwa, serta mendeskripsikan fenomena sosial yang artinya hipotesis memberikan informasi kepada peneliti tentang apa yang nyatanyata terjadi secara empirik, (Supranto, 1986).
13
2.7.1 Uji Hipotesis Parameter Faktor Tetap (fixed effect) dan Faktor Acak (random effect) Pada Model Linear Campuran
Uji hipotesis pada faktor tetap (fixed effect) dan faktor acak (random effect) menggunakan pendekatan distribusi F, yang akan disajikan pada Tabel 2, tetapi sebelum disajikannya tabel bagi model linear campuran, akan disajikan terlebih dahulu Expected Mean Square (EMS), Sebagai berikut : Tabel 2. Tabel Expected Mean Square (EMS) F a i
R b j
R r k
0
b
r
2 ๐๐2 + ๐๐๐ฝ๐พ + ๐๐ ฮฑ2 /(๐ โ 1)
๐พ๐
a
1
r
2 ๐๐2 + ๐๐๐ฝ๐พ + ๐๐๐๐พ2
(๐ฝ๐พ)๐๐
1
1
r
2 ๐๐2 + ๐๐๐ฝ๐พ
๐ ๐๐
1
1
1
๐๐2
Fixed/Random Jumlah dari level Faktor ๐ฝ๐
๐
Expected Mean Square (EMS)
Selanjutnya dapat dibentuk tabel ANOVA dari model linear campuran sebagai berikut : Tabel 3. Tabel ANOVA Faktor ๐ฝ๐
Df (a-1)
SS (Sum Square)
๐พ๐
(b-1)
Y โฒ [ Ja โ Ib โ
(๐ฝ๐พ)๐๐
(a-1)(b-1)
Y โฒ [ Ia โ J]Y
๐ ๐๐
ab(r-1)
Y โฒ [Ia โ Ib โ Ir โ
๐
Y โฒ [ Ia โ
1 a
1
1
b
r
Ja โ J โ
1
1
a
b
1
J a a
Jb โ
โ Ib โ 1 r
1 b
MS (Mean Square) SS๐ฝ๐ /(a-1)
Jr ]Y 1 r
SS๐พ๐ /(b-)1
J]Y
Jb โ
Jr ]Y
1 r
SS(๐ฝ๐พ)๐๐ /(a-1)(b-1)
F
๐๐๐ฝ๐ ๐๐(๐ฝ๐พ)๐๐ ๐๐๐พ๐ ๐๐(๐ฝ๐พ)๐๐ ๐๐(๐ฝ๐พ)๐๐ ๐๐๐ ๐๐ ๐
SS๐ ๐๐ ๐ / ab(r-1)
14
Berdasarkan Tabel 3, yaitu tabel ANOVA, uji statistik pada model linear menggunakan pendekatan uji-F. F-Hitung didefinisikan sebagai berikut : ๐๐โฒ
F-Hitung = ๐๐โฒโฒ ~ ๐น๐,๐ Dengan derajat bebas : p=
(๐๐ โฒ )2 (๐ ๐ โฒ )2 ๐ท๐
dan
q=
(๐๐ โฒโฒ )2 (๐ ๐ โฒโฒ )2 ๐ท๐
15