Penggunaan QQ plot (Quantile-Quantile Plot) adalah salah satu cara untuk membandingkan normalitas dua dataset yang berbeda. Dalam artikel ini, kita akan melihat bagaimana menggunakan QQ plot dalam Minitab dan R untuk menganalisis data yang diperoleh dari dua rumah sakit.

Menggunakan Minitab

Dalam Minitab, kita dapat menggunakan menu "Graph" > "Probability Plot" > "Simple" untuk membuat plot. Pilih kolom "Mahal2" sebagai variable grafik, lalu pilih distribusi gamma dan ukuran 2. Pada bagian "Scale", centang "Transpose Y and X". Plot tersebut akan muncul di area hasil.

Menggunakan R

Dalam R, kita dapat menggunakan paket mtbr untuk mengintegrasikan Minitab dengan R. Berikut adalah contoh script R yang dapat digunakan:

require(mtbr)
column_names <- commandArgs(trailingOnly = TRUE)

if (length(column_names) == 0) {
 current_index <- 1
 while (length(mtb_get_column(paste0("C",current_index))) > 0) {
 column_names[current_index] <- paste0("C", current_index)
 current_index <- current_index + 1
 }
}

if (length(column_names) == 0 || length(mtb_get_column(column_names[1])) == 0) {
 stop("Worksheet is empty or column data could not be found!\n\tPass columns to RSCR or move first column to C1.")
}

png("qqplot.png")
par(mfrow=c(ceiling(length(column_names)/round(sqrt(length(column_names)))), round(sqrt(length(column_names)))))

for (column_name in column_names) {
 column <- mtb_get_column(column_name)
 qqnorm(column, main=paste0("Normal Q-Q Plot of ", column_name))
 qqline(column)
}

graphics.off()
mtb_add_image("qqplot.png")

Script di atas akan membaca data dari Minitab dan membuat plot QQ untuk setiap kolom. Plot tersebut akan disimpan sebagai file PNG dengan nama "qqplot.png".

Hasil

Dalam hasil, kita dapat melihat bagaimana dua dataset yang berbeda memiliki normalitas yang sama atau tidak. Dengan menggunakan QQ plot, kita dapat mengetahui apakah data tersebut memenuhi syarat normalitas atau tidak.