已成功建立了时间序列模型，如何给模型曲线上下两侧加上95%可信区间（注意，这里不是给预测值加可信区间，我的意思是给所有模型拟合出来的值加上可信区间），如下图（这个是我看人家用metlab做的，我很想知道R怎么做，因为好多教程只是给预测值加上了可信区间），R应该用什么函数来做呢？求教大师。