require(gtools)
require(tidyverse)
require(ggforce)  
options(kableExtra.latex.load_packages = TRUE)
require(kableExtra)
require(mosaicCalc)
require(gridExtra)
require(numDeriv)
require(DT)
#Para numerar tabelas e figuras:
require(captioner)
fig_nums <- captioner(prefix = "Figura")
table_nums <- captioner(prefix = "Tabela")
quadro_nums<- captioner(prefix="Quadro")

Probabilidade: Axiomas e Propriedades

Definição Frequentista

Definição: Dado um experimento aleatório, sendo \(\Omega\) o seu espaço amostral e \(P(A)\): probabilidade de um evento A, (\(A \subset \Omega\)) é dada por: \[ P(A) = \frac{n(A)}{n(\Omega)} = \frac{\mbox{Número de casos favoráveis (A)}}{\mbox{Número total de casos}}. \]

Limitações: \(\Omega\) tem que ser finito e equiprovável (todos os eventos com a mesma probabilidade de ocorrência).

Exemplo 1: Considerando o Lançamento de uma moeda honesta e o evento A: “obter cara”, temos: \[ \small \Omega = \left\{\mbox{cara},\mbox{coroa}\right\} \Rightarrow n_\Omega = 2\\ A = \left\{\mbox{cara}\right\} \Rightarrow n_A = 1 \] Logo: \(P(A) = \frac{1}{2}\). Esse resultado mostra que se lançarmos uma moeda equilibrada, temos a probabilidade de 50% que apareça cara na face superior.

Definição Axiomática de Probabilidade

Definição de probabilidades através de axiomas e propriedades.

Leis de Morgan:

  • \(A^c \cup B^c = (A\cap B)^c\): a união dos complementares é o complementar da intersecção;
  • \(A^c \cap B^c = (A\cup B)^c\): a interseccão dos complementares é o complementar da união.
  • As leis são extendidas para três ou mais eventos.

Definições: A probabilidade é definida em uma classe de subconjuntos de \(\Omega\) com

    1. \(P(\Omega)=1\): a probabilidade do espaço amostral é igual a 1.
    1. Para \(A, B \subseteq \Omega\) disjuntos: \(P(A \cup B) = P(A)+P(B)\) (regra da aditividade).

Eventos disjuntos são quando \(A \cap B = \emptyset\). A regra da aditividade serve para três ou mais eventos disjuntos dois a dois.

Propriedades da probabilidade:

    1. \(0 \leq P(A) \leq 1, \forall A \subseteq \Omega\): a probabilidade é um número entre zero e um;
    1. \(P(\emptyset) = 0\); a probabilidade do evento vazio é igual a zero;
    1. \(P(A^c) = 1 - P(A)\); a probabilidade do evento complementar é igual a 1 menos a probabilidade deste evento;
    1. Se \(A \subseteq B,\) com \(A, B \subseteq \Omega\), então \(P(A) \leq P(B);\) se A está contido em B, então a probabilidade de A é menor ou igual à probabilidade de B.
    1. Regra da Adição: \[\small P(A \cup B) = P(A) + P(B) - P(A \cap B);\]

Exemplo 2: Um baralho de 52 cartas é subdividido em 4 naipes: copas, espadas, ouros e paus.

    1. Retirando-se uma carta ao acaso, qual é a probabilidade de que ela seja de ouros ou de copas?
    1. Retirando-se duas cartas ao acaso com reposição da 1ª carta, qual é a probabilidade de ser a 1ª de ouros e a 2ª de copas?
    1. Recalcular a probabilidade anterior se não houver reposição da 1ª carta.
    1. Havendo reposição, qual é a probabilidade de sair a 1ª carta de ouros ou então a 2ª de copas?

Resolução do exemplo 2: \[ \small \begin{array}{|llll|} \hline && a) \mbox{Como os eventos são disjuntos, basta somar as probabilidades:}\\ && \mbox{ Sejam: } \left\{ \begin{array}{ll} \mbox{O: A carta é de ouros}\\ \mbox{C: A carta é de copas} \end{array} \right. \Rightarrow \begin{array}{lll} P(\mbox{O} \cup \mbox{C}) &= P(\mbox{O}) + P(\mbox{C}) \\ &= \frac{13}{52}\times 2 = 50\% \end{array}\\ \hline && b) \mbox{Neste caso a ordem é importante, então resolve-se pelo princípio} \\ &&\mbox{fundamental da contagem:}\\ && \mbox{P(1ª carta = O, 2ª carta = C)}=\frac{13}{52}.\frac{13}{52}=\frac{169}{2704}\approx 0.0625\\ \hline && c) \mbox{Também resolve-se pelo princípio fundamental da contagem:}\\ && \mbox{P(1ª carta = O, 2ª carta = C)}=\frac{13}{52}.\frac{13}{51}=\frac{169}{2652}\approx 0.0637\\ \hline && d) \mbox{Pela regra da adição, sendo as retiradas independentes:}\\ && \mbox{P(1ª carta = O} \cup \mbox{2ª carta = C)}=\\ && \mbox{P(1ª carta = O) + P(2ª carta = C)-P(1ª carta = O, 2ª carta = C)}=\\ && \frac{13}{52}+\frac{13}{52}-\frac{169}{2704}=\frac{676.2-169}{2704}=0.4363\\ \hline \end{array} \]

Probabilidade Condicional e Independência de Eventos

