require(tidyverse)

## Carregando pacotes exigidos: tidyverse

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.0     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

Criando uma sequencia de números

a=seq(1,10,1)
a

##  [1]  1  2  3  4  5  6  7  8  9 10

b=seq(-5,5,2)
b

## [1] -5 -3 -1  1  3  5

• outro modo de criar a sequência, sendo apenas de 1 em 1 (intervalo)

a=1:10
a

##  [1]  1  2  3  4  5  6  7  8  9 10

b=-5:5
b

##  [1] -5 -4 -3 -2 -1  0  1  2  3  4  5

Importar o banco de dados:

• este banco de dados é referente a 16 pacientes submetidos a baixa e alta ingestão de sódio, e foi medida a pressão arterial nestes dois contextos.

dados do exemplo abrir um novo chunk: Ctrl + Alt + I

require(DT) #pacote datatable

## Carregando pacotes exigidos: DT

dados=read.table("dados.txt",sep=";",head=TRUE)
datatable(dados,caption="Tabela 1",class = 'cell-border order-column compact hover') 

dados

##    Paciente baixo_teor_sodio alto_teor_sodio
## 1         1              138             143
## 2         2              147             154
## 3         3              146             147
## 4         4              154             147
## 5         5              142             157
## 6         6              156             158
## 7         7              134             164
## 8         8              146             156
## 9         9              143             151
## 10       10              175             182
## 11       11              156             151
## 12       12              117             116
## 13       13              157             154
## 14       14              143             149
## 15       15              127             126
## 16       16              134             138
## 17       17              112             115
## 18       18              144             159
## 19       19              117             124
## 20       20              128             125

Para construir um boxplot com grupos, será necessário dispor o banco de dados em formato longitudinal, com uma coluna de teor (baixo x alto), e uma coluna de pressão.

require(readxl)

## Carregando pacotes exigidos: readxl

dados=readxl::read_excel("banco_modificado.xlsx")

dados$teor=as.factor(dados$teor)

Em geral, os valores de pressão arterial são inferiores com baixa ingestão de sódio, no entanto é necessário realizar um teste de comparação de médias (que no caso é pareado):

• Teste t pareado (verificar as pressuposições);
• Teste não paramétrico de Wilcoxon;
• Como melhorias (tarefas):
  • Colocar os níveis baixo e alto, nessa ordem, no eixo x do boxplot;
  • Construir um gráfico de perfis, próprio para dados pareados (em que o mesmo indivíduo aparece com duas medições de pressão (baixo e alto nível de sódio).

ggplot(dados, aes(x=teor, y=pressao, fill=teor)) + 
    geom_boxplot()