Ein Gesundheitsberater möchte die Normalverteilung der Zufriedenheit der Patienten für zwei Krankenhäuser mit einem Quantil-Quantil-Diagramm (QQ-Diagramm) vergleichen. QQ-Diagramm zeigen, wie gut die beide Gruppen der Bewertungen der Zufriedenheit von Patienten einer Normalverteilung folgen.
Das R-Beispielskript liest die Daten aus Spalten in Minitab. Das Skript berechnet die Quantile und erstellt ein QQ-Diagramm für jede Spalte. Anschließend überträgt endet das Skript die Plots an das Ausgabefenster von Minitab.
Alle Dateien, auf die in diesem Handbuch verwiesen wird, sind in dieser .ZIP-Datei verfügbar: r_guide_files.zip.
Verwenden Sie die folgenden Dateien, um die Schritte in diesem Abschnitt auszuführen:
| Datei | Beschreibung |
|---|---|
| qq_plot.R | Ein R-Skript, das Spalten aus einem Minitab-Arbeitsblatt ausliest und das QQ-Diagramm für jede Spalte anzeigt. |
Das R-Skript im folgenden Beispiel erfordert die folgenden R-Pakete:
- mtbr
- Das R-Paket, das Minitab und R integriert. Im Beispiel übertragen Funktionen aus diesem Modul R-Ergebnisse an Minitab.
- Stellen Sie sicher, dass Sie die erforderlichen Module installiert haben: mtbr.
- Um das mtbr-Paket zu installieren, führen Sie den folgenden Befehl im R-Terminal aus:
install.packages("url-to-mtbr")
- Um das mtbr-Paket zu installieren, führen Sie den folgenden Befehl im R-Terminal aus:
- Speichern Sie die R-Skriptdatei, qq_plot.R in Ihrem Minitab-Standardspeicherort für Dateien. Weitere Informationen dazu, wo Minitab nach R-Skriptdateien sucht, finden Sie unter Standardordner für R-Dateien für Minitab.
- Öffnen Sie den Beispieldatensatz HotelvergleichUngestapelt.MTW.
- Geben Sie im Minitab-Fenster Befehlszeile
RSCR "qq_plot.R" "Hospital A" "Hospital B"ein. - Klicken Sie auf Durchlauf.
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")Ergebnisse
R Script
These results are from external software.