A probabilidade do evento A ocorrer, dado que B ocorreu, é de: \[\small P(A|B) = \frac{P(A \cap B)}{P(B)}, \quad P(B) \neq 0. \] A probabilidade de A ocorrer, dado que B ocorreu, é igual a probabilidade de ocorrência simultânea de A e B dividida pela probabilidade de ocorrência de B.

  • Sempre que calculamos \(P(A|B)\), estamos essencialmente calculando \(P(A)\) em relação ao espaço amostral reduzido devido a \(B\) ter ocorrido, em lugar de faze-lo em relação ao espaço amostral original \(\Omega\).

Regra do Produto

\[ \begin{array}{l} P(A \cap B) = P(B) \times P(A|B), \mbox{ou equivalentemente,}\\ P(A \cap B) = P(A) \times P(B|A). \end{array} \] Caso especial: Para A e B independentes, \[ P(A \cap B) = P(A) \times \underbrace{P(B)}_{P(B|A)} \]

Regra da probabilidade total

A probabilidade de um evento B pode ser escrita em partes:

\[\small P(B) = \underbrace{P(A).P(B|A)}_{P(A\cap B)}+\underbrace{P(A^c).P(B|A^c)}_{P(A^c \cap B)} \]

  • Esta regra pode ser extendida para três eventos A, B, C, desde que:
    • \(A \cup B \cup C = \Omega\), e
    • os eventos A, B e C são disjuntos dois a dois (a intersecção é vazia).

Exemplo 3: (cálculo de probabilidades em tabela de dupla entrada) Um grupo de 60 pessoas apresenta a seguinte composição:

require(kableExtra)
tabela = data.frame(
                Condicao = c("Menores","Adultos","Total"),
                Homens = c(15,18,33),
                Mulheres = c(17,10,27),
                Total = c(32,28,60))
          
  colnames(tabela) <- c("", "Homens","Mulheres","Total")
  cap = table_nums("tab2","Resultados da pesquisa")
  tabela %>%
  kbl(caption = cap) %>%
  kable_classic(full_width = F, html_font = "Cambria") %>%
     kable_styling(bootstrap_options = c("striped","hold_position"),
                full_width = F) %>%
  add_header_above(c("Condição", "","Número de pessoas",""))
Tabela 1: Resultados da pesquisa
Condição
Número de pessoas
Homens Mulheres Total
Menores 15 17 32
Adultos 18 10 28
Total 33 27 60

Uma pessoa é escolhida ao acaso. Pergunta-se:

    1. qual é a probabilidade de ser homem?
    1. qual é a probabilidade de ser adulto?
    1. qual é a probabilidade de ser menor e ser mulher?
    1. sabendo-se que a pessoa escolhida é adulto, qual é a probabilidade de ser homem?
    1. dado que a escolhida é mulher, qual é a probabilidade de ser menor?

Teorema de Bayes

Vem direto da da probabilidade condicional e regra do produto: \[ P(A|B) = \frac{P(A\cap B)}{P(B)}=\frac{P(A).P(B|A)}{P(A)P(B|A)+P(A^c).P(B|A^c)} \]

  • Pode ser extendido para três ou mais eventos.

Exemplo 4: Um restaurante popular apresenta apenas dois tipos de refeições: salada completa ou um prato à base de carne. Considere que 20% dos fregueses adultos preferem a salada, 30% das crianças escolhem carne, 75% dos fregueses são adultos e os seguintes eventos:

  • A: freguês é adulto;
  • \(A^c\): freguês é criança;
  • B: freguês prefere salada;
  • \(B^c\): freguês prefere carne;

Calcule:

    1. \(P(A)\), \(P(B|A)\), \(P(B^c|A^c)\);
    1. \(P(A\cap B)\), \(P(A\cup B)\);
    1. \(P(A^c|B)\);

Exemplo 5: Em certo colégio, 5% dos homens e 2% das mulheres têm mais de 1,80m de altura. Por outro lado, 60% dos estudantes são homens. Sorteando-se um estudante aleatoriamente, qual é a probabilidade de:

    1. Ser mulher e ter mais de 1,80m?
    1. Ter mais de 1,80m?
    1. Um estudante é escolhido ao acaso e tem mais de 1,80m. Qual é a probabilidade de que o estudante seja mulher?

Resolução do exemplo 5: Considere os seguintes eventos:

  • \(A\): Estudante do sexo masculino;

  • \(A^c\): Estudante do sexo feminino;

  • \(B\): Estudante com mais de 1,80m de altura;

  • \(B^{c}\): Estudante com menos de 1,80m de altura;

  • Ser mulher \((A^c)\) e ter mais de 1,80m (\(B\))? \[\small P(A^c \cap B) = P(A^c) \times P(B|A^c)= 0,40 \times 0,02 = 0,008 \]

  • Ter mais de 1,80m? \[ P(B) = \underbrace{P(A).P(B|A)}_{P(A\cap B)}+\underbrace{P(A^c).P(B|A^c)}_{P(A^c \cap B)}\\ = 0,6.0,05+0,4.0,02=0,038 \]

  • Um estudante é escolhido ao acaso e tem mais de 1,80m. Qual é a probabilidade de que o estudante seja mulher? Fórmula da probabilidade condicional utilize as letras anteriores para resolver! \[ P(A^c|B) =\frac{P(A^c \cap B)}{P(B)}=\frac{0,008}{0,038}=0,2105 \]