(XLISP::CHECK-FSL-VERSION 1 4) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function ARGS: (STRING CHARS &optional (START 0)) Returns the position of the first occurrence starting at position START of the string CHARS in the string STRING, or NIL if CHARS doesn't appear in STRING." FUNCTION-DOCUMENTATION STRING-FIND %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(STRING-FIND #K(#(77 1 3 1 1 1 0 6 23 2 5 1 2 23 2 6 1 1 2 43 3 5 7 37 76 1 6 7 6 2 45 6 5 6 23 3 6 2 3 6 16 6 7 21 1 2 5 6 38 4 8 38 5 9 2 43 7 5 10 23 6 8 6 1 2 8 7 9 10 21 0 3 4 8 11 0 7 17 6 6 16 6 7 21 1 2 5 6 11 0 4 11 0 4) #(0 29 48 74 77 88 91) #(NIL #(0) LENGTH ISEQ :START1 :END1 STRING-EQUAL) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function ARGS: (STRING POSITION CHARS) Inserts the string CHARS after POSITION in STRING. If POSITION is negative, then CHARS is inserted at the start of STRING; if POSITION >= 1 less than the length of STRING, then CHARS is inserted at the end of STRING." FUNCTION-DOCUMENTATION INSERT-CHARACTERS %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(INSERT-CHARACTERS #K(#(80 4 2 0 1 3 37 60 1 2 2 5 25 2 0 3 4 3 1 23 3 6 1 1 61 109 6 6 37 71 3 4 2 6 25 2 0 3 4 1 3 23 4 6 2 5 2 23 5 6 2 1 6 61 112 2 7 23 3 8 1 1 61 109 8 8 23 4 7 2 7 8 23 5 7 2 1 7 25 2 0 4 4 6 3 7) #(0 12 19 34 41) #(STRING 0 CONCATENATE LENGTH ISEQ SELECT) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function ARGS: (STRING POSITION LEN) Deletes LEN characters from STRING starting at POSITION." FUNCTION-DOCUMENTATION DELETE-CHARACTERS %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(DELETE-CHARACTERS #K(#(80 4 3 0 1 2 3 37 61 1 2 2 6 23 3 7 1 1 61 109 7 7 23 4 7 2 3 7 23 5 7 2 1 7 25 6 0 2 4 7 23 3 7 1 1 2 43 2 3 8 37 61 3 4 7 8 23 3 7 1 1 23 7 7 3 7 3 5 23 4 7 2 6 7 23 5 7 2 1 7 25 6 0 2 4 7 23 3 7 1 1 61 109 7 7 2 45 2 5 8 23 4 8 2 6 8 23 5 8 2 1 8 2 43 2 3 9 23 4 9 2 9 7 23 5 9 2 1 9 25 6 0 3 4 8 9) #(0 13 40 56 86) #(STRING 1 0 LENGTH ISEQ SELECT CONCATENATE -) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function ARGS: (STRING CHARS WITH-CHARS &key (TIMES 'ALL) (POSN 0)) Starting at position POSN, replace up to TIMES occurrences of the string CHARS in STRING with the string WITH-CHARS" FUNCTION-DOCUMENTATION STRING-REPLACE %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(STRING-REPLACE #K(#(77 3 4 0 0 0 5 4 4 0 1 2 3 10 23 6 10 1 2 23 6 11 1 3 23 6 12 1 1 50 2 1 2 4 6 2 45 12 10 13 1 3 0 4 13 3 61 13 9 14 21 0 5 4 14 23 7 15 3 5 10 7 3 62 15 12 14 21 0 6 7 14 11 0 1 23 8 14 3 1 2 5 21 0 8 9 14 23 9 15 3 1 14 10 61 109 14 16 23 10 15 3 15 16 3 38 11 16 2 45 13 8 17 38 12 18 2 43 14 11 19 25 13 0 7 15 2 3 16 17 18 19 11 0 1) #(0 35 42 45 55 67 72 75 87 132) #(ALL -2 1 0 (ALL 0) (:TIMES :POSN) LENGTH + STRING-FIND DELETE-CHARACTERS INSERT-CHARACTERS :TIMES :POSN STRING-REPLACE) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :RECODE-DIALOG #K(#(80 2 14 0 1 2 3 4 5 6 7 8 9 10 11 12 13 20 38 14 16 38 15 17 23 16 16 3 16 17 15 38 17 17 23 18 18 1 14 23 16 17 3 1 17 18 38 19 18 38 15 19 23 16 18 2 18 19 38 20 19 38 15 20 23 16 19 5 19 20 17 13 18 38 14 20 38 15 21 23 16 20 3 20 21 12 38 20 21 38 15 22 23 21 23 1 11 23 16 21 5 21 22 23 10 18 38 14 22 38 15 23 23 16 22 3 22 23 9 38 22 23 38 15 24 38 23 25 23 16 23 5 23 24 11 25 8 38 14 24 38 15 25 23 16 24 3 24 25 7 38 22 25 38 15 26 38 23 27 23 16 25 5 25 26 11 27 6 38 14 26 38 24 27 23 16 26 3 26 27 5 38 25 27 38 15 28 38 26 29 51 3 30 0 23 16 27 5 27 28 4 29 30 38 25 28 38 15 29 38 26 30 51 2 31 3 25 23 18 23 16 28 5 28 29 3 30 31 38 25 29 38 15 30 38 26 31 51 1 32 3 25 23 18 23 16 29 5 29 30 2 31 32 38 27 30 38 15 31 23 21 32 1 16 23 21 33 2 20 19 23 21 34 2 22 21 23 21 35 2 24 23 23 21 33 3 33 34 35 23 21 34 2 26 25 23 21 34 1 34 23 21 35 3 28 29 27 23 21 32 4 32 33 34 35 23 16 29 3 30 31 32 38 28 30 23 16 0 3 29 30 28 38 29 30 25 16 0 2 29 30 80 4 1 30 3 38 31 5 23 16 5 2 3 5 38 24 6 23 16 6 2 2 6 38 24 7 23 16 7 2 1 7 25 21 0 4 4 5 6 7 80 4 1 32 3 38 31 5 23 16 5 2 3 5 38 24 6 23 16 6 2 2 6 38 24 7 23 16 7 2 1 7 25 21 0 4 4 5 6 7 80 1 1 32 0 11 0 1) #(0 341 381 421) #("Again" "OK" "Cancel" "Recode expression" 60 "New variable name (optional)" 20 "Variable to recode" SECONDARY "" "Candidates" PRIMARY (:VARIATE :TEXT) "Recode a Variable" TEXT-ITEM-PROTO :NEW SEND :TYPES APPEND ALLOCATE-ITEMS-PROTO SUBORDINATE-LIST-PROTO LIST EDIT-TEXT-ITEM-PROTO :TEXT-LENGTH :TEXT MODAL-BUTTON-PROTO :ACTION MODAL-DIALOG-PROTO :DEFAULT-BUTTON :MODAL-DIALOG T :VALUE NIL) 0 NIL) DATASET-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:RECODE-DIALOG) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Recode variable..." RECODE-MENU-ITEM :RECODE RC-MENU-ITEM) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #(*ARC-DATA-MENU-ITEMS* MARK-AS-SPECIAL) 0 NIL) #K(#(80 1 1 0 2 38 1 2 23 2 2 1 2 38 0 3 23 3 2 2 2 3 25 4 0 2 1 2) #(0) #(*ARC-DATA-MENU-ITEMS* RECODE-MENU-ITEM LIST APPEND SET) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(*ARC-DATA-MENU-ITEMS*) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :RECODE #K(#(80 2 3 0 1 2 10 38 3 5 23 4 5 2 1 5 21 0 1 6 5 71 2 2 5 6 71 2 2 6 6 71 3 4 5 7 50 2 2 3 2 7 16 6 7 23 5 8 1 5 38 6 9 16 6 10 23 4 9 3 1 9 10 23 7 9 2 9 8 38 8 10 38 9 11 38 10 12 38 11 13 16 6 14 23 12 14 4 4 3 14 8 23 4 0 8 1 10 11 9 12 7 13 14 16 5 5 21 0 4 5 5 38 13 5 25 4 0 2 1 5 11 0 4 11 0 4) #(0 21 42 45 112 121 124) #("" "Recode of ~a: ~a" NIL :RECODE-DIALOG SEND FOURTH :NAME-TO-DATA RECODE-VARIABLE :VARIATE :DATA :NAME :INFO FORMAT :RECODE) 0 NIL) DATASET-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:RECODE) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args (x) Returns T if all of the elements of the sequence X are T, otherwise returns nil." FUNCTION-DOCUMENTATION ALL %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(ALL #K(#(80 2 2 0 1 1 21 0 1 4 1 16 1 4 21 0 2 3 4 17 1 4 25 2 0 1 4 11 0 2 11 0 3) #(0 11 19 27 30) #(NIL T ALL) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (variable expression) Recodes VARIABLE according to the instructions in the string EXPRESSION, returning the recoded variable. Recoding instructions are of the form (old values -> new-value), (min:max -> new-value), (min: -> new-value), or (:max -> new-value), and EXPRESSION consists of a sequence of such specifications; if a value matches more than one specification, the first one is effective. A list of old values can be separated by spaces or commas. In formulating recode expressions, you can use = in place of ->; to in place of : ; lo or low for the smallest value of the variable; and hi or high for the largest value. Finally, the last specification can be of the form (else -> new-value). Examples: (recode-variable education \"(:11->1)(12->2)(13:->3)\") (recode-variable numbers \"(1 3 5 -> \"odd\")(2,4,6=\"even\")\") (recode-variable \"(low to 30=1)(else=2)\")" FUNCTION-DOCUMENTATION RECODE-VARIABLE %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(RECODE-VARIABLE #K(#(80 3 26 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 14 23 26 29 3 2 28 27 23 26 29 3 29 26 25 23 26 29 3 29 27 24 23 26 29 3 29 23 22 23 27 30 1 1 39 28 31 23 29 31 2 31 30 23 30 31 1 31 21 0 1 2 31 23 31 32 1 30 1 3 0 17 32 21 0 4 5 31 23 32 33 1 30 1 6 0 17 33 23 33 34 1 1 23 33 35 1 1 23 34 36 4 21 20 29 19 23 35 36 1 36 23 36 36 1 36 23 37 37 1 34 23 38 0 3 34 37 18 0 36 37 16 37 38 21 1 7 39 37 16 38 39 50 2 8 12 16 39 23 39 40 1 38 37 35 10 9 40 3 71 2 2 38 41 50 2 11 10 41 13 23 40 0 2 12 38 71 3 4 38 40 23 39 41 1 35 23 41 35 2 40 41 1 38 23 42 40 2 15 38 21 0 13 35 40 23 43 40 3 33 11 38 23 43 40 3 33 10 40 23 43 40 3 32 9 40 23 43 40 3 32 8 40 23 44 41 2 15 40 50 2 14 17 41 4 23 39 42 1 40 50 2 15 16 42 6 71 3 4 40 42 42 42 13 41 1 18 0 17 41 1 18 0 17 41 21 0 19 20 41 71 2 2 40 41 3 76 1 41 41 23 45 41 1 41 23 46 42 1 40 23 38 0 3 35 41 42 1 38 23 44 41 2 15 40 50 2 21 24 41 7 23 39 42 1 40 50 2 22 23 42 5 23 46 42 1 40 42 42 13 41 1 25 0 17 41 1 25 0 17 41 21 0 26 27 41 16 40 41 3 71 1 41 41 71 3 4 40 42 3 76 1 42 42 23 47 41 2 41 42 23 45 41 1 41 23 48 42 1 40 23 38 0 3 35 41 42 1 38 23 44 41 2 15 40 50 2 28 31 41 7 23 39 42 1 40 50 2 29 30 42 6 71 3 4 40 42 42 42 13 41 1 32 0 17 41 1 32 0 17 41 21 0 33 34 41 16 40 41 3 71 1 41 41 23 45 41 1 41 23 46 42 1 40 23 38 0 3 35 41 42 1 38 23 40 0 2 12 40 1 38 23 39 40 1 38 2 45 40 14 40 23 49 40 1 40 23 50 40 2 38 40 23 36 41 1 38 16 41 41 23 36 42 1 38 71 2 2 42 42 50 2 37 36 42 13 23 40 0 2 12 38 23 51 42 2 40 34 23 38 0 3 35 42 41 17 37 37 16 37 38 21 1 7 39 37 39 52 39 23 29 39 2 39 36 23 42 39 2 18 39 21 0 41 40 39 23 37 40 1 1 38 53 41 23 38 0 3 35 40 41 11 0 35) #(0 82 89 92 97 104 107 158 167 178 189 195 213 224 264 275 286 291 294 299 328 340 351 362 367 370 375 418 430 441 452 457 460 465 492 500 545 551 564 575 595 610) #(3 0 5 4 1 LOW LO HIGH HI "Recode expression not well formed: ~a" -> 2 TO ELSE NIL ? ")" "(" STRING " " "," " -> " " to " ":" "->" "=" STRING-REPLACE REMOVE-MISSING NUMBERP MAPCAR ALL MIN MAX COPY-LIST CONCATENATE READ-FROM-STRING REVERSE FIND-MISSING SET-SELECT LENGTH ERROR REPEAT MEMBER SUBST POSITION WHICH FOURTH PAND FIFTH ISEQ SELECT WHERE FIRST NOT-A-NUMBER) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (a b) Parallel AND for the sequences A and B" FUNCTION-DOCUMENTATION PAND %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(PAND #K(#(80 3 2 0 1 5 23 2 5 1 1 23 2 6 1 2 2 109 5 6 5 23 3 5 2 4 5 23 2 6 1 5 0 3 7 37 60 1 5 7 6 0 7 8 23 4 9 2 1 7 21 0 2 3 9 23 4 9 2 2 7 1 4 0 4 9 23 5 0 3 5 8 9 61 112 7 7 37 60 1 5 7 6 11 0 5) #(0 41 55 63 66 83) #(0 NIL LENGTH REPEAT SELECT SET-SELECT) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (list1 list2) Finds the positions of all elements of LIST1 in LIST2 (or NIL). Example: (where '(1 2) '(1 2 3 4 1 2 3 4))" FUNCTION-DOCUMENTATION WHERE %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(WHERE #K(#(80 3 1 0 3 39 1 4 20 3 1 5 51 1 6 1 2 23 2 5 2 6 5 25 3 0 2 4 5 80 3 0 1 39 4 3 23 2 3 3 3 2 1 25 5 0 1 3) #(0 29) #(NIL APPEND MAP-ELEMENTS APPLY EQUAL WHICH) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (term) Returns the degrees of freedom for a term in the model." :DEGREES-OF-FREEDOM-TERM #K(#(80 3 0 2 38 0 3 23 1 3 2 1 3 38 2 4 23 1 3 3 3 4 2 25 3 0 1 3) #(0) #(:GET-DATASET SEND :DATA LENGTH) 0 NIL) GLIM-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:DEGREES-OF-FREEDOM-TERM) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Args: (hyp &key (right-hand-side 0) (quiet nil)) Method for computing the F-statistic for a linear hypothesis in a regression model. hyp is a hypothesis matrix, each row specifying a linear combination of coefficients, with as many columns as coefficients in the model. right-hand-side is a scalar or vector giving the right-hand-side of the hypothesis (default 0). If quiet is t, then a vector of the sum of squares, F, numerator and denominator df, and p-value for the hypothesis is returned. If quiet is nil (the default), then the method prints a report." :TEST-LINEAR-HYPOTHESIS #K(#(77 3 3 0 0 0 5 4 5 0 1 2 0 3 10 23 6 10 1 2 21 0 1 2 10 0 2 10 1 3 23 7 11 1 2 23 8 11 1 11 20 5 11 11 23 9 10 2 11 2 23 10 11 1 10 16 11 11 38 11 12 23 12 12 2 1 12 38 13 13 23 12 13 2 1 13 38 14 14 23 12 14 2 1 14 38 15 15 23 12 15 2 1 15 38 16 16 23 12 16 2 1 16 23 17 16 2 16 9 2 42 16 14 16 23 18 17 1 10 23 19 17 2 15 17 2 45 17 3 17 23 18 18 1 10 23 19 18 3 10 16 18 23 20 18 1 18 23 19 19 2 10 15 2 45 19 3 19 23 19 14 3 17 18 19 2 47 14 11 15 2 42 14 13 16 2 47 16 12 16 23 21 17 3 15 11 12 2 45 8 17 17 21 0 4 5 4 25 8 0 5 16 15 11 12 17 25 22 0 7 7 6 16 11 12 15 17) #(0 25 30 50 197 206) #(1 "sum of squares = ~a F[~a,~a] = ~a p = ~,4f~%" T 2 (0 NIL) (:RIGHT-HAND-SIDE :QUIET) MATRIXP LENGTH LIST MATRIX ARRAY-DIMENSIONS :DF SEND :SUM-OF-SQUARES :XTXINV :COEF-ESTIMATES :SIGMA-HAT ^ TRANSPOSE %* INVERSE F-CDF FORMAT) 0 NIL) NORMALREG-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:TEST-LINEAR-HYPOTHESIS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Args: (hyp &key (right-hand-side 0) (quiet nil)) Method for computing the Wald test for a linear hypothesis in a generalized linear model. hyp is a hypothesis matrix, each row specifying a linear combination of coefficients, with as many columns as coefficients in the model. right-hand-side is a scalar or vector giving the right-hand-side of the hypothesis (default 0). If quiet is t, then a vector of the chisquare, df, and p-value for the hypothesis is returned. If quiet is nil (the default), then the method prints a report." :TEST-LINEAR-HYPOTHESIS #K(#(77 3 3 0 0 0 4 3 4 0 1 2 0 7 23 5 9 1 2 21 0 1 2 9 0 2 9 1 3 23 6 10 1 2 23 7 10 1 10 20 5 10 10 23 8 9 2 10 2 23 9 10 1 9 16 10 10 38 10 11 23 11 11 2 1 11 38 12 12 23 11 12 2 1 12 23 13 13 1 9 23 14 13 2 12 13 2 45 13 3 13 23 13 14 1 9 23 14 14 3 9 11 14 23 15 14 1 14 23 14 15 2 9 12 2 45 15 3 15 23 14 11 3 13 14 15 23 16 12 2 11 10 2 45 8 12 12 21 0 4 5 4 25 7 0 3 11 10 12 25 17 0 5 7 6 11 10 12) #(0 24 29 49 142 149) #(1 "Wald chisquare statistic = ~a df = ~a p = ~,4f~%" T (0 NIL) (:RIGHT-HAND-SIDE :QUIET) MATRIXP LENGTH LIST MATRIX ARRAY-DIMENSIONS :XTXINV SEND :COEF-ESTIMATES TRANSPOSE %* INVERSE CHISQ-CDF FORMAT) 0 NIL) GLIM-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:TEST-LINEAR-HYPOTHESIS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 0 1 39 0 1 49 1 1 0 57 0) #(0) #(INVERSE INV) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(INVERSE #K(#(80 2 3 0 1 2 2 39 3 5 23 4 6 1 1 23 5 5 2 5 6 37 61 1 2 4 5 68 1 3 3 5 2 47 4 5 5 23 6 5 1 5 25 7 0 2 2 5 38 8 5 23 6 6 1 1 25 9 0 2 5 6) #(0 27 48) #((1 1) 0 1 * ARRAY-DIMENSIONS REDUCE LIST MATRIX INV APPLY) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Prints an analysis-of-variance table with an F-test for each term in the model, using the sum of squares for that term 'after' all of the other terms." :F-TESTS #K(#(80 2 13 0 1 2 3 4 5 6 7 8 9 10 11 12 13 38 13 15 23 14 15 2 1 15 23 15 16 1 15 51 12 17 1 1 23 16 17 2 17 15 23 17 18 1 17 38 18 19 23 14 19 2 1 19 23 15 19 1 19 37 35 1 2 18 19 23 19 0 2 14 2 11 0 9 23 20 18 1 16 23 21 18 2 18 17 38 22 19 23 14 19 2 1 19 21 0 3 4 19 23 23 18 2 3 18 23 15 19 1 18 23 24 19 1 19 23 15 20 1 18 23 20 20 1 20 38 22 21 23 14 21 2 1 21 21 0 5 6 21 38 25 22 23 26 23 3 19 10 20 38 27 24 23 14 21 5 1 22 23 24 14 1 7 0 9 21 38 28 22 23 14 22 2 1 22 23 19 0 3 14 13 22 23 19 0 2 14 12 23 19 0 2 14 11 38 22 22 23 14 22 2 1 22 21 0 8 9 22 23 26 23 2 21 10 23 26 24 2 21 8 23 26 25 2 21 7 23 26 26 2 21 6 23 19 0 7 14 5 23 24 25 26 4 0 10 21 37 60 10 11 21 16 23 26 0 2 17 21 3 61 21 18 22 23 29 22 1 22 23 26 22 3 19 22 20 38 25 23 38 27 24 23 14 22 5 1 23 22 24 14 23 26 23 2 22 10 23 26 24 2 22 8 23 26 25 2 22 7 23 26 26 2 22 6 23 26 27 2 15 21 23 19 0 7 14 5 23 24 25 26 27 61 112 21 21 37 60 10 11 21 16 11 0 9 80 3 0 1 38 30 3 25 14 0 3 1 3 2) #(0 67 76 101 107 141 165 168 210 245 254 343 346) #("~%Coefficients aliased, no tests calculated.~%" -1 "Constant" "~15g ~8d ~10g ~10,5f ~a~%" 4 1 2 NIL 0 "----------------------------------------------------------~%" "~%Sum of Squares df F p Source~%" "~%~% F-Tests for terms in Model ~a" T :PREDICTORS SEND LENGTH MAPCAR SUM :BASIS FORMAT ISEQ REPEAT :INTERCEPT ADJOIN IDENTITY-MATRIX :TEST-LINEAR-HYPOTHESIS SELECT :QUIET :NAME WHICH :DEGREES-OF-FREEDOM-TERM) 0 NIL) NORMALREG-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:F-TESTS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Prints a Wald chisquare test for each term in the model, testing the hypothesis that the term (i.e., set of coefficients) is zero." :WALD-TESTS #K(#(80 2 12 0 1 2 3 4 5 6 7 8 9 10 11 12 38 12 14 23 13 14 2 1 14 23 14 15 1 14 51 12 16 1 1 23 15 16 2 16 14 23 16 17 1 16 38 17 18 23 13 18 2 1 18 23 14 18 1 18 37 35 1 2 17 18 23 18 0 2 13 2 11 0 8 23 19 17 1 15 23 20 17 2 17 16 38 21 18 23 13 18 2 1 18 21 0 3 4 18 23 22 17 2 3 17 23 14 18 1 17 23 23 18 1 18 23 14 19 1 17 23 19 19 1 19 38 21 20 23 13 20 2 1 20 21 0 5 6 20 38 24 21 23 25 22 3 18 9 19 38 26 23 23 13 20 5 1 21 22 23 13 1 7 0 8 20 38 27 21 23 13 21 2 1 21 23 18 0 3 13 12 21 23 18 0 2 13 11 23 18 0 2 13 10 38 21 21 23 13 21 2 1 21 21 0 8 9 21 23 25 22 2 20 7 23 25 23 2 20 9 23 25 24 2 20 6 23 18 0 6 13 5 22 23 24 4 0 9 20 37 60 10 11 20 15 23 25 0 2 16 20 3 61 20 17 21 23 28 21 1 21 23 25 21 3 18 21 19 38 24 22 38 26 23 23 13 21 5 1 22 21 23 13 23 25 22 2 21 7 23 25 23 2 21 9 23 25 24 2 21 6 23 25 25 2 14 20 23 18 0 6 13 5 22 23 24 25 61 112 20 20 37 60 10 11 20 15 11 0 8 80 3 0 1 38 29 3 25 13 0 3 1 3 2) #(0 66 75 100 106 140 164 167 209 237 246 328 331) #("~%Coefficients aliased, no tests calculated.~%" -1 "Constant" "~8d ~10g ~10,5f ~a~%" 2 1 NIL 0 "------------------------------------------------~%" "~% df Wald Chisq. p Source~%" "~%~% Wald-Tests for terms in Model ~a" T :PREDICTORS SEND LENGTH MAPCAR SUM :BASIS FORMAT ISEQ REPEAT :INTERCEPT ADJOIN IDENTITY-MATRIX :TEST-LINEAR-HYPOTHESIS SELECT :QUIET :NAME WHICH :DEGREES-OF-FREEDOM-TERM) 0 NIL) GLIM-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:WALD-TESTS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :PCA #K(#(77 3 2 0 0 0 11 10 10 0 1 2 3 4 5 6 7 8 9 17 38 12 15 23 13 15 2 1 15 38 14 16 23 15 17 1 15 23 13 16 4 1 16 14 17 38 14 17 23 15 18 1 15 23 13 17 4 1 17 13 18 38 14 18 23 15 19 1 15 23 13 18 4 1 18 12 19 61 109 15 19 23 16 20 3 16 17 18 2 42 19 20 19 23 17 20 1 16 23 17 21 1 17 23 17 22 1 18 23 18 20 3 20 21 22 23 19 21 1 19 71 2 2 21 22 16 21 23 23 20 24 1 23 23 21 25 1 23 21 0 1 2 2 23 18 26 3 7 6 5 1 3 23 18 26 3 9 9 9 21 0 4 5 4 38 22 27 23 13 27 2 1 27 23 23 26 2 27 8 38 24 27 16 22 28 9 25 14 29 2 47 29 12 29 38 25 30 23 26 31 2 26 14 23 13 0 7 1 27 20 28 29 30 31 21 0 6 7 3 38 24 27 71 2 2 22 28 9 25 13 29 2 47 29 12 29 38 25 30 23 26 31 2 26 13 23 13 0 7 1 27 20 28 29 30 31 21 0 8 9 3 38 24 27 71 3 4 22 28 9 25 12 29 2 47 29 12 29 38 25 30 23 26 31 2 26 12 23 13 0 7 1 27 20 28 29 30 31 16 22 25 9 25 14 25 16 22 26 9 26 13 26 16 22 27 9 27 12 27 23 27 0 5 11 10 25 26 27 21 0 10 11 3 71 2 2 22 25 9 25 14 25 71 2 2 22 26 9 26 13 26 71 2 2 22 27 9 27 12 27 23 27 0 5 11 10 25 26 27 21 0 12 13 3 71 3 4 22 25 9 25 14 25 71 3 4 22 26 9 26 13 26 71 3 4 22 27 9 27 12 27 23 27 0 5 11 10 25 26 27 9 23 14 25 2 47 25 24 25 9 23 13 26 2 47 26 24 26 9 23 12 27 2 47 27 24 27 25 27 0 5 11 10 25 26 27) #(0 143 152 159 164 179 219 256 261 298 333 369 374 410) #(BLUE YELLOW RED 3 NIL "~10,6f ~10,6f ~10,6f~%" T 2 1 0 (T T NIL) (:COLORS :ALLAXES :ERASE) :NUM-POINTS SEND :POINT-COORDINATE ISEQ COVARIANCE-MATRIX MEAN LIST EIGEN SUM SQRT :BACK-COLOR REPEAT :DIRCOS :COLOR SELECT FORMAT) 0 NIL) SPIN-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:PCA) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :DIRCOS #K(#(77 3 5 0 0 0 1 0 0 4 38 2 6 2 42 4 3 7 2 43 2 7 7 2 42 4 3 8 2 45 2 8 8 38 3 9 25 4 0 6 1 6 7 8 9 5) #(0) #((NIL) (:COLOR) :ABLINE :COLOR SEND) 0 NIL) SPIN-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:DIRCOS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :ABLINE #K(#(77 3 4 0 0 0 1 0 0 4 21 0 1 2 4 0 4 5 1 3 38 2 5 23 3 5 2 1 5 38 4 6 23 5 7 2 2 3 38 6 8 25 3 0 5 1 6 7 8 5) #(0 15 20 29) #((NIL) (:COLOR) :DRAW-COLOR SEND :ADD-LINES MAKE-PAIRS :COLOR) 0 NIL) SPIN-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:ABLINE) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(MAKE-PAIRS #K(#(80 3 0 2 23 0 3 1 1 23 1 3 1 3 51 1 4 2 2 1 25 2 0 2 4 3 80 4 0 2 9 2 3 4 9 1 3 5 25 3 0 2 4 5) #(0 26) #(LENGTH ISEQ MAPCAR LIST) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (variables &key (miss not-a-number) (test #'eq)) VARIABLES is a list of sequences, containing the variables; MISS gives the missing data code (default not-a-number); TEST gives the equality test (default #'eq). The function returns a list a sequences containing only observations with valid data for all variables. Example: (delete-missing (list x y z) :miss '*)" FUNCTION-DOCUMENTATION DELETE-MISSING %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(DELETE-MISSING #K(#(77 3 2 0 0 0 1 0 0 2 21 2 1 2 2 0 2 4 1 3 38 2 4 21 2 4 5 3 0 3 5 1 6 39 3 5 51 8 4 2 4 5 23 4 4 2 4 1 23 5 4 1 4 23 6 4 1 4 16 1 5 23 7 5 1 5 23 8 5 1 5 23 9 4 2 5 4 23 10 4 1 4 51 7 4 1 4 25 4 0 2 4 1 80 3 0 0 25 11 0 2 2 1 80 4 0 1 23 12 4 3 2 1 3 25 13 0 1 4) #(0 15 20 23 28 33 36 93 103) #((%NOT-SUPPLIED %NOT-SUPPLIED) (:MISS :TEST) NOT-A-NUMBER EQ MAPCAR COMBINE REMOVE-DUPLICATES LENGTH ISEQ SET-DIFFERENCE SORT-DATA SELECT MAP-ELEMENTS WHICH) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (x) X is a data matrix; correlations are calculated and returned for the columns of X, removing rows with missing data, represented by not-a-number." FUNCTION-DOCUMENTATION CORRELATION-MATRIX %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(CORRELATION-MATRIX #K(#(80 2 1 0 2 23 1 3 1 1 23 2 3 1 3 39 3 4 23 4 3 2 4 3 23 5 3 1 3 23 6 4 1 3 23 7 4 2 4 2 23 8 4 2 4 4 2 47 3 4 0 57 0) #(0) #(0.5 COLUMN-LIST DELETE-MISSING BIND-COLUMNS APPLY COVARIANCE-MATRIX DIAGONAL ^ OUTER-PRODUCT) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((XLISP::ADD-METHOD NORMALREG-PROTO (QUOTE :VIFS) (FUNCTION (LAMBDA (SELF) (BLOCK :VIFS (WHEN (NOT (SEND SELF :INTERCEPT)) (FORMAT T "~%No intercept, no VIFs calculated.~%") (RETURN FROM :VIFS)) (LET* ((R (CORRELATION-MATRIX (SEND SELF :X))) (R-INV (INVERSE R)) (DET-R (DETERMINANT R)) (TERMS (SEND SELF :PREDICTORS)) (NUMBER-TERMS (LENGTH TERMS)) (DF (MAPCAR (FUNCTION (LAMBDA (TERM) (SEND SELF :DEGREES-OF-FREEDOM-TERM TERM))) TERMS))) (WHEN (/= (SUM DF) (LENGTH (SEND SELF :BASIS))) (FORMAT T "~%Coefficients aliased, no VIFs calculated.~%") (RETURN-FROM :VIFS)) (LET ((POSITIONS (REPEAT (ISEQ NUMBER-TERMS) DF))) (FORMAT T "~%Generalized Variance-Inflation Factors for Model ~a" (SEND SELF :NAME)) (FORMAT T "~% df GVIF GVIF^(1/2df) term") (FORMAT T "~%----------------------------------------~%") (DOTIMES (TERM NUMBER-TERMS) (LET* ((DF-TERM (SELECT DF TERM)) (POS-TERM (WHICH (= TERM POSITIONS))) (POS-OTHERS (WHICH (/= TERM POSITIONS))) (VIF (IF (= DF-TERM 1) (AREF R-INV (FIRST POS-TERM) (FIRST POS-TERM)) (/ (* (DETERMINANT (SELECT R POS-TERM POS-TERM)) (DETERMINANT (SELECT R POS-OTHERS POS-OTHERS))) DET-R))) (ROOT-VIF (^ VIF (/ 1 (* 2 DF-TERM))))) (FORMAT T "~4d ~10g ~10g ~a~%" DF-TERM VIF ROOT-VIF (SELECT TERMS TERM))))))))) "Message args: () Prints the generalized variance-inflation factor for each term in the model. For 1 df effects, this is just the usual VIF.") EVAL) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:VIFS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("F-tests" F-TESTS-MENU-ITEM :F-TESTS RC-MENU-ITEM) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Variance-Inflation Factors" VIFS-MENU-ITEM :VIFS RC-MENU-ITEM) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #(*ARC-LINEAR-REGRESSION-MENU-ITEMS* MARK-AS-SPECIAL) 0 NIL) #K(#(80 1 4 0 1 2 3 2 38 0 5 23 4 6 3 4 3 2 23 5 5 2 5 6 25 6 0 2 1 5) #(0) #(*ARC-LINEAR-REGRESSION-MENU-ITEMS* DASH VIFS-MENU-ITEM F-TESTS-MENU-ITEM LIST APPEND SET) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(*ARC-LINEAR-REGRESSION-MENU-ITEMS*) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Wald tests" WALD-TESTS-MENU-ITEM :WALD-TESTS RC-MENU-ITEM) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #(*ARC-GLM-REGRESSION-MENU-ITEMS2* MARK-AS-SPECIAL) 0 NIL) #K(#(80 1 3 0 1 2 2 38 0 4 23 3 5 2 3 2 23 4 4 2 4 5 25 5 0 2 1 4) #(0) #(*ARC-GLM-REGRESSION-MENU-ITEMS2* DASH WALD-TESTS-MENU-ITEM LIST APPEND SET) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(*ARC-GLM-REGRESSION-MENU-ITEMS2*) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (x y &key (other-x nil) (tol .01) (score-test t)) Calculates and returns the maximum-likelihood estimates of the transformation parameters for the Box-Tidwell regression model. x a list of lists or vectors containing the independent variables that are candidates for transformation y a list or vector of dependent-variable values :other-x a list of lists or vectors containing independent variables that are not candidates for transformation; default NIL. :tol convergence criterion. When the maximum change in the transformation estimates is less than tol iterations cease. :quiet if T, iteration history is suppressed; default NIL. Example: (box-tidwell (list income education) prestige :other-x (list percent-women)" FUNCTION-DOCUMENTATION BOX-TIDWELL %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(BOX-TIDWELL #K(#(77 3 3 0 0 0 7 6 6 0 1 2 3 4 5 16 39 8 12 23 9 12 2 12 1 39 10 13 23 9 12 3 13 1 12 23 11 13 1 1 23 12 13 2 11 13 23 13 14 2 12 3 38 14 15 23 15 14 4 14 2 15 10 38 16 15 23 17 15 2 14 15 23 18 15 2 15 13 23 19 16 1 12 2 42 12 16 16 23 13 16 3 16 12 3 38 14 17 23 15 16 4 16 2 17 10 38 16 17 23 17 17 2 16 17 23 18 17 2 17 13 2 47 17 15 18 61 112 18 0 38 20 18 23 17 18 2 16 18 23 18 18 2 18 13 2 47 17 18 17 23 21 17 1 17 23 22 18 1 17 2 45 11 18 18 2 42 9 18 18 23 23 0 4 8 7 17 18 0 12 17 23 19 18 1 12 2 42 12 18 18 0 14 19 38 16 20 23 17 20 2 19 20 23 18 20 2 20 13 0 16 21 38 16 22 23 17 22 2 21 22 23 18 22 2 22 13 2 47 22 20 23 61 112 23 23 0 11 24 0 23 25 2 45 25 24 26 23 21 26 1 26 23 8 26 1 26 37 62 4 1 4 26 21 0 3 2 5 23 23 0 3 8 6 25 23 24 17 2 12 25 23 19 26 1 17 2 42 17 26 18 23 13 26 2 17 3 38 14 27 23 15 19 4 26 2 27 10 38 16 26 23 17 26 2 19 26 23 18 20 2 26 13 23 13 26 3 18 17 3 38 14 27 23 15 21 4 26 2 27 10 38 16 26 23 17 26 2 21 26 23 18 22 2 26 13 2 47 22 20 26 61 112 26 23 0 25 24 2 42 25 23 25 2 45 25 24 26 23 21 26 1 26 23 8 26 1 26 37 62 4 1 4 26 11 0 25) #(0 263 268 275 394) #("~% powers = ~a " "score statistics ~a~%p-values ~a" T 2 NIL 1 (NIL 0.01 NIL) (:OTHER-X :TOL :QUIET) MAX MAPCAR / LENGTH ISEQ APPEND :PRINT REGRESSION-MODEL :COEF-ESTIMATES SEND SELECT LOG :COEF-STANDARD-ERRORS ABS NORMAL-CDF FORMAT ^) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args: (y) Compute constructed variable for Box-Cox transformation of y" FUNCTION-DOCUMENTATION BOX-COX-CONSTRUCTED-VARIABLE %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(BOX-COX-CONSTRUCTED-VARIABLE #K(#(80 2 1 0 1 23 1 3 1 1 2 47 1 3 3 23 2 3 1 3 2 45 3 2 3 2 42 1 3 0 57 0) #(0) #(1 GEOMETRIC-MEAN LOG) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 3 0 1 2 4 38 3 4 39 4 5 23 5 6 1 4 21 0 1 2 6 0 4 6 1 3 23 6 6 1 4 23 7 4 2 5 6 0 3 5 16 5 6 21 1 4 5 5 38 8 7 23 9 0 3 4 7 6 17 5 5 16 5 6 21 1 4 5 5 38 10 7 23 9 0 4 4 7 2 1 23 11 0 2 2 4 11 0 2) #(0 23 28 33 50 71) #((LOGL INFO GRADIENT COEF-ESTIMATES WEIGHTS-NAME RESPONSE-CATEGORIES W MISSING NEEDS-COMPUTING COMPUTE-INFO-MATRIX MENU) MLOGIT-MODEL-PROTO NIL GLIM-PROTO MAKE-OBJECT LISTP LIST APPLY :ADD-SLOT SEND :MAKE-PROTOTYPE SET) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :ISNEW #K(#(77 2 2 0 2 0 1 0 63 4 3 0 25 2 0 0) #(0) #(T NEEDS-COMPUTING CALL-NEXT-METHOD) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:ISNEW) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional (compute nil set)) If T, compute information matrix. If nil, use scores to approximate the information matrix." :COMPUTE-INFO-MATRIX #K(#(77 1 2 1 2 2 0 1 2 64 2 5 21 2 1 2 2 0 2 6 1 3 0 3 6 21 0 4 5 5 63 4 6 0 62 4 0 57 0) #(0 17 22 25 30 34) #(NIL COMPUTE-INFO-MATRIX #(%NOT-SUPPLIED)) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:COMPUTE-INFO-MATRIX) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional (needs nil set)) If T, model needs re-computing." :NEEDS-COMPUTING #K(#(77 1 2 1 2 2 0 1 2 64 2 5 21 2 1 2 2 0 2 6 1 3 0 3 6 21 0 4 5 5 63 4 6 0 62 4 0 57 0) #(0 17 22 25 30 34) #(NIL NEEDS-COMPUTING #(%NOT-SUPPLIED)) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:NEEDS-COMPUTING) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional cats) Set or retrieve response variable categories." :RESPONSE-CATEGORIES #K(#(77 1 2 1 2 2 0 1 0 21 0 1 2 2 63 4 2 0 62 4 0 57 0) #(0 14 18) #(NIL RESPONSE-CATEGORIES #(NIL)) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:RESPONSE-CATEGORIES) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional miss) Set or retrieve logical list of missing-values indicators." :MISSING #K(#(77 1 2 1 2 2 0 1 0 21 0 1 2 2 63 4 2 0 62 4 0 57 0) #(0 14 18) #(NIL MISSING #(NIL)) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:MISSING) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional name) Set or retrieve name of case-weight variable." :WEIGHTS-NAME #K(#(77 1 2 1 2 2 0 1 0 21 0 1 2 2 63 4 2 0 62 4 0 57 0) #(0 14 18) #(NIL WEIGHTS-NAME #(NIL)) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:WEIGHTS-NAME) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional w) Set or retrieve matrix of response indicators, filtering out missing and excluded data." :W #K(#(77 1 2 1 2 2 0 1 3 21 0 1 2 2 63 4 2 0 62 4 5 38 3 6 23 4 6 2 1 6 23 5 6 1 6 38 6 7 23 4 7 2 1 7 23 7 0 3 6 7 3 25 8 0 2 5 6) #(0 14 18) #(NIL W #(NIL) :INCLUDED SEND COPY-LIST :MISSING SET-SELECT INCLUDED) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:W) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional info) Set or retrieve information matrix, recomputing model if necessary." :INFO #K(#(77 1 2 1 3 3 0 1 2 2 21 0 1 2 2 63 5 2 0 1 4 38 4 6 23 5 6 2 1 6 21 0 3 4 6 38 6 6 38 7 7 23 5 0 4 1 6 7 3 62 5 0 57 0) #(0 15 21 35 49) #(T NIL INFO #(NIL) :NEEDS-COMPUTING SEND :COMPUTE :QUIET) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:INFO) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional logL) Set or retrieve log-likelihood, recomputing model if necessary." :LOGL #K(#(77 1 2 1 3 3 0 1 2 2 21 0 1 2 2 63 5 2 0 1 4 38 4 6 23 5 6 2 1 6 21 0 3 4 6 38 6 6 38 7 7 23 5 0 4 1 6 7 3 62 5 0 57 0) #(0 15 21 35 49) #(T NIL LOGL #(NIL) :NEEDS-COMPUTING SEND :COMPUTE :QUIET) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:LOGL) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional gradient) Set or retrieve gradient vector. If estimates fail to converge, contains list of log-likelihood, gradient, XtX matrix, and estimates at last iteration." :GRADIENT #K(#(77 1 2 1 3 3 0 1 2 2 21 0 1 2 2 63 5 2 0 1 4 38 4 6 23 5 6 2 1 6 21 0 3 4 6 38 6 6 38 7 7 23 5 0 4 1 6 7 3 62 5 0 57 0) #(0 15 21 35 49) #(T NIL GRADIENT #(NIL) :NEEDS-COMPUTING SEND :COMPUTE :QUIET) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:GRADIENT) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional (coefs nil set)) Set or retrieve coefficients, recomputing model if necessary. Current coefficients are used as start values; if nil, then zeroes are used as start values." :COEF-ESTIMATES #K(#(77 1 2 1 3 3 0 1 2 4 64 2 6 21 2 1 2 2 0 2 7 1 3 0 4 7 21 0 4 5 6 63 5 7 0 1 7 38 4 8 23 5 8 2 1 8 21 0 6 7 8 38 6 8 38 7 9 23 5 0 4 1 8 9 3 62 5 0 57 0) #(0 18 23 26 31 37 51 65) #(T NIL COEF-ESTIMATES #(%NOT-SUPPLIED) :NEEDS-COMPUTING SEND :COMPUTE :QUIET) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:COEF-ESTIMATES) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns coefficient standard errors, recomputing the model if necessary." :COEF-STANDARD-ERRORS #K(#(80 2 2 0 1 3 38 2 4 23 3 4 2 1 4 21 0 1 2 4 38 4 5 38 5 6 23 3 0 4 1 5 6 2 38 6 4 23 3 4 2 1 4 23 7 4 1 4 23 8 4 1 4 25 9 0 1 4) #(0 20 34) #(T NIL :NEEDS-COMPUTING SEND :COMPUTE :QUIET :INFO INVERSE DIAGONAL SQRT) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:COEF-STANDARD-ERRORS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns the covariance matrix of the coefficients, by inverting the contents of the information-matrix slot, recomputing the model if necessary." :XTXINV #K(#(80 2 2 0 1 3 38 2 4 23 3 4 2 1 4 21 0 1 2 4 38 4 5 38 5 6 23 3 0 4 1 5 6 2 38 6 4 23 3 4 2 1 4 25 7 0 1 4) #(0 20 34) #(T NIL :NEEDS-COMPUTING SEND :COMPUTE :QUIET :INFO INVERSE) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:XTXINV) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Args (list) Returns list rotated one position to the right, with last element at the front." FUNCTION-DOCUMENTATION ROTATE %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(ROTATE #K(#(80 2 0 2 23 0 2 1 1 23 1 3 1 1 17 3 3 23 1 3 1 3 25 2 0 2 2 3) #(0) #(LAST REVERSE COMBINE) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&optional new-y) Sets or retrieves the response variable y, also setting the response indicators w and response-categories slots." :Y #K(#(77 1 2 1 3 3 0 1 2 8 21 0 1 2 2 63 5 2 6 38 4 7 23 5 8 1 2 38 6 9 39 7 10 23 8 8 3 8 9 10 23 9 8 1 8 23 10 7 3 1 7 8 38 11 8 23 10 8 3 1 8 3 38 12 9 38 4 10 23 10 10 2 1 10 23 13 11 1 10 23 14 11 1 11 23 15 11 1 11 23 16 10 3 2 10 11 23 5 11 1 10 23 17 11 1 11 23 14 11 1 11 39 18 12 51 6 13 1 10 23 19 13 2 13 11 23 20 10 2 12 13 23 10 9 3 1 9 10 23 21 0 4 6 7 8 9 62 5 0 57 0 80 3 2 22 23 0 37 61 4 5 2 1 11 0 3 11 0 4 80 3 0 1 51 3 3 1 2 25 19 0 2 3 1) #(0 15 146 151 163 166 169) #(T NIL Y #(NIL) :RESPONSE-CATEGORIES REMOVE-MISSING :TEST EQUAL REMOVE-DUPLICATES SORT-DATA SEND :NEEDS-COMPUTING :W LENGTH ISEQ ROTATE RECODE-VALUES MAX BIND-COLUMNS MAPCAR APPLY LIST 1 0) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:Y) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns the model-matrix, filtering for missing and excluded cases." :X-MATRIX #K(#(80 2 3 0 1 2 6 38 3 5 23 4 5 2 1 5 38 5 6 23 4 6 2 1 6 23 6 7 2 6 4 38 7 8 23 4 8 2 1 8 23 8 8 1 8 38 9 9 23 4 9 2 1 9 21 0 1 2 5 23 10 10 2 2 7 23 11 10 2 10 6 1 5 23 12 10 1 6 21 0 3 4 10 0 6 10 1 5 23 13 10 2 7 2 23 14 10 2 10 6 23 15 0 3 8 9 3 25 16 0 2 10 8) #(0 59 73 83 88 100) #(1 NIL 0 :INTERCEPT SEND :X ARRAY-DIMENSION :INCLUDED COPY-LIST :MISSING REPEAT BIND-COLUMNS MATRIXP LIST MATRIX SET-SELECT INCLUDED) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:X-MATRIX) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Prints a summary of the model, recomputing the model if necessary." :SUMMARY #K(#(80 2 19 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 27 38 19 21 23 20 21 2 1 21 21 0 1 2 21 38 21 22 23 20 0 2 1 22 38 22 21 23 20 21 2 1 21 38 23 22 23 20 21 2 21 22 38 23 22 23 20 22 2 1 22 38 24 23 23 20 23 2 1 23 23 25 23 1 23 38 26 24 23 20 24 2 1 24 38 27 25 23 20 25 2 1 25 38 28 26 23 20 26 2 1 26 2 45 26 25 26 38 29 27 23 20 27 2 1 27 38 30 28 23 20 28 2 1 28 38 31 29 23 20 29 2 1 29 38 32 30 23 20 30 2 1 30 21 0 3 4 27 23 33 31 2 2 28 1 5 0 28 31 23 34 32 1 31 23 35 33 1 32 38 36 34 23 20 34 2 1 34 23 34 35 1 34 2 45 35 20 35 38 37 36 23 20 36 2 1 36 38 38 37 23 20 37 2 1 37 23 39 38 2 35 32 23 40 38 2 38 36 23 39 39 2 35 32 23 40 39 2 39 37 38 41 40 23 20 40 2 1 40 23 42 0 3 23 24 19 38 43 41 23 20 41 2 1 41 39 44 42 23 45 42 2 42 23 23 46 42 1 42 23 47 41 2 41 42 23 48 0 4 18 17 21 22 23 48 0 2 18 16 23 48 0 3 18 15 29 23 48 0 3 18 14 28 23 48 0 3 18 13 30 23 48 0 3 18 12 25 37 35 6 9 26 6 37 60 7 8 26 5 23 48 0 3 18 3 41 1 9 23 48 0 3 18 4 26 23 35 41 1 35 16 41 42 21 1 10 13 41 2 43 42 20 43 23 47 43 2 34 43 16 34 44 23 48 0 6 18 9 29 43 29 44 23 48 0 2 18 8 0 33 43 16 43 44 21 1 11 12 43 23 47 45 3 38 42 44 23 47 46 3 39 42 44 23 47 47 2 31 44 23 48 0 5 18 7 45 46 47 17 43 43 16 43 44 21 1 11 12 43 17 41 41 16 41 42 21 1 10 13 41 2 42 11 40 43 61 45 43 43 25 48 0 3 18 10 43) #(0 37 46 160 168 171 335 341 350 357 370 411 451 462) #("Constant" "~%Cases excluded: ~a" "~%Number of cases excluded = ~a" 20 0 "~%~10f ~10f ~a" "~%Coefficient Std. Error" "~%~%Logit [~a=~a/~a=~a]" "~%~%-2 log-likelihood = ~10f~%" 2 "~%Number of cases used = ~a" "~%Case Weights = ~a" "~%Predictors = ~a" "~%Response = ~a" "~%Multinomial Logit Model" "~%Data set = ~a, Name of Model = ~a" T NIL 1 :NEEDS-COMPUTING SEND :COMPUTE :GET-DATASET :NAME :INCLUDED COPY-LIST :MISSING :NUM-INCLUDED :NUM-CASES :INTERCEPT :PREDICTOR-NAMES :RESPONSE-NAME :WEIGHTS-NAME COMBINE LENGTH ISEQ :RESPONSE-CATEGORIES :COEF-ESTIMATES :COEF-STANDARD-ERRORS LIST MATRIX :LOGL SET-SELECT :CASE-LABELS NOT MAPCAR WHICH SELECT FORMAT) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:SUMMARY) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns the log-likelihood, gradient vector, and XtX matrix at the current parameter estimates; these are left in the gradient slot if there is no convergence." :DERIV #K(#(80 2 3 0 1 2 23 38 3 5 38 4 6 23 5 6 2 1 6 38 6 7 23 5 7 2 1 7 38 7 8 23 5 8 2 1 8 38 8 9 23 5 9 2 1 9 23 9 10 2 8 4 23 9 11 2 8 3 23 9 12 2 7 3 0 4 13 23 10 14 2 12 11 23 11 14 2 14 6 23 12 14 1 14 23 13 15 1 11 23 14 16 2 4 12 2 42 11 12 17 23 10 18 2 17 17 2 42 17 17 19 23 14 19 2 4 19 23 11 17 2 18 19 0 4 18 37 60 1 6 18 10 23 15 19 3 8 18 15 0 4 20 0 3 21 0 2 22 0 4 23 37 60 2 3 23 12 23 15 24 3 14 15 23 23 16 24 2 19 24 68 7 18 23 25 2 42 25 24 25 2 43 20 25 20 23 17 24 1 24 2 43 21 24 21 61 112 23 23 37 60 2 3 23 12 23 15 24 2 9 18 23 18 25 1 21 61 45 25 25 2 43 20 25 25 2 42 24 25 24 2 43 13 24 13 0 4 24 37 60 4 5 24 12 23 15 25 3 14 15 24 23 16 25 2 19 25 23 17 25 1 25 2 47 25 21 25 68 7 18 24 26 2 42 26 19 26 23 16 27 2 25 19 2 45 26 27 25 20 25 22 22 23 15 26 2 16 24 23 15 27 2 9 18 2 42 27 25 27 2 43 26 27 25 23 19 0 3 16 24 25 61 112 24 24 37 60 4 5 24 12 39 20 25 23 21 26 1 22 23 22 25 2 25 26 23 15 26 2 9 18 23 23 27 2 25 25 2 42 26 27 25 2 43 17 25 17 61 112 18 18 37 60 1 6 18 10 68 13 4 4 19 39 20 20 23 22 20 2 20 16 23 10 19 3 19 20 17 25 5 0 3 1 5 19) #(0 132 157 205 244 331 377) #(NIL 1 0 :GRADIENT :COEF-ESTIMATES SEND :W :X-MATRIX :WEIGHTS ARRAY-DIMENSION LIST MATRIX TRANSPOSE ISEQ REPEAT SELECT MATMULT EXP LOG SET-SELECT COMBINE REVERSE APPLY OUTER-PRODUCT) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:DERIV) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Computes and returns the information matrix." :INFO-MATRIX #K(#(80 2 2 0 1 22 38 2 4 38 3 5 23 4 5 2 1 5 38 5 6 23 4 6 2 1 6 38 6 7 23 4 7 2 1 7 38 7 8 23 4 8 2 1 8 23 8 9 2 7 3 23 8 10 2 7 2 23 8 11 2 6 2 23 9 12 2 11 10 23 10 12 2 12 5 23 11 12 1 12 23 12 13 1 10 23 13 14 1 5 23 12 15 1 11 23 9 16 2 14 14 2 42 14 14 17 23 14 17 2 3 17 23 10 16 2 16 17 0 3 17 37 60 1 4 17 9 23 15 18 3 7 17 13 0 16 19 23 15 20 2 8 17 0 3 21 0 3 22 37 60 2 3 22 11 23 15 23 3 12 13 22 23 16 23 2 18 23 23 17 23 1 23 2 43 21 23 21 61 112 22 22 37 60 2 3 22 11 2 43 2 21 23 23 15 23 3 23 3 3 39 18 24 51 8 25 5 15 23 13 12 18 23 19 25 2 25 15 23 20 23 2 24 25 2 42 20 23 23 2 43 19 23 16 61 112 17 17 37 60 1 4 17 9 25 4 0 3 1 4 16 80 7 2 21 1 4 37 61 6 7 6 5 23 15 9 3 3 2 6 23 16 9 2 4 9 23 15 9 3 9 8 8 23 11 10 1 4 23 16 10 2 10 4 23 17 11 1 9 23 17 12 1 9 2 45 1 12 12 23 22 9 3 10 11 12 23 23 10 2 1 7 2 47 9 10 0 57 0 23 11 9 1 4 23 16 9 2 9 4 23 15 10 3 3 2 5 23 15 11 3 3 2 6 2 43 10 11 10 23 16 10 2 4 10 23 17 10 1 10 23 15 10 3 10 8 8 2 42 9 10 9 23 23 10 2 1 7 2 47 9 10 9 61 45 9 0 57 0 80 7 0 2 39 24 7 51 5 8 5 2 3 4 5 6 23 19 8 2 8 1 25 20 0 2 7 8) #(0 127 155 188 244 251 263 329 399) #(1 0 :INFO :COEF-ESTIMATES SEND :W :X-MATRIX :WEIGHTS ARRAY-DIMENSION LIST MATRIX TRANSPOSE ISEQ LENGTH REPEAT SELECT MATMULT EXP BIND-ROWS MAPCAR APPLY 2 * ^ BIND-COLUMNS) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:INFO-MATRIX) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function args: (x y) A convergence check from Tierney (1996)." FUNCTION-DOCUMENTATION MAX-RELATIVE-ERROR %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(MAX-RELATIVE-ERROR #K(#(80 3 0 3 2 45 1 2 3 23 0 3 1 3 38 1 4 23 2 4 1 4 23 0 5 1 2 23 3 4 2 4 5 2 47 3 4 3 25 4 0 1 3) #(0) #(ABS MACHINE-EPSILON SQRT PMAX MAX) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #(0.001 T *TOLERANCE* MARK-AS-SPECIAL) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(*TOLERANCE*) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&key quiet (tol *tolerance*) (maxiter 20) Computes maximum-likelihood estimates by the modified method of scoring. If :quiet is T, then the iteration history is suppressed; :tot (default *tolerance*) is the maximum relative change in the coefficients for convergence; :maxiter (default 20) is the maximum number of iterations." :COMPUTE #K(#(77 3 2 0 0 0 9 8 8 0 1 2 3 4 5 6 7 9 21 2 1 2 3 0 3 13 1 3 38 10 13 38 11 14 23 12 0 3 1 14 12 38 13 14 23 12 14 2 1 14 38 14 15 23 12 15 2 1 15 23 15 0 2 15 11 23 15 15 2 15 10 23 15 16 2 14 10 38 16 17 23 12 17 2 1 17 21 0 4 5 17 0 17 14 1 6 2 42 15 16 18 23 17 14 2 11 18 0 14 15 38 16 16 23 12 0 3 1 16 14 0 11 16 37 60 7 12 16 4 0 14 15 38 18 17 23 12 17 2 1 17 16 17 18 71 2 2 17 19 71 3 4 17 20 21 0 9 8 2 23 19 0 4 6 5 16 18 23 20 21 1 20 23 21 21 2 21 19 2 43 14 21 14 38 16 21 23 12 0 3 1 21 14 23 22 21 2 14 15 37 62 10 11 13 21 23 23 13 4 18 19 20 14 1 13 61 112 16 16 37 60 7 12 16 4 0 12 13 21 0 15 14 13 23 24 0 1 7 38 25 14 23 26 15 2 13 11 23 12 0 3 1 14 15 38 27 14 23 26 15 2 13 10 23 12 0 3 1 14 15 38 28 14 38 29 15 23 12 15 2 1 15 21 0 16 17 15 38 30 16 23 12 15 2 1 16 1 18 23 26 15 2 13 8 23 12 0 3 1 14 15 38 16 14 23 26 15 2 13 9 25 12 0 3 1 14 15) #(0 23 28 31 91 96 107 129 159 167 205 215 225 228 233 238 287 298 304) #("~&Iteration ~a: logL = ~a~%" T "Convergence failed, maximum iterations exceeded." 2 3 1 0 NIL (NIL %NOT-SUPPLIED 20) (:QUIET :TOL :MAXITER) *TOLERANCE* :NEEDS-COMPUTING SEND :W :X-MATRIX ARRAY-DIMENSION :COEF-ESTIMATES REPEAT :DERIV FORMAT INVERSE MATMULT MAX-RELATIVE-ERROR LIST ERROR :LOGL SELECT :GRADIENT :INFO :COMPUTE-INFO-MATRIX :INFO-MATRIX) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:COMPUTE) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Fit multinomial logit model..." MLOGIT-MENU-ITEM :FIT-MLOGIT RC-MENU-ITEM) 0 NIL) #K(#(80 1 3 0 1 2 2 38 1 4 23 3 4 2 3 4 21 0 1 2 4 11 0 1 23 4 0 1 2 38 1 4 23 5 5 1 3 23 6 4 2 4 5 23 7 0 2 2 4 11 0 2) #(0 21 24) #(NIL *ARC-FIT-MENU-ITEMS* MLOGIT-MENU-ITEM POSITION MARK-AS-SPECIAL LIST APPEND SET) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Arc fit-method for multinomial logit models." :FIT-MLOGIT #K(#(80 2 3 0 1 2 4 38 3 5 23 4 0 2 1 5 38 5 5 23 4 5 2 1 5 23 6 5 3 4 3 5 38 7 6 23 8 7 1 2 23 4 6 3 1 6 7 38 9 7 38 10 8 23 11 5 5 6 7 5 8 6 21 0 1 2 5 39 4 6 38 12 7 25 13 0 4 6 1 7 5 11 0 4) #(0 67 81) #((:ONES :VARIATE :FACTOR :FACTOR1 :EFFECT :INTERACTION :TEXT) "M~a" NIL :ONES SEND :MODEL-NUMBER FORMAT :TYPES APPEND :NAME :CANDIDATES MLOGIT-DIALOG :MAKE-MLOGIT APPLY) 0 NIL) DATASET-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:FIT-MLOGIT) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function args: (names &key reg-name type) (mod-type \"Multinomial Logit Model\") (name \"mlogit\") (weights \"\") (response \"\") (candidates names) (predictors (repeat \"\" (length names))) (intercept t) (compute-info-matrix t)) Opens an Arc dialog to define a multinomial logistic regression." FUNCTION-DOCUMENTATION MLOGIT-DIALOG %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(MLOGIT-DIALOG #K(#(77 3 2 0 0 0 18 17 17 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 24 21 2 1 2 6 0 6 27 1 3 0 1 27 21 2 4 5 7 0 7 28 1 6 23 19 28 1 1 23 20 28 2 10 28 38 21 29 38 22 30 38 23 31 23 24 31 3 26 25 31 23 25 29 3 29 30 31 38 26 30 38 22 31 23 25 30 2 30 31 38 27 31 38 22 32 23 19 33 1 1 23 28 33 2 27 33 23 25 31 5 31 32 33 24 30 38 21 32 38 22 33 23 25 32 3 32 33 23 38 27 33 38 22 34 23 19 35 1 1 23 28 35 2 28 35 23 25 33 5 33 34 35 22 30 38 21 34 38 22 35 23 25 34 3 34 35 21 38 27 35 38 22 36 23 29 37 1 5 23 25 35 5 35 36 37 20 30 38 21 36 38 22 37 23 25 36 3 36 37 19 38 27 37 38 22 38 21 0 7 8 4 23 29 39 1 4 1 9 0 11 39 23 25 37 5 37 38 39 20 30 38 21 38 38 22 39 23 25 38 3 38 39 18 38 30 39 38 22 40 38 31 41 23 25 39 5 39 40 17 41 8 38 30 40 38 22 41 38 31 42 23 25 40 5 40 41 16 42 9 38 32 41 38 22 42 38 33 43 23 25 41 5 41 42 3 43 15 38 21 42 38 22 43 23 24 44 3 26 14 2 23 25 42 3 42 43 44 38 34 43 38 22 44 23 25 43 3 43 44 13 38 34 44 38 22 45 38 35 46 51 10 47 6 2 1 40 39 41 30 23 25 44 5 44 45 12 46 47 38 36 45 38 22 46 23 29 47 3 29 42 41 23 29 48 4 32 31 36 38 23 29 49 4 34 33 35 37 23 29 50 4 39 40 44 43 23 29 48 3 48 49 50 23 29 47 2 47 48 23 25 45 3 45 46 47 38 37 46 23 25 0 3 45 46 44 38 38 46 25 25 0 2 45 46 80 7 5 39 40 41 42 43 17 38 31 12 23 25 12 2 6 12 38 44 13 23 25 13 2 6 13 23 45 14 2 12 11 21 0 11 14 14 23 45 15 2 12 8 21 0 12 14 15 23 45 15 2 12 7 21 0 13 14 15 38 46 15 23 45 16 2 12 11 38 47 17 23 45 18 2 12 9 38 48 19 23 45 20 2 12 8 16 20 20 38 49 21 23 45 22 2 12 7 16 22 22 38 50 23 38 31 24 23 25 24 2 4 24 38 51 25 38 31 26 23 25 26 2 3 26 38 52 27 38 53 28 23 25 28 2 5 28 25 29 0 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 23 54 0 1 10 38 55 14 38 52 15 38 53 16 23 25 16 2 5 16 38 47 17 23 45 18 2 13 9 38 46 19 23 45 20 2 13 11 38 48 21 23 45 22 2 13 8 16 22 22 38 49 23 23 45 24 2 13 7 16 24 24 38 50 25 38 31 26 23 25 26 2 4 26 38 51 27 38 31 28 23 25 28 2 3 28 25 56 0 17 2 14 1 15 16 17 18 19 20 21 22 23 24 25 26 27 28) #(0 32 37 40 45 50 61 218 225 228 441 479 490 501 597) #("" ("") "Done" "Cancel" "Name for ~%~a" 14 "Calculate Info Matrix" "Fit Intercept" "Case Weights... " "Response... " ONEITEM "Predictors" SECONDARY "Candidates" PRIMARY "Arc ~a" NIL ("Multinomial Logit Model" "mlogit" "" "" %NOT-SUPPLIED %NOT-SUPPLIED T T) (:MOD-TYPE :NAME :WEIGHTS :RESPONSE :CANDIDATES :PREDICTORS :INTERCEPT :COMPUTE-INFO-MATRIX) LENGTH REPEAT TEXT-ITEM-PROTO :NEW *ARC-VERSION* FORMAT SEND ALLOCATE-ITEMS-PROTO SUBORDINATE-LIST-PROTO PAD-LIST LIST TOGGLE-ITEM-PROTO :VALUE EDIT-TEXT-ITEM-PROTO :TEXT-LENGTH MODAL-BUTTON-PROTO :ACTION MODAL-DIALOG-PROTO :DEFAULT-BUTTON :MODAL-DIALOG 3 2 0 "Predictor, response, or case weights missing" 1 :ALL-VALUES SELECT :PREDICTORS :CANDIDATES :RESPONSE :WEIGHTS :INTERCEPT :COMPUTE-INFO-MATRIX :NAME :TEXT MESSAGE-DIALOG :MOD-TYPE MLOGIT-DIALOG) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (&rest args &key predictors response intercept compute-info-matrix weights name) Sets up multinomial logit model from Arc." :MAKE-MLOGIT #K(#(77 3 2 0 1 0 6 5 5 0 1 2 3 4 8 39 7 14 51 7 15 1 1 23 8 15 2 15 3 23 9 14 2 14 15 38 0 15 23 10 15 3 1 15 4 16 15 15 39 11 16 20 15 14 17 23 8 16 2 16 17 23 12 16 1 16 23 13 16 1 16 23 14 17 1 16 23 15 16 2 16 17 38 0 17 23 10 17 3 1 17 7 16 17 17 38 16 18 38 17 19 38 18 20 23 10 18 4 18 19 20 8 23 19 19 1 8 23 20 19 1 19 23 21 0 2 19 18 38 22 19 23 10 19 2 18 19 38 23 20 23 10 0 2 19 20 38 22 19 23 10 19 2 18 19 38 24 20 23 10 0 3 19 20 8 38 22 19 23 10 19 2 18 19 38 25 20 23 10 0 2 19 20 38 26 19 23 10 0 3 1 19 18 38 27 19 23 10 0 3 18 19 14 38 28 19 23 10 0 3 18 19 15 38 29 19 23 10 0 3 18 19 7 38 30 19 23 10 0 3 18 19 17 38 31 19 23 10 0 3 18 19 16 38 32 19 23 10 0 3 18 19 5 38 33 19 23 10 0 3 18 19 6 38 34 19 51 6 20 1 1 23 8 20 2 20 3 51 5 21 0 23 8 20 2 21 20 23 12 20 1 20 23 10 0 3 18 19 20 38 35 19 23 10 0 3 18 19 4 38 36 19 38 37 20 23 10 0 3 18 19 20 38 38 19 23 10 0 4 18 19 13 2 38 38 19 23 10 0 4 18 19 12 1 51 4 19 1 1 23 39 0 4 18 11 19 10 51 3 19 0 23 39 0 4 18 9 19 10 38 40 19 23 10 0 2 18 19 38 41 19 21 0 1 2 19 38 42 20 23 10 0 2 18 20 11 0 18 80 2 1 3 0 62 2 0 57 0 77 2 3 0 0 2 39 10 4 38 2 5 25 9 0 4 4 1 5 3 80 2 0 1 38 43 2 25 10 0 2 1 2 80 3 0 1 38 44 3 25 10 0 3 1 3 2 80 3 0 1 38 0 3 25 10 0 3 1 3 2) #(0 386 395 398 408 428 441 455) #(:DATA NIL :INCLUDED DATA ARGS (NIL NIL NIL NIL NIL NIL) (:PREDICTORS :RESPONSE :INTERCEPT :COMPUTE-INFO-MATRIX :WEIGHTS :NAME) APPEND MAPCAR APPLY SEND FIND-MISSING COMBINE REMOVE-DUPLICATES WHICH SELECT MLOGIT-MODEL-PROTO :NEW :NAME STRING-UPCASE INTERN SET :MENU :REMOVE :TITLE :INSTALL :MODELS :X :Y :WEIGHTS-NAME :WEIGHTS :MISSING :INTERCEPT :COMPUTE-INFO-MATRIX :PREDICTOR-NAMES :RESPONSE-NAME :TYPE :MLOGIT-MODEL :ADD-SLOT XLISP::ADD-METHOD :COMPUTE *SHOW-DISPLAY-FIT* :SUMMARY :LABELS :FIND-DATALIST) 0 NIL) DATASET-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:MAKE-MLOGIT) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 0 25 3 0 3 3 2 1) #(0) #("Function args: (data included) Filters list, vector, or matrix data for logical list of included observations." FUNCTION-DOCUMENTATION INCLUDED %SET-GET) 0 NIL) #K(#(80 1 2 0 1 1 23 2 3 1 2 25 3 0 2 1 3) #(0) #(INCLUDED #K(#(80 3 2 0 1 2 23 2 5 1 2 23 3 6 1 1 21 0 1 2 6 25 4 0 2 1 5 23 5 6 1 1 21 0 3 4 6 25 4 0 2 1 5 23 6 6 1 1 21 0 5 6 6 23 7 6 2 1 3 23 8 6 1 6 25 4 0 3 1 5 6 25 9 0 1 4) #(0 21 27 37 43 53 71) #(1 "~%Not a list, vector, or matrix~%" WHICH LISTP SELECT VECTORP MATRIXP ARRAY-DIMENSION ISEQ ERROR) 0 NIL) BYTE-CODE-CLOSE XLISP::INSTALL-FUNCTION) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns number of included observations." :NUM-INCLUDED #K(#(80 2 1 0 2 38 1 3 23 2 3 2 1 3 23 3 3 1 3 38 4 4 23 2 4 2 1 4 23 5 0 3 3 4 2 23 6 3 1 3 25 7 0 1 3) #(0) #(NIL :INCLUDED SEND COPY-LIST :MISSING SET-SELECT WHICH LENGTH) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:NUM-INCLUDED) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :MENU #K(#(77 3 2 0 0 0 3 2 2 0 1 3 62 4 5 21 0 1 2 2 0 2 6 1 3 0 3 6 21 0 4 5 5 11 0 5 38 4 7 23 5 5 3 1 7 6 63 4 5 0 57 0) #(0 20 25 28 33 36) #("MLogit" MENU (NIL) (:TITLE) :MAKE-MENU SEND) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:MENU) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Arc menu items for multinomial logit model." :MENU-ITEMS #K(#(80 2 0 2 38 0 2 51 1 3 1 1 25 1 0 2 3 2 80 3 0 0 25 2 0 2 2 1) #(0 18) #(*ARC-MULTINOMIAL-LOGIT-MODEL-MENU-ITEMS* MAPCAR MAKE-MENU-ITEM) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:MENU-ITEMS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Arc list of model-related variables for multinomial logit model." :PLOTLIST #K(#(80 2 0 1 38 0 2 25 1 0 1 2) #(0) #(:DEV-COMPONENTS LIST) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:PLOTLIST) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #(NIL :REMOVE-MENU #K(#(80 2 0 2 38 0 2 23 1 2 2 1 2 38 2 3 23 1 0 2 2 3 38 3 2 23 1 2 2 1 2 38 4 3 25 1 0 3 2 3 1) #(0) #(:MENU SEND :REMOVE :DATA :DELETE-MODEL) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:REMOVE-MENU) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Summary" MLOGIT-SUMMARY-ITEM :SUMMARY RC-MENU-ITEM) 0 NIL) #K(#(80 1 2 0 1 1 38 2 3 25 3 0 3 2 1 3) #(0) #("Remove model" MLOGIT-REMOVE-ITEM :REMOVE-MENU RC-MENU-ITEM) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #(*ARC-MULTINOMIAL-LOGIT-MODEL-MENU-ITEMS* MARK-AS-SPECIAL) 0 NIL) #K(#(80 1 4 0 1 2 3 1 23 4 5 3 4 3 2 25 5 0 2 1 5) #(0) #(*ARC-MULTINOMIAL-LOGIT-MODEL-MENU-ITEMS* MLOGIT-REMOVE-ITEM WALD-TESTS-MENU-ITEM MLOGIT-SUMMARY-ITEM LIST SET) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(*ARC-MULTINOMIAL-LOGIT-MODEL-MENU-ITEMS*) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Returns a list of case-wise components of the deviance for the model." :DEV-COMPONENTS #K(#(80 2 4 0 1 2 3 19 38 4 6 23 5 6 2 1 6 38 6 7 23 5 7 2 1 7 38 7 8 23 5 8 2 1 8 38 8 9 23 5 9 2 1 9 0 5 10 23 9 11 2 8 4 23 9 12 2 8 3 23 9 13 2 7 3 23 10 14 2 13 12 23 11 14 2 14 6 23 12 14 1 14 23 13 15 1 12 0 4 16 37 60 1 6 16 11 23 14 17 2 9 16 0 4 18 0 4 19 37 60 2 3 19 13 68 7 16 19 20 23 14 21 3 8 16 15 23 14 22 3 14 15 19 23 15 21 2 21 22 2 42 20 21 20 2 43 18 20 18 61 112 19 19 37 60 2 3 19 13 0 18 20 0 4 21 0 4 22 37 60 4 5 22 13 23 14 23 3 8 16 15 23 14 24 3 14 15 22 23 15 23 2 23 24 23 16 23 1 23 2 43 21 23 21 61 112 22 22 37 60 4 5 22 13 2 43 3 21 23 23 17 23 1 23 61 45 23 23 2 43 20 23 23 23 18 23 3 2 17 23 68 23 4 4 23 20 23 10 10 61 112 16 16 37 60 1 6 16 11 38 19 17 25 5 0 3 1 17 10) #(0 96 114 159 174 214 259) #(-2 1 0 NIL :COEF-ESTIMATES SEND :W :X-MATRIX :WEIGHTS ARRAY-DIMENSION LIST MATRIX TRANSPOSE ISEQ SELECT MATMULT EXP LOG * :EXPAND-TO-MISSING) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:DEV-COMPONENTS) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: (x) Pads x with missing data corresponding to excluded observations." :EXPAND-TO-MISSING #K(#(80 3 0 3 38 0 3 23 1 3 2 1 3 38 2 4 23 3 5 1 3 23 4 4 2 4 5 23 5 5 1 3 23 6 0 3 4 5 2 11 0 4) #(0) #(:INCLUDED SEND NOT-A-NUMBER LENGTH REPEAT WHICH SET-SELECT) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:EXPAND-TO-MISSING) 0 NIL)) EVAL) 0 NIL) #K(#(80 1 1 0 0 25 1 0 1 1) #(0) #((PROGN #K(#(80 1 3 0 1 2 2 38 3 4 23 4 5 1 3 25 5 0 4 4 2 5 1) #(0) #("Message args: () Prints a Wald chisquare test for each term in the model, testing the hypothesis that the term (i.e., set of coefficients) is zero." :WALD-TESTS #K(#(80 2 12 0 1 2 3 4 5 6 7 8 9 10 11 14 38 12 14 23 13 14 2 1 14 38 14 15 23 13 14 2 14 15 38 15 15 23 13 15 2 1 15 23 16 16 1 15 51 10 17 1 1 23 17 17 2 17 15 38 18 18 23 13 18 2 1 18 23 16 18 1 18 2 45 18 13 18 23 19 19 1 16 23 20 19 2 19 17 38 21 20 23 13 20 2 1 20 21 0 1 2 20 23 22 19 2 2 19 23 16 20 1 19 23 23 20 1 20 23 16 21 1 19 23 19 21 1 21 38 21 22 23 13 22 2 1 22 21 0 3 4 22 38 24 23 23 25 24 3 20 7 21 27 15 24 2 24 18 38 26 25 23 13 22 5 1 23 24 25 12 1 5 0 6 22 38 14 23 23 13 23 2 1 23 23 27 0 3 12 11 23 23 27 0 3 12 10 14 23 27 0 2 12 9 23 27 0 2 12 8 38 21 23 23 13 23 2 1 23 21 0 6 7 23 23 25 24 2 22 13 23 25 25 2 22 7 23 25 26 2 22 5 23 27 0 6 12 4 24 25 26 3 0 7 22 37 60 8 9 22 16 23 25 0 2 17 22 3 61 22 19 23 23 28 23 1 23 23 25 23 3 20 23 21 27 15 23 2 23 18 38 24 24 38 26 25 23 13 23 5 1 24 23 25 12 23 25 24 2 23 13 23 25 25 2 23 7 23 25 26 2 23 5 23 25 27 2 15 22 23 27 0 6 12 4 24 25 26 27 61 112 22 22 37 60 8 9 22 16 11 0 6 80 3 0 1 38 29 3 25 13 0 3 1 3 2 80 5 0 0 37 61 12 13 4 3 11 0 1 11 0 2 80 5 0 3 39 30 5 23 19 6 1 3 51 11 7 3 1 2 4 23 17 6 2 7 6 25 31 0 2 5 6 80 3 2 5 11 4 23 32 5 2 1 4 23 32 6 2 1 3 23 33 7 2 6 5 2 42 6 5 8 23 20 8 2 3 8 23 34 5 2 7 8 39 35 6 23 19 7 1 2 51 14 8 3 1 5 2 23 17 7 2 8 7 25 31 0 2 6 7) #(0 103 109 143 173 176 225 253 262 350 353 367 377 380 383 414) #(-1 "Constant" "~8d ~10g ~10,5f ~a~%" 2 NIL 0 "------------------------------------------------~%" "~% df Wald Chisq. p Source~%" ", Dataset = ~a" "~%~% Wald-Tests for terms in Model ~a" T 1 :GET-DATASET SEND :NAME :PREDICTORS LENGTH MAPCAR :RESPONSE-CATEGORIES ISEQ REPEAT :INTERCEPT ADJOIN IDENTITY-MATRIX :TEST-LINEAR-HYPOTHESIS SELECT :QUIET FORMAT WHICH :DEGREES-OF-FREEDOM-TERM BIND-COLUMNS APPLY ARRAY-DIMENSION LIST MATRIX BIND-ROWS) 0 NIL) MLOGIT-MODEL-PROTO BYTE-CODE-CLOSE XLISP::ADD-METHOD) 0 NIL) #K(#(80 1 1 0 0 11 0 1) #(0) #(:WALD-TESTS) 0 NIL)) EVAL) 0 NIL)