使用arima.sim模拟带有R漂移的arima 1,1,1

我试图使用ARIMA sim包模拟带有漂移的ARIMA模拟。我的问题是我似乎无法让它工作

我需要这样的东西：

但我的代码产生了以下结果：

> mean(datatime)
[1] 15881.56
> sd(datatime)
[1] 8726.893
> length(datatime)
[1] 123

#原始序列的平均值和方差
原始经验指数=15881.56
原始经验变量=8726.893*8726.893
原始经验n=123
#使用arima.sim进行模拟
ts.simARIMA过程根据定义没有任何漂移/趋势。受交叉验证答案的启发，并考虑到您想要的价值：

set.seed(123)
intercept <- 4500
b <- (32000 - intercept) / 123
x <- 1:123
y <- b * x + arima.sim(model=list(c(1, 0, 1)),
                    n = 123, mean=intercept, sd=2000)


> sd(y)
[1] 8020
> mean(y)
[1] 18370

set.seed（123）
拦截
set.seed(123)
intercept <- 4500
b <- (32000 - intercept) / 123
x <- 1:123
y <- b * x + arima.sim(model=list(c(1, 0, 1)),
                    n = 123, mean=intercept, sd=2000)


> sd(y)
[1] 8020
> mean(y)
[1] 18370