(def x(g-normal 5 100))

;; rational function

(def y (/ (nth 0 x) ( + .5 (** (+ (nth 1 x) 1.5) 2))))
 

;; apply sir

(def out (sir-i-ii-model x y))

;; input  an integer when asked, suggested values : between 5 and 20,
 

;; now  remove the spin-plot window for sir-ii

;; to compare what is seen from sir-i window with the best view

(spin-plot  (list (nth 0 x) y (nth 1 x)))
 

;; to compare with pca model

(def out-pca (pca-model  x))

;; find  all 5 principal components

  (def newx(send out-pca :scores))

;; plot y  against the first two principal components

 (spin-plot (list (nth 0 newx) y (nth 1 newx)))