Sergio, boa tarde!
Você teria um conjunto de dados pra testar? No teste que fiz usei o logaritmo dos valores de dados para representar no gráfico, substituindo a notação do eixo (-1, 0, 1, 2) pelos valores na escala original (.1, 1, 10, 100).
### <code r>
dados <- structure(c(1, 3, 4, 5, 8, 13, 25, 32, 36, 45, 55, 35, 12, 10,
8, 7, 3, 2, 1, NA, 66, 29, 15, 15, 14, 12, 3, 0.5, 0.4, 0.3,
72, 50, 35, 33, 23, 16, 4, 1, NA, NA), .Dim = c(10L, 4L), .Dimnames = list(
NULL, c("rank", "alto", "medio", "baixo")))
dados.log <- cbind(dados[,1], log10(dados[,2:4]))
log10(c(.1, 1, 10, 100)) # [1] -1 0 1 2
plot(1:50, seq(-1, 2, len=50), type="n", xlab="Ranking", ylab="Abundância relativa %", yaxt="n")
axis(side=2, at=c(-1, 0, 1, 2), labels=c(.1,1,10,100), tick = TRUE)
lines(dados.log[,c(1,2)], lty=1, lwd=2, col=2)
lines(dados.log[,c(1,3)], lty=2, lwd=2, col=3)
lines(dados.log[,c(1,4)], lty=3, lwd=2, col=4)
abline(h=c(-1,0,1,2), lty=3, col=1)
legend(35,1.5, legend = c("alto", "médio", "baixo"), col = 2:4, lty=1:3, cex = .75, lwd = 2)
### </code>