r - 为数百万年的数据指定时间序列参数

Question

我有超过 600 万年的 foram delta 18O 数据，每 0.003 Myrs 有一个数据点，如下所示：

Age (Myr) | d18O
  0       | 3.43
  0.003   | 3.37
  0.006   | 3.54
  0.009   | 3.87
  0.012   | 4.36
  0.015   | 4.90
  0.018   | 5.01
  0.021   | 4.96
  0.024   | 4.87
  0.027   | 4.67
  0.03    | 4.58

ts_d18O <- ts(d18O[,2])
plot(ts_d18O)

然而，当我告诉 R 将其视为一个时间序列并绘制它时，我在 x 轴上得到了一直到 2000 的值，即规模不是以百万年为单位。这是情节

我该如何解决？我需要它是一个时间序列，因为我必须对其进行光谱分析，例如周期图

Answer 1

您需要指定观察的频率。您每 3000 年有 1 次观察，因此如果您想要以年为单位的答案，频率将为 1/3000。我假设你在数百万年后想要你的答案。这意味着您的频率是 1/0.003 或 333.333（即这是每百万年采集的样本数量）。

您还应该指定开始时间（本例中的 -5 表示 5 Mya）。

最后，您可以根据需要标记您的 x 轴。

ts_d18O <- ts(d18O[,2], start = c(0.003, -5 / 0.003), frequency = 1/0.003)
ts
#> Time Series:
#> Start = -5 
#> End = -4.97 
#> Frequency = 333.333333333333 
#>  [1] 3.43 3.37 3.54 3.87 4.36 4.90 5.01 4.96 4.87 4.67 4.58
plot(ts_d18O, xlab = "Million years ago")