(defun assign_strata (g8ctrl g8regon g8urban g8minor)    
(case g8ctrl
      (1 (case g8regon (1 (setf h 1))
                       (2 (setf h 2))
                       (3 (setf h 3))
                       (4 (setf h 4))    
                       (8 (setf h 20))))
      (2 (case g8regon (1 (setf h 5))
                       (2 (setf h 6))
                       (3 (setf h 7))
                       (4 (setf h 8))    
                       (8 (setf h 20))))
      (3 (case g8regon (1 (setf h 9))
                       (2 (setf h 10))
                       (3 (setf h 11))
                       (4 (setf h 12))    
                       (8 (setf h 20))))
      (4 (case g8regon (1 (setf h 13))
                       (2 (setf h 14))
                       (3 (setf h 15))
                       (4 (setf h 16))    
                       (8 (setf h 20))))
)
)

(setf h_list(repeat 0 1035))
(dotimes (i 1035)
(assign_strata (nth i g8ctrl) (nth i g8regon) (nth i g8urban) (nth i g8minor))
(setf (nth i h_list) h)
)

(def school_id(nth 15 data-nonmiss))
(def g8ctrl(nth 9 data-nonmiss))
(def bysc47b(nth 0 data-nonmiss))
 (def bysc47c(nth 1 data-nonmiss))
(def bysc47d(nth 2 data-nonmiss))
 (def bysc47g(nth 3 data-nonmiss))
(def bysc47j(nth 4 data-nonmiss))
(def bysc47m(nth 5 data-nonmiss))
(def bysc47n(nth 6 data-nonmiss))
(def bysc47o(nth 7 data-nonmiss))
(def g8type(nth 13  data-nonmiss))
(def byscenrl(nth 8  data-nonmiss))
(def g8enrol(nth 10  data-nonmiss))
(def g8urban(nth 14  data-nonmiss))
(def g8regon(nth 12  data-nonmiss))
(def g8minor(nth 11 data-nonmiss))

(length (select h_list (which (= 1 h_list))))


(setf h1 (select g8enrol (which (= 1  h_list))))
(setf h2 (select g8enrol (which (= 2  h_list))))
(setf h3 (select g8enrol (which (= 3  h_list))))
(setf h4 (select g8enrol (which (= 4  h_list))))
(setf h5 (select g8enrol (which (= 5  h_list))))
(setf h6 (select g8enrol (which (= 6  h_list))))
(setf h7 (select g8enrol (which (= 7  h_list))))
(setf h8 (select g8enrol (which (= 8  h_list))))
(setf h9 (select g8enrol (which (= 9  h_list))))
(setf h10 (select g8enrol (which (= 10  h_list))))
(setf h11 (select g8enrol (which (= 11  h_list))))
(setf h12 (select g8enrol (which (= 12  h_list))))
(setf h13 (select g8enrol (which (= 13  h_list))))
(setf h14 (select g8enrol (which (= 14  h_list))))
(setf h15 (select g8enrol (which (= 15  h_list))))
(setf h16 (select g8enrol (which (= 16  h_list))))


(setf g8enrol-prop-h1 (/ (sum h1)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h2 (/ (sum h2)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h3 (/ (sum h3)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h4 (/ (sum h4)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h5 (/ (sum h5)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h6 (/ (sum h6)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h7 (/ (sum h7)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h8 (/ (sum h8)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h9 (/ (sum h9)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h10 (/ (sum h10)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h11 (/ (sum h11)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h12 (/ (sum h12)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h13 (/ (sum h13)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h14 (/ (sum h14)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h15 (/ (sum h15)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)
(setf g8enrol-prop-h16 (/ (sum h16)
                         (sum         
                         (sum h1) (sum h2) (sum h3) (sum h4) 
                         (sum h5) (sum h6) (sum h7) (sum h8)
                         (sum h9) (sum h10) (sum h11)  (sum h12)
                         (sum h13) (sum h14) (sum h15) (sum h16)
                         )
                      )
)

(load "g8enrol-prop-h1")
(load "g8enrol-prop-h2")
(load "g8enrol-prop-h3")
(load "g8enrol-prop-h4")
(load "g8enrol-prop-h5")
(load "g8enrol-prop-h6")
(load "g8enrol-prop-h7")
(load "g8enrol-prop-h8")
(load "g8enrol-prop-h9")
(load "g8enrol-prop-h10")
(load "g8enrol-prop-h11")
(load "g8enrol-prop-h12")
(load "g8enrol-prop-h13")
(load "g8enrol-prop-h14")
(load "g8enrol-prop-h15")
(load "g8enrol-prop-h16")

(mapcar #'(lambda (x)
              (* x r (/ 1 (elt alpha 0)))) alpha)