Com os dados da pergunta aqui estão duas maneiras de traçar o gráfico com a curva prevista.
#R base
R base
O problema das linhas emaranhadas resolve-se ordenando a variável do eixo dos x e, depois, usar essa variável ordenada para calcular os valores de y.
b <- sort(df$body)
new <- data.frame(body = b)
with(df, plot(body, n.kill))
with(df, lines(b, predict(fit, newdata = new)))
Em vez de ordenar body
, outra maneira, talvez melhor, é criar um vetor para os valores do eixo dos x com seq
.
b <- seq(min(df$body), max(df$body), length.out = 500)
O resto do código é exatamente igual.
#Pacote ggplot2
Pacote ggplot2
Com o pacote ggplot2
, não é necessário ajustar o modelo antes de traçar o gráfico, a geom_smooth
faz isso desde que também se passem os argumentos usados no fit.
Mas cuidado, a geom_smooth
só usa variáveis com os nomes x
e y
, as mesmas de aes()
. É portanto necessário modificar a fórmula para usar esses nomes.
library(ggplot2)
ggplot(df, aes(x = body, y = n.kill)) +
geom_point() +
geom_smooth(method = 'nls',
formula = y ~ a*(1 - exp(-b*x))^c,
se = FALSE,
method.args = list(start = list(a = ca, b = 0.05, c = 1),
algorithm = 'port',
control = nls.control(maxiter = 100, tol = eps)))