Cara Membuat Error Bar (Standar Error) di SPSS 21

Diagram error bar seringkali ditemukan di berbagai jurnal sains terutama yang berkaitan dengan jurnal Biologi. Diagram error bar dibangun dari peubah numerik (peubah pertama). Pada kasus percobaan satu faktor, peubah kedua membagi peubah pertama menjadi kelompok (misalnya, Jenis pupuk atau varietas tanaman). Sedangkan dalam kasus percobaan dua faktor, peubah ketiga dapat digunakan untuk membagi sub kelompok.

Mengapa diagram error bar penting?

Karena Error bar dapat membantu menjelaskan dan menginterpretasikan data dengan benar. Umumnya Error Bar menunjukkan range, interval kepercayaan, standar error (galat baku), standar deviasi (Geoff Cumming, dkk, 2007).

Error Bar	Rumus	Deskripsi
Range	Data terbesar dikurangi data terkecil	Range (jangkauan) merupakan ukuran penyebaran data yang mencakup nilai terendah dan tertinggi.
Standar deviasi		Simpangan rata-rata antara data observasi dan rataannya.
Standar Error		Mencerminkan keakuratan contoh yang dipilih terhadap populasinya.
Interval Kepercayaan (95%)		Rentang nilai pada selang kepercayaan 95% diyakini memuat nilai rataan sebenarnya

(Geoff Cumming, dkk, 2007).

Jenis error bar yang mana sebaiknya digunakan?

Pada penelitian rancangan percobaan biasanya mencoba membandingkan hasil eksperimen dengan kontrol, biasanya tepat digunakan diagram error bar, seperti SE atau Interval Kepercayaan, daripada SD. Namun, jika n sangat kecil (misalnya n = 3), daripada menampilkan error bar dan statistik, lebih baik menggunakan plot sederhana.

Langkah Membuat Bagan Error Bar (SE) di SPSS 21

1.       Graphs -> klik Chat Builder -> OK

2.     Ikuti langkah di bawah sampai selesai

