Some Nutrient Plots by Substrate Texture
Ezra Citron
Monday, May 18, 2015
Read in the data
#elwha = read.csv("ElwhaMicronutrientStudy_Data_Citron_150518.csv", header=TRUE)
elwha = read.csv(file.choose(), header=TRUE)
head(elwha)## sample lat long site plant subtype subwater
## 1 141106-Aldwell-01-01- 48.07300 -123.5803 1 1 course 4.558490566
## 2 141106-Aldwell-01-02- 48.07300 -123.5803 1 2 course 3.01734104
## 3 141106-Aldwell-01-03- 48.07300 -123.5803 1 3 course 2.9376
## 4 141106-Aldwell-01-04- 48.07300 -123.5803 1 4 course 3.23723229
## 5 141107-Aldwell-02-01- 48.07737 -123.5766 2 1 fine 34.141527
## 6 141107-Aldwell-02-02- 48.07737 -123.5766 2 2 fine 36.19723866
## ht diam denscwd denstot N C Na Mg
## 1 1.97 18.0 1.909859317 1.973521294 1.456195 46.67445 28.70920 2941.645
## 2 1.54 16.0 0.671987537 0.742723068 1.142410 47.67860 30.86497 2366.566
## 3 1.67 22.5 0.101859164 0.101859164 1.383349 46.11044 27.56143 2756.648
## 4 2.13 25.0 0.381971863 0.458366236 1.491973 47.10913 22.43267 2683.669
## 5 3.79 13.0 N/A 15.59718442 2.082052 47.22429 24.90818 3376.432
## 6 3.16 11.0 N/A 12.73239545 1.570839 47.89363 29.98803 3177.777
## Al K Ca V Cr Mn Fe
## 1 82.33288 6009.130 13496.154 0.15305722 1.3676623 179.7861 92.52694
## 2 50.86306 3966.413 15253.888 0.07682285 0.9996529 144.7912 49.44936
## 3 37.84200 3418.769 19899.017 0.06611103 0.8835022 183.4733 53.52568
## 4 51.79446 3131.030 9662.997 0.10174977 0.8127058 184.5613 68.29711
## 5 37.44854 12540.366 8646.277 0.04920077 0.6216552 499.0030 45.01753
## 6 94.85604 6353.565 9565.148 0.06388365 1.0872650 499.1377 56.37670
## Co Ni Cu Zn Mo Ag Cd
## 1 1.374003 58.96946 4.684390 172.2747 0.25558312 11.315894 0.6594435
## 2 1.144042 40.49927 4.508751 234.2418 0.18081291 4.739673 0.6459760
## 3 1.779135 39.47961 4.109195 201.8206 0.12964434 3.030595 1.0391765
## 4 6.516558 41.50019 4.451739 454.0714 0.22169222 3.771253 1.0438423
## 5 19.489099 43.29201 5.610286 166.3499 0.46458580 24.567737 0.6561657
## 6 32.130744 70.74431 6.547453 185.6355 0.08685194 5.036120 0.9172400
## Ba
## 1 5.140928
## 2 5.931509
## 3 7.434964
## 4 7.127952
## 5 2.620932
## 6 5.344352str(elwha) #take a look at variables in data.frame## 'data.frame': 30 obs. of 30 variables:
## $ sample : Factor w/ 30 levels "141106-Aldwell-01-01- ",..: 1 2 3 4 9 10 11 12 13 14 ...
## $ lat : num 48.1 48.1 48.1 48.1 48.1 ...
## $ long : num -124 -124 -124 -124 -124 ...
## $ site : int 1 1 1 1 2 2 2 2 3 3 ...
## $ plant : int 1 2 3 4 1 2 3 4 1 2 ...
## $ subtype : Factor w/ 4 levels "course","fine",..: 1 1 1 1 2 2 2 2 4 4 ...
## $ subwater: Factor w/ 29 levels "14.96494157",..: 24 5 2 6 16 20 17 15 1 12 ...
## $ ht : num 1.97 1.54 1.67 2.13 3.79 3.16 3.93 4.3 5.05 4.54 ...
## $ diam : num 18 16 22.5 25 13 11 22.5 22.5 26.5 22.5 ...
## $ denscwd : Factor w/ 13 levels "0.101859164",..: 8 7 1 4 13 13 13 13 13 13 ...
## $ denstot : Factor w/ 25 levels "0.101859164",..: 6 3 1 2 16 10 8 11 19 16 ...
## $ N : num 1.46 1.14 1.38 1.49 2.08 ...
## $ C : num 46.7 47.7 46.1 47.1 47.2 ...
## $ Na : num 28.7 30.9 27.6 22.4 24.9 ...
## $ Mg : num 2942 2367 2757 2684 3376 ...
## $ Al : num 82.3 50.9 37.8 51.8 37.4 ...
## $ K : num 6009 3966 3419 3131 12540 ...
## $ Ca : num 13496 15254 19899 9663 8646 ...
## $ V : num 0.1531 0.0768 0.0661 0.1017 0.0492 ...
## $ Cr : num 1.368 1 0.884 0.813 0.622 ...
## $ Mn : num 180 145 183 185 499 ...
## $ Fe : num 92.5 49.4 53.5 68.3 45 ...
## $ Co : num 1.37 1.14 1.78 6.52 19.49 ...
## $ Ni : num 59 40.5 39.5 41.5 43.3 ...
## $ Cu : num 4.68 4.51 4.11 4.45 5.61 ...
## $ Zn : num 172 234 202 454 166 ...
## $ Mo : num 0.256 0.181 0.13 0.222 0.465 ...
## $ Ag : num 11.32 4.74 3.03 3.77 24.57 ...
## $ Cd : num 0.659 0.646 1.039 1.044 0.656 ...
## $ Ba : num 5.14 5.93 7.43 7.13 2.62 ...Some plots
attach(elwha)
plot.default(C,N, ylab="leaf %N", xlab="leaf %C")
CBox <- boxplot(C~subtype, ylab = "Leaf %C", xlab="Substrate Texture", plot=0)
boxplot(C~subtype, ylab = "Leaf %C", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(N~subtype, ylab = "Leaf %N", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(C~subtype, ylab = "Leaf %C", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Na~subtype, ylab = "Leaf Na (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Mg~subtype, ylab = "Leaf Mg (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Al~subtype, ylab = "Leaf Al (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(K~subtype, ylab = "Leaf K (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Ca~subtype, ylab = "Leaf Ca (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(V~subtype, ylab = "Leaf V (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Cr~subtype, ylab = "Leaf Cr (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
MnBox <- boxplot(Mn~subtype, ylab = "Leaf Fe (ppm)", xlab="Substrate Texture", plot=0)
boxplot(Mn~subtype, ylab = "Leaf Fe (ppm)", xlab="Substrate Texture", names=paste(MnBox$names, "(n=", MnBox$n, ")"))
boxplot(Co~subtype, ylab = "Leaf Co (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Ni~subtype, ylab = "Leaf Mg (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Cu~subtype, ylab = "Leaf Cu (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Zn~subtype, ylab = "Leaf Zn (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Cd~subtype, ylab = "Leaf Cd (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))
boxplot(Ba~subtype, elwha, ylab = "Leaf Ba (ppm)", xlab="Substrate Texture", names=paste(CBox$names, "(n=", CBox$n, ")"))