"Achse" wird keine X-Achse zu Boxplot hinzufügen

Question

Ich versuche, ein Boxplot mit benutzerdefinierten Achsenbeschriftungen zu machen, aber ich kann nicht scheinen, eine X-Achse zum Plot hinzuzufügen.

Zum Beispiel:

test <- data.frame(year=as.integer(rep(1963:2014, each=10)), 
     response=rnorm(520)) 
boxplot(response~year, data=test, ylim=c(-3,3), xlab="", ylab="", 
    range=0, xaxt="n", yaxt="n") 
responselabs <- as.numeric(c(-3:3, by=1)) 
yearlabs <- as.integer(seq(1965,2015, by=5)) 
axis(2, at=responselabs, tck=0.03, las=1) 
axis(1, at=yearlabs, tck=0.03) 

gibt den boxplot, aber keine x-Achsenbeschriftungen:

Den Versuch, es in die andere Richtung zu hacken, indem zunächst ein leeres Grundstück zu machen, kann ich Holen Sie sich die Achsen, aber es wird nicht das Boxplot hinzugefügt:

plot(NA, ylim=c(-3, 3), xlim=c(1962, 2015), xaxt="n", yaxt="n", ylab="", xlab="") 
axis(2, at=responselabs, tck=0.03, las=1) 
axis(1, at=yearlabs, tck=0.03) 
boxplot(response~year, data=test, ylim=c(-3,3), xlab="", ylab="", 
    range=0, xaxt="n", yaxt="n", add=T) 

Was ist hier los?

Answer 1

Ich denke, was passiert ist, dass boxplotyear zu einem factor umwandelt. Wir können mithilfe der labels Argument in axis dieses Problem umgehen:

boxplot(response~year, data=test, ylim=c(-3,3), xlab="", ylab="", 
    range=0, xaxt="n", yaxt="n") 
responselabs <- as.numeric(c(-3:3, by=1)) 
yearlabs <- as.integer(seq(1965,2015, by=5)) 
axis(2, at=responselabs, tck=0.03, las=1) 
axis(1, at = yearlabs - 1962, labels = yearlabs) 

Answer 2

Warum nicht ggplot2 verwenden?

library(ggplot2) 
p<-ggplot(test,aes(x=year,y=response,group=year))+ 
geom_boxplot()+ 
scale_x_continuous(breaks=round(seq(min(test$year),max(test$year),by=5),1)) 

Wenn Sie im scale_x_continuous() Argument anzupassen der Code ist ziemlich einfach auf die nächsten 5 dann gerundet werden.

p<-ggplot(test,aes(x=year,y=response,group=year))+ 
    geom_boxplot()+ 
    scale_x_continuous(breaks = round(seq(round(min(test$year)/5,0)*5,round(max(test$year)/5,0)*5, by = 5),1)) 

Answer 3

Wie andere gesagt haben, die Werte verwechselt haben die x-Achse zu Grunde liegen. Herumspielen mit abline() zeigt, dass abline(v=2010) nicht angezeigt wird, aber so etwas wie abline(v=50) wird angezeigt.

Hier ist ein Plot Sie nach, glaube ich (mit einigen zusätzlichen Änderungen):

year_vals <- as.integer(rep(1963:2014, each=10)) # pulling out of next line for easy reference 
test <- data.frame(year=year_vals,response=rnorm(520)) 
boxplot(response~year, data=test, ylim=c(-3,3), xlab="", ylab="", range=0, xaxt="n", yaxt="n") 
responselabs <- as.numeric(c(-3:3, by=1)) 
yearlabs <- year_vals[year_vals%%5==0] # alternate, more general definition; note that 2015 wasn't in the original 'test' data.frame #as.integer(seq(1965,2015, by=5)) 
axis(2, at=responselabs, tck=0.03, las=1) 
# axis(1, at=1:length(yearlabs), tck=0.03) 
axis(1, at=which(unique(year_vals)%in%yearlabs), labels=unique(yearlabs), tck=0.03)