Jika ingin menampilkan error bar (standar deviasi atau interval kepercayaan cukup ganti pada langkah no.5 pada gambar di atas

3.       Selesai.  

Diagram error bar (standar error) dapat digunakan untuk mengetahui apakah perbedaan antar perlakuan signifikan atau tidak. Kita bisa melihat bagan error bar (standar error) yang saling tumpang tindih (overlap) atar bagan error bar.

Gambar 1. Jika standar error bar tumpang tindih cukup banyak, itu menandakan bahwa perbedaannya tidak signifikan secara statistik. Namun, teman-teman tatap harus melakukan uji statistik untuk menarik kesimpulan.

Gambar 2. Jika standar error bar tumpang tindih sedikit, itu menandakan bahwa perbedaannya mungkin tidak signifikan secara statistik. Karena gambar 2 sulit ditentukan apakah signifikan atau tidak maka teman-teman harus melakukan uji statistik untuk menarik kesimpulan.

Gambar 3. Jika standar error bar tidak tumpang tindih, itu menandakan bahwa perbedaannya mungkin signifikan, tetapi jika teman-teman tidak yakin. Teman-teman bisa melakukan uji statistik untuk menarik kesimpulan.

Jadi, diagram standar error bar hanya memberikan petunjuk bahwa perbedaan suatu perlakuan signifikan atau tidak. Sehingga tetap diperlukan uji statistik untuk menarik kesimpulan yang valid. 

Sumber:

Geoff Cumming, dkk, 2007. Error bars in experimental biology. The Journal of Cell Biology. 177: 7–11.

https://www.biologyforlife.com/interpreting-error-bars.html

RANCANGAN ACAK LENGKAP(RAL) MENGGUNAKAN MINITAB : UJI LANJUT TUKEY DAN INTERPRETASI PLOT

Pada tulisan sebelumnya, saya telah membahas RAL menggunakan SPSS 21. Kali ini masih topik bahasan yang sama namun menggunakan soal dan software yang berbeda. Dan tentu saja saya memodifikasi soal di bawah agar menghasilkan hasil analisis yang berbeda dari soal artikel sebelumnya. Kalau pada link RAL SPSS 21 lebih banyak menggunakan uji exact, artikel ini lebih banyak  menampilkan grafik dan bagaimana menginterpretasikan grafik.

Soal di bawah, saya peroleh dari modul Dr. Jianbiao (John) Pan yang berjudul Minitab Tutorial for Design and Analysis of Experiment. Namun saya telah mengubah sebagian data pada tingkat etsa yang memiliki tingkat daya RF 180W

Masalah: Dalam industri mekanik dan elektronik banyak menggunakan proses Etsa Plasma (plasma etching). Seorang peneliti tertarik menyelidiki hubungan antara pengaturan daya RF dan tingkat etsa. Dia tertarik pada gas khusus (C₂F₆) dan celah (0,80 cm) dan ingin menguji empat tingkat daya RF: 160W, 180W, 200W, dan 220W. Percobaan dilakukan dengan lima kali ulangan.
                                                      Berikut data hasil percobaan:

Ulangan	Tingkat Daya RF
	160 W	180 W	200 W	220 W
1	54.2	57.9	65.1	71.9
2	53	55	62.9	72.5
3	57	57	60	71.5
4	57.5	56.3	63.7	68.5
5	53.9	56.5	61	70

Langkah-langkah menggunakan Minitab untuk uji asumsi, ANOVA dan uji Lanjut

Menampilkan grafik untuk uji asumsi

Untuk Menampilkan output Uji lanjut 

Langkah terakhir klik OK.

Sebelum melakukan analisis varian dan uji lanjut, kita perlu mengecek terlebih dahulu apakah asumsi-asumsi pada RAL telah terpenuhi. Mengapa harus mengecek asumsi? Mungkin teman-teman pernah menemukan kasus seperti hasil uji ANOVA tidak signifikan tetapi ternyata hasil diuji lanjut menunjukkan ada yang signifikan. Padahal seharusnya jika hasil uji ANOVA tidak signifikan maka hasil uji lanjut tidak ada yang signifikan.  Mengapa hal tersebut bisa terjadi? Hal itu terjadi kemungkinan disebabkan ada asumsi tidak terpenuhi. Oleh karena itu, kita akan mengecek terlebih dahulu asumsi-asumsi berikut:

1. Model bersifat Aditif

2. Sisaan saling bebas satu dengan lainnya.

3. Sisaan menyebar normal

3. Sisaan memiliki varian yang sama

• Untuk mengecek apakah sisaan saling bebas maka kita dapat melihat plot sisaan vs urutan amatan seperti berikut:

Plot di atas menunjukkan bahwa sebaran sisaan tidak membentuk pola tertentu sehingga dapat disimpulkan bahwa sisaan saling bebas.

• Untuk mengecek apakah sisaan menyebar normal

Plot Normal Probability merupakan normal plot sisaan yang bertujuan untuk mengecek apakah asumsi sisaan menyebar normal. Plot di atas menunjukkan bahwa hampir semua titik-titik sisaan pada plot mendekati dan mengikut garis lurus. Oleh karena itu , dapat dikatakan bahwa sisaan menyebar normal. Selain itu, untuk lebih meyakinkan dapat dicek secara exact yaitu menggunakan uji Anderson Darling, Ryan-Joiner, Kolmogorov-smirnov ataupun shapiro-wilk.

Langka-langkah Uji Normalitas di Minitib

Stat > Basic Statistics > Normality test > maka akan muncul seperti berikut:

Maka akan muncul plot di bawah lengkap dengan hasil uji exact.

Interpretasi output Uji Normalitas

Hipotesis:

Ho : Sisaan menyebar normal

H1 : Sisaan tidak menyebar normal

Berdasarkan uji normalitas menggunakan Anderson-Darling diperoleh p- value 0.332 > 0.05 sehingga dapat disimpulkan bahwa Tolak H1, artinya dengan tingkat kepercayaan 95% belum cukup bukti untuk menyatakan bahwa sisaan tidak menyebar normal. (versi orang statistik )

Berdasarkan uji normalitas menggunakan Anderson-Darling diperoleh p- value 0.332 > 0.05 sehingga dapat disimpulkan bahwa Terima Ho, artinya sisaan menyebar normal. (versi orang awam).

Mengapa ada seperti itu? Bagikan jawaban anda di kolom komentar.

• Untuk mengecek apakah sisaan memiliki varian yang sama

Plot di atas bertujuan untuk mengecek apakah sisaan memiliki varian yang sama atau konstan. salah satu cara mengeceknya adalah dengan melihat lebar pitanya. Plot di atas menunjukkan lebar pita yang sama sehingga dapat disimpulkan bahwa sisaan memiliki varian yang sama.

Output ANOVA

Analysis of Variance

Source   DF  Adj SS   Adj MS  F-Value  P-Value

Daya RF   3  773.65  257.885    85.83    0.000

Error    16   48.08    3.005

Total    19  821.73

Hipotesis:

Ho : Tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 

H1:  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 

Berdasarkan tabel ANOVA di atas menunjukkan bahwa p-value 0.000 < 0.05. Oleh karena itu, disimpulkan bahwa Terima H1 artinya  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF.

Uji ANOVA merupakan syarat untuk melakukan uji lanjut, jika tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF maka uji lanjut tidak perlu dilakukan. 

Model Summary

      S    R-sq  R-sq(adj)  R-sq(pred)

1.73342  94.15%     93.05%      90.86%

Dalam RAL R-square ataupun R-square Adjusted tidak diinterperetasikan. Mengapa? Karena konteks soal ini bukan untuk memprediksi.

Means

Daya RF  N    Mean  StDev       95% CI

160      5  55.120  2.002  (53.477, 56.763)

180      5  56.480  1.085  (54.837, 58.123)

200      5  62.540  2.053  (60.897, 64.183)

220      5  70.880  1.619  (69.237, 72.523)

Pooled StDev = 1.73342

Tukey Pairwise Comparisons

Grouping Information Using the Tukey Method and 95% Confidence

Daya RF  N    Mean  Grouping

220      5  70.880  A

200      5  62.540    B

180      5  56.480      C

160      5  55.120      C

Means that do not share a letter are significantly different.

  

Uji lanjut dalam pada kasus ini digunakan uji tukey atau dikenal dengan uji Beda Nyata Jujur (BNJ). Hasil uji Tukey pada tabel Tukey Pairwise Comparisons sama saja dengan plot di atas. Oleh karena itu, saya akan menjelaskan salah satunya saja yaitu dengan menginterpretasikan plot hasil uji tukey.

Cara membaca plot di atas adalah

·        Garis yang menyentuh titik nol artinya tidak ada perbedaan

·        Garis yang tidak menyentuh titik nol artinya ada perbedaan  

Hipotesis 1:

Ho : Tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 180W

H1:  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 180W

·         Berdasarkan plot di atas bahwa garis 160W dan 180W menyentuh titik nol sehingga dapat disimpulkan bahwa Tolak H1, artinya dengan tingkat kepercayaan 95% belum cukup bukti untuk menyatakan  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 180W (versi orang statistik hehehe)

·         Berdasarkan plot di atas bahwa garis 160W dan 180W menyentuh titik nol sehingga dapat disimpulkan bahwa terima Ho, artinya tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 180W (versi orang awam)

Hipotesis 2:

Ho : Tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 200W

H1:  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 200W

Berdasarkan plot di atas bahwa garis 160W dan 200W tidak menyentuh titik nol sehingga dapat disimpulkan bahwa terima H1, artinya ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 200W

Hipotesis 3:

Ho : Tidak ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 220W

H1:  ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 220W

Berdasarkan plot di atas bahwa garis 160W dan 2200W tidak menyentuh titik nol sehingga dapat disimpulkan bahwa terima H1, artinya ada perbedaan rata-rata tingkat etsa dengan pada pengaturan daya RF 160W dan 220W

Silahkan teman-teman membuat Hipotesis ke 4 sampai ke 6 caranya masih sama dengan yang di atas. Jika ada pertanyaan dan koreksi silahkan disampaikan di kolom komentar. Mari saling belajar.

Sumber:

1.https://support.minitab.com/en-us/minitab/18/help-and-how-to/modeling-statistics/regression/supporting-topics/residuals-and-residual-plots/residual-plots-in-minitab/
2.http://www.philender.com/courses/linearmodels/cnotes1/cr1a.html

RANCANGAN ACAK LENGKAP (RAL) ULANGAN TIDAK SAMA

Pak Tono seorang peternak itik, dia ingin melakukan sebuah percobaan untuk membandingkan 4 jenis ransum terhadap pertumbuhan bobot itik. Pak Tono memiliki beberapa kandang itik. Setiap kandang berisi 12 ekor itik dan diberi satu jenis ransum. Oleh karena itu, ulangannya adalah satu kandang itik.

Berikut adalah data rata-rata bobot itik. (Data Karangan)

Ulangan	Perlakuan				Total
	A	B	C	D
1	6.4	4	5.5	7.2	23.1
2	6.9	4.2	5.7	7.4	24.2
3	6.8		6	7.3	20.1
4				7.8	7.8
Total	20.1	8.2	17.2	29.7	75.2
Rata-rata	6.70	4.10	5.73	7.43

Data hasil percobaan dengan ulangan tidak sama pada tabel di atas juga dapat dikerjakan di SPSS dan tentu saja hasilnya akan sama yang diperoleh dengan langkah manual di bawah. Berikut cara menginput data di atas pada SPSS 21. Langkah analisisnya sama halnya untuk data dengan ulangan sama seperti di link RAL Ulangan Sama Pada SPSS 21

Tabel ANOVA

Sumber keragaman	Derajat bebas	Jumah Kuadrat	Kuadrat Tengah	F Hitung
Perlakuan	3	16,173	5,391	87,374
Galat	8	0,494	0,0617
Total	11	16,667

     Hipotesis:

Karena F hitung > Ftabel (87,374>4,07) maka tolak Ho. Artinya pada taraf  alpa 0,05, cukup bukti untuk menyatakan bahwa ada jenis ransum berpengaruh terhadap bobot ayam.

Kemudian kita ingin melihat apakah ransum A berbeda dengan Ransum B menggunakan uji BNT?

Hipotesis:

Statistik Uji BNT

Statistik uji dibandingkan dengan selisih

Kesimpulan:

Tolak Ho artinya sudah cukup bukti untuk menyatakan bahwa rata-rata bobot itik dengan pemberian ransum A dan ransum B berbeda , dalam hal ini nilai rata-rata bobot itik yang diberi ransum A lebih berat daripada yang diberi ransum B.

Pada uji BNT perlu diperhatikan pasangan yang dibandingkan. Uji BNT hanya benar jika satu perlakuan yang dibandingkan. Karena semakin banyak kombinasi pasangan yang bandingkan maka alpa akan membesar sehingga berakibat pada kesalahan jenis 1. Resiko kesalahan pada saat menolak Ho menjadi sangat besar. Oleh karena itu, untuk menghindari kesalahan tersebut bisa digunakan uji turkey atau BNJ. 

Teman-teman boleh mencoba sendiri untuk melihat misalnya ransum A dengan Ransum C, Ransum D. langkahnya sama seperti di atas.