Я думаю, что вы должны поменять порядок logTL
и logMass
на вашем графике, т.е.
plot(logTL, logMass, xlab="log (base 10) total length", ylab="log (base 10) mass")
, так как вы сделали регрессию logMass
относительно logTL
, т.е. lma<-lm(logMass~logTL)
В противном случае вы можете адаптировать его по-другому, например, abline(1/coefficients(lma))
.