This project is a scala/Spark port of Arthur Chapman's reverse jackknifing approach to finding bioclimatic outliers within occurrence data. It is intended for internal usage within GBIF. It runs in 25 minutes on the current cluster setup.
Currently it uses bioclim data from 19 bioclimatic surfaces at a 0.1 degree resolution. Bioclim is only available for land, so no aquatic species could be run.
Run bioclim_extract.r to create bioclim_0.1_extract.tsv
Copy the bioclim_0.1_extract.tsv
onto the server and put into hdfs.
scp -r /cygdrive/c/Users/ftw712/Desktop/gbif_reverse_jackknife/data/bioclim_0.1_extract.tsv jwaller@c5gateway-vh.gbif.org:/home/jwaller/
hdfs dfs -put bioclim_0.1_extract.tsv bioclim_0.1_extract.tsv
hdfs dfs -ls
- Build project
sbt package
- Copy packaged jar onto cluster
scp -r /cygdrive/c/Users/ftw712/Desktop/gbif_reverse_jackknife/target/scala-2.11/gbif_reverse_jackknife_2.11-0.1.jar jwaller@c5gateway-vh.gbif.org:/home/jwaller/
- Run with spark2-submit
spark2-submit --num-executors 40 --executor-cores 5 --driver-memory 8g --driver-cores 4 --executor-memory 16g gbif_reverse_jackknife_2.11-0.1.jar
Output is a parquet file in hdfs called rjack_outliers_export.parquet
Here I plot an example of the results of a primate species with one outlier in Oregon State.
The point(s) occur(s) at the Oregon Regional Primate Research Center.
Cumulative frequency plots are to the right. Bioclimatic surface values are on the x-axis and frequency on the y-axis. Black-colored points very much to the right or left or right are outliers. I also write the number of other surfaces that flagged this point as an outlier as a little number next to the point.
Because of multiple-comparisons, a point should not be considered a "true outlier" unless it an outlier in > 5 surfaces.
Around 2M occurrences get flagged with at least 1 surface (of 19) as an outlier. Around 50K occurrences have >5 five surfaces as outliers.
The current implementation runs on all terrestial Fungi, Animals, Plants in around 25 minutes.
- Only includes occurrences with hasgeospatialissues = false
- Excludes basis record = FOSSIL_SPECIMEN, UNKNOWN, LIVING_SPECIMEN
- Only run on Kingdoms = Fungi, Animals, Plants
- decimallatitude & decimallongitude rounded to nearest 0.1 degrees
- Only unique bioclimatic values for each surface used
- Only specieskeys with > 10 unique bioclimatic values for each surface were run
This implementation is a translation from R to scala of biogeo::rjack
source. It is also similar to what is available in DivaGIS.
rjack <-
function (d)
{
xx <- d
d <- unique(d)
rng <- diff(range(d))
mx <- mean(d)
n <- length(d)
n1 <- n - 1
t1 <- (0.95 * sqrt(n)) + 0.2
x <- sort(d)
y <- rep(0, n1)
for (i in 1:n1) {
x1 <- x[i + 1]
if (x[i] < mx) {
y[i] <- (x1 - x[i]) * (mx - x[i])
}
else {
y[i] <- (x1 - x[i]) * (x1 - mx)
}
}
my <- mean(y)
z <- y/(sqrt(sum((y - my)^2)/n1))
out <- rep(0, length(xx))
if (any(z > t1)) {
f <- which(z > t1)
v <- x[f]
if (v < median(x)) {
xa <- (xx <= v) * 1
out <- out + xa
}
if (v > median(x)) {
xb <- (xx >= v) * 1
out <- out + xb
}
}
else {
out <- out
}
f <- which(out == 1)
}
Around 6M outliers are flagged, so it is possible to export locally.
In the end, a file is created in hdfs called rjack_outliers_export.parquet
.
rm -r /mnt/auto/misc/download.gbif.org/custom_download/jwaller/rjack_outliers_export.parquet.zip
hadoop fs -get /user/jwaller/rjack_outliers_export.parquet /mnt/auto/misc/download.gbif.org/custom_download/jwaller/rjack_outliers_export.parquet
cd /mnt/auto/misc/download.gbif.org/custom_download/jwaller/
zip -r rjack_outliers_export.parquet.zip rjack_outliers_export.parquet
rm -r /mnt/auto/misc/download.gbif.org/custom_download/jwaller/rjack_outliers_export.parquet
Read multi-part parquet into R.
install.packages("arrow")
devtools::install_github("jhnwllr/parqr")
path = "C:/Users/ftw712/Desktop/"
parqr::parquet_readr(paste0(path,"rjack_outliers_export.parquet"))
Chapman, A.D. (2005) Principles and Methods of Data Cleaning - Primary Species and Species- Occurrence Data, version 1.0. Report for the Global Biodiversity Information Facility, Copenhagen.