[R] Interaction plots (six on one page)

[Richard M. Heiberger]

I would do this as a lattice plot.  Continuing with your data:

tmp <- data.frame(sapply(data, tapply, data[1:2], mean))
tmp$time <- factor(tmp$time)
xyplot(thanaa+thalcho+thalino+ponaa+pocho+poino ~ time,
       group=BMIakt, data=tmp, type="l", scales=list(relation="free"),
       auto.key=list(title="BMIakt", border=TRUE))

Rich