# \donttest{
require(minpack.lm)
require(ggplot2)
## Temperature response example
data(maizeleafext)
## Fit model
fit <- nlsLM(rate ~ SSbeta5(temp, mu, tb, a, tc, b), data = maizeleafext)
## Visualize
ndat <- data.frame(temp = 0:45)
ndat$rate <- predict(fit, newdata = ndat)
ggplot() +
geom_point(data = maizeleafext, aes(temp, rate)) +
geom_line(data = ndat, aes(x = temp, y = rate))
# }
Run the code above in your browser using DataLab