Suplemen Praktikum
Pertemuan
ANALISIS DATA KATEGORIK (STK351)
10+
Departemen Statistika – FMIPA IPB Pokok Bahasan
Sub Pokok Bahasan
Referensi
Analisis Nonparameterik dan Data Kategorik dengan SAS dan MINITAB
Menggunakan SAS dan MINITAB untuk memecahkan beberapa persoalan dalam analisis nonparametrik dan data kategorik
SAS Institute (2009)
Waktu
Jumat 7 Januari 2011 05.00 – 18.00
Data : IPK * No.
Nama Peubah
Tipe
Label
Nilai
1 2 3 4 5 6 7 8 9
NoResp Tahun Jkelamin JMasuk Org UgBulan IpkTPB IPK KelIPK
Char Char Char Char Char Num Num Num Char
No Responden Tahun Masuk Jenis Kelamin Jalur Masuk Aktivitas Organisasi Uang Bulanan (Rp.000) IPK TPB IPK Terakhir Kategori IPK
MHS 2003, 2004, 2005,2006 Putra, Putri BUD, SPMB, USMI Aktif, Tidak 200 s.d. 1600 2.58 s.d. 4.00 2.26 s.d. 3.88 Baik (>=3), Kurang (<3)
Data : Ramal * No.
Nama Peubah
Tipe
Label
Nilai
1 2 3
Warna Ramal Frek
Char Char Num
Warna Sebenarnya Warna Ramalan
Hitam, Putih Hitam, Putih 2, 3
*) Data dapat diunduh melalui tautan : http://pendek.in/02yub atau http://bit.ly/h9NKku
Sebelum melanjutkan buatlah SAS-Library baru, misalnya dengan nama ADK, yang merujuk pada direktori data Anda. Sebagai contoh, data Anda di simpan pada D:\Data, maka library ADK dapat dibuat dengan perintah : libname ADK "D:\Data";
UJI TANDA SATU CONTOH Data : IPK Peubah : IPK SAS ods select TestsForLocation; proc univariate data=adk.ipk mu0=3.175; var ipk; run;
MINITAB Stat Nonparametrics 1-Sample Sign…
Catatan : Dapat juga menggunakan PROC CAPABILITY
UJI PERINGKAT BERTANDA WILCOXON Data : IPK Peubah : IPK SAS ods select TestsForLocation; proc univariate data=adk.ipk mu0=3.175; var ipk; run;
MINITAB Stat Nonparametrics 1-Sample Wilcoxon…
Catatan : Dapat juga menggunakan PROC CAPABILITY
2 / 10
UJI PENGAMATAN BERPASANGAN Data : IPK Peubah : diffIPK = IPK – ipkTPB SAS data ipk; set adk.ipk; difIPK = ipk - ipkTPB; run;
MINITAB Stat Nonparametrics Pairwise Differences…
Catatan : Perintah di atas digunakan untuk membuat peubah baru (diffIPK) yang merupakan selisih dua peubah. Jalankan PROC UNIVARIATE atau PROC CAPABILITY untuk melakukan analisis bagi peubah diffIPK
Lakukan prosedur 1-Sampe Sign… atau 1-Sample Wilcoxon… menggunakan peubah pada C9
UJI MANN-WHITNEY Data : IPK Peubah respon : IPK
Faktor : JKelamin
SAS proc npar1way data=adk.ipk wilcoxon; class jkelamin; var ipk; run;
MINITAB Pisahkan terlebih dahulu nilai peubah menjadi 2 kolom berdasarkan taraf faktor, misalnya IPKPutra dan IPK-Putri. Stat Nonparametrics Mann-Whitney…
Catatan : Kata kunci WILCOXON digunakan untuk analisis-analisis berbasis peringkat Wilcoxon. Jika faktor terdiri dari dua taraf, kata kunci ini akan melakukan uji Mann-Withney
3 / 10
UJI MEDIAN Data : IPK Peubah respon : IPK
Faktor : JKelamin
SAS proc npar1way data=adk.ipk median; class jkelamin; var ipk; run;
MINITAB Stat Nonparametrics Mood’s Median Test…
Catatan : Kata kunci MEDIAN digunakan untuk melakukan uji Median (Moods)
UJI KRUSKAL-WALLIS Data : IPK Peubah respon : IPK
Faktor : JMasuk
SAS proc npar1way data=adk.ipk wilcoxon; class jmasuk; var ipk; run;
MINITAB Stat Nonparametrics Kruskal-Wallis…
Catatan : Kata kunci WILCOXON digunakan untuk analisis-analisis berbasis peringkat Wilcoxon. Jika faktor terdiri dari tiga atau lebih taraf, kata kunci ini akan melakukan uji Kruskal-Wallis
4 / 10
UJI KOLMOGOROV FRIEDMAN Data : IPK Peubah : IPK
Faktor : JMasuk
Kelompok : JKelamin
MINITAB Stat Nonparametrics Friedman…
UJI KEBAIKAN SUAI KHI-KUADRAT (SEBARAN KONTINU) Data : IPK Peubah : IPK SAS ods select ParameterEstimates GoodnessOfFit; proc capability data=adk.ipk; var ipk; histogram / normal (mu=3.175 sigma=0.5); run;
Catatan : Digunakan apabila mean dan sigma diketahui. Apabila tidak diketahui, hilangkan kata kunci MU= dan SIGMA=. Sebagai alternatif dapat menggunakan PROC UNIVARIATE Sebaran yang dapat diuji : NORMAL, BETA, EXPONENTIAL, GAMMA, LOGNORMAL, WEIBULL
5 / 10
UJI KOLMOGOROV SMIRNOV SATU CONTOH Data : IPK Peubah : IPK SAS ods select ParameterEstimates GoodnessOfFit; proc capability data=adk.ipk; var ipk; histogram / normal (mu=3.175 sigma=0.5); run;
MINITAB Untuk menguji kenormalan data, nilai tengah dan simpangan baku tidak diketahui Stat Basic Statistic Normality Test…
Catatan : Apabila mean dan sigma tidak diketahui, hilangkan kata kunci MU= dan SIGMA=. Alternatif gunakan PROC UNIVARIATE Sebaran yang dapat diuji : NORMAL, BETA, EXPONENTIAL, GAMMA, LOGNORMAL, WEIBULL
UJI KOLMOGOROV SMIRNOV DUA CONTOH Data : IPK Peubah : IPK
Faktor : JKelamin
SAS proc npar1way data=adk.ipk edf; class jkelamin; var ipk; run;
6 / 10
UJI KORELASI SPEARMAN DAN TAU-KENDALL Data : IPK Peubah 1 : IPK SAS proc corr data=adk.ipk spearman kendall; var ipk ipkTPB; run;
Peubah 2 : IpkTPB MINITAB Stat Tables Cross Tabulation and Chi Square… Other Stat
Catatan : Dapat pula dengan perintah : proc corr data=adk.ipk spearman kendall; var ipk; with ipkTPB; run;
UJI KORELASI SPEARMAN DAN TAU-KENDALL PARSIAL Data : IPK Peubah 1 : IPK SAS proc corr data=adk.ipk spearman kendall; var IPK UgBulan; partial IpkTPB; run;
Peubah 2 : UgBulan
Konstan : IpkTPB
Catatan : Dapat pula dengan perintah : proc corr data=adk.ipk spearman kendall; var IPK; with UgBulan; partial IpkTPB; run;
7 / 10
UJI PROPORSI BINOMIAL DAN SELANG KEPERCAYAAN Data : IPK Peubah 1 : KelIPK
Peubah 2 : UgBulan
SAS proc freq data=adk.ipk; tables KelIPK / binomial (p=0.8 level='Baik (>=3)' all); run;
Konstan : IpkTPB
MINITAB Stat Basic Statistic 1 Proportion… Option : cek Use test and CI base on normal distribution
Catatan : Dapat pula dengan perintah : Level : taraf faktor yang akan diuji All : menampilkan semua jenis selang kepercayaan
TABEL KONTINGENSI Data : IPK Baris : JKelamin
Kolom : KelIPK
SAS proc freq data=adk.ipk; tables JKelamin*KelIPK; run;
Layer : Org MINITAB Stat Tables Cross Tabulation and Chi Square…
Catatan : Untuk menambahkan layers : proc freq data=adk.ipk; tables Org*JKelamin*KelIPK ; run; Gunakan kata kunci EXPECTED dan CELLCHI2 untuk menampilkan nilai frekuensi harapan dan khi-kuadrat pada setiap sel. Atau dapat juga gunakan kata kunci NOCOL atau NOROW untuk menghilangkan persentase kolom dan baris pada tabel.
Untuk menambahkan layer, masukkan Org pada input For layers
8 / 10
TABEL KONTINGENSI (ODDS RATIO DAN RELATIVE RISK) Data : IPK Baris : JKelamin
Kolom : KelIPK
SAS proc freq data=adk.ipk; tables JKelamin*KelIPK / relrisk; run;
TABEL KONTINGENSI (UJI KHI KUADRAT DAN LIKELIHOOD RATIO) Data : IPK Baris : JMasuk
Kolom : KelIPK
SAS proc freq data=adk.ipk; tables JMasuk*KelIPK / chisq; run;
MINITAB Stat Tables Cross Tabulation and Chi Square… Chi Square…
TABEL KONTINGENSI (UJI POLA LINEAR) Data : IPK Baris : Tahun
Kolom : KelIPK
SAS proc freq data=adk.ipk; tables Tahun*KelIPK / cmh1 score=ridit; run;
Catatan CMH1 merupakan kata kunci untuk pengujian pola linear dengan metode Cochran-Mantel-Haenszel
SCORE=RIDIT merupakan kata kunci agar pengujian menggunakan skor mid-rank.
9 / 10
TABEL KONTINGENSI (FISHER EXACT TEST) Data : RAMAL Baris : Warna
Kolom : Ramal
SAS proc freq data=adk.ramal; table Warna*Ramal / chisq norow; weight Frek; run;
Catatan WEIGHT digunakan karena data harus diboboti terlebih dahulu Untuk dimensi tabel yang lebih dari 2 x 2, gunakan pilihan FISHER
TABEL KONTINGENSI TIGA ARAH Data : IPK Baris : JKelamin
Kolom : KelIPK
SAS proc tabulate data=adk.ipk; class Org JKelamin KelIPK; tables (Org*JKelamin all), (KelIPK all); run;
Layer : Org Catatan Gunakan CLASS untuk mendefinisikan faktor (peubah kategorik) Kata kunci ALL berfungsi untuk menampilkan frekuensi pada setiap baris dan/atau kolom
Agar output SAS Anda dapat disimpan secara otomatis dalam file *.rtf, misalnya pada direktori D:\Data dengan nama SAS-Out.rtf, apitlah perintah-perintah SAS yang anda tulis dengan pernyataan berikut : ods rtf file="D:\Data\SAS-Out.rtf"; ---- perintah SAS ---ods rtf close;
CMIIW : Correct me if I am wrong …
10 / 10