2VAProva

VADeaths

    library(plotly)
    library(dplyr)
    library(tidyr)
    library(tibble)

    df_VADeaths <- as.data.frame(VADeaths) %>%
      rownames_to_column(var = "Faixa Etária") %>%
      pivot_longer(cols = -`Faixa Etária`, names_to = "Categoria", values_to = "Taxa de Mortalidade")

    graficoVADeaths <- plot_ly(
      data = df_VADeaths,
      x = ~`Faixa Etária`, 
      y = ~`Taxa de Mortalidade`, 
      color = ~Categoria,
      type = 'bar',
      text = ~sprintf("%.1f", `Taxa de Mortalidade`),
      textposition = 'inside',
      insidetextfont = list(color = 'white', size = 12, family = 'Arial')
    )

    graficoVADeaths <- graficoVADeaths %>% layout(
      barmode = 'stack',
      title = list(text = "<b> Taxas de Mortalidade na Virgínia (1940)</b>", x = 0.5),
      xaxis = list(title = "<b> Faixa Etária</b>"),
      yaxis = list(title = "<b> Taxa de Mortalidade</b>"),
      
      legend = list(
        title = list(text = '<b> Categoria</b>'),
        orientation = "v",                      
        x = 1.02,                               
        y = 1,                                  
        bgcolor = 'rgba(255, 255, 255, 0.8)',   
        bordercolor = 'Gainsboro',              
        borderwidth = 2                         
      )
    )
graficoVADeaths

ClassificaçãoDoença

dados <- factor(c("Moderado", "Leve", "Leve", "Severo", "Leve", "Moderado", 
                      "Moderado", "Moderado", "Leve", "Leve", "Severo", "Leve", 
                      "Moderado", "Moderado", "Leve", "Severo", "Moderado", 
                      "Moderado", "Moderado", "Leve"),
                    levels = c("Leve", "Moderado", "Severo"))

    df_contagem <- as.data.frame(table(dados))
    names(df_contagem) <- c("Estagio", "Contagem")

    cores <- c("#56B4E9", "#E69F00", "#D55E00")

    graficoDoenca <- plot_ly(
      data = df_contagem,
      labels = ~Estagio,
      values = ~Contagem,
      type = 'pie',
      textinfo = 'label+percent',
      marker = list(colors = cores,
                    line = list(color = '#FFFFFF', width = 2))
    )
    graficoDoenca <- graficoDoenca %>% layout(
      title = list(text = "<b>Distribuição dos Estágios da Doença em 20 Pacientes</b>", x = 0.5),
      legend = list(
        title = list(text = '<b>Estágios</b>'),
        orientation = "v",
        x = 1, y = 0.9,
        bordercolor = 'Gainsboro',
        borderwidth = 2            
      )
    )
graficoDoenca

Teorema

Histograma e a curva de densidade do conjunto de dados “flu”.

    library(gss)
    library(plotly)
    
    flu <- read.csv("https://www.dropbox.com/scl/fi/bvf1mhw33x4h6lvtty3ks/flu.csv?authuser=0&rlkey=e9kreupfbwrfhc3425tm3dq32&e=1&dl=1")
    populacao <- flu$age

    densidade_pop <- density(populacao)

    plot_ly(x = ~populacao, type = 'histogram', histnorm = 'probability density', 
            name = 'Histograma da População', alpha = 0.7) %>%
      add_trace(x = ~densidade_pop$x, y = ~densidade_pop$y, type = 'scatter', mode = 'lines', fill = 'tozeroy', name = 'Curva de Densidade', yaxis = 'y1') %>%
      layout(
        title = "<b>Distribuição da População: Idade das Mortes por Gripe (1918)</b>",
        xaxis = list(title = "Idade"),
        yaxis = list(title = "Densidade",tickformat = ".3f")
      )

Medias de amostras da população

    n_amostras <- 200
    tamanho_amostra <- 35
    
    medias_amostrais <- replicate(n_amostras, mean(sample(populacao, size = tamanho_amostra, replace = TRUE)))
    head(medias_amostrais)

## [1] 30.85714 47.94286 37.68571 42.25714 43.82857 45.02857

Histograma com a curva de densidade para a médias das amostras

    densidade_amostras <- density(medias_amostrais)

    media_pop <- mean(populacao)
    dp_pop <- sd(populacao)
    dp_amostras_teorico <- dp_pop / sqrt(tamanho_amostra)

    x_teorico <- seq(min(medias_amostrais), max(medias_amostrais), length.out = 200)
    y_teorico <- dnorm(x_teorico, mean = media_pop, sd = dp_amostras_teorico)

    plot_ly(x = ~medias_amostrais, type = 'histogram', histnorm = 'probability density', name = 'Histograma das Médias', alpha = 0.7) %>%
      add_trace(x = ~densidade_amostras$x, y = ~densidade_amostras$y, type = 'scatter', mode = 'lines', fill = 'tozeroy', name = 'Curva de Densidade (Dados)')  %>%
      layout(
        title = "<b>Distribuição das Médias Amostrais (n=35)</b>",
        xaxis = list(title = "Média da Idade por Amostra"),
        yaxis = list(title = "Densidade"),
        legend = list(orientation = 'h', xanchor = "center", x = 0.5, y = -0.2)
      )