Resultados

ClassificaçãoDoença

dados_pacientes <- c("moderado", "leve", "leve", "severo", "leve", "moderado", 
                     "moderado", "moderado", "leve", "leve", "severo", "leve", 
                     "moderado", "moderado", "leve", "severo", "moderado", 
                     "moderado", "moderado", "leve")

tabela <- table(dados_pacientes)

porcentagens <- round(100 * tabela / sum(tabela), 1)
rotulos <- paste(porcentagens, "%", sep="")

cores <- c("green", "blue", "salmon")

pie(tabela, 
    labels = rotulos, 
    col = cores, 
    main = "Classificação dos Estágios da Doença, ou seja n=20")

legend("topright", 
       legend = names(tabela), 
       fill = cores, 
       title = "Estágio", 
       cex = 0.8) 

Teorema

# Carrega o dataset "flu" (Gripe Espanhola na Suíça, 1918)
flu <- read.csv('C:/Users/Merti/Downloads/flu.csv')

table(flu$age)
## 
##    0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15 
## 6409 1307  723  504  445  368  368  300  248  271  233  270  238  281  335  396 
##   16   17   18   19   20   21   22   23   24   25   26   27   28   29   30   31 
##  544  620  696  840  923  944  986 1062 1076 1117 1105 1144 1166 1100 1075 1027 
##   32   33   34   35   36   37   38   39   40   41   42   43   44   45   46   47 
##  992  997  832  877  872  785  733  757  751  689  727  738  598  622  629  617 
##   48   49   50   51   52   53   54   55   56   57   58   59   60   61   62   63 
##  697  574  591  661  656  745  794  791  792  814  799  849  848  947  909  839 
##   64   65   66   67   68   69   70   71   72   73   74   75   76   77   78   79 
##  891  994  966 1039 1135 1154 1054 1190 1229 1221 1171 1189 1112 1030 1039  922 
##   80   81   82   83   84   85   86   87   88   89   90   91   92   93   94   95 
##  753  727  677  579  501  392  354  299  224  144  133  100   56   43   30   16 
##   96   97   98   99  100 
##   11    6    6    3    1
# Histograma e curva de densidade
hist(flu$age, 
     freq = FALSE, 
     col = "lavender", 
     main = "Idades das Mortes", 
     xlab = "Idade",
     breaks = 30)
lines(density(flu$age, na.rm=TRUE), col = "red", lwd = 2)

# Média por Amostragem - 200 médias com n = 35
n_amostras <- 200
tam_amostra <- 35
medias_amostrais <- numeric(n_amostras) 

set.seed(72)

for(i in 1:n_amostras) {
  amostra <- sample(flu$age, size = tam_amostra, replace = FALSE)
  medias_amostrais[i] <- mean(amostra)
}

# Histograma e curva de densidade das médias das amostras
hist(medias_amostrais, 
     probability = TRUE, 
     col = "green", 
     main = "Médias Amostras n=35", 
     xlab = "Média da Idade",
     breaks = 20)
lines(density(medias_amostrais), col = "black", lwd = 2)

VADeaths

head(VADeaths)
##       Rural Male Rural Female Urban Male Urban Female
## 50-54       11.7          8.7       15.4          8.4
## 55-59       18.1         11.7       24.3         13.6
## 60-64       26.9         20.3       37.0         19.3
## 65-69       41.0         30.9       54.6         35.1
## 70-74       66.0         54.3       71.1         50.0
# Definindo as cores dos grupos
cores <- c("pink", "orange", "lavender", "purple", "cyan")

barplot(height = VADeaths, 
        beside = TRUE,                             
        col = cores,                               
        legend.text = TRUE,                       
        main = "Taxas de Mortalidade na Virgínia (1940)", 
        xlab = "Grupos Demográficos",              
        ylab = "Taxa de Mortalidade (por 1000)",   
        ylim = c(0, 100),                          
        args.legend = list(title = "Faixa Etária", x = "topright") 
)