# NOT RUN {
library(ggtrendline)
x<-1:5
y<-c(2,4,8,20,25)
xy<-data.frame(x,y)
getInitial(y ~ SSexp2P(x,a,b), data = xy)
## Initial values are in fact the converged values
fitexp2P <- nls(y~SSexp2P(x,a,b), data=xy)
summary(fitexp2P)
prediction <- predFit(fitexp2P , data.frame(x=x), se.fit = TRUE,
level = 0.95, interval = "confidence")
yfitexp2P <- prediction$fit
yfitexp2P # output a matrix of predictions and bounds with column names fit, lwr, and upr.
# }
Run the code above in your browser using DataLab