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...

Boxplot
Modificando as cores
Formatando a figura
Boxplot
Para os exemplos a seguir, clicar aqui para abrir o banco de dados em formato Excel.
# Função boxplot()
boxplot(Peso)
# A função boxplot.stats mostra os detalhes
boxplot.stats(Peso)
$stats
[1] 47.9 62.9 69.3 74.1 89.6
$n
[1] 40
$conf
[1] 66.50202 72.09798
$out
[1] 111.5
# Valores atípicos (outliers)
boxplot.stats(Peso)$out
[1] 111.5
# Comparando dois grupos
fet <- ifelse(Idade<=50,0,1)
boxplot(Peso~fet)
boxplot(Peso~fet,names=c("Até 50 anos","Mais de 50 anos"),main="(a) Peso", ylab="Peso (kg)",xlab="Idade")

# Cor do boxplot
boxplot(Peso~fet,names=c("Até 50 anos","Mais de 50 anos"),main="(a) Peso",ylab="Peso (kg)",xlab="Idade",
las=1,col="bisque",pch=19)
# staplewex e whisklty
boxplot(Peso~fet,names=c("Até 50 anos","Mais de 50 anos"),main="(a) Peso",ylab="Peso (kg)",xlab="Idade",
las=1, col="bisque",pch=19,staplewex=0.3,whisklty=1)

# Retirando o "contorno" do gráfico
boxplot(Peso~fet, main="(a) Peso",ylab="Peso (kg)",xlab="Idade", las=1,col="bisque",pch=19,
staplewex=0.3, whisklty=1,boxwex=0.6,axes=F)
axis(2,las=1)
axis(1,at=1:2,label=c("Até 50 anos","Mais de 50 anos"))

# Boxplot com as estatísticas descritivas
boxplot(Peso~fet, main="(a) Peso",ylab="Peso (kg)", las=1,col="bisque",pch=19,staplewex=0.3,
whisklty=1,boxwex=0.6,axes=F)
axis(2,las=1)
axis(1,at=1:2,label=c("Até 50 anos","Mais de 50 anos"))
mtext(side=1,at=0.3,line=1,"Idade")
mtext(side=1,at=0.3,line=2,"Média")
mtext(side=1,at=0.3,line=3,"DP")
mtext(side=1,at=1,line=2,round(mean(Peso[fet==0]),2))
mtext(side=1,at=2,line=2,round(mean(Peso[fet==1]),2))
mtext(side=1,at=1,line=3,round(sd(Peso[fet==0]),2))
mtext(side=1,at=2,line=3,round(sd(Peso[fet==1]),2))
