第一次发帖请教，因为是在西班牙 所以变量名称用的西班牙语。措辞也不是很准确
`
source("fitModels.R")
edad<-20:90
periodo<-1950:2009
#女士数据拟合李卡特模型
ajuste_m<-fit701(edad,periodo,t(esp_m[edad+1,]),t(muerto_m[edad+1,]),matrix(1,length(periodo),length(edad)))
plot(ajuste_m$kappa,col="red",type='l',ylab="kappas",main = "kappas de mujeres")
kappa_mujer<-ajuste_m$kappa
#男士数据拟合李卡特模型
ajuste_h<-fit701(edad,periodo,t(esp_h[edad+1,]),t(muerto_h[edad+1,]),matrix(1,length(periodo),length(edad)))
kappa_hombre<-ajuste_h$kappa
##转化为时间序列
kappa_t_mujer<-ts(kappa_mujer,frequency=1, start=1950)
kappa_t_hombre<-ts(kappa_hombre,frequency=1, start=1950)
预测
require(forecast)
#女士
ARIMAfit_m<- auto.arima(kappa_t_mujer, approximation=FALSE,trace=FALSE)
summary(ARIMAfit_m)
#显示为 ARIMA(1,1,0)
#男士kappa的拟合
ARIMAfit_h<- auto.arima(kappa_t_hombre, approximation=FALSE,trace=FALSE)
summary(ARIMAfit_h)
#显示为ARIMA(1,2,1)
#下面开始预测 奇怪的是 男士的可以预测 女士的报错 报错内容为数据不是向量(vector )
(pred_m <-predict(ARIMAfit_m, n.ahead = 3))
(pred_h <-predict(ARIMAfit_h, n.ahead = 3))
`
在其他地方也有查到说有出现过这种情况 建议该用forecast()函数 但是我想知道出现这种情况的愿意。
赐教
不知道插入图片成功没。。。