From e12150797ca75f1c35d55e16c73848c30937ff77 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Tue, 21 May 2024 07:49:27 -0500
Subject: Writing json-ld

---
 scripts/precompute/list-traits-to-compute.scm | 28 +++++++++++++++++----------
 1 file changed, 18 insertions(+), 10 deletions(-)

(limited to 'scripts/precompute/list-traits-to-compute.scm')

diff --git a/scripts/precompute/list-traits-to-compute.scm b/scripts/precompute/list-traits-to-compute.scm
index 7e0660d..2c48d83 100755
--- a/scripts/precompute/list-traits-to-compute.scm
+++ b/scripts/precompute/list-traits-to-compute.scm
@@ -38,7 +38,7 @@ you should see
 
 Now list the next 1000 trait IDs:
 
-    . .guix-shell -- guile -L . -e main -s ./scripts/precompute/list-traits-to-compute.scm --first 0 --next 1000
+    . .guix-shell -- guile -L . -e main -s ./scripts/precompute/list-traits-to-compute.scm --next 1000
 
 The current logic is to list all datasets that contain a
 BXD. (bxd-strain-id-names #:used-for-mapping? #t) fetches all ids and
@@ -71,27 +71,35 @@ When that is the case we might as well write the phenotype file because we have
              (gn util convert)
              (gn runner gemma)
              ; (rnrs base)
+             (hashing sha-2)
              (ice-9 getopt-long)
              (ice-9 match)
              (ice-9 textual-ports)
              (json)
+             (rnrs bytevectors)
              (srfi srfi-1)
              (srfi srfi-19) ; time
              )
 
 
-(define (write-json id recs)
+(define (write-json-ld id recs)
+  ;; see also https://www.w3.org/2018/jsonld-cg-reports/json-ld/
   (display id)
   (newline)
   (let* [(traits (map (lambda (r)
-                                                (match r
-                                                  [(strain-id . value) (cons (bxd-name strain-id) value)]
-                                                  ))
+                        (match r
+                          [(strain-id . value) (cons (bxd-name strain-id) value)]
+                          ))
                       (reverse recs)))
-         (json-data `((meta .
-                            ((id . ,id)
-                             (time . ,(date->string (time-utc->date (current-time))))))
-                      (traits . ,traits)))]
+         (sha256 (sha-256->string (sha-256 (string->utf8 (scm->json-string traits)))))
+         (json-data `(("@context" . "https://genenetwork.org/resource")
+                      (type . traits)
+                      (meta . ((version . "0.1")
+                               (steps . ())
+                               (sha256 . ((input-traits . ,sha256)))
+                               (time . ,(date->string (time-utc->date (current-time))))))
+                      (traits .
+                      ((,id . ,traits)))))]
     (call-with-output-file (string-append (number->string id) ".json")
       (lambda (port)
         (put-string port (scm->json-string json-data))))
@@ -135,7 +143,7 @@ When that is the case we might as well write the phenotype file because we have
              (for-each (lambda (r)
                          (match r
                            ((id . recs) (if (has-bxd? recs)
-                                              (write-json id recs)
+                                              (write-json-ld id recs)
                                             ))
                             )) nrecs)
              (if (> rest 0)
-- 
cgit v1.2.3