| 所在主题: | |
| 文件名: electricity.txt | |
| 资料下载链接地址: https://bbs.pinggu.org/a-2106835.html | |
| 附件大小: | |
|
大神们好:
小弟学习R语言不久,现在在看时间序列方面的书,其中有3个问题吧,一直想了好久不知道怎么解决,现在向大神们求助: 1, 为发电量时间序列,我想对其进行预测分析,现在代码如下: data("electricity") plot(log(electricity)) month. <- season(electricity) trend <- time(electricity) lm.sol <- lm(log(electricity) ~ month. + trend - 1) summary(lm.sol) lines(x = as.vector(time(electricity)), y = fitted(lm.sol), col = 'red', lty = 2, type = 'l', lwd = 2) 从图中看模拟结果还算不错, 黑线为实际发电量,红虚线为拟合值。 现在的疑问是不知道怎么进行预测,predict函数不会写。。。 2,针对这个题目,有这样的参考答案: data(electricity); month.=season(electricity); trend=time(electricity) model=arima(log(electricity),order=c(0,0,0), xreg=as.matrix(model.matrix(~month.+trend-1))[,-1]) newmonth.=season(ts(rep(1,24),start=c(2006,1),freq=12)) newtrend=time(electricity)[length(electricity)]+(1:24)*deltat(electricity) result=plot(model,n.ahead=24,n1=c(2001,1),ylab='Electricity',xlab='Year', newxreg=as.matrix(model.matrix(~newmonth.+ newtrend-1))[,-1],pch=19,transform=exp) 请问标红的[,-1]是怎么回事?为什么要减掉这一列呢? 3,关于acf函数,为什么有的时间序列就能直接用acf函数,而有的却需要as.vector一下呢? 例arima为时间序列,为什么有的可以直接acf(arima),而有的却需要acf(as.vector(arima))呢? 问的有点乱,请大神解答,谢过。。。 |
|
熟悉论坛请点击新手指南
|
|
| 下载说明 | |
|
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。 2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。 3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明