Gecko [ 122.154.134.11 ]

Name	Size	Permission
README	7.89 KB	-rw-r--r--
README.prev	31.02 KB	-rw-r--r--
calc-aent.elc	26.85 KB	-rw-r--r--
calc-alg.elc	46.35 KB	-rw-r--r--
calc-arith.elc	81.62 KB	-rw-r--r--
calc-bin.elc	21.53 KB	-rw-r--r--
calc-comb.elc	30.12 KB	-rw-r--r--
calc-cplx.elc	8.5 KB	-rw-r--r--
calc-embed.elc	29.83 KB	-rw-r--r--
calc-ext.elc	108.83 KB	-rw-r--r--
calc-fin.elc	11.84 KB	-rw-r--r--
calc-forms.elc	48.05 KB	-rw-r--r--
calc-frac.elc	5.58 KB	-rw-r--r--
calc-funcs.elc	25.78 KB	-rw-r--r--
calc-graph.elc	35.09 KB	-rw-r--r--
calc-help.elc	19.55 KB	-rw-r--r--
calc-incom.elc	5.29 KB	-rw-r--r--
calc-keypd.elc	15.83 KB	-rw-r--r--
calc-lang.elc	57.72 KB	-rw-r--r--
calc-loaddefs.el	10.61 KB	-rw-r--r--
calc-macs.elc	6.09 KB	-rw-r--r--
calc-map.elc	29.56 KB	-rw-r--r--
calc-math.elc	60.12 KB	-rw-r--r--
calc-menu.elc	32.85 KB	-rw-r--r--
calc-misc.elc	20.54 KB	-rw-r--r--
calc-mode.elc	19.45 KB	-rw-r--r--
calc-mtx.elc	6.77 KB	-rw-r--r--
calc-nlfit.elc	13.59 KB	-rw-r--r--
calc-poly.elc	24.53 KB	-rw-r--r--
calc-prog.elc	51.27 KB	-rw-r--r--
calc-rewr.elc	40.01 KB	-rw-r--r--
calc-rules.elc	17.6 KB	-rw-r--r--
calc-sel.elc	19.83 KB	-rw-r--r--
calc-stat.elc	13.97 KB	-rw-r--r--
calc-store.elc	17.31 KB	-rw-r--r--
calc-stuff.elc	8.07 KB	-rw-r--r--
calc-trail.elc	4.67 KB	-rw-r--r--
calc-undo.elc	3.35 KB	-rw-r--r--
calc-units.elc	37.86 KB	-rw-r--r--
calc-vec.elc	37.25 KB	-rw-r--r--
calc-yank.elc	16.66 KB	-rw-r--r--
calc.elc	101.08 KB	-rw-r--r--
calcalg2.elc	80.01 KB	-rw-r--r--
calcalg3.elc	42.38 KB	-rw-r--r--
calccomp.elc	32.86 KB	-rw-r--r--
calcsel2.elc	6.04 KB	-rw-r--r--

Code Editor : calc-nlfit.elc

;ELC���
;;; Compiled by mockbuild@c6b8.bsys.dev.centos.org on Thu Feb 19 13:02:38 2015
;;; from file /builddir/build/BUILD/emacs-23.1/lisp/calc/calc-nlfit.el
;;; in Emacs version 23.1.1
;;; with all optimizations.

