Un consultor de salud quiere comparar la normalidad de las calificaciones de satisfacción del paciente de dos hospitales utilizando una parcela cuantilo-cuantilo (QQ). Las gráficas de QQ muestran lo bien que cada conjunto de calificaciones de satisfacción del paciente se ajusta a una distribución normal.
El script R de ejemplo lee los datos de las columnas de Minitab. El script calcula los cuantilos y crea un trazado de QQ para cada columna. A continuación, el script envía los trazados al panel Salida de Minitab.
Todos los archivos a los que se hace referencia en esta guía están disponibles en este archivo .ZIP: .
Utilice el siguiente archivo para llevar a cabo los pasos de esta sección:
| Archivo | Descripción |
|---|---|
| qq_plot.R | Un script R que toma columnas de una hoja de cálculo minitab y muestra el trazado QQ para cada columna. |
El script R en el ejemplo siguiente requiere los siguientes paquetes R
- mtbr
- El paquete R que integra Minitab y R. En el ejemplo, las funciones de este módulo envían resultados R a Minitab.
- Asegúrese de haber instalado los módulos necesarios: mtbr.
- Para obtener información sobre cómo instalar el paquete R de Minitab, vaya al paso 2: Instale mtbr.
- Guarde el archivo R script qq_plot.RR, en la ubicación predeterminada del archivo De Minitab. Para obtener más información sobre dónde Minitab busca archivos de script R vaya a Carpetas predeterminadas para archivos R para Minitab.
- Abra el conjunto de datos de la muestra
- En el Línea de
comandos panel Minitab, escriba
RSCR "qq_plot.R" "Hospital A" "Hospital B" - Haga clic en Corrida.
qq_plot. R
library(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")
Resultados
R Script
These results are from external software.
