3 (asdf:operate
'asdf
:compile-op
'cffi
)
4 (asdf:operate
'asdf
:compile-op
'lispstat
)
6 ;; (asdf:operate 'asdf:load-op 'cffi)
7 ;; (asdf:operate 'asdf:load-op 'rclg)
9 (asdf:oos
'asdf
:load-op
:lispstat
)
12 (setf *my-base-directory
*
13 #p
"/home/tony/sandbox/CLS.git/"
14 #p
"/Users/ungil/lisp/CommonLispStat/"
17 (load "/Users/ungil/lisp/CommonLispStat/init.lisp")
18 (load "/Users/ungil/lisp/CommonLispStat/cffiglue.lsp")
23 (lisp-stat-basics:binomial-quant
0.95 3 0.4) ;;; DOESN'T RETURN
24 (lisp-stat-basics:binomial-quant
0 3 0.4) ;;; -2147483648
28 (binomial-quant 0.95 3 0.4) ;;; DOESN'T RETURN
29 (binomial-quant 0 3 0.4) ;;; -2147483648
33 ;; optimization-related functions not tested
34 ;;(lisp-stat-optimize:numgrad
35 ;;(lisp-stat-optimize:numhess
36 ;;(lisp-stat-optimize::minfo-maximize
38 ;; the rest return reasonable values, except BINOMIAL-QUANT:
39 ;; (lisp-stat-basics:binomial-quant 0.95 3 0.4) ;;; DOESN'T RETURN
40 ;; (lisp-stat-basics:binomial-quant 0 3 0.4) ;;; -2147483648
42 (lisp-stat-linalg:chol-decomp
#2A
((2 3 4) (1 2 4) (2 4 5)))
43 ;; (#2A((1.7888543819998317 0.0 0.0)
44 ;; (1.6770509831248424 0.11180339887498929 0.0)
45 ;; (2.23606797749979 2.23606797749979 3.332000937312528e-8))
48 (lisp-stat-linalg:lu-decomp
#2A
((2 3 4) (1 2 4) (2 4 5)))
49 ;; (#2A((2.0 3.0 4.0) (1.0 1.0 1.0) (0.5 0.5 1.5)) #(0 2 2) -1.0 NIL)
51 (lisp-stat-linalg:lu-solve
52 (lisp-stat-linalg:lu-decomp
#2A
((2 3 4) (1 2 4) (2 4 5)))
54 ;; #(-2.333333333333333 1.3333333333333335 0.6666666666666666)
56 (lisp-stat-linalg:inverse
#2A
((2 3 4) (1 2 4) (2 4 5)))
57 ;; #2A((2.0 -0.33333333333333326 -1.3333333333333335)
58 ;; (-1.0 -0.6666666666666666 1.3333333333333333)
59 ;; (0.0 0.6666666666666666 -0.3333333333333333))
61 (lisp-stat-linalg:sv-decomp
#2A
((2 3 4) (1 2 4) (2 4 5)))
62 ;; (#2A((-0.5536537653489974 0.34181191712789266 -0.7593629708013371)
63 ;; (-0.4653437312661058 -0.8832095891230851 -0.05827549615722014)
64 ;; (-0.6905959164998124 0.3211003503429828 0.6480523475178517))
65 ;; #(9.699290438141343 0.8971681569301373 0.3447525123483081)
66 ;; #2A((-0.30454218417339873 0.49334669582252344 -0.8147779426198863)
67 ;; (-0.5520024849987308 0.6057035911404464 0.5730762743603965)
68 ;; (-0.7762392122368734 -0.6242853493399995 -0.08786630745236332))
71 (lisp-stat-linalg:qr-decomp
#2A
((2 3 4) (1 2 4) (2 4 5)))
72 ;; (#2A((-0.6666666666666665 0.7453559924999298 5.551115123125783e-17)
73 ;; (-0.3333333333333333 -0.2981423969999719 -0.894427190999916)
74 ;; (-0.6666666666666666 -0.5962847939999439 0.44721359549995787))
75 ;; #2A((-3.0 -5.333333333333334 -7.333333333333332)
76 ;; (0.0 -0.7453559924999292 -1.1925695879998877)
77 ;; (0.0 0.0 -1.3416407864998738)))
79 (lisp-stat-linalg:rcondest
#2A
((2 3 4) (1 2 4) (2 4 5)))
82 (lisp-stat-linalg:eigen
#2A
((2 3 4) (1 2 4) (2 4 5)))
83 ;; (#(10.656854249492381 -0.6568542494923802 -0.9999999999999996)
84 ;; (#(0.4999999999999998 0.4999999999999997 0.7071067811865475)
85 ;; #(-0.49999999999999856 -0.5000000000000011 0.7071067811865474)
86 ;; #(0.7071067811865483 -0.7071067811865466 -1.2560739669470215e-15))
89 (lisp-stat-linalg:spline
#(1.0
1.2 1.3 1.8 2.1 2.5)
90 #(1.2
2.0 2.1 2.0 1.1 2.8) :xvals
6)
91 ;; ((1.0 1.3 1.6 1.9 2.2 2.5)
92 ;; (1.2 2.1 2.2750696543866313 1.6465231041904045 1.2186576148879609 2.8))
94 ;;; using KERNEL-SMOOTH-FRONT, not KERNEL-SMOOTH-CPORT
95 (lisp-stat-linalg:kernel-smooth
#(1.0
1.2 1.3 1.8 2.1 2.5)
96 #(1.2
2.0 2.1 2.0 1.1 2.8) :xvals
5)
97 ;; ((1.0 1.375 1.75 2.125 2.5)
98 ;; (1.6603277642110226 1.9471748095239771 1.7938127405752287
99 ;; 1.5871511322219498 2.518194783156392))
101 (lisp-stat-linalg:kernel-dens
#(1.0
1.2 2.5 2.1 1.8 1.2) :xvals
5)
102 ;; ((1.0 1.375 1.75 2.125 2.5)
103 ;; (0.7224150453621405 0.5820045548233707 0.38216411702854214
104 ;; 0.4829822708587095 0.3485939156929503))
106 (lisp-stat-linalg:fft
#(1.0
1.2 2.5 2.1 1.8))
107 ;; #(#C(1.0 0.0) #C(1.2 0.0) #C(2.5 0.0) #C(2.1 0.0) #C(1.8 0.0))
109 (lisp-stat-descriptive-statistics:lowess
110 #(1.0
1.2 2.5 2.1 1.8 1.2) #(1.2
2.0 2.1 2.0 1.1 2.8))
111 ;; (#(1.0 1.2 1.2 1.8 2.1 2.5))
113 ;;;; Log-gamma function
115 (lisp-stat-basics:log-gamma
3.4) ;;1.0923280596789584
117 ;;;; Normal distribution
119 (lisp-stat-basics:normal-quant
0.95) ;;1.6448536279366268
120 (lisp-stat-basics:normal-cdf
1.3) ;;0.9031995154143897
121 (lisp-stat-basics:normal-dens
1.3) ;;0.17136859204780736
122 (lisp-stat-basics:normal-rand
2) ;;(-0.40502015f0 -0.8091404f0)
124 (lisp-stat-basics:bivnorm-cdf
0.2 0.4 0.6) ;;0.4736873734160288
126 ;;;; Cauchy distribution
128 (lisp-stat-basics:cauchy-quant
0.95) ;;6.313751514675031
129 (lisp-stat-basics:cauchy-cdf
1.3) ;;0.7912855998398473
130 (lisp-stat-basics:cauchy-dens
1.3) ;;0.1183308127104695
131 (lisp-stat-basics:cauchy-rand
2) ;;(-1.06224644160405 -0.4524695943939537)
133 ;;;; Gamma distribution
135 (lisp-stat-basics:gamma-quant
0.95 4.3) ;;8.178692439291645
136 (lisp-stat-basics:gamma-cdf
1.3 4.3) ;;0.028895150986674906
137 (lisp-stat-basics:gamma-dens
1.3 4.3) ;;0.0731517686447374
138 (lisp-stat-basics:gamma-rand
2 4.3) ;;(2.454918912880936 4.081365384357454)
140 ;;;; Chi-square distribution
142 (lisp-stat-basics:chisq-quant
0.95 3) ;;7.814727903379012
143 (lisp-stat-basics:chisq-cdf
1 5) ;;0.03743422675631789
144 (lisp-stat-basics:chisq-dens
1 5) ;;0.08065690818083521
145 (lisp-stat-basics:chisq-rand
2 4) ;;(1.968535826180572 2.9988646156942997)
147 ;;;; Beta distribution
149 (lisp-stat-basics:beta-quant
0.95 3 2) ;;0.9023885371149876
150 (lisp-stat-basics:beta-cdf
0.4 2 2.4) ;;0.4247997418541529
151 (lisp-stat-basics:beta-dens
0.4 2 2.4) ;;1.5964741858913518
152 (lisp-stat-basics:beta-rand
2 2 2.4) ;;(0.8014897077282279 0.6516371997922659)
156 (lisp-stat-basics:t-quant
0.95 3) ;;2.35336343484194
157 (lisp-stat-basics:t-cdf
1 2.3) ;;0.794733624298342
158 (lisp-stat-basics:t-dens
1 2.3) ;;0.1978163816318102
159 (lisp-stat-basics:t-rand
2 2.3) ;;(-0.34303672776089306 -1.142505872436518)
163 (lisp-stat-basics:f-quant
0.95 3 5) ;;5.409451318117459
164 (lisp-stat-basics:f-cdf
1 3.2 5.4) ;;0.5347130905510765
165 (lisp-stat-basics:f-dens
1 3.2 5.4) ;;0.37551128864591415
166 (lisp-stat-basics:f-rand
2 3 2) ;;(0.7939093442091963 0.07442694152491144)
168 ;;;; Poisson distribution
170 (lisp-stat-basics:poisson-quant
0.95 3.2) ;;6
171 (lisp-stat-basics:poisson-cdf
1 3.2) ;;0.17120125672252395
172 (lisp-stat-basics:poisson-pmf
1 3.2) ;;0.13043905274097067
173 (lisp-stat-basics:poisson-rand
5 3.2) ;;(2 1 2 0 3)
175 ;;;; Binomial distribution
177 (lisp-stat-basics:binomial-quant
0.95 3 0.4) ;;; DOESN'T RETURN
178 (lisp-stat-basics:binomial-quant
0 3 0.4) ;;; -2147483648
179 (lisp-stat-basics:binomial-cdf
1 3 0.4) ;;0.6479999999965776
180 (lisp-stat-basics:binomial-pmf
1 3 0.4) ;;0.4320000000226171
181 (lisp-stat-basics:binomial-rand
5 3 0.4) ;;(2 2 0 1 2)