R Markdown

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