From 5f8dedd0af383ab14582f538443bf8ddd207df33 Mon Sep 17 00:00:00 2001
From: Munyoki Kilyungi
Date: Fri, 6 Feb 2026 06:10:12 +0300
Subject: Add marker count.

Signed-off-by: Munyoki Kilyungi <me@bonfacemunyoki.com>
---
 examples/genotype-datasets.scm | 16 +++++++++++++++-
 1 file changed, 15 insertions(+), 1 deletion(-)

(limited to 'examples/genotype-datasets.scm')

diff --git a/examples/genotype-datasets.scm b/examples/genotype-datasets.scm
index ba31088..ebe2349 100755
--- a/examples/genotype-datasets.scm
+++ b/examples/genotype-datasets.scm
@@ -36,6 +36,19 @@
     (set gnt:has_strain (string->identifier "set" (field InbredSet Name InbredSetName) #:separator "_"))
     (set dct:created (annotate-field (field GenoFreeze CreateTime) '^^xsd:datetime))))
 
+(define-transformer gn:dataset->marker/snp-count
+  (tables (GenoFreeze
+           (inner-join InbredSet "ON InbredSet.Id = GenoFreeze.InbredSetId")
+           (inner-join Species "ON InbredSet.SpeciesId = Species.Id")
+           (inner-join Geno "ON Geno.SpeciesId = Species.Id"))
+          "WHERE GenoFreeze.public > 0 AND Species.Name != 'monkey' GROUP BY GenoFreeze.Name")
+  (triples (string->identifier "dataset" (field GenoFreeze Name) #:separator "_")
+    (set gnt:has_marker_count
+         (string->symbol
+          (format #f "'~s'^^xsd:integer"
+                  (field
+                   ("COUNT(DISTINCT Geno.Marker_Name)" MarkerCount)))))))
+
 
 (let* ((option-spec
         '((settings (single-char #\s) (value #t))
@@ -67,7 +80,8 @@
    (inputs
     (list
      gn:set->gn:dataset
-     gn:dataset->metadata))
+     gn:dataset->metadata
+     gn:dataset->marker/snp-count))
    (outputs
     `(#:documentation ,documentation
       #:rdf ,output))))
-- 
cgit 1.4.1