-
Notifications
You must be signed in to change notification settings - Fork 0
/
R_timewise_iso_plots.R
133 lines (114 loc) · 6.04 KB
/
R_timewise_iso_plots.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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
#March 2019
#Joy Buongiorno
#Analysis of ndividual oncolite carbonate isotope values
library(ggplot2)
iso_cross_cols<-c( "d18O" = "blue", "d13C" = "red")
########### LN-13-1 #################
LN.13.1<-read.csv("Isotopes_LN.13.1.txt", header=TRUE, sep="\t")
carbon_iso_1<- ggplot(data= LN.13.1, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,16), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_1<- ggplot(data= LN.13.1, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(-3,6), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()
########### LN-13-5 #################
LN.13.5<-read.csv("Isotopes_LN.13.5.txt", header=TRUE, sep="\t")
carbon_iso_5<- ggplot(data= LN.13.5, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,16), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_5<- ggplot(data= LN.13.5, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(-3,6), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()
########### LN-13-10 #################
LN.13.10<-read.csv("Isotopes_LN.13.10.txt", header=TRUE, sep="\t")
carbon_iso_10<- ggplot(data= LN.13.10, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,16), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_10<- ggplot(data= LN.13.10, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(-3,6), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()
########### LN-13-13 #################
LN.13.13<-read.csv("Isotopes_LN.13.13.txt", header=TRUE, sep="\t")
carbon_iso_13<- ggplot(data= LN.13.13, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,16), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_13<- ggplot(data= LN.13.13, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(-3,6), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()
########### LN-13-20 #################
LN.13.20<-read.csv("Isotopes_LN.13.20.txt", header=TRUE, sep="\t")
carbon_iso_20<- ggplot(data= LN.13.20, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,16), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_20<- ggplot(data= LN.13.20, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(-3,6), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme(text = element_text(size=30)) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()
########### LN-13-28 #################
LN.13.28<-read.csv("Isotopes_LN.13.28.txt", header=TRUE, sep="\t")
carbon_iso_28<- ggplot(data= LN.13.28, aes(x=Layer, y=d13C)) +
geom_point(aes(fill=d13C), pch=23, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{13}, "C (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(6,20),labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme_bw(base_size=30) +
scale_fill_gradient(low = "lightblue1", high= "midnightblue") +
scale_x_reverse() +
coord_flip()
oxygen_iso_28<- ggplot(data= LN.13.28, aes(x=Layer, y=d18O)) +
geom_point(aes(fill=d18O), pch=21, colour = "black", size=5) +
labs(x="Layer from outside", y=expression(paste(delta^{18}, "O (\u2030 vs. PDB)"))) +
scale_y_continuous(limits=c(0,10), labels = function(x) ifelse(x>0, paste0("+", x), x)) +
theme_bw(base_size=30) +
scale_fill_gradient(low = "mistyrose", high= "orangered4") +
scale_x_reverse() +
coord_flip()