## Data 
library(AnalystHelper)
library(reshape2)
library(plyr)
library(dataRetrieval)

## Spatial
library(sf)
library(terra)
library(EVERSpatDat)

library(downloadthis)
library(rvest)


data(LOK)

# renv::snapshot(update = TRUE,force = TRUE)

Updated: 2025-07-14 09:43 EDT

Including Plots

# set.seed(123)
n.val <- 500
dat <- rnorm(n.val)


plot(dat,ann=F,type="n",axes=F,ylim=c(-3,3))
points(dat,pch=21,bg="pink",lwd=0.1,col="red",cex=2)
abline(h=mean(dat))
abline(h=c(mean(dat)+sd(dat),mean(dat)-sd(dat)),lty=2,col="blue")
axis_fun(1,seq(0,n.val,n.val/5),seq(0,n.val,n.val/50),seq(0,n.val,n.val/5))
axis_fun(2,seq(-3,3,1),seq(-3,3,0.5),seq(-3,3,1))
box(lwd=1)

summary(dat)
##     Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
## -3.17566 -0.57909  0.08748  0.06533  0.74674  2.82984
samp.site <- st_sample(LOK,10)

plot(st_geometry(LOK),col="lightblue",bg="cornsilk",border="grey")
plot(st_geometry(samp.site),add=T,pch=21,bg="lightgreen",col="red",cex=1.5,lwd=0.5)
box(lwd=1)