r11162: support for map/smap and ucon ancestors
authorKevin M. Rosenberg <kevin@rosenberg.net>
Sat, 16 Sep 2006 13:58:04 +0000 (13:58 +0000)
committerKevin M. Rosenberg <kevin@rosenberg.net>
Sat, 16 Sep 2006 13:58:04 +0000 (13:58 +0000)
class-support.lisp
classes.lisp
sql-classes.lisp

index 92b39aab617f3a7549c0ed5f587cc1eff0078377..7948dec091732de7c0337687ef34d0f2ca2aa400 100644 (file)
          "Other"))))
 
 
-(defun ucon-parents (con &key sab include-rb)
-  (ucon-ancestors con :sab sab :include-rb include-rb 
-                  :single-level t))
-
-(defun ucon-ancestors (ucon &key sab single-level include-rb)
+(defun ucon-parents (ucon &key sab include-rb)
+  (ucon-ancestors ucon :sab sab :include-rb include-rb
+                  :ancestors nil))
+
+(defun is-ucon-in-ancestors (ucon ancestors)
+  (cond
+   ((null ancestors) nil)
+   ((atom ancestors) (eql (cui ucon) (cui ancestors)))
+   ((listp (car ancestors))
+    (or (is-ucon-in-ancestors ucon (car ancestors))
+        (is-ucon-in-ancestors ucon (cdr ancestors))))
+   (t
+    (or (eql (cui ucon) (cui (car ancestors)))
+        (is-ucon-in-ancestors ucon (cdr ancestors))))))
+
+
+(defun ucon-ancestors (ucon &key sab include-rb ancestors)
   "Returns a list of ancestor lists for a concept"
   (let* ((parent-rels (append (filter-urels-by-rel (s#rel ucon) "par")
-                              (when include-rb 
+                              (when include-rb
                                 (filter-urels-by-rel (s#rel ucon) "rb"))))
-        (anc nil))
+         (parents nil))
     (when sab
       (setq parent-rels (delete-if-not
                         (lambda (rel) (string-equal sab (sab rel)))
                         parent-rels)))
-    (dolist (rel parent-rels (nreverse anc))
+    (dolist (rel parent-rels)
       (let ((parent (find-ucon-cui (cui2 rel))))
-       (push
-        (if single-level
-            (list parent)
-          (list* parent (car (ucon-ancestors parent (sab rel) nil))))
-        anc)))))
+        ;; (format t "~S ~S ~S ~S~%" rel ucon parent ancestors)
+        (unless (is-ucon-in-ancestors parent ancestors)
+          (push
+           (list*
+            parent
+            (ucon-ancestors parent :sab (sab rel) :ancestors (append (list parent) ancestors)))
+           parents))))
+    (nreverse parents)))
+
+
+(defun find-minimum-distance (anc1 anc2 &key (distance 0) (minimum 0) path)
+  (cond
+    ((or (null anc1) (null anc2) nil)
+    ((and (atom (car anc1)) (atom (car anc2)))
+     (when (eql (cui ucon1) (cui ucon2))
+       (return-from find-minimum-distance (values distance path)))
+     (incf distance)
+     (when (> distance minimum)
+       (return-from find-minimum-distance nil))
+     (multiple-value-bind (dist1 path1)
+         (find-minimum-distance anc1 (cdr anc2)
+                                :distance distance :minimum minimum
+                                :path (list* anc1 path))
+       (when (and dist1 (< dist1 minimum))
+         (setf minimum dist1
+               path path1)))
+     (multiple-value-bind (dist2 path2)
+         (find-minimum-distance (cdr anc1) anc2
+                                :distance distance :minimum minimum
+                                :path (list* anc2 path))
+       (when (and dist2 (< dist2 minimum))
+         (setf minimum dist2
+               path path2)))
+       (return-from find-minimum-distance (values distance path))
+
+
+           (min2 (find-minimum-distance
+                  (cdr anc1) anc2
+                  :distance distance :minimum minimum
+                  :path (list* anc2 path))))
+
+       (when (and min2 (< min2 minimum))
+         (setf minimum min2
+               path (list* anc2 path)))))
+
+
+  )
+
+(defun ucon-cdist (ucon1 ucon2 &key include-rb sab)
+  "Compute James Cimino's CDist metric"
+  (let ((anc1 (ucon-ancestors ucon1 :include-rb include-rb :sab sab))
+        (anc2 (ucon-ancestors ucon2 :include-rb include-rb :sab sab)))
+    (find-minimum-distance (list ucon1 anc1) (list ucon2 anc2))))
+
 
 (defgeneric cxt-ancestors (obj))
 (defmethod cxt-ancestors ((con ucon))
 (defun ucon-has-sab (ucon sab)
   (and (find-if (lambda (uso) (string-equal sab (sab uso))) (s#so ucon)) t))
 
-    
+
 #+scl
 (dolist (c '(urank udef usat uso ucxt ustr ulo uterm usty urel ucoc uatx ucon uxw uxnw uxns lexterm labr lagr lcmp lmod lnom lprn lprp lspl ltrm ltyp lwd sdef sstr sstre1 sstre2 usrl))
     (let ((cl (find-class c)))
index cee2fcd3b8db1569c77bf39a97d5713f336a58c8..cdc4d561648b62b234ad293f6ba58b6b446086f2 100644 (file)
@@ -32,7 +32,7 @@
   (:default-print-slots sab srl)
   (:description "Custom Table: Source Restriction Level"))
 
-  
+
 (defclass urank (umlsclass)
   ((rank :value-type fixnum :initarg :rank :reader rank)
    (sab :value-type string :initarg :sab :reader sab)
   (:user-name "String")
   (:default-print-slots sui stt lrl str suppress))
 
-       
+
 (defclass uso (umlsclass)
   ((aui :value-type fixnum :initarg :aui :reader aui :print-formatter fmt-aui
         :hyperlink find-ucon-aui)
   (:metaclass hyperobject-class)
   (:user-name "Relationship")
   (:default-print-slots stype1 rel cui2 aui2 stype2 rela rui srui sab sl rg dir suppress pfstr2))
-       
+
 (defclass ucoc (umlsclass)
   ((cui1 :value-type fixnum :initarg :cui1 :reader cui1 :print-formatter fmt-cui)
    (aui1 :value-type fixnum :initarg :aui1 :reader aui1 :print-formatter fmt-aui)
    (pfstr :value-type cdata :initarg :pfstr :reader pfstr)
    (s#def :reader s#def :subobject (find-udef-cui cui))
    (s#so :reader s#so :subobject (find-uso-cui cui))
+   (s#hier :reader s#hier :subobject (find-uhier-cui cui))
+   (s#map :reader s#map :subobject (find-umap-cui cui))
+   (s#smap :reader s#smap :subobject (find-usmap-cui cui))
    (s#sty :reader s#sty :subobject (find-usty-cui cui))
    (s#lo :reader s#lo :subobject (find-ulo-cui cui))
    (s#term :reader s#term :subobject (find-uterm-cui cui))
index 1187d5931bd9bc52a0791c43752dbaa434b2e865..594103fcd95c0d9a7b3303f5119b1f3512854136 100644 (file)
@@ -229,19 +229,13 @@ is OBJNAME from TABLE where WHERE-NAME field = WHERE-VALUE with FIELDS"
   (ensure-cui-integer cui)
   (unless cui (return-from find-ucon-cui nil))
 
-  (let ((tuple (car (mutex-sql-query
-                     (format nil
-                             "SELECT kcuilrl,str FROM MRCONSO WHERE CUI=~D AND KPFENG=1~A"
-                             cui (if srl (format nil " AND SRL<=~D" srl) ""))))))
-    (unless tuple
-      (setq tuple (car (mutex-sql-query
-                        (format nil
-                                "SELECT kcuilrl,str FROM MRCONSO WHERE CUI=~D"
-                                cui (if srl (format nil " AND SRL<=~D" srl) nil))))))
-    (unless tuple
-      (return-from find-ucon-cui nil))
-    (make-instance 'ucon :cui cui :pfstr (second tuple)
-                   :lrl (ensure-integer (first tuple)))))
+  (or
+   (collect-umlisp-query (mrconso (kcuilrl str) srl cui cui :single t :filter "KPFENG=1")
+    (make-instance 'ucon :cui cui :pfstr str
+                   :lrl kcuilrl))
+   (collect-umlisp-query (mrconso (kcuilrl str) srl cui cui :single t)
+    (make-instance 'ucon :cui cui :pfstr str
+                   :lrl kcuilrl))))
 
 (defun find-ucon-cui-sans-pfstr (cui &key (srl *current-srl*))
   "Find ucon for a cui"
@@ -572,11 +566,11 @@ is OBJNAME from TABLE where WHERE-NAME field = WHERE-VALUE with FIELDS"
 
 (defun find-uso-aui (aui &key (srl *current-srl*))
   (ensure-sui-integer aui)
-  (collect-umlisp-query (mrconso (sab cui sui code srl tty saui sdui scui lat str) srl aui
-                                aui :lrl srl :single t)
+  (collect-umlisp-query (mrconso (sab cui sui code srl tty saui sdui scui lat
+                                      str) srl aui aui :lrl srl :single t)
     (make-instance 'uso :aui aui :cui cui :sab sab :code code :srl srl :tty tty
-                  :sui sui :saui saui :sdui sdui :scui scui
-                   :lat lat :str str)))
+                  :sui sui :saui saui :sdui sdui :scui scui :lat lat
+                   :str str)))
 
 (defun find-ucxt-cuisui (cui sui &key (srl *current-srl*))
   (ensure-cui-integer cui)
@@ -696,6 +690,15 @@ is OBJNAME from TABLE where WHERE-NAME field = WHERE-VALUE with FIELDS"
                   :torule torule :tores tores :maprule maprule :maptype maptype :mapatn mapatn
                   :mapatv mapatv :cvf cvf)))
 
+(defun find-usmap-cui (cui)
+  (ensure-cui-integer cui)
+  (collect-umlisp-query (mrsmap (mapsetsab fromexpr fromtype rel rela toexpr totype cvf)
+                              nil mapsetcui cui)
+    (make-instance 'usmap :mapsetcui cui :mapsetsab mapsetsab
+                  :fromexpr fromexpr :fromtype fromtype
+                  :rel rel :rela rela :toexpr toexpr :totype totype
+                   :cvf cvf)))
+
 ;;;; Cross table find functions
 
 (defun find-ucon-tui (tui &key (srl *current-srl*))