R 中第一个差分记录系列的持续集成

Question

我有一个第一个差异记录系列，我需要将其转换回原始级别单位。我如何在 R 中做到这一点？

以下是我的数据系列和尝试的代码：

原创系列：

1.1, 2.6, 3.6, 4.8, 5.1, 6.0, 7.3, 8.8, 9.4, 10.5

尝试的代码：

o <- c(1.1, 2.6, 3.6, 4.8, 5.1, 6.0, 7.3, 8.8, 9.4, 10.5)
l <- log(o)
dl <- diff(l)

exp(diffinv(dl, differences = 1)) # Attempt to recover o

电流输出：

[1] 1.000000 2.363636 3.272727 4.363636 4.636364 5.454545 6.636364 8.000000 
[9] 8.545455 9.545455

期望的输出：

[1] 1.1 2.6 3.6 4.8 5.1 6.0 7.3 8.8
[9] 9.4 10.5

Answer 1

在不知道第一个元素的情况下，不可能完全恢复原始系列：

exp(diffinv(dl) + log(1.1))
# [1]  1.1  2.6  3.6  4.8  5.1  6.0  7.3  8.8  9.4 10.5

或

exp(cumsum(c(log(1.1), dl)))
# [1]  1.1  2.6  3.6  4.8  5.1  6.0  7.3  8.8  9.4 10.5

---

设X₁, X₂, X₃, X₄ 为原序列。那么最终你有

Z₁ = ln(X₂) - ln(X₁),

Z₂ = ln(X₃) - ln(X₂)

Z₃ = ln(X₄) - ln(X₃)

给定 Z₁, Z₂, Z₃，你希望恢复 X₁, X2、X₃、X₄。然后注意

Z₁ = ln(X₂) - ln(X₁),

Z₂ + Z₁ = ln(X₃) - ln(X₁),

Z₃ + Z₂ + Z₁ = ln(X₄) - ln(X1),

这样

Z₁ + ln(X₁) = ln(X₂),

Z₂ + Z₁ + ln(X₁) = ln(X₃),

Z₃ + Z₂ + Z₁ + ln(X₁) = ln(X4),

和

exp(Z₁ + ln(X₁)) = X₂,

exp(Z₂ + Z₁ + ln(X₁)) = X₃,

exp(Z₃ + Z₂ + Z₁ + ln(X₁)) = X ₄,

这正是（比第一个解决方案更透明）

exp(cumsum(c(log(1.1), dl)))

确实如此。因此，

 exp(diffinv(dl))

单独仅在原始系列从 1 = exp(0) 开始时才有效。

因此，您现在必须将系列的初始级别从任何类型的差异中恢复，这是非常自然的。想象一下，我们只知道您每天将赚取和花费多少。在不知道初始金额的情况下，绝不可能说出您总共有多少钱。