;;; This file uses dynamic docstrings, first added in Emacs 19.29.
(if (and (boundp 'emacs-version)
	 (< (aref emacs-version (1- (length emacs-version))) ?A)
	 (or (and (boundp 'epoch::version) epoch::version)
	     (string-lessp emacs-version "19.29")))
    (error "`calc-nlfit.el' was compiled for Emacs 19.29 or later"))

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

(byte-code "\300\301!\210\300\302!\207" [require calc-arith calcalg3] 2)
#@70 Return the parameters A and B for the best least squares fit y=a+bx.
(defalias 'math-nlfit-least-squares #[(xdata ydata &optional sdata sigmas) "G\n\203�\306\307\n\"\202�\310	\311\"\n\203�\312\202�	\312\211\312\211\313\203\274�@@@\314
\203L�\315\"\202N�\"\314\203a�\315\"\202c�\"\314\203z�\315\316\211\"\"\202�\316\211\"\"\314\203\227�\315\316\"\"\202\235�\316\"\"\n\203\255�\314\f\315\311\"\"+AAA\202+�\317\316\f\"\316
\211\"\"\315\317\316\"\316
\"\"\"\315\317\316\f\"\316
\"\"\"\203.\315\"\320\315
\"!\315\f\" !\321\322!!E\321\322!E\323\211! E\323 EEE+\2023D.\n\207" [xdata n sdata s2data S Sx mapcar calcFunc-sqr make-list 1 0 nil math-add math-div math-mul math-sub math-neg sdev calcFunc-sqrt vec Sy Sxx Sxy D ydata s y x B A sigmas C22 C12 C11] 8 (#$ . 689)])
(defalias 'math-nlfit-get-rates-from-cumul #[(tdata qdata) "\303\304A@@\"\304	A@	@\"\"CG\305V\203J�\306\307\310\303\304\3058A@\"\304\305	8	A@\"\"\303\304A@@\"\304	A@	@\"\"\"\"\nBA\211\202�\303\304A@@\"\304	A@	@\"\"\nB\311\n!)\207" [qdata tdata pdata math-div math-sub 2 math-mul (float 5 -1) math-add reverse] 10])
(defalias 'math-nlfit-get-cumul-from-rates #[(tdata pdata q0) "C\nA\203+�\304	@\305\211\306\304\nA@\n@\"\"\307A@@\"\"\"	B\nAA\202�\310	!)\207" [q0 qdata pdata tdata math-add math-mul (float 5 -1) math-sub reverse] 8])
(defalias 'math-nlfit-find-qmax #[(qdata pdata tdata) "\306\307	#\310\n\"\311	@	A\"\312\307\fA@\f@\"!\313
\"\203,�\314\315
\"\202.�,\207" [pdata qdata ratios tdata lsdata qmax math-map-binop math-div math-nlfit-least-squares math-max-list math-neg math-lessp math-add (float 5 -1) a] 4])
(defalias 'math-nlfit-find-logistic-parameters #[(qdata pdata tdata) "\306	\n#\307\310\"\311\n\f\"\312
@!
A@E+\207" [qdata pdata tdata a newqdata bandc math-nlfit-find-qmax mapcar #[(q) "\302\303\304	\"\305\"!\207" [a q calcFunc-ln math-sub math-div 1] 5] math-nlfit-least-squares calcFunc-exp] 4])
(defalias 'math-nlfit-find-qmaxhalf #[(pdata tdata) "\304@A\"\305\306@\n\"\2031�\307	\310\211\311\307A@@\"\"\312A@@\"\"\"AA\202	�	*\207" [pdata qmh pmax tdata math-max-list 0 math-lessp math-add math-mul (float 5 -1) math-sub] 8])
(defalias 'math-nlfit-find-q0 #[(pdata tdata) "\306	\"\307\310\n\"\311	#\312\313\314\315\f\"	#\307\316\317\n\"\"\"\203,�\317\n\"\202�\320\n\"\313\314\321\f\"	#\322\211\323W\203u�\307\324\317
\"\"\312\313\314\325\f\"	#\307\326\317\n\"\"\"\203i�\202l�T\211\202B�\307\327\317
\"\".\207" [pdata tdata qhalf q0 qdata qmin math-nlfit-find-qmaxhalf math-mul 2 math-nlfit-get-cumul-from-rates math-lessp math-nlfit-find-qmax mapcar #[(q) "\302	\"\207" [q0 q math-add] 3] (float 5 -1) math-add math-sub #[(q) "\302	\"\207" [q0 q math-add] 3] 0 10 (float 5 -1) #[(q) "\302	\"\207" [q0 q math-add] 3] (float 5 -1) (float 5 -1) qmax qt i] 8])
(byte-code "\301B\302\301!\204�\303\301\304\"\210\305B\302\305!\204�\303\305\306\"\210\307B\302\307!\204-�\303\307\310\"\210\301\207" [current-load-list math-nlfit-delta default-boundp set-default (float 1 -8) math-nlfit-epsilon (float 1 -5) math-nlfit-max-its 100] 3)
(defalias 'math-nlfit-set-elt #[(vec i x) "S	\233\n\240\207" [i vec x] 2])
(defalias 'math-nlfit-get-elt #[(vec i) "S	8\207" [i vec] 2])
(defalias 'math-nlfit-make-matrix #[(i j) "\305\306\"\307\306	\fW\203�\310!\nB	T\211\202\n�\n+\207" [j k mat row i make-list 0 nil copy-sequence] 4])
(defalias 'math-nlfit-set-matx-elt #[(mat i j x) "S	S\n8\233\240\207" [j i mat x] 3])
(defalias 'math-nlfit-get-matx-elt #[(mat i j) "S	S\n88\207" [j i mat] 3])
(defalias 'math-nlfit-givens #[(C d) "\306!\306	!@GG\307\310\311\"\310\311\"\312\211\312\211\f\nX\203\fT\211X\203w\313\f#\313\f\211# \314\311 \"\204m\315\316!\317!\316 !\"\"\203r�\320 !\311\307\211\202\266�\317\321!\322\323\317\211\"\317 \211\"\"!\"\324\"\320\324 \"!\315\316!\"\203\254�\202\266�\324\321!\" \325\f\211$\210\325\f$\210\fT\"\"\nX\203+\313\f\"##\313\"#$\326\317#\"\317$\"\"%\323\317#\"\317$\"\"$%#\325\"$$\210\325\f\"#$\210\"T\"+\202\323�)\327	\"&\327	\f\"'\326\317'\"\317&\"\"%\323\317'\"\317&\"\"&%'\330	&#\210\330	\f'#\210+*T\211\2023�)\fT\211\202)�\n\312(\211\307Y\203\337\330
\311#\210\327	\"(T\211\"\nX\203\302\323(\317\313\"#\327\"\"\"\"(\"T\211\"\202\237)\330\320\324(\313\211#\"!#\210S\211\202\206*\nT\211X\203\376\330
\327	\"#\210T\211\202\345)\n\211\307Y\203\267\211\fTY\203\257\313\f#\314\307\"\203*\311\307\202b\315\316!\307\"\203G\322\326\307\317\211\"\"!\202b\324\307\316!\"\317\321!\322\326\307\317\211\"\"!\"\327
\"\327
\f\"\312%)*\323\317)\"\317*\"\"%\326\317*\"\317)\"\"*%)\330
*#\210\330
\f)#\210+S\211\202)\fS\211\202).\207" [C d n N j r 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 x w gamma sigma rho i cjj cij math-nlfit-delta k cjk cik h di dj s rj ri] 11])
(defalias 'math-nlfit-jacobian #[(grad xlist parms &optional slist) "\306	\203'�\307\n	@#
\203�\310\311\f\"\202�\fB)
A	A\211\204�\312!)\207" [j xlist grad parms row slist nil apply mapcar #[(x) "\302	@\"\207" [x slist math-div] 3] reverse] 5])
(defalias 'math-nlfit-make-ident #[(l n) "\304\211\"\305	X\203�\306\n	\211$\210	T\211\202�\n*\207" [n i m l math-nlfit-make-matrix 1 math-nlfit-set-matx-elt] 6])
(defalias 'math-nlfit-chi-sq #[(xlist ylist parms fn &optional slist) "\306	\2034�\307\310\n	@#\f@\"\f\203�\311
\f@\"\312\313
\211\"\")	A\fA\fA\f\202�)\207" [cs xlist fn parms ylist c 0 math-sub apply math-div math-add math-mul slist] 5])
(defalias 'math-nlfit-init-lambda #[(C) "\305@GG\203-�@\211\203%�\306\307\f@\f@\"\"\fA\211\204�)A\211\204
�\310\311\307\n	\"\"!+\207" [C N n l row 0 math-add math-mul calcFunc-sqrt math-div] 7])
(defalias 'math-nlfit-make-Ctilda #[(C l) "@G\304\n	\"\305\"*\207" [C n l bot math-nlfit-make-ident append] 3])
(defalias 'math-nlfit-make-d #[(fn xdata ydata parms &optional sdata) "\306	\2032�\307\310\n	@#\f@\"\203�\311
@\"\202 �
)B	A\fAA\202�\312!)\207" [d xdata fn parms ydata dd nil math-sub apply math-div reverse sdata] 5])
(defalias 'math-nlfit-make-dtilda #[(d n) "\302\303	\304\"\"\207" [d n append make-list 0] 5])
(defalias 'math-nlfit-fit #[(xlist ylist parms fn grad &optional slist) "\306	\n$\307
	\n%\310	\n
%\311\f!\312\313\204\276�W\203\276�T\312\204\251�\314\f\" \315\f@G\"!\316 !\"\"\317\320\321\n!\"##\310	#
%$\322$\"\203\236�\322\323\324$\"$\"%\"\203\212�\325\323\326\"$#\325\202\245�\327\326\"-\202;�\306	\n$\307
	\n%)\202&�\n.D\207" [grad xlist parms slist C fn 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 ylist d chisq lambda really-done iters math-nlfit-max-its done Ctilda dtilda zeta newparms newchisq math-nlfit-epsilon] 6])
(defalias 'math-nlfit-s-logistic-fn #[(x a b c) "\304\305\306\307	\310\307\n\"!\"\"\"\207" [a b c x math-div math-add 1 math-mul calcFunc-exp] 10])
(defalias 'math-nlfit-s-logistic-grad #[(x a b c) "\306\307	\"!\310\311\307\n\"\"\307\f\211\"\312\311\f\"\313\312\307\f\n\"
\"!\313\312\307\f\307\307	\n\"\"\"
\"!E+\207" [c x ep b d d2 calcFunc-exp math-mul math-add 1 math-div math-neg a] 11])
(defalias 'math-nlfit-b-logistic-fn #[(x a c d) "\305\306\307	\n\"\"!\310\306\f\"\311\312\"\306\f\211\")\")\207" [c x d ex a calcFunc-exp math-mul math-sub math-div math-add 1] 6])
(defalias 'math-nlfit-b-logistic-grad #[(x a c d) "\306\307\310	\n\"\"!\311\312\"\310	\n\"\313\f\307\211\")\"\310\313\307\307
\"\"\f\307\211\")\"\313\307\314\307\307
\307\211\")\"\"\"\315\f\316\"\"\"\310\313\307\314\307\307\307\211\")\"\"\"\315\f\316\"\"\313\307\307\"\"\f\307\211\")\"\"E+\207" [c x d ex ex1 xd calcFunc-exp math-mul math-sub math-add 1 math-div 2 math-pow 3 a] 13])
(defalias 'math-nlfit-find-covar #[(grad xlist pparms) "\304	\203�\305\306\n	@#BB	A\211\204�\305\307!B\310\311!\"\312!)\207" [j xlist grad pparms nil vec apply reverse math-mul calcFunc-trn calcFunc-inv] 6])
(defalias 'math-nlfit-get-sigmas #[(grad xlist pparms chisq) "\306\307	\n#\211GS\nG\310\f
V\2039�\f
X\2035�\311\f\211\f88!B\fT\211\f\202�\312!\f-D\207" [sgs grad xlist pparms covar n nil math-nlfit-find-covar 1 calcFunc-sqrt reverse N i] 5])
(defalias 'math-nlfit-s-logistic-params #[(xdata ydata) "\303	\"\304	\n#)\207" [xdata ydata pdata math-nlfit-get-rates-from-cumul math-nlfit-find-logistic-parameters] 4])
(defalias 'math-nlfit-b-logistic-params #[(xdata ydata) "\306	\"\307	\n#\310	#\211A@\311\f8\312\313\f@\313
\"\"!\312\314\315
!\"!\314
\"\211E.\207" [ydata xdata q0 qdata abc B math-nlfit-find-q0 math-nlfit-get-cumul-from-rates math-nlfit-find-logistic-parameters 2 math-neg math-mul math-div calcFunc-ln C A D] 7])
(defalias 'math-nlfit-s-logistic-solnexpr #[(pms var) "@A@\3058\306\307\310\311\n\312\311	\fE!EEE+\207" [pms c b a var 2 / + 1 * calcFunc-exp] 10])
(defalias 'math-nlfit-b-logistic-solnexpr #[(pms var) "@A@\3058\306\307\310\307\n\311\f	EE!E\312\313\314\310\307\n\311\f	EE!E\305EE+\207" [pms d c a var 2 / * calcFunc-exp - ^ + 1] 11])
(defalias 'math-nlfit-enter-result #[(n prefix vals) "\304\n#\210\305 \207" [prefix calc-aborted-prefix n vals calc-pop-push-record-list calc-handle-whys] 4])
(defalias 'math-nlfit-fit-curve #[(fn grad solnexpr initparms &optional sdv) "\300\301`\"\207" [calc-do #[nil "\306=\206	�\307=\310\311\312!\211A@AAA@A\313
!\205%�\314\315
\"\"\314\316
\"\310\211#$\312%\317\312\320\"&#A@'$A()\f
\"*\321\f
*+,\"&\211-A@.	\205o�\322,\f.-@$/	\205x�/@\2110\203\210�\323\324.0#\202\212�..1.'\"2\325\31034(\203\273�\326(@.@E3B3.A.(A\211(\204\240�\327\3303\237B!3\306=\203\324�\331\312\3322#\210\202%\307=\203\3105\33023/A@\333-@\fG.G67\"\203
76V\203
\334-@76Z\"\202\335*\2575\331\312\3365#\210)\202%\331\312\3372#\210\3403\341\".\207" [sdv sdevv calc-display-working-message data xdata ydata calcFunc-efit calcFunc-xfit nil calc-top 1 math-contains-sdev-p mapcar #[(x) "\301\302\"\207" [x math-get-sdev t] 3] #[(x) "\301!\207" [x math-get-value] 2] calc-get-fit-variables 3 math-nlfit-fit math-nlfit-get-sigmas math-map-binop #[(x y) "\302	E\207" [x y sdev] 3] t calcFunc-eq calc-normalize vec math-nlfit-enter-result "efit" (vec) calcFunc-utpc (var nan var-nan) "xfit" "fit" calc-record "parm" sdata calc-curve-varnames calc-curve-coefnames calc-curve-nvars fitvars var parms initparms parmguess fn grad fit finalparms sigmacovar sigmas solnexpr soln traillist calc-fit-to-trail sln m n] 11]] 3])
(defalias 'calc-fit-s-shaped-logistic-curve #[(arg) "\301\302\303\304\305%\207" [arg math-nlfit-fit-curve math-nlfit-s-logistic-fn math-nlfit-s-logistic-grad math-nlfit-s-logistic-solnexpr math-nlfit-s-logistic-params] 6 nil "P"])
(defalias 'calc-fit-bell-shaped-logistic-curve #[(arg) "\301\302\303\304\305%\207" [arg math-nlfit-fit-curve math-nlfit-b-logistic-fn math-nlfit-b-logistic-grad math-nlfit-b-logistic-solnexpr math-nlfit-b-logistic-params] 6 nil "P"])
(defalias 'calc-fit-hubbert-linear-curve #[(&optional sdv) "\300\301`\"\207" [calc-do #[nil "\306=\206	�\307=\310\311\312!\211A@AAA@A\313
!\205%�\314\315
\")\314\316
\"\317\320
\f#*\321\f*)	$\211+@\322\320+@+A@\"!D,\310\211-.\312/\323\312\324\"0-A@1.A2\325,@\326\312\3271,A@EEE3\330\31045\331\3322@,@E\3322A@,A@EE4\306=\203\245�\333\312\3343#\210\202\307=\203\310\335\f*,@A@,A@A@D\336)%67\33134\324+8\331\337\340\327\341\342EDE6\fG8)\203\363�8\324V\203\363�\34368\324Z\"\202\364�\344)\2577\333\312\3457#\210*\202\333\312\3463#\210\3474\350\".\207" [sdv sdevv calc-display-working-message data qdata pdata calcFunc-efit calcFunc-xfit nil calc-top 1 math-contains-sdev-p mapcar #[(x) "\301\302\"\207" [x math-get-sdev t] 3] #[(x) "\301!\207" [x math-get-value] 2] 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 #[(x a b) "\303\304\305\306	\n\"\"\"\207" [a x b math-mul math-sub 1 math-div] 7] (calcFunc-fitdummy 1) calcFunc-neg (calcFunc-fitdummy 1) (calcFunc-fitdummy 2) calcFunc-utpc (var nan var-nan) "xfit" "fit" calc-record "parm" sdata poverqdata parmvals finalparms calc-curve-varnames calc-curve-coefnames calc-curve-nvars fitvars var parms soln traillist calc-fit-to-trail chisq sln n] 11]] 3])
(provide 'calc-nlfit)

${this.title}

Code Editor : calc-nlfit.elc