Определение границ бокса вокруг барплота в R

Question

Я создал барплот с 24 барами по оси X (0-23) и фоновым затенением, используя следующий код:

#Data
Hours = seq(from=0, to=23)
Mean = rnorm(24, mean=5, sd=2)

#Create number seq for tick mark locations
at_tick = seq_len(length(Hours)+1) 

#Plot with background rectangle shading
x=barplot(Mean,names.arg=Hours, border="white", ylab="Freq", xlab="Hour", 
          ylim=c(0,10), axes=FALSE, space=0, col="grey50")
X = c(0,5)
Y = c(0,10)
rect(X[1], Y[1], X[2], Y[2], border = "gray80", col = "gray80")
X2 = c(19,24)
Y2 = c(0,10)
rect(X2[1], Y2[1], X2[2], Y2[2], border = "gray80", col = "gray80")
barplot(Mean,names.arg=Hours, ylim=c(0,10), border="white", ylab="", xlab="", axes=FALSE, space=0, col="gray50", add=TRUE) 
axis(2, las=2, pos=0)
axis(1, at = at_tick -1, pos=0, labels = FALSE)

box(which="plot", bty="]") #add a box around the plot

Это создает график с окружающим прямоугольником, который выходит за пределы оси X в обоих направлениях. Вместо этого я хотел бы добавить рамку вокруг графика, которая выравнивается по границам оси (т. Е. Ось x: 0-23, ось y: 0-10). Я потратил целую вечность, пытаясь найти способ сделать это без удачи. Любая помощь могла бы быть полезна. Спасибо!

Answer 1

Как насчет рисования отдельных линий? Вы можете использовать segment функция вместо box сделать это:

segments(24,10, 24,0)
segments(0,10, 24,10)

Полный код:

#Data
Hours = seq(from=0, to=23)
Mean = rnorm(24, mean=5, sd=2)

#Create number seq for tick mark locations
at_tick = seq_len(length(Hours)+1) 

#Plot with background rectangle shading
x=barplot(Mean,names.arg=Hours, border="white", ylab="Freq", xlab="Hour", 
          ylim=c(0,10), axes=FALSE, space=0, col="grey50")
X = c(0,5)
Y = c(0,10)
rect(X[1], Y[1], X[2], Y[2], border = "gray80", col = "gray80")
X2 = c(19,24)
Y2 = c(0,10)
rect(X2[1], Y2[1], X2[2], Y2[2], border = "gray80", col = "gray80")
barplot(Mean,names.arg=Hours, ylim=c(0,10), border="white", ylab="", xlab="", axes=FALSE, space=0, col="gray50", add=TRUE) 
axis(2, las=2, pos=0)
axis(1, at = at_tick -1, pos=0, labels = FALSE)

segments(24,10, 24,0)
segments(0,10, 24,10)