Dear all, I found really difficult with the time series questions, please help me with this monthly airline series! I have run the following r code, and there is an error appeared at the end. The data files was enclosed in the email. I'm sorry the errors message appeared in chinese, but it says "plot.xy(xy.coords(x, y), type = type, ...) : errors in argument has more than 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 and regard, Cathy Wong _________________________________________________________________ Got a cool Hotmail story? Tell us now -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: spirits Rdata.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100118/c1be065b/attachment.txt> -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: SPIRITS.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100118/c1be065b/attachment-0001.txt>
That line: lines(yfit1,c=3) does not appear in the code that you sent. I assume that you meant 'col=3'. The message may be saying that there are multiple partial matches for 'c' as a parameter. The code you included worked fine on my system. Do you manually input this line? 2010/1/18 snowcat wong <cathy_white725@hotmail.com>> > Dear all, > > > > I found really difficult with the time series questions, please help me > with this monthly airline series! > > > > I have run the following r code, and there is an error appeared at the end. > The data files was enclosed in the email. > > I'm sorry the errors message appeared in chinese, but it says > "plot.xy(xy.coords(x, y), type = type, ...) : > > errors in argument has more than 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 and regard, > > Cathy Wong > > _________________________________________________________________ > Got a cool Hotmail story? Tell us now > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html> > and provide commented, minimal, self-contained, reproducible code. > >-- Jim Holtman Cincinnati, OH +1 513 646 9390 What is the problem that you are trying to solve? [[alternative HTML version deleted]]
Dear all, I am having difficulty to built a model of quarter sales of spirits data, and deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not appeared right at the end. The data and script is enclosed with this email. I am using the harmonic regression model to exam the periodogram of the residuals. I am not sure this step is right or not? Which is the best formula to analysis the trend and seasonality for this data set? AR, MA or ARMA and how to decide? Please help!! many thanks and regard, Cathy Wong _________________________________________________________________ Got a cool Hotmail story? Tell us now -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: SPIRITS.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100119/1be36fc0/attachment.txt> -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: spirit script.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100119/1be36fc0/attachment-0001.txt>
Dear all, I am having difficulty to built a model of quarter sales of spirits data, and deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not appeared right at the end. The data and script is enclosed with this email. I am using the harmonic regression model to exam the periodogram of the residuals. I am not sure this step is right or not? Which is the best formula to analysis the trend and seasonality for this data set? AR, MA or ARMA and how to decide? Please help!! many thanks and regard, Cathy Wong Got a cool Hotmail story? Tell us now _________________________________________________________________ Got a cool Hotmail story? Tell us now -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: SPIRITS.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100119/dc1589df/attachment.txt> -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: spirit script.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100119/dc1589df/attachment-0001.txt>
Dear all, I am having difficulty to built a model of quarter sales of spirits data, and deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not appeared right at the end. The data and script is enclosed with this email. I am using the harmonic regression model to exam the periodogram of the residuals. I am not sure this step is right or not? Which is the best formula to analysis the trend and seasonality for this data set? AR, MA or ARMA and how to decide? Please help!! many thanks and regard, Cathy Wong _________________________________________________________________ Send us your Hotmail stories and be featured in our newsletter -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: SPIRITS.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100121/1570c47a/attachment.txt> -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: spirit script.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100121/1570c47a/attachment-0001.txt>
Cathy,
How does this model look?
[(1-B**4)]Y(T) = 20.767
+[X1(T)][(1-B**4)][(+ 56.1962)] :PULSE 7/
4 I~P00028test
+[X2(T)][(1-B**4)][(+ 74.4301)] :PULSE 9/
4 I~P00036test
+[X3(T)][(1-B**4)][(- 59.9872)] :PULSE 6/
3 I~P00023test
+[X4(T)][(1-B**4)][(- 27.2187)] :PULSE 7/
1 I~P00025test
+ [(1- .435B** 1)]**-1 [A(T)]
--
View this message in context:
http://r.789695.n4.nabble.com/errors-appears-in-my-time-Series-regression-fomula-tp1016593p3054417.html
Sent from the R help mailing list archive at Nabble.com.
Reasonably Related Threads
- RGL 3D curvilinear shapes
- Beginners Question: Make nlm work
- question about the aggregate function with respect to order of levels of grouping elements
- Multiple series plot with different 'type' argument
- Need something like multiple plots in ts.plot but then change x values