ISCTE - IUL - Probabilidades - e - Estatística - (Cap1)
ISCTE - IUL - Probabilidades - e - Estatística - (Cap1)
ISCTE - IUL - Probabilidades - e - Estatística - (Cap1)
(2024)
Probabilidades e Estatı́stica
Probabilidades e Estatı́stica ISCTE - IUL
Conteúdo
1 Introdução à Estatı́stica Descritiva 4
1.1 Tipos de variáveis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Tabela de frequências e representações gráficas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Medidas de localização . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3.1 Tendência central . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3.2 Tendência não central . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3.3 Relação entre quartis, decis e percentis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4 Medidas de dispersão, assimetria e curtose . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.1 Medidas de dispersão . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.2 Medidas de assimetria e curtose . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4.3 Box Plot e outliers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5 Medidas de associação entre variáveis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.5.1 Correlação de Pearson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.5.2 Coeficiente V de Cramer e coeficiente Phi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.6 Estatı́stica Descritiva: Exemplos de outputs com o R/RSTUDIO . . . . . . . . . . . . . . . . . . . . . . . . 12
1.7 Estatı́stica Descritiva: Exercı́cios propostos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Lista de Figuras
4 Distribuições bimodais e unimodal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
5 Distribuições (as)simétricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
6 Skewness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
7 Kurtosis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
8 Boxplots . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
9 Histogramas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
10 Outliers moderados e severos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
11 Representação gráfico dos dados: Preferência de férias. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
12 Representação gráfica dos dados: Notas (amostra2). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
13 Boxplot dos dados: Notas (amostra2). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
14 Layout do Excel: Agregado familiar. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
15 Representação gráfica das idades dos habitantes de uma aldeia. . . . . . . . . . . . . . . . . . . . . . . . . . 16
16 Layout do Excel: Ogiva. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
17 Layout do Excel: Boxplot. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
18 Layout do Excel: Colunas empilhadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
19 Layout do Excel: Gráfico de dispersão com reta de regressão . . . . . . . . . . . . . . . . . . . . . . . . . . 17
20 Layout do Excel: Gráfico de dispersão. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
21 Diagrama de Dispersão TV/GPA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
22 Boxplot e Histograma: Dados (amostra3). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Lista de Tabelas
1 Tabela de frequências . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2 Nı́vel de associação (proposto) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3 Tabela de contingência . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4 Nı́vel de associação (Cohen, 2013) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
5 Agregado familiar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
6 Idades dos habitantes da aldeia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
7 Index - HEI. Fonte healthy-eating-index-hei. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
8 TV hours week - Grade point average . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
9 Número de sinistros. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
10 Pesos dos estudantes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
11 Classificação. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
12 Dados para o coeficiente de correlação. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
13 Disposição dos alunos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
14 O tabagismo e a presença de doença pulmonar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
15 Tipo de sangue - Sintomas (Esp.) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
Livros recomendados
https://silabo.pt/categoria-produto/ciencias-exatas-e-naturais/estatistica/
https://diascurto.wixsite.com/sitedc/livros
A Estatı́stica Descritiva consite na recolha, apresentação, análise e interpretação de dados numéricos através da criação de
instrumentos adequados: quadros, gráficos e indicadores numéricos (Reis, 2008).
Na Inferência Estatı́stica é (primeiro) necessário introduzir os conceitos de população e amostra.
Uma amostra é um subconjunto retirado de uma população. A dimensão/tamanho da amostra é representado por n.
Uma amostra é um subconjunto (representativo) da população . É selecionada de forma aleatória ou não aleatória, com o
objetivo de extrair informações e fazer inferências sobre a população. O tamanho da amostra pode variar e é geralmente muito
menor do que o tamanho da população, mas deve ser grande o suficiente para ser estatisticamente significativa e representativa.
A Inferência Estatı́stica é uma metodologia que permite fazer estimativas e tirar conclusões sobre uma população a partir da
informação contida numa amostra (representativa da população de onde foi retirada).
Definição 3. Parâmetro e Estatı́stica
Um parâmetro é uma medida numérica que descreve uma caracterı́stica especı́fica de uma população.
Uma estatı́stica é uma medida numérica que descreve uma caracterı́stica especı́fica de uma amostra.
k Fi k
• fi =
• ∑ Fi = n n • ∑ fi = 1
i=1 i=1
Xi Fi fi CumFi Cum fi
x1 F1 f1 CumF1 = F1 Cum f1
x2 F2 f2 CumF2 = F1 + F2 Cum f2 = f1 + f2
x3 F3 f3 CumF3 = F1 + F2 + F3 Cum f3 = f1 + f2 + f3
··· ··· ··· ··· ···
xk Fk fk n 1
Total n 1
Dois do tipos de gráficos mais usados são os gráficos de barras e os gráficos circulares. Nos dados qualitativos e quantita-
tivos discretos a largura das barras (deve deixar-se um espaço entre as barras adjacentes) é igual para todas as categorias e a
altura é igual à frequência. As categorias são representadas no eixo horizontal e as frequências no eixo vertical.
Os gráficos circulares (normalmente usados nos dados qualitativos) estão divididos em setores circulares cuja área (e ângulo
ao centro correspondente) é proporcional à frequência da categoria que representam.
Quando a variável é qualitativa, numa escala ordinal, ou quantitativa, faz sentido falar em frequências acumuladas. O CumFi
dá-nos o número que observações com valor inferior ou igual à caracterı́stica xi [F(xi ) = CumFi ].
O gráfico da distribuição de frequências acumuladas de dados quantitativos discretos é uma ‘escada de graus’, onde a altura
de cada grau é a frequência acumulada e o ‘salto’ entre classes é a frequência (simples). Os pontos de descontinuidade existem
nos ‘saltos’, pelo que a função cumulativa é contı́nua à direita desses pontos.
Quanto trabalhamos com dados quantitativos contı́nuos é, quase sempre, necessário agrupá-los em classes para os represen-
tar graficamente (os métodos para agrupar os dados em classes não será objeto de estudo).
Definição 4. Histograma
É uma representação gráfica de dados quantitativos onde as classes são representadas no eixo horizontal e as frequências
no eixo vertical. Se as classes tiverem todas a mesma amplitude, a altura representa a frequência.
No caso das classes não terem todas a mesma amplitude, a área das barras é proporcional à sua frequência.
A média aritmética (ou só média) é a soma de todos os valores numéricos observados a dividir pelo número de observações
(se os valores forem todos iguais a uma constante a, isto é, se não existir dispersão, a média é a):
∑ni=1 xi
x̄ = .
n
Quando os dados estão organizados numa tabela de frequências para determinarmos a média usamos:
∑ki=1 Fi xi k
x̄ = = ∑ fi xi .
n i=1
A mediana divide ao meio o conjunto de valores observados (ordenados por ordem crescente ou decrescente). É menor valor
da variável que contém pelo menos 50% das observações.
O conjunto de dados é dividido de forma a que o número de observações inferiores ou iguais à mediana seja igual ao número
de observações superiores ou igual à mediana (o cálculo da mediana pressupõe que os dados já foram previamente ordenados).
Consideremos x1 , x2 , ..., xn um conjunto de observações. Usando a notação x(i) para representar o i-ésimo valor observado
mais pequeno, passamos a ter a nossa amostra ordenada, pelo que se verifica que x(1) ≤ x(2) ≤ · · · ≤ x(n) .
Em particular x(1) e x(n) são respetivamente o mı́nimo e o máximo da amostra.
x( ) + x( 2n +1)
n
2 se n par,
2
Me = x̃ =
x( n+1 ) se n ı́mpar.
2
A moda é o valor que ocorre com maior frequência num conjunto de observações. A moda pode não existir (conjunto amodal)
e se existir pode não ser única. Se for única a distribuição diz-se unimodal, se tiver duas modas diz-se bimodal, se tiver mais
de duas diz-se multimodal. Se os dados estiverem agrupados em classes numa tabela de frequências o cálculo da moda é mais
complexo.
(a) Distribuição simétrica (b) Distribuição assimétrica positiva (c) Distribuição assimétrica negativa
(np) + x(np+1)
x
se np inteiro,
Qp = 2
x[np+1] se np não inteiro.
As medidas de assimetria indicam se uma distribuição é assimétrica ou não. O coeficiente de assimetria é dado por uma
medida normalizada (para poderem serem feitas comparações) que se obtém dividindo m3 (momento centrado na média de
ordem três) pelo cubo do desvio padrão.
Existem dois coeficientes. Os dois coeficientes tomam valores muito próximos quando o n é grande (n → +∞).
Figura 6: Skewness
As medidas de curtose indicam sobre o ‘peso’ das caudas da distribuição. coeficiente de curtose é dado por uma medida
normalizada (para comparações) que se obtém dividindo m4 (momento centrado na média de ordem quatro) pelo desvio
padrão elevado a 4.
O R/RSTUDIO também dá o excesso de curtose (Dias, 2021). Os valores dos dois coeficientes diferem aproximadamente 3
unidades com n grande (n → +∞), k2 = k1 − 3 (excesso de curtose):
m4 n2 (n + 1) 3(n − 1)2
k1 = 4
−→ k2 = kurt = × k1 − .
s (n − 1)(n − 2)(n − 3) (n − 2)(n − 3)
(a) Mesocúrtica: k1 ≈ 3 (k2 ≈ 0) (b) Leptocúrtica: k1 > 3 (k2 > 0) (c) Platicúrtica: k1 < 3 (k2 < 0)
Figura 7: Kurtosis
A interpretação não é fácil por isso costuma-se comparar o seu valor com a da curva Normal que, com k1 = 3 (k2 = 0), é
designada de Mesocúrtica (a). Assim se k1 > 3 (k2 > 0) a curva da distribuição é mais esguia, designada Leptocúrtica (b),
com caudas mais ‘pesadas’ do que a distribuição Normal (‘caudas pesadas e um pico alto’). Se k1 < 3 (k2 < 0) a curva
da distribuição é mais achatada, designada por Platicúrtica (c), com caudas menos ‘pesadas’ do que a distribuição Normal
(‘caudas leves e um pico baixo’).
Figura 8: Boxplots
Figura 9: Histogramas
A covariância é uma medida estatı́stica que reflete o grau de variação conjunta de duas variáveis aleatórias, indicando se
tendem a variar na mesma direção (covariância positiva), ou seja se uma aumenta a outra também aumenta, ou em direções
opostas (covariância negativa), se uma aumenta a outra diminui. O valor da covariância pode assumir qualquer valor real
[−∞ < Cov(X,Y ) < +∞], por isso apenas podemos verificar a direção da relação mas não a intensidade da relação (depende
das unidades de medida das variáveis). A covariância fornece uma indicação inicial da relação entre variáveis, diferenciando-
se da correlação, que normaliza essa relação, permitindo a interpretação por que o valor está limitado ente -1 e 1.
A existência de correlação não implica necessariamente uma relação de causa e efeito entre as variáveis.
Quando uma das variáveis é qualitativa ordinal e a outra quantitativa, ou as duas variáveis são qualitativas ordinais, o coefici-
ente de correlação de Spearman é o mais usado (no entanto não o iremos calcular este ano).
Y Totais |a × d − b × c|
0≤φ = p ≤ 1(perfeita associação)
a b a+b (a + b)(a + c)(b + d)(c + d)
X
c d c+d
Totais a+c b+d a+b+c+d
Preferência de férias
Xi Fi fi
Casa 8 0.08
Praia 46 0.46
Campo 32 0.32
Outra 14 0.14
Total 100 1
Xi Fi fi
1 8 0.08
2 20 0.20
3 30 0.30
4 26 0.26
5 12 0.12
6 4 0.04
[li , Li [ Fi fi
[15, 25[ 1 0.02
[25, 35[ 4 0.08
[35, 45[ 9 0.18
[45, 55[ 14 0.28
[55, 65[ 12 0.24
[65, 75[ 9 0.18
[75, 85[ 1 0.06
Figura 15: Representação gráfica das idades dos habitantes de uma aldeia.
8. Considere a seguinte tabela, de dupla entrada, relativa ao nı́vel de atividade de 4460 participantes. Represente os dados
fazendo Inserir Gráfico-Colunas 2D-Colunas Empilhadas. Deve formatar o gráfico: Trocar linhas/colunas.
Homens Mulheres
Sedentário 957 1 226
Ativo 340 417
Muito ativo 842 678
Total 2 139 2 321
9. Comece por criar aleatoriamente duas colunas com 30 observações cada. Calcule o coeficiente de regressão de Pearson
e interprete o resultado. Represente graficamente os dados fazendo Inserir Gráfico-Dispersão. Selecionar o gráfico,
clicar no botão direito e Adicionar linha de tendência (mostrar equação no gráfico).
10. Represente o gráfico de dispersão dos dados, calcule o coeficiente de correlação e interprete o resultado:
(a) Excel (Ferramentas-Suplementos-Análise de Dados).
(b) R/RSTUDIO (Recorrendo aos dados introduzidos no Excel).
TV GPA
20 2.35
5 3.8
8 3.5
10 2.75
13 3.25
7 3.4
13 2.9
5 3.5
25 2.25
14 2.75
Figura 20: Layout do Excel: Gráfico de dispersão.
1 rm(list=ls())
2 setwd("/Users/user/Desktop")
3 TV <- c(20, 5, 8, 10, 13, 7, 13, 5, 25, 14)
4 GPA <- c(2.35, 3.8, 3.5, 2.75, 3.25, 3.4,
5 2.9, 3.5, 2.25, 2.75)
6 coeficiente_correlacao <- cor(TV, GPA)
7 plot(TV, GPA, main="Diagrama de Dispersão TV/GPA",
8 xlab="Horas de TV", ylab="GPA",
9 pch=19, col="red", cex=0.5)
Estatı́sticas descritivas
Dados em bruto
1 rm(list=ls())
2 setwd("/Users/user/Desktop")
3 getwd()
4 dir()
5 install.packages("readxl")
6 library(readxl) #Para ler ficheiros de Excel
7 amostra3 <- read_excel("PE1.xlsx", sheet="Exe11")
8 attach(amostra3) #"diz" ao R que a primeira linha tem o nome das variáveis
9 head(amostra3, 6)
10 # Application: Boxplot and Histogram
11 par(mfrow=c(1,2))
12 boxplot(amostra3$Dados, horizontal = TRUE, col = "blue", border ="green", main = "Dados (amostra3)")
13 hist(amostra3$Dados, freq = TRUE, col = "red", border ="yellow", breaks = "Sturges",
14 main = "Dados (amostra3)", xlab = "Classes")
{69, 85, 75, 89, 73, 61, 62, 75, 98, 63}.
Determine e interprete:
(a) A média, a mediana, o 1.º quartil e o 3.º quartil.
(b) A variância e o desvio padrão.
(c) O coeficiente de assimetria e o coeficiente (excesso) de curtose.
{7.3, 12.7, 7.4, 8.5, 9.0, 9.5, 9.2, 12.4, 10.3, 12.3, 11.3, 10.4, 7.3, 9.4, 8.2, 12.9, 8.9, 8.9}.
Determine e interprete:
(a) A média, a mediana, a moda, o 1.º quartil e o 3.º quartil.
(b) A variância e o desvio padrão.
(c) O coeficiente de assimetria e o coeficiente (excesso) de curtose.
1 install.packages("DescTools")
2 library(DescTools)
3 install.packages("moments")
4 library(moments)
5 dados <- c(7.3, 12.7, 7.4, 8.5, 9.0, 9.5, 9.2, 12.4, 10.3,
6 12.3, 11.3, 10.4, 7.3, 9.4, 8.2, 12.9, 8.9, 8.9)
7 Kurt(dados, method = 1)#k2
8 Skew(dados, method = 2)#G1
9. O número de sinistros participados por cada apólice do ramo automóvel é uma informação de grande importância para
que as seguradoras possam determinar o valor do prémio do seguro. A tabela seguinte apresenta a informação relativa
à sinistralidade em Responsabilidade Civil de uma seguradora francesa.
(a) Complete a tabela de frequências, calculando as frequências absolutas e relativas simples e acumuladas.
(b) Represente a informação por meio de um gráfico de barras.
(c) Calcule a média, a mediana, a moda, a variância e o desvio padrão do número de sinistros por apólice.
(d) Determine e interprete o coeficiente de variação.
10. Considere a tabela seguinte, que mostra a distribuição dos pesos de 40 estudantes do ISCTE.
11. Considere a tabela seguinte, que mostra a classificação de 100 estudantes do ISCTE no exame de Estatı́stica.
Classificação Frequência
[2, 4[ 4
[4, 6[ 12
[6, 8[ 16
[8, 10[ 23
[10, 12[ 19
[12, 14[ 15
[14, 16[ 6
[16, 18[ 4
[18, 20[ 1
• Turma A: {0,0,0,0,0,20,20,20,20,20};
• Turma B: {0,1,2,3,8,9,11,12,17,18,19,20};
• Turma C: {10,10,10,10,10,10,10,10,10}.
A média de cada uma das turmas é igual a 10. Isso quer dizer que as turmas são semelhantes?
15. Os ordenados (em euros) da Empresa X são: {500, 500, 600, 1300, 1400, 2200, 60000}. Comente:
16. Considere os seguintes dados (pesos em kg) e comente a afirmação (verdadeira ou falsa):
• Homem: Média é de 80 kg e o desvio padrão de 11 kg;
• Elefante: Média é de 4 200 kg e o desvio padrão de 294 kg.
”A dispersão nos pesos é mais elevada nos Elefantes do que nos Homens, porque o desvio padrão é maior.”
17. Uma companhia de transportes de carga pretende comparar a variabilidade do peso com a variabilidade do volume de
caixas que transporta. Com os dados da tabela sugira um método para efetuar essa comparação.
18. Observe os histogramas apresentados em baixo (com base em três amostras aleatórias de 1000 observações cada), diga
qual das distribuições tem maior (e menor) coeficiente, de excesso, de curtose (kurtosis):
20. Represente o gráfico de dispersão e calcule o coeficiente de correlação de Pearson dos dados na Tabela 12 :
(a) Excel (Ferramentas-Suplementos-Análise de Dados).
(b) R/RSTUDIO (Recorrendo aos dados introduzidos no Excel).
1 library(readxl)#Carrega o package necessário
2 # Ler o ficheiro do Excel como um dataframe
3 dados <- read_excel("caminho/.../dados.xlsx")
4 #Calcular a correlação entre duas colunas do dataframe
5 correlacao <- cor(dados$X, dados$Y)
6 print(correlacao)#Mostrar o resultado
X 1 3 4 6 8 9 11 14
Y 1 2 4 4 5 7 8 9
22. Um investigador efetuou uma experiência para determinar se o género do professor, influência a disposição para os
alunos (homens) assistirem à aula. Para isso foi recolhida uma amostra, e os resultados encontram-se na tabela seguinte.
(a) Calcule o coeficiente Phi para avaliar a força da associação entre as variáveis em estudo.
(b) Interprete o resultado obtido.
23. Foi realizado um estudo para investigar a relação entre o hábito de fumar e a ocorrência de doenças pulmonares. Os
dados foram recolhidos numa amostra de indivı́duos e categorizados conforme apresentado na Tabela 14:
(a) Calcule o coeficiente Phi para avaliar a força da associação entre as variáveis em estudo.
(b) Que recomendações poderiam ser feitas com base na sua análise?
24. A Covid-19 é uma doença infecciosa causada pelo vı́rus SARS-CoV-2, declarada pela OMS como uma pandemia.
Recentes estudos sugeriram um efeito protetor de anticorpos anti-A, associando a presença destes no soro (grupos O e
B) às apresentações brandas da doença, enquanto a ausência de circulação de anticorpos anti-A (grupos A e AB) esteve
relacionada às formas graves. O dados encontram-se registados na Tabela 15:
(a) Calcule o coeficiente Phi para avaliar a força da associação entre as variáveis em estudo.
(b) Que análise faz do valor obtido?
Referências
Cohen, J. (2013). Statistical power analysis for the behavioral sciences. Academic press.
Dias, C. (2021). Estatı́stica com R: Aprenda Fazendo. 1.ª Edição.
Reis, E. (2008). Estatı́stica descritiva. Lisboa: Edições Sı́labo.
Reis, E., Melo, P., Andrade, R., and Calapez, T. (2015). Estatı́stica aplicada. Lisboa: Edições Sı́labo.