R - boxplot con più etichette di fattore

Question

Sto lavorando per provare a creare un boxplot in R-cran che è classificato da due fattori diversi sull'asse x. Il mio problema sta nel creare etichette per un fattore con livelli di +20 che copre l'intero grafico in modo appropriato mentre si utilizza una legenda per etichettare il secondo fattore che ha solo 2 o 3 livelli.

Ecco uno script di test che circa imita il mio attuale set di dati:

d<-data.frame(x=rnorm(1500),f1=rep(seq(1:20),75),f2=rep(letters[1:3],500)) 
# first factor has 20+ levels 
d$f1<-factor(d$f1) 
# second factor a,b,c 
d$f2<-factor(d$f2) 

boxplot(x~f2*f1,data=d,col=c("red","blue","green"),frame.plot=TRUE,axes=FALSE) 

# y axis is numeric and works fine 
yts=pretty(d$x,n=5) 
axis(2,yts) 

# I know this doesn't work; what I'd like is to spread the factors out 
# so the each group of three(a,b,c) is labeled correctly 
axis(1,at=seq(1:20)) 

# Use the legend to handle the f2 factor labels 
legend(1, max(d$x), c("a", "b","c"),fill = c("red", "blue","green")) 

Grazie per qualsiasi aiuto

Answer 1

FWIW, una soluzione ggplot2:

library(ggplot2) 
ggplot(data = d, aes(x = f1, y = x)) + 
    geom_boxplot(aes(fill = f2), width = 0.8) + theme_bw() 

Answer 2

Se si vuole un'etichetta a metà di ogni gruppo di 3 scatole, provare qualcosa di simile :

axis(1,at=seq(2,60,3),labels=1:20,cex.axis=0.7) 

Per generalizzare, questo sarebbe:

groups <- 20 
numbox <- 3 
total <- groups * numbox 
xpoints <- seq(median(1:numbox),total,numbox)