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

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 年，即尺度不是以百万年为单位。 Here's the plot

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

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")

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

Specify time series parameters for data over millions of years

r

time-series

spectral