Anova rancangan acak kelompok (RAK) dan uji lanjut Duncan

Pada postingan sebelumnya, saya menjelaskan secara ringkas dan mudah cara menganalisis split plot dan uji lanjutan Tukey. Penelitian saat skripsi tersebut memang menggunakan rancangan splitplot. Sekedar untuk latihan, data (bisa didownload disini) tersebut kita analisis menggunakan rancangan acak kelompok. Seperti biasa, dengan asumsi jendela Rcmdr sudah terbuka, masukkan script berikut.

library(agricolae)

#Rancangan Acak Kelompok/ Randomized Complete Block Design (RCBD)
model <- aov(Batang ~ Rambatan*Pupuk, bobot)
summary(model)
cv.model(model)

#Uji Duncan
duncan.test(model,”Rambatan”,alpha=0.05,console=TRUE, main=”Bobot Batang Saga Manis”)
duncan.test(model,”Pupuk”,alpha=0.05,console=TRUE, main=”Bobot Batang Saga Manis”)

Hasilnya sebagai berikut:

> library(agricolae)
> model <- aov(Batang ~ Rambatan*Pupuk, bobot)
> summary(model)
Df Sum Sq Mean Sq F value Pr(>F)
Rambatan 1 8.99E+07 89924161 2.182 0.1526
Pupuk 3 1.40E+09 4.65E+08 11.281 8.24E-05 ***
Rambatan:Pupuk 3 4.53E+08 1.51E+08 3.662 0.0264 *
Residuals 24 9.89E+08 41210765
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ‘ 1
> cv.model(model)
[1] 66.8057
> duncan.test(model,”Rambatan”,alpha=0.05,console=TRUE, main=”Bobot Batang Saga Manis”)
Study: Bobot Batang Saga Manis
Duncan’s new multiple range test
for Batang
Mean Square Error:  41210765
Rambatan,  means
Batang std r Min Max
R0 7932.958 9525.989 16 0.762 24958
R1 11285.65 9917.613 16 0.928 30642
alpha: 0.05 ; Df Error: 24
Critical Range
2
4684.343
Means with the same letter are not significantly different.
Groups, Treatments and means
a  R1 11290
a  R0 7933
> duncan.test(model,”Pupuk”,alpha=0.05,console=TRUE, main=”Bobot Batang Saga Manis”)
Study: Bobot Batang Saga Manis
Duncan’s new multiple range test
for Batang
Mean Square Error:  41210765
Pupuk,  means
Batang std r Min Max
P0 5164.154 4769.845 8 7.89 13832
P1 20417.51 9607.771 8 24.08 30642
P2 9387.583 9749.089 8 11.25 23756
P3 3467.961 2952.076 8 0.762 7744
alpha: 0.05 ; Df Error: 24
Critical Range
       2        3        4
6624.661 6957.886 7171.831
Means with the same letter are not significantly different.
Groups, Treatments and means
a  P1 20420
b  P2 9388
b  P0 5164
b  P3 3468

Hasil sidik ragam neunjukkan perbedaan sangat nyata pada perlakuan Pupuk serta terdapat interaksi antara Rambatan dan Pupuk. Uji Duncan dianggap lebih sensitif untuk uji nilai tengah rancangan 2 factorial atau lebih akrena memiliki beberapa nilai kritikal sebagai pembanding, dibanding hanya satu pada uji lain. Ini berguna iji nilai terendah dan tertinggi rada jauh. Default ji ini adalah alpha 5% (0.05), tetapi bisa disesuaikan misal 1% untuk uji lebih sensitif pada bidang kedokteran, atau 10% pada bidang ilmu sosial, sesuaikan dengan keperluan penelitian. Ganti saja nilai alpha sesuai yang dikehendaki pada script. Uji Duncan, sama seperti Uji Tukey atau LSD hanya bisa untuk uji lanjut faktor tunggal. Pada kasus ini terdapat interaksi antara Rambatan dan Pupuk. Sementara ini dahulu yang dipelajari untuk anova rancangan acak lengkap dan uji Duncan. Uji nilai tengah interaksi dibahas pada lain kesempatan.

Location: Göttingen, German

Time: May 2017

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s