O_NIVEL_Intervalos

install.packages("classInt")

## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.5'
## (as 'lib' is unspecified)

library(classInt)
install.packages("BAMMtools")

## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.5'
## (as 'lib' is unspecified)

library(BAMMtools)

## Loading required package: ape

Cargar datos

setwd("/cloud/project")
sondeo<-read.csv("ONivel_collar.txt",header=T, sep=",",dec=".")
muestra<-read.csv("ONivel_muestra.txt",header=T, sep=",",dec=".")
compCristina<-read.csv("O_NIVEL_Composito_Cristina.txt",header=T, sep = "",dec=".")
compJane<-read.csv("O_NIVEL_Composito_Jane.txt",header=T, sep= "",dec=".")
compX<-read.csv("O_NIVEL_Composito_X.txt",header=T, sep= "",dec=".")

MUESTRAS Extraer la variable de interes

Oro <- muestra$AU_PPM
Oro <- na.omit(Oro)
CajaOro <- boxplot(Oro, horizontal =  T)

OroOut <- CajaOro$out

Intervalos iguales

IntIgual <- classIntervals(OroOut, n = 5, style = "equal")
breaks_igual <- IntIgual$brks
cat("Breaks - Intervalos iguales:\n", breaks_igual, "\n\n")

## Breaks - Intervalos iguales:
##  6.89 72.952 139.014 205.076 271.138 337.2

Percentiles

IntPerc <- classIntervals(OroOut, n = 5, style = "quantile")
breaks_perc <- IntPerc$brks
cat("Breaks - Percentiles:\n", breaks_perc, "\n\n")

## Breaks - Percentiles:
##  6.89 8.48 10.36 20.162 33.4 337.2

Jenks

IntJenks <- classIntervals(OroOut, n = 5, style = "jenks")
breaks_jenks <- IntJenks$brks
cat("Breaks - Jenks:\n", breaks_jenks, "\n\n")

## Breaks - Jenks:
##  6.89 17.8 32.7 66.32 123.4 337.2

VETA CRISTINA Extraer la variable de interes

OroCristina <- compCristina$Au.ppm.
OroCristina <- na.omit(OroCristina)
CajaOroCristina <- boxplot(OroCristina, horizontal =  T)

OroOutCristina <- CajaOroCristina$out

Intervalos iguales

IntIgualCrist <- classIntervals(OroOutCristina, n = 5, style = "equal")
breaks_igualCrist <- IntIgualCrist$brks
cat("Breaks - Intervalos iguales:\n", breaks_igualCrist, "\n\n")

## Breaks - Intervalos iguales:
##  5.78 10.266 14.752 19.238 23.724 28.21

Percentiles

IntPercCrist <- classIntervals(OroOutCristina, n = 5, style = "quantile")
breaks_percCrist <- IntPercCrist$brks
cat("Breaks - Percentiles:\n", breaks_percCrist, "\n\n")

## Breaks - Percentiles:
##  5.78 6.144 16.386 23.628 27.97 28.21

Jenks

IntJenksCrist <- classIntervals(OroOutCristina, n = 5, style = "jenks")
breaks_jenksCrist <- IntJenksCrist$brks
cat("Breaks - Jenks:\n", breaks_jenksCrist, "\n\n")

## Breaks - Jenks:
##  5.78 6.18 12.42 19.5 25.15 28.21

VETA JANE Extraer la variable de interes

OroJane <- compJane$Au.ppm.
OroJane <- na.omit(OroJane)
CajaOroJane <- boxplot(OroJane, horizontal =  T)

OroOutJane <- CajaOroJane$out

Intervalos iguales

IntIgualJane <- classIntervals(OroOutJane, n = 5, style = "equal")
breaks_igualJane <- IntIgualJane$brks
cat("Breaks - Intervalos iguales:\n", breaks_igualJane, "\n\n")

## Breaks - Intervalos iguales:
##  23.81 86.488 149.166 211.844 274.522 337.2

Percentiles

IntPercJane <- classIntervals(OroOutJane, n = 5, style = "quantile")
breaks_percJane <- IntPercJane$brks
cat("Breaks - Percentiles:\n", breaks_percJane, "\n\n")

## Breaks - Percentiles:
##  23.81 29.48 34.234 40.11 46.72 337.2

Jenks

IntJenksJane <- classIntervals(OroOutJane, n = 5, style = "jenks")
breaks_jenksJane <- IntJenksJane$brks
cat("Breaks - Jenks:\n", breaks_jenksJane, "\n\n")

## Breaks - Jenks:
##  23.81 40.5 66.32 123.4 260.32 337.2

VETA X Extraer la variable de interes

OroX <- compX$Au.ppm.
OroX <- na.omit(OroX)
CajaOroX <- boxplot(OroX, horizontal =  T)

Intervalos iguales

IntIgualX <- classIntervals(OroX, n = 5, style = "equal")
breaks_igualX <- IntIgualX$brks
cat("Breaks - Intervalos iguales:\n", breaks_igualX, "\n\n")

## Breaks - Intervalos iguales:
##  0.95 1.402 1.854 2.306 2.758 3.21

Percentiles

IntPercX <- classIntervals(OroX, n = 5, style = "quantile")
breaks_percX <- IntPercX$brks
cat("Breaks - Percentiles:\n", breaks_percX, "\n\n")

## Breaks - Percentiles:
##  0.95 1.6 1.908 2.586 3.01 3.21

Jenks

IntJenksX <- classIntervals(OroX, n = 5, style = "jenks")
breaks_jenksX <- IntJenksX$brks
cat("Breaks - Jenks:\n", breaks_jenksX, "\n\n")

## Breaks - Jenks:
##  0.95 1.22 1.86 2.41 2.81 3.21