Укажите параметры временных рядов для данных за миллионы лет

У меня есть данные о дельте 18O за 6 миллионов лет, с одной точкой данных каждые 0,003 млн лет, например:

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, то есть масштаб не в миллионах лет. Вот сюжет

Как это исправить? Мне нужно, чтобы это был временной ряд, потому что мне нужно провести спектральный анализ на нем, например периодограммы


r time-series spectral
person Nicola Gambaro    schedule 02.01.2020    source источник
comment
Мое первое предположение состоит в том, что проблема в том, что тип данных, который R использует для хранения даты / времени, не может возвращаться более чем на несколько тысяч лет или что-то в этом роде. Возможно, вам придется рассматривать ось x просто как числовое значение, а не как дату / время. Если вы покажете свой код, это поможет другим понять, что происходит, и, следовательно, поможет другим помочь вам.   -  person Robert Dodier    schedule 02.01.2020

Ответы (1)


arrow_upward
2
arrow_downward

Вам необходимо указать периодичность наблюдений. У вас будет 1 наблюдение каждые 3000 лет, поэтому частота будет 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")

введите здесь описание изображения

person Allan Cameron    schedule 02.01.2020