;ELC
;;; Compiled
;;; in Emacs version 28.2
;;; with all optimizations.



(byte-code "\300\301!\210\300\302!\207" [require calc-arith calcalg3] 2)
#@111 Return the parameters A and B for the best least squares fit y=a+bx.

(fn XDATA YDATA &optional SDATA SIGMAS)
(defalias 'math-nlfit-least-squares #[1026 "G\203
\300\301\"\202\302
\303\"\203\304\202
\304\211\211\211\305\203\246@@@\306\203;\307\"\202<\"\262\306\203M\307\"\202N\"\262\306\203c\307\310\211\"\"\202g\310\211\"\"\262\306\203|\307\310\"\"\202\200\310\"\"\262\f\203\222\306	\307\303\"\"\262	\266A\262\f\nA\262A\262\202\311\310\"\310\211\"\"\262
\307\311\310\"\310\"\"\"\307\311\310	\"\310		\"\"\"\n\203
\307\"\312\307	\"!\307\n\"\313\314!E\313\314!E\315\211E\315EEE\207

D\207" [mapcar calcFunc-sqr make-list 1 0 nil math-add math-div math-mul math-sub math-neg sdev calcFunc-sqrt vec] 24 (#$ . 158)])
#@20 

(fn TDATA QDATA)
(defalias 'math-nlfit-get-rates-from-cumul #[514 "\300\301A@@\"\301A@@\"\"C
G\302V\203R\303\304\305\300\301\3028A@\"\301\302	8	A@\"\"\300\301A@@\"\301	A@\n@\"\"\"\"
B\262

A\262\202\300\301A@@\"\301A@@\"\"
B\262
\306
!\207" [math-div math-sub 2 math-mul (float 5 -1) math-add reverse] 12 (#$ . 990)])
#@23 

(fn TDATA PDATA Q0)
(defalias 'math-nlfit-get-cumul-from-rates #[771 "\211CA\2031\300
@\301\211\302\300A@	@\"\"\303A@	@\"\"\"
B\262
A\262A\262\202\304
!\207" [math-add math-mul (float 5 -1) math-sub reverse] 12 (#$ . 1341)])
#@26 

(fn QDATA PDATA TDATA)
(defalias 'math-nlfit-find-qmax #[771 "\300\301#\302
\"\303@A\"\304\301A@@\"!\305
\"\203%\306\307\"\207\207" [math-map-binop math-div math-nlfit-least-squares math-max-list math-neg math-lessp math-add (float 5 -1)] 10 (#$ . 1588)])
#@26 

(fn QDATA PDATA TDATA)
(defalias 'math-nlfit-find-logistic-parameters #[771 "\300#\301\302\303\"\"\304\"\305@!A@E\207" [math-nlfit-find-qmax mapcar make-closure #[257 "\301\302\303\300\"\304\"!\207" [V0 calcFunc-ln math-sub math-div 1] 6 "\n\n(fn Q)"] math-nlfit-least-squares calcFunc-exp] 9 (#$ . 1863)])
#@20 

(fn PDATA TDATA)
(defalias 'math-nlfit-find-qmaxhalf #[514 "\300@A\"\301\302@\"\2036\303
\304\211\305\303	A@\n@\"\"\306A@@\"\"\"\262
A\262A\262\202\207" [math-max-list 0 math-lessp math-add math-mul (float 5 -1) math-sub] 12 (#$ . 2189)])
#@20 

(fn PDATA TDATA)
(defalias 'math-nlfit-find-q0 #[514 "\300\"\301\302\"C\303\242#\304\305\306\307\310\"\"#\301\311\312\242\"\"\"\2037
\312\242\"\240\210\202\313\242\"\242\305\306\307\314\"\"#\210\315\316\211\317W\203\221\301\320\312\"\"\240\210\304\305\306\307\321\n\"\"#\301\322\312\242\"\"\"\203\206\242\262\202\212\242\262\211T\262
\202O\301\323\312\"\"\207" [math-nlfit-find-qmaxhalf math-mul 2 math-nlfit-get-cumul-from-rates math-lessp math-nlfit-find-qmax mapcar make-closure #[257 "\301\300\242\"\207" [V0 math-add] 4 "\n\n(fn Q)"] (float 5 -1) math-add math-sub #[257 "\301\300\242\"\207" [V0 math-add] 4 "\n\n(fn Q)"] nil 0 10 (float 5 -1) #[257 "\301\300\242\"\207" [V0 math-add] 4 "\n\n(fn Q)"] (float 5 -1) (float 5 -1)] 16 (#$ . 2453)])
(defvar math-nlfit-delta '(float 1 -8))
(defvar math-nlfit-epsilon '(float 1 -5))
(defvar math-nlfit-max-its 100)
#@16 

(fn VEC I X)
(defalias 'math-nlfit-set-elt #[771 "
S\233
\240\207" [] 5 (#$ . 3388)])
#@14 

(fn VEC I)
(defalias 'math-nlfit-get-elt #[514 "\211S8\207" [] 4 (#$ . 3483)])
#@12 

(fn I J)
(defalias 'math-nlfit-make-matrix #[514 "\300
\301\"\302\301\211W\203\303!B\262\211T\262
\202
\207" [make-list 0 nil copy-sequence] 7 (#$ . 3571)])
#@18 

(fn MAT I J X)
(defalias 'math-nlfit-set-matx-elt #[1028 "
SS8\233
\240\207" [] 7 (#$ . 3744)])
#@16 

(fn MAT I J)
(defalias 'math-nlfit-get-matx-elt #[771 "\211SS88\207" [] 6 (#$ . 3850)])
#@12 

(fn C D)
(defalias 'math-nlfit-givens #[514 "\301!\301!
@GG\302\303\304\"\303\304\"\305\211\211\211	X\203W
T\211	X\203N
\306\f\n#\306
\n\211#\307\304\"\204E
\310\311!\312\311!\"\"\203_\313!\262\304\262\302\262\302\262\202\236\312\314!\315\316\312\211\"\312\211\"\"!\"\262\317
\"\262\313\317	\"!\262\310\311!\"\203\224\262\202\236\317\314!\"\262\262
\262\320\211\n$\210\320\f$\210	T\211
X\203\n
\306\f#\306#\321\312\n\"\312\n\"\"\316\312\n\"\312\f\"\"\262\211\262\320$\210\320$\210T\266\204\202\274\210\322
\"\322\f\"\321\312	\"\312	\"\"\316\312	\"\312\"\"\262\211\262\323#\210\323#\266\266\211T\262
\202#\210T\262\202\305
\302Y\203\255
\323\304#\210\322\f\"\262

T\211\fX\203\223
\316\312\306#\322\f\"\"\"\262\211T\262
\202p
\210\323\313\317\306	\211#\"!#\210
S\262\202Z
\266T\211	X\203\313
\323\322\"#\210\211T\262
\202\262
\210\211\302Y\203s\211TY\203k\306
#\262\307\302\"\203\365
\304\262\302\262\202)\310\311!\302\"\203\262\315\321\302\312\211\"\"!\262\202)\317\302\311!\"\262\312\314!\315\321\302\312\n\211\"\"!\"\262\322\"\322	\"\305\316\312	\"\312	\"\"\262
\321\312	\"\312	\"\"\262\211\262\323#\210\323#\266\211S\262
\202\326
\210\211S\262
\202\316
\210\207" [math-nlfit-delta copy-tree 1 make-list 0 nil math-nlfit-get-matx-elt math-equal math-lessp calcFunc-abs math-mul math-neg calcFunc-sign calcFunc-sqrt math-add math-div math-nlfit-set-matx-elt math-sub math-nlfit-get-elt math-nlfit-set-elt] 25 (#$ . 3948)])
#@41 

(fn GRAD XLIST PARMS &optional SLIST)
(defalias 'math-nlfit-jacobian #[1027 "\211C\300\2031\301@#\242\203\302\303\304\"\"\202 \211B\266\202
\211\242A\240\210A\262\202\305
!\207" [nil apply mapcar make-closure #[257 "\301
\300\242@\"\207" [V0 math-div] 4 "\n\n(fn X)"] reverse] 11 (#$ . 5613)])
#@12 

(fn L N)
(defalias 'math-nlfit-make-ident #[514 "\300
\211\"\301\211X\203\302\211$\210\211T\262
\202
\207" [math-nlfit-make-matrix 1 math-nlfit-set-matx-elt] 9 (#$ . 5934)])
#@45 

(fn XLIST YLIST PARMS FN &optional SLIST)
(defalias 'math-nlfit-chi-sq #[1284 "\300\2035\301\302@#@\"\203\303
@\"\262
\304\305\211\"\"\266\202A\262A\262
A\262\202
\207" [0 math-sub apply math-div math-add math-mul] 12 (#$ . 6125)])
#@10 

(fn C)
(defalias 'math-nlfit-init-lambda #[257 "\300
@GG\203*@\211\203\"\301\302@@\"\"\262\211A\262
\202\f\210A\262\202\303\304\302\"\"!\207" [0 math-add math-mul calcFunc-sqrt math-div] 10 (#$ . 6386)])
#@12 

(fn C L)
(defalias 'math-nlfit-make-Ctilda #[514 "
@G\300\"\301\"\207" [math-nlfit-make-ident append] 7 (#$ . 6615)])
#@45 

(fn FN XDATA YDATA PARMS &optional SDATA)
(defalias 'math-nlfit-make-d #[1284 "\300\2034\301\302@#@\"\203\303
@\"\202\211\262

B\262
A\262A\262
A\262\202
\304
!\207" [nil math-sub apply math-div reverse] 11 (#$ . 6745)])
#@12 

(fn D N)
(defalias 'math-nlfit-make-dtilda #[514 "\300\301\302\"\"\207" [append make-list 0] 7 (#$ . 6995)])
#@50 

(fn XLIST YLIST PARMS FN GRAD &optional SLIST)
(defalias 'math-nlfit-fit #[1541 "\302$\303%\304%\305!\306\307
\204\261\211W\203\261\211T\262
\306\211\204\223\310\"\311	@G\"\312\"\313\314\315!#\304%\316
\"\203\207\316\317\320\f\"\"	\"\203t\321\262\317	\322\"\262	\211\262\n
\262\266\321\2023\323	\322\"\262	\266\2023\302	
$\262\303\n\f%\262\210\202$\nD\207" [math-nlfit-max-its math-nlfit-epsilon math-nlfit-jacobian math-nlfit-make-d math-nlfit-chi-sq math-nlfit-init-lambda nil 0 math-nlfit-make-Ctilda math-nlfit-make-dtilda math-nlfit-givens math-map-binop math-add copy-tree math-lessp math-div math-sub t 10 math-mul] 23 (#$ . 7114)])
#@16 

(fn X A B C)
(defalias 'math-nlfit-s-logistic-fn #[1028 "\300\301\302\303\304\303\f\"!\"\"\"\207" [math-div math-add 1 math-mul calcFunc-exp] 14 (#$ . 7846)])
#@16 

(fn X A B C)
(defalias 'math-nlfit-s-logistic-grad #[1028 "\300\301\"!\302\303\301\"\"\301
\211\"\304\303\"\305\304\301	\"\"!\305\304\301\n\301\301\f\"\"\"\"!E\207" [calcFunc-exp math-mul math-add 1 math-div math-neg] 18 (#$ . 8019)])
#@16 

(fn X A C D)
(defalias 'math-nlfit-b-logistic-fn #[1028 "\300\301\302\"\"!\303\301\"\304\305\"\301
\211\"\262
\"\207" [calcFunc-exp math-mul math-sub math-div math-add 1] 11 (#$ . 8279)])
#@16 

(fn X A C D)
(defalias 'math-nlfit-b-logistic-grad #[1028 "\300\301\302\"\"!\303\304\"\302\"\305\301
\211\"\262
\"\302\305\301	\301	\"\"\301
\211\"\262
\"\305\301\306\301\f\301	\f\301
\211\"\262
\"\"\"\307\310\"\"\"\302\305\301\306\301\f\301
\f\301
\211\"\262
\"\"\"\307\310\"\"\305\301\301\f\"\"\301
\211\"\262
\"\"E\207" [calcFunc-exp math-mul math-sub math-add 1 math-div 2 math-pow 3] 21 (#$ . 8482)])
#@26 

(fn GRAD XLIST PPARMS)
(defalias 'math-nlfit-find-covar #[771 "\300\203\301\302@#B
B\262
A\262\202
\301\303!B\262
\304\305!\"\262
\306
!\207" [nil vec apply reverse math-mul calcFunc-trn calcFunc-inv] 9 (#$ . 8928)])
#@32 

(fn GRAD XLIST PPARMS CHISQ)
(defalias 'math-nlfit-get-sigmas #[1028 "\300\301#\211GSG\302
V\2030\211X\203+\303
\21188!B\262\211T\262
\202\304!\262D\207" [nil math-nlfit-find-covar 1 calcFunc-sqrt reverse] 13 (#$ . 9165)])
#@20 

(fn XDATA YDATA)
(defalias 'math-nlfit-s-logistic-params #[514 "\300\"\301#\207" [math-nlfit-get-rates-from-cumul math-nlfit-find-logistic-parameters] 7 (#$ . 9414)])
#@20 

(fn XDATA YDATA)
(defalias 'math-nlfit-b-logistic-params #[514 "\300
\"\301#\302
#\211A@\3038\304\305@\305\"\"!\304\306\307!\"!\306\"\211E\207" [math-nlfit-find-q0 math-nlfit-get-cumul-from-rates math-nlfit-find-logistic-parameters 2 math-neg math-mul math-div calcFunc-ln] 13 (#$ . 9594)])
#@16 

(fn PMS VAR)
(defalias 'math-nlfit-s-logistic-solnexpr #[514 "
@A@\3008\301\302\303\304\305\304\fE!EEE\207" [2 / + 1 * calcFunc-exp] 15 (#$ . 9910)])
#@16 

(fn PMS VAR)
(defalias 'math-nlfit-b-logistic-solnexpr #[514 "
@A@\3008\301\302\303\302\304\nEE!E\305\306\307\303\302\304\f\nEE!E\300EE\207" [2 / * calcFunc-exp - ^ + 1] 16 (#$ . 10076)])
#@22 

(fn N PREFIX VALS)
(defalias 'math-nlfit-enter-result #[771 "
\301#\210\302 \207" [calc-aborted-prefix calc-pop-push-record-list calc-handle-whys] 7 (#$ . 10284)])
#@49 

(fn FN GRAD SOLNEXPR INITPARMS &optional SDV)
(defalias 'math-nlfit-fit-curve #[1284 "\300\301\302&`\"\207" [calc-do make-closure #[0 "\304\306=\206	\304\307=\310\311\312!\211A@A
AA@A\313
!\205!\314\315\"\314\316\"\310\211\"#\312$\317\312\320\"\210\310\"A@#A\303\"\321\300\301&\211A@\205\\\322\301@$\f\205c\211@\211\203o\323\324#\202p\302
	\"\325\310
%	\203\226\326\n@@E
B\262
A\262	A\262\n\202z\327\330\237B!\262
\304\306=\203\255\331\312\332#\210\202\364\304\307=\203\356\310\330	A@\333
@GG\203\335

V\203\335\334@Z\"\266\202\202\340\266\335\257\262
\331\312\336#\266\202\364\331\312\337#\210\340
\341\")\266\202,\207" [V0 V1 V2 V3 V4 calc-display-working-message calcFunc-efit calcFunc-xfit nil calc-top 1 math-contains-sdev-p mapcar #[257 "\300
\301\"\207" [math-get-sdev t] 4 "\n\n(fn X)"] #[257 "\300
!\207" [math-get-value] 3 "\n\n(fn X)"] calc-get-fit-variables 3 math-nlfit-fit math-nlfit-get-sigmas math-map-binop #[514 "\300E\207" [sdev] 5 "\n\n(fn X Y)"] t calcFunc-eq calc-normalize vec math-nlfit-enter-result "efit" (vec) calcFunc-utpc (var nan var-nan) "xfit" "fit" calc-record "parm" calc-curve-varnames calc-curve-coefnames calc-curve-nvars calc-fit-to-trail] 32]] 13 (#$ . 10460)])
#@12 

(fn ARG)
(defalias 'calc-fit-s-shaped-logistic-curve #[257 "\300\301\302\303\304%\207" [math-nlfit-fit-curve math-nlfit-s-logistic-fn math-nlfit-s-logistic-grad math-nlfit-s-logistic-solnexpr math-nlfit-s-logistic-params] 7 (#$ . 11755) "P"])
#@12 

(fn ARG)
(defalias 'calc-fit-bell-shaped-logistic-curve #[257 "\300\301\302\303\304%\207" [math-nlfit-fit-curve math-nlfit-b-logistic-fn math-nlfit-b-logistic-grad math-nlfit-b-logistic-solnexpr math-nlfit-b-logistic-params] 7 (#$ . 12007) "P"])
#@22 

(fn &optional SDV)
(defalias 'calc-fit-hubbert-linear-curve #[256 "\300\301\302\"`\"\207" [calc-do make-closure #[0 "\300\306=\206	\300\307=\310\311\312!\211A@A
AA@A\313
!\205!\314\315\"\314\316\"\317\320#\321\n$\211@\322\320@A@\"!D\310\211\312\323\312\324\"\210\310\nA@A\325@\326\312\327	A@EEE\330\310
\331\332@	@E\332A@\nA@EE\262
\300\306=\203\206\333\312\334#\210\202\342\300\307=\203\334\310\335\n@A@A@A@D\336%\331\324
8\331\337\340\327\341\342EDEG\203\313\211\324V\203\313\343\324Z\"\202\314\344\262
\257\262\333\312\345#\266\202\342\333\312\346#\210\347
\350\")\266\202,\207" [V0 calc-display-working-message calc-curve-varnames calc-curve-coefnames calc-curve-nvars calc-fit-to-trail calcFunc-efit calcFunc-xfit nil calc-top 1 math-contains-sdev-p mapcar #[257 "\300
\301\"\207" [math-get-sdev t] 4 "\n\n(fn X)"] #[257 "\300
!\207" [math-get-value] 3 "\n\n(fn X)"] math-map-binop math-div math-nlfit-least-squares math-neg calc-get-fit-variables 2 * - / t vec calcFunc-eq math-nlfit-enter-result "efit" math-nlfit-chi-sq #[771 "\300\301\302\303\"\"\"\207" [math-mul math-sub 1 math-div] 10 "\n\n(fn X A B)"] (calcFunc-fitdummy 1) calcFunc-neg (calcFunc-fitdummy 1) (calcFunc-fitdummy 2) calcFunc-utpc (var nan var-nan) "xfit" "fit" calc-record "parm"] 29]] 5 (#$ . 12262)])
(provide 'calc-nlfit)