# 12-01-2016: lists parameters from Tab. 2 using the summary(lm*) function in R # "*" denotes not significant at 0.01 significance level infile <- "data_dx500_dz10_1.1.dat" source("read_mapped_data.R") # CA n=907 lm.deltanat.co3.temp.depAll <- lm(c13sednatb - c13wcnatb ~ co3b + tempb + depb) # r = 0.17, sigma = 0.29, intercept = 0.51+-0.04, slope_co3 = -0.0022+-0.0003, slope_temp = -0.019+-0.005, slope_dep = (-7.7+-0.8)e-5 lm.delta.co3.temp.depAll <- lm(c13sedb - c13wcb ~ co3b + tempb + depb) # r = 0.16, sigma = 0.29, intercept = 0.42+-0.04, slope_co3 = -0.0011+-0.0003, slope_temp = 0.003+-0.005*, slope_dep = (-9.2+-0.8)e-5 lm.deltanat.co3.depAll <- lm(c13sednatb - c13wcnatb ~ co3b + depb) # CA1 r = 0.15, sigma = 0.29, intercept = 0.49+-0.04, slope_co3 = -2.4+-0.3, slope_dep = (-7.1+-0.8)e-5 lm.delta.co3.depAll <- lm(c13sedb - c13wcb ~ co3b + depb) # CA2 r = 0.16, sigma = 0.29, intercept = 0.42+-0.04, slope_co3 = -1.0+-0.3, slope_dep = -9.3+-0.8 lm.deltanat.co3All <- lm(c13sednatb - c13wcnatb ~ co3b) # CA3 r = 0.07, sigma = 0.31, intercept = 0.35+-0.04, slope_co3 = -0.0027+-0.0003 lm.deltanat.tempAll <- lm(c13sednatb - c13wcnatb ~ tempb) # r = 0.01, sigma = 0.32, intercept = 0.07+-0.01, slope_temp = -0.02+-0.004 lm.deltanat.depAll <- lm(c13sednatb - c13wcnatb ~ depb) # CA4 r = 0.10, sigma = 0.30, intercept = 0.23+-0.02, slope_dep = (-7.7+-0.8)e-5 lm.deltanat.1All <- lm(c13sednatb - c13wcnatb ~ 1) # sigma = 0.32, intercept = 0.04+-0.01 lm.deltanat.0All <- lm(c13sednatb - c13wcnatb ~ 0) # CA5 sigma = 0.32 lm.c13wcnat.All <- lm(c13sednatb ~ c13wcnatb) # CA6 r = 0.60, sigma = 0.31, intercept = 0.14+-0.02, slope_c13wc = 0.86+-0.02 lm.c13wc.All <- lm(c13sedb ~ c13wcb) # CA7 r = 0.53, sigma = 0.31, intercept = 0.14+-0.02, slope_c13wc = 0.90+-0.03 # CL n=507 lm.deltanat.co3.temp.depLH <- lm(c13sednatLHb - c13wcnatLHb ~ co3LHb + tempLHb + depLHb) # r = 0.19, sigma = 0.25, intercept = 0.33+-0.06, slope_co3 = -0.0030+-0.0003, slope_temp = 0.017+-0.008*, slope_dep = -1.6+-1.4* lm.delta.co3.temp.depLH <- lm(c13sedLHb - c13wcLHb ~ co3LHb + tempLHb + depLHb) # r = 0.21, sigma = 0.25, intercept = 0.37+-0.05, slope_co3 = -0.0013+-0.0003, slope_temp = 0.02+-0.008*, slope_dep = -8.0+-1.4* lm.deltanat.co3.depLH <- lm(c13sednatLHb - c13wcnatLHb ~ co3LHb + depLHb) # CL1 r = 0.19, sigma = 0.25, intercept = 0.41+-0.04, slope_co3 = -0.0028+-0.0003, slope_dep = -3.8+-1.0 lm.delta.co3.depLH <- lm(c13sedLHb - c13wcLHb ~ co3LHb + depLHb) # CL2 r = 0.21, sigma = 0.25, intercept = 0.46+-0.04, slope_co3 = -0.0010+-0.0003, slope_dep = -10.4+-1.0 lm.deltanat.co3LH <- lm(c13sednatLHb - c13wcnatLHb ~ co3LHb) # CL3 r = 0.17, sigma = 0.26, intercept = 0.34+-0.04, slope_co3 = -0.0030+-0.0003 lm.deltanat.tempLH <- lm(c13sednatLHb - c13wcnatLHb ~ tempLHb) # r = 0.01, sigma = 0.28, intercept = -0.05+-0.02, slope_temp = 0.017+-0.006* lm.deltanat.depLH <- lm(c13sednatLHb - c13wcnatLHb ~ depLHb) # CL4 r = 0.05, sigma = 0.27, intercept = 0.14+-0.03, slope_dep = -5.5+-1.1 lm.deltanat.1LH <- lm(c13sednatLHb - c13wcnatLHb ~ 1) # sigma = 0.28, intercept = 0.005+-0.01* lm.deltanat.0LH <- lm(c13sednatLHb - c13wcnatLHb ~ 0) # CL5 sigma = 0.28 lm.c13wcnatLH <- lm(c13sednatLHb ~ c13wcnatLHb) # CL6 r = 0.69, sigma = 0.25, intercept = 0.14+-0.02, slope_c13wc = 0.76+-0.02 lm.c13wcLH <- lm(c13sedLHb ~ c13wcLHb) # CL7 r = 0.65, sigma = 0.27, intercept = 0.16+-0.02, slope_c13wc = 0.84+-0.03 # CW n = 247 lm.deltanat.co3.depWUE <- lm(c13sednatWUEb - c13wcnatWUEb ~ co3WUEb + depWUEb) # CW1 r = 0.25, sigma = 0.22, intercept = 0.44+-0.05, slope_co3 = -0.0029+-0.0004, slope_dep = -3.1+-1.4 lm.deltanat.0WUE <- lm(c13sednatWUEb - c13wcnatWUEb ~ 0) # CW5 sigma = 0.27 lm.deltanat.1WUE <- lm(c13sednatWUEb ~ c13wcnatWUEb) # CW6 r = 0.76, sigma = 0.22, intercept = 0.19+-0.02, slope_c13wc = 0.74+-0.03