Rancangan Acak Lengkap
(Completely Randomized Design)
1
Pendahuluan
Misalkan kita mempunyai t perlakuan dan N = tr satuan percobaan. Bagilah satuan-satuan percobaan tersebut secara acak menjadi t kelompok yang masing-masing terdiri atas r satuan percobaan. Alokasikan t perlakuan ke dalam kelompok tersebut sedemikian rupa sehingga perlakuan ke-i diterapkan sebanyak r kalidi setiap satuan percobaan dalam kelompok ke-i (i = 1, 2, …, t). Prosedur ini mendefinisikan rancangan acak lengkap dengan ulangan yang sama untuk t perlakuan.
2
Proses Pengacakan
Beri nomor satuan percobaan dengan 1, 2, …, N. Buat gulungan kertas dan beri nomor k = 1, 2, … N. Aduk gulungan kertas dan ambil satu, tandai gulungan tersebut dengan label (11). Sisihkan gulungan tersebut, dan ambil gulungan kertas berikutnya lalu beri label (12). Lakukan berkali-kali sampai seluruh gulungan terambil dan yang terakhir beri label (tr). Setiap gulungan sekarang telah mendapat nomor dan label. Misalnya suatu gulungan bernomor k dan berlabel (ij), maka ini merepresentasikan satuan percobaan ke k merupakan ulangan ke j dari perlakuan ke i.
3
Proses Pengacakan …
Asosiasikan bilangan 1, 2, …, N = tr secara acak dengan satu set label 11,12, …, 1r, …, tr. Misalnya N = 12 dan t = 3 dan r = 4. Dengan dua digit bilangan random, kita dapat menggunakan 01, 02, 03, …, 12 dan mengabaikan bilangan yang lain. Alternatif lain adalah dengan mengasosiasikan bilangan acak 00 s.d 07 dengan satuan percobaan 1, bilangan 08 s.d 15 dengan satuan percobaan 2, dan seterusnya sampai bilangan 88 s.d 95 dengan satuan percobaan 12 dan mengabaikan bilangan 96, 97, 98, 99.
Misalnya, bilangan acak yang terambil sebagai sampel adalah sebagai berikut:
4
63 | 07 | 84 | 70 | 89 | 78 | 48 | 41 | 48 | 51 | 70 | 65 | 62 | 52 | 16 | 53 | 34 | 19 | 21 | 98 | 18 | 21 | 15 | 31 |
bilangan-bilangan tersebut berasosiasi dengan satuan percobaan berikut:
8 | 1 | 11 | 9 | 12 | 10 | 7 | 6 | 7 | 7 | 9 | 9 | 8 | 7 | 3 | 7 | 5 | 3 | 3 | × | 3 | 3 | 2 | 4 |
Masing-masing satuan percobaan berasosiasi dengan label berikut: | |||||||||||||||||||||||
11 | 12 | 13 | 14 | 21 | 22 | 23 | 24 | | | | | | | 31 | | 32 | | | | | | 33 | 34 |
Artinya, satuan percobaan 8, 1, 11 dan 9 mendapat perlakuan 1; satuan percobaan 12, 10, 7 dan 6 mendapat perlakuan 2, dan seterusnya
Proses Pengacakan …
Bangkitkan bilangan1, 2, …, N = tr di kolom C1 untuk merepresentasikan satuan percobaan. Bangkitkan di kolom lainnya, misalnya di C2, label 11, 12, …, 1r, …, tr untuk merepresentasikan perlakuan. Misalnya t = 3 dan r = 4.
5
Pilih menu
Calc > Random Data > Sample From Columns …
6
Hasil pengacakan menunjukkan bahwa satuan percobaan 8, 9, 3 dan 10 mendapat perlakuan 1; satuan percobaan 2, 11, 1 dan 5 mendapat perlakuan 2; dan satuan percobaan 12, 7, 4 dan 6 mendapat perlakuan 3.
Model
7
Tabel Analisis Varians
8
Sumber | d.b. | JK | KT | E(KT) |
Perlakuan | t – 1 | JK(Perl) | | |
Galat | t(r – 1) | JK(Galat) | | |
Total | tr – 1 | JK(Total) | | |
Perhatikan bahwa
JK(Total) = JK(Perl) + JK(Galat)
Nilai Harapan bagi Kuadrat Tengah
9
Nilai Harapan bagi Kuadrat Tengah …
10
Nilai Harapan bagi Kuadrat Tengah …
11
Pengujian Hipotesis
H0: τ1 = τ2 = …= τt = 0
H1: paling tidak ada satu τi ≠ 0
dengan derajat bebas t – 1 dan t(r – 1)
12
Keuntungan dan Kekurangan RAL
Keuntungan
Kekurangan
13
Kadar nitrogen pada red clover yang diinokulasi kultur strain Rhizobium trifolii dan Rhizobium meliloti (steel & Torries, 1993 hal 171)
14
3DOk1 | 3DOk5 | 3DOk4 | 3DOK7 | 3DOk13 | Komposit |
19.4 | 17.7 | 17.0 | 20.7 | 14.3 | 17.3 |
32.6 | 24.8 | 19.4 | 21.0 | 14.4 | 19.4 |
27.0 | 27.9 | 9.1 | 20.5 | 11.8 | 19.1 |
32.1 | 25.2 | 11.9 | 18.8 | 11.6 | 16.9 |
33.0 | 24.3 | 15.8 | 18.6 | 14.2 | 20.8 |
Percobaan dilakukan dalam rumah kaca dengan menggunakan rancangan acak lengkap dengan lima pot per perlakuan
Tabel Analisis Varians
15
Sumber | d.b. | JK | KT | F |
Perlakuan | 6 – 1 = 5 | | | |
Galat | 6(5 – 1)= 24 | | | |
Total | 30 – 1 = 29 | | | |
Percobaan Sederhana
16
Output Minitab
One-way ANOVA: setelah versus treatment
Source DF SS MS F P
treatment 2 6447 3224 7.68 0.004
Error 18 7557 420
Total 20 14004
S = 20.49 R-Sq = 46.04% R-Sq(adj) = 40.04%
Individual 95% CIs For Mean Based on
Pooled StDev
Level N Mean StDev +---------+---------+---------+---------
1 7 96.00 6.93 (-------*-------)
2 7 119.43 21.09 (-------*-------)
3 7 138.86 27.69 (-------*--------)
+---------+---------+---------+---------
80 100 120 140
Pooled StDev = 20.49
17
RAL dengan ulangan tidak sama
18
Model & Pengujian hipotesis
dimana i = 1, 2, …, t sedangkan j = 1, 2, …, ri
dengan derajat bebas t – 1 dan
19
Tabel Analisis Varians
20
Sumber | d.b. | JK | KT | E(KT) |
Perlakuan | t – 1 | JK(Perl) | | |
Galat | | JK(Galat) | | |
Total | | JK(Total) | | |
Suatu percobaan ‘tanpa rancangan’
21
5 M
A = A. aulacocarpa
Ar = A. auriculiformis
C = A. crassicarpa
M = A. mangium
4 A
3 A
9 A
10 M
1 C
8 Ar
7 M
6 Ar
2 C
12 M
11 M
A. aulacocarpa | A. auriculiformis | A. crassicarpa | A. mangium |
28,32 | 27,25 | 30,10 | 49,18 |
33,40 | 32,57 | 32,70 | 47,68 |
29,70 | | | 41,62 |
| | | 38,25 |
| | | 44,82 |
Tabel Analisis Varians
Sumber keragaman | db | JK | KT | F | P |
Perlakuan | 3 | 552,33 | 184,11 | 13,31 | 0,002 |
Galat | 8 | 110,63 | 13,83 | | |
Total | 11 | 662,96 | | | |
Subsampling dalam RAL
dengan i = 1, 2, …, t sedangkan j = 1, 2, …, r’ dan k = 1, 2, …, n. Dimana εij adalah galat percobaan dan ηijk adalah galat pengamatan. Diasumsikan bahwa εij i.i.d. dan ηijk i.i.d.
Perhatikan bahwa
23
Tabel Analisis Varians
sehingga
24
Tabel Analisis Varians …
25
Sumber | d.b. | JK | KT | E(KT) |
Perlakuan | t – 1 | JK(Perl) | KT(Perl) | |
Galat Percobaan | t(r’ – 1) | JK(GalatPercobaan) | KT(GalatPercobaan) | |
Galat Sampling | tr'(n – 1) | JK(GalatSampling) | KT(GalatSampling) | |
Total | tr‘n – 1 | JK(Total) | | |
Inferens dengan Subsampling
dengan derajat bebas t – 1 dan t(r’ – 1).
dan
26
Pertumbuhan Batang Tanaman Mint�(Steel & Torries, 1993 hal 189)
27
Sampling | Perlakuan 1 | Perlakuan 2 | Perlakuan 3 | ||||||
Pot 1 | Pot 2 | Pot 3 | Pot 1 | Pot 2 | Pot 3 | Pot 1 | Pot 2 | Pot 3 | |
1 | 3.5 | 2.5 | 3.0 | 5.0 | 3.5 | 4.5 | 5.0 | 5.5 | 5.5 |
2 | 4.0 | 4.5 | 3.0 | 5.5 | 3.5 | 4.0 | 4.5 | 6.0 | 4.5 |
3 | 3.0 | 5.5 | 2.5 | 4.0 | 3.0 | 4.0 | 5.0 | 5.0 | 6.5 |
4 | 4.5 | 5.0 | 3.0 | 3.5 | 4.0 | 5.0 | 4.5 | 5.0 | 5.5 |
yij• | 15.0 | 17.5 | 11.5 | 18.0 | 14.0 | 17.5 | 19.0 | 21.5 | 22.0 |
yi•• | 44.0 | 49.5 | 62.5 | ||||||
y••• | 156 | ||||||||
Penghitungan Jumlah Kuadrat
28
Analisis dengan Minitab
29
Stat>Anova>General Linear Model
30
Output Minitab
General Linear Model: Data versus Perl, Pot
Factor Type Levels Values
Perl fixed 3 1, 2, 3
Pot(Perl) fixed 9 C1, C2, C3, C4, C5, C6, C7, C8, C9
Analysis of Variance for Data, using Adjusted SS for Tests
Source DF Seq SS Adj SS Adj MS F P
Perl 2 15.0417 15.0417 7.5208 15.33 0.000
Pot(Perl) 6 8.2083 8.2083 1.3681 2.79 0.031
Error 27 13.2500 13.2500 0.4907
Total 35 36.5000
S = 0.700529 R-Sq = 63.70% R-Sq(adj) = 52.94%
Unusual Observations for Data
Obs Data Fit SE Fit Residual St Resid
5 2.50000 4.37500 0.35026 -1.87500 -3.09 R
R denotes an observation with a large standardized residual
31
Asumsi yang Mendasari Analisis Varians
32
Transformasi
33
Transformasi …
34
Transformasi …
35
Transformasi …
36