#c1-c3(Oe horizon)
library(ggplot2)
summer.2017<-read.csv("summersmrdata.csv",header=TRUE)
summer.2017
##    Sample   Forest   Age Stand Plot Treatment Horizon    T4
## 1       1 Bartlett Young    C1    1         P      Oe 11.62
## 2       2 Bartlett Young    C1    2         N      Oe  9.08
## 3       3 Bartlett Young    C1    3       Con      Oe 10.08
## 4       4 Bartlett Young    C1    4        NP      Oe  5.09
## 5       5 Bartlett Young    C2    1        NP      Oe  6.69
## 6       6 Bartlett Young    C2    2       Con      Oe 12.45
## 7       7 Bartlett Young    C2    3         P      Oe  8.60
## 8       8 Bartlett Young    C2    4         N      Oe  7.55
## 9       9 Bartlett Young    C3    1        NP      Oe  6.18
## 10     10 Bartlett Young    C3    2         P      Oe  9.82
## 11     11 Bartlett Young    C3    3         N      Oe  5.65
## 12     12 Bartlett Young    C3    4       Con      Oe 11.32
summer.2017.plot<- ggplot(summer.2017,aes(x= Treatment,y=T4,colour=Stand))+geom_point()+theme_bw()
summer.2017.plot

#c4toc6 (Oe horizon)
summer.2017b<-read.csv("c4toc6.csv",header=TRUE)
summer.2017b
##    Sample   Forest Age Stand Plot Treatment Horizon    T4
## 1      13 Bartlett Mid    C4    1        NP      Oe  9.51
## 2      14 Bartlett Mid    C4    2         N      Oe 11.11
## 3      15 Bartlett Mid    C4    3       Con      Oe 14.26
## 4      16 Bartlett Mid    C4    4         P      Oe 17.44
## 5      17 Bartlett Mid    C5    1       Con      Oe 10.04
## 6      18 Bartlett Mid    C5    2        NP      Oe  8.66
## 7      19 Bartlett Mid    C5    3         N      Oe  9.17
## 8      20 Bartlett Mid    C5    4         P      Oe 14.66
## 9      21 Bartlett Mid    C6    1        NP      Oe  7.39
## 10     22 Bartlett Mid    C6    2       Con      Oe 13.54
## 11     23 Bartlett Mid    C6    3         N      Oe 15.03
## 12     24 Bartlett Mid    C6    4         P      Oe 12.16
summer.2017.plotb<- ggplot(summer.2017b,aes(x= Treatment,y=T4,colour=Stand))+geom_point()+theme_bw()
summer.2017.plotb

#c7tohb (Oe horizon)
summer.2017c<-read.csv("c7tohb.csv",header=TRUE)
summer.2017c
##    ï..Sample        Forest    Age Stand Plot Treatment Horizon    T4
## 1         25      Bartlett Mature    C7    1         N      Oe  9.74
## 2         26      Bartlett Mature    C7    2        NP      Oe 12.39
## 3         27      Bartlett Mature    C7    3         P      Oe 15.29
## 4         28      Bartlett Mature    C7    4       Con      Oe  9.20
## 5         29      Bartlett Mature    C8    1         P      Oe  9.14
## 6         30      Bartlett Mature    C8    2       Con      Oe  4.45
## 7         31      Bartlett Mature    C8    3         N      Oe  9.53
## 8         32      Bartlett Mature    C8    4        NP      Oe  6.90
## 9         33      Bartlett Mature    C9    1       Con      Oe 11.13
## 10        34      Bartlett Mature    C9    2         P      Oe 12.12
## 11        35      Bartlett Mature    C9    3        NP      Oe  7.75
## 12        36      Bartlett Mature    C9    4         N      Oe 11.03
## 13        37 Hubbard Brook    Mid   HBM    1        NP      Oe  7.49
## 14        38 Hubbard Brook    Mid   HBM    2         N      Oe  7.50
## 15        39 Hubbard Brook    Mid   HBM    3       Con      Oe  9.00
## 16        40 Hubbard Brook    Mid   HBM    4         P      Oe  7.44
## 17        41 Hubbard Brook Mature   HBO    1         P      Oe  7.63
## 18        42 Hubbard Brook Mature   HBO    2         N      Oe  5.08
## 19        43 Hubbard Brook Mature   HBO    3        NP      Oe  6.58
## 20        44 Hubbard Brook Mature   HBO    4       Con      Oe 11.80
summer.2017.plotc<- ggplot(summer.2017c,aes(x= Treatment,y=T4,colour=Stand))+geom_point()+theme_bw()
summer.2017.plotc