R语言绘制plot图形出错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

#加载程序包
library(tseries)
library(forecast)


#加的可能是错的
library(Metrics)
#加的可能是错的



#读取数据
data <-read.table("yichang.txt",header=FALSE)
data <- as.matrix(data)
#画出data变量图
plot(data)

n<-length(data)
m<-floor(n*0.8)
L<-n-m

#将数据拆分为训练集
#取data的80%进行训练

train<-data[1:m,1]

#tsdisplay用来看ACF和PACF

tsdisplay(train)
auto.arima(train)

#diff():求差和求导

###  differential   ###
s1<-diff(train,12)
adf.test(s1)
tsdisplay(s1)

ndiffs(train)
#######################

###  forecast ###
#nrow:行 ncol：列 matrix:矩阵
fore<- matrix(NA, nrow=L, ncol=1)
#循环拟合，for(var in seq),var:变量名，seq：向量表达式
for (i in 1:L) {

fit1<-arima(data[1:(m+i-1),1],order=c(1,0,1))
tsdiag(fit1)
f.p1<-forecast(fit1,h=3,level=c(99.5))
fore<-Vectorize(fore)
fore[i,1]<-f.p1$mean
}


f.p1$mean
# 7/25 #

plot(fore)

lines(fore,col="green")
lines(data[m+1:n],col="red")

write.table(fore,"fore_mean_h=3.csv",sep=",")
obs<-data[m+1:n]
write.table(obs,"obs_h=3.csv",sep=",")

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：plot R语言 Forecast Library tseries

这个是提示

Error in plot.window(...) : 'ylim'值不能是无限的
此外: Warning messages:
1: In min(x) : min里所有的参数都不存在; 回覆Inf
2: In max(x) : max里所有的参数都不存在；回覆-Inf
>

看看原始数据加载过程中，是否读入R后出现了错误，或者存在缺失值。

xjg1983 发表于 2019-7-26 09:00
看看原始数据加载过程中，是否读入R后出现了错误，或者存在缺失值。

找到了。前面有个Error
Error in fore[i, 1] <- f.p1$mean : 被替换的项目不是替换值长度的倍数
这个是什么意思呢？