vignettes/directional_stats_vig.R

In maxgav13/GMisc: Miscellaneous Functions for Geology (and its areas)

## ----setup, include=FALSE-----------------------------------------------------
library(GMisc)

knitr::opts_chunk$set(
  echo = TRUE,
  message = FALSE,
  warning = FALSE,
  error = FALSE,
  fig.align = "center",
  out.width = "80%"
)

## -----------------------------------------------------------------------------
x <- c(255, 239, 222, 231, 199, 271, 222, 274, 228, 
       246, 177, 199, 257, 201, 237, 209, 216)
dir_stats_2D(x, dir = 1)

## -----------------------------------------------------------------------------
dir_stats_2D(carolina, dir = 0)

## -----------------------------------------------------------------------------
rose_diag_2D(x, width = 30, dir = 1)

## -----------------------------------------------------------------------------
rose_diag_2D(carolina, width = 10, dir = 0)

## -----------------------------------------------------------------------------
dir_unif_test_2D(x, dir = 1)

## -----------------------------------------------------------------------------
y = c(225, 208, 172, 198, 204, 183, 190, 212, 247, 127, 167, 234, 217, 192, 212, 171, 169, 210, 245, 222, 185, 227, 193, 178, 187, 182, 194, 217, 168, 211, 234, 204, 221, 198, 261, 228, 146, 201, 146, 231)
dir_2sample_test_2D(x, y, dir = 1, conf.level = 0.95)

## -----------------------------------------------------------------------------
dir = c(12,18,22,15,10,20)
dip = c(42,40,48,30,42,30)
dir_stats_3D(dir, dip, conf.level = 0.95)

## -----------------------------------------------------------------------------
dip.dir = c(12,18,22,15,10,20)
strike = d2s(dip.dir)
dip = c(42,40,48,30,42,30)
rose_diag_3D(az = dip.dir, dip = dip, width = 30, dir = 1, conf.level = 0.95)
rose_diag_3D(az = strike, dip = dip, width = 30, dir = 0, conf.level = 0.95)

## -----------------------------------------------------------------------------
dir_unif_test_3D(dir, dip)

## -----------------------------------------------------------------------------
d1 = c(12,18,22,15,10,20)
i1 = c(42,40,48,30,42,30)
d2 = c(111.7,109.2,185.1,83.1,235.2,226.7)
i2 = c(40.4,31.4,25.9,26.2,27,31.5)
dir_2sample_test_3D(d1, i1, d2, i2)

## -----------------------------------------------------------------------------
set.seed(4101)
dir <- runif(min = 30, max = 80, n = 20)
dip <- runif(min = 10, max = 60, n = 20)
stereo_lines(dir, dip, points.shape = 3, points.col = "blue", add = FALSE)

## -----------------------------------------------------------------------------
set.seed(4101)
strike <- runif(min = 30, max = 80, n = 20)
dip.dir <- s2d(strike)
dip <- runif(min = 10, max = 60, n = 20)
strike2 <- runif(min = 240, max = 290, n = 20)
dip2 <- runif(min = 10, max = 60, n = 20)
stereo_plane(strike, dip)
stereo_plane(dip.dir, dip, dir = "dip")
stereo_plane(strike2, dip2, col = "red", add = TRUE)

## -----------------------------------------------------------------------------
set.seed(4101)
dir <- runif(min = 30, max = 80, n = 20)
dip <- runif(min = 10, max = 60, n = 20)
stereo_mean(dir, dip, conf.level = 0.95, mean.col = "red", add = FALSE)

## -----------------------------------------------------------------------------
set.seed(4101)
dip.dir <- runif(min = 120, max = 170, n = 20)
dip <- runif(min = 10, max = 60, n = 20)
stereo_pole(dip.dir, dip)
stereo_plane(dip.dir, dip, dir = "dip")
stereo_pole(dip.dir, dip, add = TRUE)

## -----------------------------------------------------------------------------
set.seed(4101)
dir <- runif(min = 30, max = 80, n = 20)
dip <- runif(min = 10, max = 60, n = 20)
stereo_mean(dir, dip)
stereo_ps(dir, dip)