Вы можете подойти очень близко, вопрос в том, зачем вам это нужно? В любом случае, попробуйте что-нибудь вроде этого:
x <- rnorm(100)
# calculate the axis limits
limitsX <- c(-1, length(x))
limitsY <- range(x) * 1.1
# get pretty ticks
tickMarksX <- pretty(0:length(x))
tickMarksY <- pretty(x)
# plot
plot(NULL, xlim=limitsX, ylim=limitsY, bty = 'n', xaxs = 'i', xaxt = 'n', yaxt = 'n', yaxs = 'i')
rect(limitsX[1], limitsY[1], limitsX[2], limitsY[2], col = "#E0E0E0", border = NA)
axis(1, at = tickMarksX, lwd = 0, ticks.lwd = 1)
axis(2, at = tickMarksY, lwd = 0, ticks.lwd = 1)
# major x grid lines
abline(v = tickMarksX, col = "white", lwd=2)
# minor x grid lines
abline(v = tickMarksX-(tickMarksX[2]-tickMarksX[1])/2, col = "white")
# major y grid lines
abline(h = tickMarksY, col = "white")
# minor y grid lines
abline(h = tickMarksY-(tickMarksY[2]-tickMarksY[1])/2, col = "white")
lines(x = 1:length(x), y = x)
введите описание изображения здесь