Displaying 1 result from an estimated 1 matches for "ymod1".
Did you mean:
mod1
2010 Jan 18
5
errors appears in my time Series regression fomula
...3 matches"
> y=scan('F:/Math355/AIRLINE.txt',skip=2,nlines=18)
Read 144 items
> yts=ts(log(y))
> plot(yts,main="LOG AIRLINE TOTALS",type='p',col=2)
> n=length(y)
> time=seq(1:n)
>
> month=c(rep(seq(1:12),12))
> fmonth=as.factor(month)
> ymod1=lm(yts~time+fmonth)
> summary(ymod1)
> yfit1=ymod1$fitted
> lines(yfit1,c=3)
??? plot.xy(xy.coords(x, y), type = type, ...) :
??3?????????????
Also, for a similar data set quarterly sales of spirits series, what is the best model to use?? Should I use ARIMA??
Many thanks...