한 의료 컨설턴트는 QQ(quantile-quantile) Plot을 사용하여 두 병원의 환자 만족도 등급의 정규성을 비교하고자 합니다. QQ Plot은 각 환자 만족도 등급 집합이 일반 분포를 얼마나 잘 적합하는지 보여줍니다.
예제 R 스크립트는 Minitab의 열에서 데이터를 읽습니다. 스크립트는 분위수를 계산하고 각 열에 대한 QQ Plot을 만듭니다. 그런 다음 스크립트가 이 플롯을 Minitab 출력 창으로 보냅니다.
이 가이드에서 참조하는 모든 파일은 .ZIP 파일인 r_guide_files.zip 파일로 제공됩니다.
다음 파일을 사용하여 이 섹션의 단계를 수행합니다.
| 파일 | 설명 |
|---|---|
| qq_plot.R | Minitab 워크시트에서 열을 가져와 각 열에 대한 QQ Plot을 표시하는 R 스크립트입니다. |
아래 예제의 R 스크립트에는 다음 R 패키지가 필요합니다.
- mtbr
- Minitab과 R을 통합하는 R 패키지입니다. 예제에서는 이 모듈의 함수가 R 결과를 Minitab으로 보냅니다.
- 필요한 모듈인 mtbr 파일로 제공됩니다.
- mtbr 패키지를 설치하려면 R 터미널에서 다음 명령을 실행합니다.
install.packages("url-to-mtbr")
- mtbr 패키지를 설치하려면 R 터미널에서 다음 명령을 실행합니다.
- R 스크립트 파일 qq_plot.R를 Minitab 기본 파일 위치에 저장합니다. Minitab이 R 스크립트 파일을 찾는 위치에 대한 자세한 내용은 Minitab용 R 파일의 기본 폴더로 이동하십시오.
- 표본 데이터 세트를 엽니다 호텔비교(분할된데이터).MTW.
- Minitab 명령줄 창에서
RSCR "qq_plot.R" "Hospital A" "Hospital B"을 입력합니다. - 실행을 클릭합니다.
qq_plot. R
require(mtbr, quietly=TRUE) column_names <- commandArgs(trailingOnly = TRUE) if (length(column_names) == 0) { current_index <- 1 while (length(mtbr::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(mtbr::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 <- mtbr::mtb_get_column(column_name) qqnorm(column, main=paste0("Normal Q-Q Plot of ", column_name)) qqline(column) } graphics.off() mtbr::mtb_add_image("qqplot.png")결과
R Script
These results are from external software.
