为数百万年的数据指定时间序列参数答案

【问题标题】：Specify time series parameters for data over millions of years为数百万年的数据指定时间序列参数
【发布时间】：2020-04-21 10:03:31
【问题描述】：

我有超过 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 的值，即比例不是以百万年为单位。 Here's the plot

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

【问题讨论】：

我的第一个猜测是问题在于 R 用来存储日期/时间的数据类型不能追溯到几千年或类似的时间。您可能不得不将 x 轴视为一个数值，而不是具体的日期/时间。如果您展示您的代码，它将帮助其他人了解正在发生的事情，从而帮助其他人帮助您。

标签： r time-series spectral

【解决方案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")

【讨论】：