Eu recomendo usar os pacotes ggplot2
e Rmisc
para fazer este gráfico. O primeiro pacote faz o gráfico em si, enquanto o segundo prepara os dados para a análise. Abaixo vou explicar o passo a passo de como construí o gráfico desejado.
Em primeiro lugar, utilizo a função summarySE
para obter as médias e erros padrão do conjunto de dados desejado. Veja que bastou informar o nome da data frame, a variável resposta e a variável de agrupamento para obter o que desejamos.
library(ggplot2)
library(Rmisc)
malaria.plot <- summarySE(malaria, measurevar="peso", groupvars="grupo", na.rm=TRUE)
Em seguida, para o gráfico ficar com os labels do eixo x com os nomes corretos, sem o uso de números, converti a coluna grupo
em fator. Abaixo mostro o resultado final desta preparação dos dados.
malaria.plot$grupo <- factor(c("Controle", "Vivax", "Falciparum", "Mista"),
levels=c("Controle", "Vivax", "Falciparum", "Mista"))
malaria.plot
grupo N peso sd se ci
1 Controle 206 3225.830 510.6585 35.57927 70.14821
2 Vivax 173 3134.098 508.4377 38.65580 76.30084
3 Falciparum 100 3122.550 512.6327 51.26327 101.71744
4 Mista 56 3144.696 489.7987 65.45211 131.16896
Com os dados prontos, basta fazer o gráfico. Utilizei a função ggplot
juntamente com geom_errorbar
, considerando os cálculos presentes em malaria.plot
. Perceba que não estou plotando o conjunto de dados original, mas sim a transformação que fiz nos dados. Utilizo também geom_line
e geom_point
para fazer os pontos e as linhas unindo-os. Por fim, labs
acrescenta os nomes aos eixos.
ggplot(malaria.plot, aes(x=grupo, y=peso, group=1)) +
geom_errorbar(aes(ymin=peso-se, ymax=peso+se), width=.1) +
geom_line() +
geom_point() +
labs(x="Grupo", y="Peso (kg)")

Caso não queira o fundo cinza na imagem, adicione + theme_bw()
ao comando acima. Outros detalhes do gráfico podem ser ajustados procurando o help do ggplot2
.