Um consultor de saúde quer comparar a normalidade dos índices de satisfação dos pacientes de dois hospitais usando um gráfico quantil-quantil (QQ). Os gráficos QQ mostram em que medida cada conjunto de índices de satisfação do paciente ajusta bem uma distribuição normal.
O script R de exemplo lê os dados das colunas no Minitab. O script calcula os quantis e cria um gráfico QQ para cada coluna. Em seguida, o script envia os gráficos para o painel Saída do Minitab.
Todos os arquivos mencionados neste guia estão disponíveis neste arquivo .ZIP : r_guide_files.zip.
Use o seguinte arquivo para executar as etapas nesta seção:
| Arquivo | Descrição |
|---|---|
| qq_plot.R | Um script R que usa colunas de uma worksheet do Minitab e exibe o Gráfico QQ para cada coluna. |
O script R no exemplo abaixo requer os seguintes pacotes R:
- mtbr
- O pacote R que integra o Minitab e o R. No exemplo, as funções deste módulo enviam resultados do R para o Minitab.
- Instale os módulos necessários: mtbr.
- Para instalar o pacote mtbr, execute o seguinte comando no terminal R:
install.packages("url-to-mtbr")
- Para instalar o pacote mtbr, execute o seguinte comando no terminal R:
- Salve o arquivo de script do R, qq_plot.R, no local do arquivo padrão do Minitab. Para obter mais informações sobre onde o Minitab procura os arquivos de script do R, acesse Pastas padrão de arquivos R para o Minitab.
- Abra o conjunto de dados de amostra ComparaçãoDesempilhadaDeHotéis.MTW.
- No painel Linha de
comando do Minitab, digite
RSCR "qq_plot.R" "Hospital A" "Hospital B". - Clique em Ensaio.
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")Resultados
R Script
These results are from external software.
