From 66cb2a3e2e704c7d00fd9156b710c2dbd539c5fa Mon Sep 17 00:00:00 2001
From: Dylan Beaudette <dylanbeaudette@users.noreply.github.com>
Date: Mon, 3 Feb 2025 09:56:39 -0800
Subject: [PATCH] notes

---
 R/accuracy_uncertainty.R                          | 2 ++
 misc/sandbox/soil-texture-separate-sim-by-class.R | 6 +++++-
 2 files changed, 7 insertions(+), 1 deletion(-)

diff --git a/R/accuracy_uncertainty.R b/R/accuracy_uncertainty.R
index 8627d0218..58855dddd 100644
--- a/R/accuracy_uncertainty.R
+++ b/R/accuracy_uncertainty.R
@@ -32,6 +32,8 @@
 ## TODO: test that sum(x) == 1
 shannonEntropy <- function(x, b = 2) {
   # 0s automatically removed by na.rm=TRUE (0 * log(0) = Nan)
+  # consider: ifelse(x == 0, 0, x * log(x))
+  
   res <- -1 * sum(x * log(x, base = b), na.rm = TRUE)
   return(res)
 }
diff --git a/misc/sandbox/soil-texture-separate-sim-by-class.R b/misc/sandbox/soil-texture-separate-sim-by-class.R
index 2c0ae5c76..c65791f99 100644
--- a/misc/sandbox/soil-texture-separate-sim-by-class.R
+++ b/misc/sandbox/soil-texture-separate-sim-by-class.R
@@ -18,7 +18,7 @@ data("soiltexture", package = 'aqp')
 ## get all soil texture measurements from latest KSSL snapshot
 
 # connect
-db <- dbConnect(RSQLite::SQLite(), 'E:/NASIS-KSSL-LDM/ncss_labdata.sqlite')
+db <- dbConnect(RSQLite::SQLite(), 'E:/NASIS-KSSL-LDM/ncss_labdata-2023.sqlite')
 
 qq <- "
 SELECT
@@ -124,8 +124,12 @@ prepareCompositionalSummary <- function(i, var = 'class') {
 }
 
 
+## TODO: use rnorm.acomp()
 
 
+##
+##
+##
 sampleComposition <- function(i, var = 'class', n = 100) {
   # simulate from Dirichlet
   s <- rDirichlet.acomp(n, i$D.alpha)