-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathload.R
38 lines (34 loc) · 2.23 KB
/
load.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
library(tidyverse)
library(rjags)
D14_Mackenzie <- data.frame(
year_cal = c(2006L, 2007L, 2005L, 2002L, 1997L, 1989L, 1978L, 1976L, 1974L,
1970L, 1968L, 1964L, 1958L, 1949L, 1933L, 1918L, 1906L, 1894L),
D14C = c(-509, -636, -589, -489, -455, -471, -529, -505, -467, -453,
-436, -481, -508, -556, -536, -519, -657, -499),
sd = c(18, 104, 7, 7, 9, 8, 12, 14, 11, 12, 16, 7, 14, 8, 7, 7, 27, 14)) %>%
mutate(year_BP = 1950 - year_cal)
D14_atm <- data.frame(
year_cal = c(2000L, 1999L, 1998L, 1997L, 1996L, 1995L,
1994L, 1993L, 1992L, 1991L, 1990L, 1989L, 1988L, 1987L, 1986L,
1985L, 1984L, 1983L, 1982L, 1981L, 1980L, 1979L, 1978L, 1977L,
1976L, 1975L, 1974L, 1973L, 1972L, 1971L, 1970L, 1969L, 1968L,
1967L, 1966L, 1965L, 1964L, 1963L, 1962L, 1961L, 1960L, 1959L,
1958L, 1957L, 1956L, 1955L, 1954L, 1953L, 1952L, 1951L, 1950L,
1949L, 1948L, 1947L, 1946L, 1945L, 1944L, 1943L, 1942L, 1941L,
1940L, 1939L, 1938L, 1937L, 1936L, 1935L, 1934L, 1933L, 1932L,
1931L, 1930L, 1929L, 1928L, 1927L, 1926L, 1925L, 1924L, 1923L,
1922L, 1921L, 1920L, 1919L, 1918L, 1917L, 1916L, 1915L, 1914L,
1913L, 1912L, 1911L, 1910L, 1909L, 1908L, 1907L, 1906L, 1905L,
1904L, 1903L, 1902L, 1901L, 1900L),
D14C = c(86.25, 89.24, 97, 99.74, 102.45, 111.07, 117.39, 123.33, 131.03, 136.21, 148.28,
160.26, 167.54, 175.32, 183.16, 194.43, 204.86, 225.61, 239.93,
256.41, 267.13, 297.19, 323.21, 332.03, 350.46, 370.14, 392.2,
420.27, 465.53, 498.72, 529.57, 545.25, 565.33, 623.34, 691.59,
754.5, 835.33, 713.08, 361.17, 221.29, 212.58, 228.33, 148, 45,
20, 15, -22.5, -24.1, -25.8, -25.5, -25.8, -26, -22.1, -21.6,
-22, -22.4, -23.1, -24.3, -20.4, -19.9, -23, -20.1, -16.2, -17.2,
-16.5, -16.7, -15.6, -18.6, -20.7, -16.6, -14.5, -18.2, -15.5,
-15.7, -14.7, -12.3, -11.4, -14.1, -12.2, -12.9, -14.1, -11.1,
-11.3, -9.7, -11, -6.3, -7.1, -7.3, -8.1, -8.3, -7.5, -6.5, -8.4,
-6.2, -4.4, -5.5, -4.3, -4.2, -2.9, 1.2, -2.9)) %>%
mutate(year_BP = 1950 - year_cal)