arfima模型“'x'是一个列表,但没有组件'x'和'y'”

arfima模型“'x'是一个列表,但没有组件'x'和'y'”,r,plot,R,Plot,R软件 set.seed(1000) sim <- arfima.sim(1000, model = list(phi = c(0.2, 0.1), dfrac = 0.4, theta = 0.9)) fit <- arfima(sim, order = c(2, 0, 1), back=T) plot(fit) xy.coordsx、y、xlabel、ylabel、log中出错: “x”是一个列表,但没有组件“x”和“y” 我怎样才能使曲线拟合?我认为问题在于拟合的值应该是y值

R软件

set.seed(1000)
sim <- arfima.sim(1000, model = list(phi = c(0.2, 0.1), dfrac = 0.4, theta = 0.9))
fit <- arfima(sim, order = c(2, 0, 1), back=T)
plot(fit)
xy.coordsx、y、xlabel、ylabel、log中出错:

“x”是一个列表,但没有组件“x”和“y”

我怎样才能使曲线拟合?我认为问题在于拟合的值应该是y值,而x是表示观测值的数字1,2,3,4。有人能解决这个问题吗

我想绘制我从arfima模型得到的值。就像我试图绘制实际值和拟合值一样,问题是:您想绘制什么?如果是包含模型预测的时间序列数据,请尝试:

set.seed(1000)
sim <- arfima.sim(1000, model = list(phi = c(0.2, 0.1), dfrac = 0.4, theta = 0.9))
fit <- arfima(sim, order = c(2, 0, 1), back=T)
plot(predict(fit, n.ahead = 5))

不,我想根据我得到的arfima模型绘制拟合值,这样我就可以将其与原始模型进行比较values@korzuen可以使用fittedfit从ARFIMA fit对象中为每个模式提取拟合值。这些都是来自的有很好文档记录的函数。不客气@korzuen。请考虑通过在答案旁边设置复选标记来结束问题。
library(tidyverse);
cbind.data.frame(Data = sim, fitted(fit)) %>%
    rowid_to_column("idx") %>%
    gather(what, value, -idx) %>%
    ggplot(aes(idx, value, colour = what)) +
    facet_wrap(~ what, ncol = 1, scales = "free_y") +
    geom_line()