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