Boxplots

Lectura de datos

tabla2.7<-read.table("LEPPIK.DAT")
tabla2.7$V1<-factor(tabla2.7$V1, labels=c("Placebo","Progabide"))
tabla2.7$V2<-factor(tabla2.7$V2)
dimnames(tabla2.7)[[2]]<-c("Grupo","Indiv",paste("Semana",c(2,4,6,8) ,sep="_") )
semanas<-c(2,4,6,8)
datos<-melt(tabla2.7,id=c("Grupo","Indiv"))
head(datos)

##       Grupo Indiv variable value
## 1 Progabide   101 Semana_2    11
## 2 Progabide   102 Semana_2     8
## 3 Progabide   103 Semana_2     0
## 4 Progabide   108 Semana_2     3
## 5 Progabide   110 Semana_2     2
## 6 Progabide   111 Semana_2     4

Boxplots con gráficos tradicionales

par(mfrow=c(2,2))
for(i in levels(datos$variable)){
  dat<-subset(datos,subset = variable==i)
  boxplot(dat$value~dat$Grupo,xlab=i)
}

Usando ggplot2

Esta es una forma elegante y moderna de hacer el gráfico

library(ggplot2)
p1 <- ggplot(datos, aes(x=Grupo, y=value, fill=Grupo)) + geom_boxplot() + facet_wrap(~variable)
print(p1)

Con Lattice

Una forma no tan moderna … pero más elegante que la primera

library(lattice)
#trellis.device(color =FALSE )
print( bwplot(value~Grupo |variable,data=datos,layout = c(2,2)) )