This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
#This file includes additional calculations carried out for the Warton (Carnforth) and #Hodbarrow data, specifically to use in relation to discussion of Renyi analysis results.
#Should not need to load any packages for this data, seeing as such simple calculations are #being carried out here.
#Load the data…
FerrousGeochemData <- read.csv("CARNFORTH_AND_HODBARROW_CHEMISTRY_WITH_SITE_NAMES_MG_KG.csv", header = TRUE)
#Checking the data…
head(FerrousGeochemData)
## X pH_level SiO2 Al2O3 Fe2O3 CaO MgO Na2O
## 1 Carnforth_1 9.008 182299.6 47632.53 153873.611 164378.9 6030.359 2967.428
## 2 Carnforth_2 10.175 159862.7 52131.16 7903.508 262291.5 8924.931 2151.385
## 3 Carnforth_3 10.589 170613.7 44139.48 20283.340 204401.6 8140.985 2299.756
## 4 Carnforth_4 9.516 181832.2 35936.10 7833.566 200828.1 8080.681 2596.499
## 5 Carnforth_5 10.133 142567.7 43927.78 4826.036 250141.8 9769.182 1186.971
## 6 Carnforth_6 8.790 200997.0 31755.02 181850.631 121497.4 6030.359 3709.284
## K2O Cr2O3 TiO2 MnO P2O5 SrO BaO LOI Ag
## 1 10791.921 342.10121 3596.057 5421.202 999.8709 287.5021 1791.3027 22.86 0.00
## 2 3486.621 27.36810 4495.072 1858.698 399.9484 338.2378 1433.0422 13.90 0.15
## 3 5561.990 20.52607 2637.109 1936.144 599.9225 253.6784 716.5211 19.00 1.30
## 4 4648.828 41.05215 2457.306 1781.252 499.9354 253.6784 985.2165 20.70 0.30
## 5 2407.429 13.68405 3356.320 929.349 599.9225 253.6784 716.5211 22.20 0.30
## 6 13282.364 68.42024 1798.029 2323.372 999.8709 211.3986 806.0862 20.61 0.00
## As B Be Bi Cd Co Cu Ga Hg La Li Mo Ni Pb S Sb Sc Th Tl U V
## 1 0 0 0.0 0.0 0.00 0 0 0 0.0 0 0 0.0 0 0 10000 0.0 0 0 0 0 100
## 2 3 50 12.8 1.5 0.25 1 14 5 0.5 30 50 0.5 7 11 6000 2.0 11 10 5 10 10
## 3 11 60 8.4 2.0 2.70 3 53 10 0.5 30 40 0.5 8 212 4300 16.0 7 10 5 5 12
## 4 7 40 8.6 1.5 0.60 36 30 5 0.5 20 40 0.5 6 71 2700 2.0 5 10 5 5 11
## 5 3 40 11.4 1.5 0.25 4 17 5 0.5 30 40 0.5 3 28 8100 1.5 8 10 5 5 5
## 6 0 0 0.0 0.0 0.00 0 0 0 0.0 0 0 0.0 0 0 5000 0.0 0 0 0 0 100
## W Zn Akermanite Albite Aluminium_oxide_hydroxide Anhydrite Aragonite Augite
## 1 0 6000 1 0 0 0 0 0
## 2 5 36 0 0 0 0 0 0
## 3 5 503 0 0 0 0 0 0
## 4 5 120 0 0 0 0 0 0
## 5 5 46 1 0 0 0 0 0
## 6 0 2000 0 0 0 0 0 0
## Biotite Birnessite Calcite Clinochlore Cuspidine Diaspore Dickite Gehlenite
## 1 1 0 1 0 0 0 0 0
## 2 0 0 1 0 1 0 0 1
## 3 0 0 1 0 1 0 0 1
## 4 0 0 1 0 0 0 0 1
## 5 0 0 1 0 1 0 0 0
## 6 0 0 1 0 0 0 0 0
## Goethite Haematite Illite Kaolinite Langite Linnaeite Magnesioferrite
## 1 0 0 0 0 0 0 0
## 2 0 0 0 0 0 0 0
## 3 0 0 0 0 0 0 0
## 4 0 0 0 0 0 0 0
## 5 0 0 0 0 0 0 0
## 6 0 0 0 0 0 0 1
## Melilite Merwinite Microcline Mullite Muscovite Nitratine Orthoclase
## 1 0 0 0 0 0 0 0
## 2 0 0 0 0 0 1 0
## 3 0 0 0 0 0 1 0
## 4 0 0 0 0 0 0 0
## 5 0 0 0 0 1 1 0
## 6 0 1 0 0 0 0 0
## Orthopyroxene Periclase Phengite Pigeonite Pseudowollastonite Quartz Spinel
## 1 0 0 0 1 1 1 1
## 2 0 1 0 0 0 1 0
## 3 0 0 0 0 0 1 0
## 4 0 0 0 0 0 1 0
## 5 0 1 0 0 0 0 0
## 6 0 0 0 0 0 1 0
## Staurolite Valentinite
## 1 0 0
## 2 0 0
## 3 0 0
## 4 0 0
## 5 0 0
## 6 0 0
#Looking Cd data across the two sites…
min(FerrousGeochemData$Cd)
## [1] 0
mean(FerrousGeochemData$Cd)
## [1] 0.4633333
median(FerrousGeochemData$Cd)
## [1] 0.25
max(FerrousGeochemData$Cd)
## [1] 2.7
sd(FerrousGeochemData$Cd)
## [1] 0.6575025
#Looking B data across the two sites…
min(FerrousGeochemData$B)
## [1] 0
mean(FerrousGeochemData$B)
## [1] 31.33333
median(FerrousGeochemData$B)
## [1] 40
max(FerrousGeochemData$B)
## [1] 70
sd(FerrousGeochemData$B)
## [1] 23.25838
#Looking at Zn data across the two sites…
min(FerrousGeochemData$Zn)
## [1] 36
mean(FerrousGeochemData$Zn)
## [1] 625.6667
median(FerrousGeochemData$Zn)
## [1] 73
max(FerrousGeochemData$Zn)
## [1] 6000
sd(FerrousGeochemData$Zn)
## [1] 1568.947
#Looking at Mo data across the two sites…
min(FerrousGeochemData$Mo)
## [1] 0
mean(FerrousGeochemData$Mo)
## [1] 0.5
median(FerrousGeochemData$Mo)
## [1] 0.5
max(FerrousGeochemData$Mo)
## [1] 1
sd(FerrousGeochemData$Mo)
## [1] 0.2672612
#Looking at Hg data across the two sites…
min(FerrousGeochemData$Hg)
## [1] 0
mean(FerrousGeochemData$Hg)
## [1] 0.4666667
median(FerrousGeochemData$Hg)
## [1] 0.5
max(FerrousGeochemData$Hg)
## [1] 1
sd(FerrousGeochemData$Hg)
## [1] 0.2288689
#Looking at Cu data across the two sites…
min(FerrousGeochemData$Cu)
## [1] 0
mean(FerrousGeochemData$Cu)
## [1] 18.66667
median(FerrousGeochemData$Cu)
## [1] 18
max(FerrousGeochemData$Cu)
## [1] 53
sd(FerrousGeochemData$Cu)
## [1] 13.06395
#Looking at Pb data across the two sites…
min(FerrousGeochemData$Pb)
## [1] 0
mean(FerrousGeochemData$Pb)
## [1] 55.86667
median(FerrousGeochemData$Pb)
## [1] 48
max(FerrousGeochemData$Pb)
## [1] 212
sd(FerrousGeochemData$Pb)
## [1] 53.00252
#Looking at SrO data across the two sites…
min(FerrousGeochemData$SrO)
## [1] 42.27973
mean(FerrousGeochemData$SrO)
## [1] 216.4722
median(FerrousGeochemData$SrO)
## [1] 253.6784
max(FerrousGeochemData$SrO)
## [1] 338.2378
sd(FerrousGeochemData$SrO)
## [1] 89.70593
#Looking at Ag data across the two sites…
min(FerrousGeochemData$Ag)
## [1] 0
mean(FerrousGeochemData$Ag)
## [1] 0.3433333
median(FerrousGeochemData$Ag)
## [1] 0.3
max(FerrousGeochemData$Ag)
## [1] 1.3
sd(FerrousGeochemData$Ag)
## [1] 0.3348063
#Looking at As data across the two sites…
min(FerrousGeochemData$As)
## [1] 0
mean(FerrousGeochemData$As)
## [1] 7.066667
median(FerrousGeochemData$As)
## [1] 6
max(FerrousGeochemData$As)
## [1] 19
sd(FerrousGeochemData$As)
## [1] 5.338093
#Looking at Co data across the two sites…
min(FerrousGeochemData$Co)
## [1] 0
mean(FerrousGeochemData$Co)
## [1] 5.066667
median(FerrousGeochemData$Co)
## [1] 2
max(FerrousGeochemData$Co)
## [1] 36
sd(FerrousGeochemData$Co)
## [1] 8.947998