Bioestatística
Prof. Dr. Edson Zangiacomi Martinez
Faculdade de Medicina de Ribeirão Preto
Universidade de São Paulo (USP)
Esta página está em construção!
Todo seu conteúdo não é definitivo...
Histogramas
Tabela de frequências
Modificando cores
Histograma
Histogramas são gerados com a função hist()
Para os exemplos a seguir, clicar aqui para abrir o banco de dados em formato Excel.
hist(dados$Altura)
# Argumentos
hist(dados$Altura,
main="Distribuição das alturas", # Modifica o título do gráfico
xlab="Alturas (m)", # Label do eixo x
ylab="Frequências", # Label do eixo y
border="blue", # Cor das bordas dos retângulos
col="cyan", # Cor das bordas dos retângulos
ylim=c(0,20), # Limites do eixo y
las=1) # Rotaciona os valores nos eixos
# Atribuindo o resultado da função hist() a um objeto:
figura <- hist(dados$Altura, main="Distribuição das alturas", xlab="Alturas (m)",
ylab="Frequências", border="blue", col="cyan", ylim=c(0,20), las=1)
figura
$breaks
[1] 1.45 1.50 1.55 1.60 1.65 1.70 1.75
$counts
[1] 3 12 16 6 2 1
$density
[1] 1.5 6.0 8.0 3.0 1.0 0.5
$mids
[1] 1.475 1.525 1.575 1.625 1.675 1.725
$xname
[1] "dados$Altura"
$equidist
[1] TRUE
attr(,"class")
[1] "histogram"
# Exibindo a tabela de frequências:
k <- length(figura$counts)
intclass <- paste(figura$breaks[1:(k-1)]," - ",figura$breaks[2:(k+1)])
tabela <- data.frame(intclass,figura$counts,figura$counts/length(dados$Altura),figura$mids)
names(tabela) <- c("Intervalo","n","freq rel","Ponto médio")
tabela
Intervalo n freq rel Ponto médio
1 1.45 - 1.5 3 0.075 1.475
2 1.5 - 1.55 12 0.300 1.525
3 1.55 - 1.6 16 0.400 1.575
4 1.6 - 1.65 6 0.150 1.625
5 1.65 - 1.7 2 0.050 1.675
6 1.45 - 1.75 1 0.025 1.725
# Argumento prob
hist(dados$Altura,
main="Distribuição das alturas", # Modifica o título do gráfico
xlab="Alturas (m)", # Label do eixo x
ylab="Densidade", # Label do eixo y
border="blue", # Cor das bordas dos retângulos
col="cyan", # Cor das bordas dos retângulos
ylim=c(0,10), # Limites do eixo y
las=1, # Rotaciona os valores nos eixos
prob=T) # Densidades de probabilidade, e não frequências
# Mostrando as contagens em cada intervalo de classe:
hist(dados$Altura, main="Distribuição das alturas", xlab="Alturas (m)",
ylab="Frequências", border="brown", col="bisque", ylim=c(0,20), las=1, labels=T)
# Argumento breaks
hist(dados$Altura, main="Distribuição das alturas", xlab="Alturas (m)",
ylab="Frequências", border="brown", col="bisque", ylim=c(0,8), las=1,
xlim=c(1.45,1.75), breaks=seq(1.45,1.75,0.015), labels=T)
hist(dados$Altura, main="Distribuição das alturas", xlab="Alturas (m)",
ylab="Frequências", border="brown", col="bisque", ylim=c(0,8), las=1,
xlim=c(1.45,1.75), breaks=seq(1.45,1.75,0.015),labels=T, axes=F)
axis(2)
axis(1,at=seq(1.45,1.75,0.015))