data(eaglecol)
plot(eaglecol)
## Not run: #confidence ellipse
# library("ellipse") #needs this package!
# ansTFD <- artfima(eaglecol)
# v <- ansTFD$varbeta
# bHat <- c(ansTFD$dHat, ansTFD$lambdaHat)
# xy <- ellipse(v, centre=bHat, level=0.9)
# plot(xy, type="l", lwd=2, xlab=expression(delta), ylab=expression(lambda))
# points(matrix(bHat,ncol=2), pch=16, cex=3, col="blue")
# #setwd("D:/DropBox/R/2016/artfima/Explore_ts_data/eaglecol")
# #postscript(file="eaglecolCI.eps")
# #plot(xy, type="l", lwd=2, xlab=expression(delta), ylab=expression(lambda))
# #points(matrix(bHat,ncol=2), pch=16, cex=3, col="blue")
# #graphics.off()
# ## End(Not run)
## Not run: #forecast comparison
#
# ## End(Not run)
Run the code above in your browser using DataLab