Plot Utilities

74d0f691-a5a4-40e4-ae55-be0d0e82fe68 34 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

LMSA robust initial ... slow

Plot Utilities

vx vy char size clr plotG n vx length : plot vx 1 el vy 1 el char size clr augment : i 2 n range plot plot vx i el vy i el char size clr augment stack : for plot 41 line :

Bar plot unicolor style, bar width user. data dt BarSize v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

ListCollect

data list ListCollect collects index list select list : ct 1 : i select a ct el data i 1 el : b ct el data i 2 el : ct ct 1 + : 31 line for a b augment 51 line :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

Levenberg Marquardt's second algorithm . . . Carlos 20200115

x y β Tol LMSA Eps Tol : eta1 Eps sqrt : eta2 eta1 : h eta1 : stop 0 : 15 mat m Clear 1 : m x rows : n β rows : 12 mat β fx i 1 m range ff i el x i el ft y i el eval - : for ff 21 line : β derfh i 1 m range j 1 n range β j el β j el h + : jjh1 β fx i el : β j el β j el 2 h * - : jjh2 β fx i el : jjh0 i j el jjh1 jjh2 - 2 h * / : β j el β j el h + : 61 line for 11 line for jjh0 21 line : A n n Zeros : g n Zeros : J m n Zero : 13 mat Fn 1 Zero : dL 1 Zero : dF 1 Zero : F 1 Zero : 14 mat f β fx : J β derfh : 12 mat A J transpose J * : g J transpose f * : ng g normi : 13 mat F f transpose f * 2 / : mu eta1 A Diag max * : nu 2 : 13 mat stop ¬ ng eta2 ≤ stop 1 : p A mu n identity * + invert g - (* : break try np β normi : nx eta2 β normi + : np eta2 nx * ≤ stop 2 : stop ¬ xnew β p + : fn xnew fx : Jn xnew derfh : Fn fn transpose fn * 2 / : dL p transpose mu p * g - (* 2 / : dF F Fn - : dL 1 el 0 > (dF 1 el 0 > (& β xnew : F Fn : J Jn : f fn : A J transpose J * : g J transpose f * : ng g normi : mu mu 13 / 12 dF 1 el dL 1 el / * 1 - (3^- Max * : nu 2 : 91 line mu mu nu * : nu 2 nu * : 21 line if 71 line 11 line if if 41 line if 11 line while β 121 line :

Javier data XY 30507011014116319132122242925262618278293330360400440172 mat :

β x y Clear 1

t0 1 time :

Original data set U XY :

XY XY 117 range 12 range el :

Isolate/store the conserved range ... x XY 1 col eval : y XY 2 col eval : 12 mat

support U . 10 black plot :

shifted symmetric Gaussian x ft β 1 el x β 2 el - (2^β 3 el / - exp * :

vector of guesses ... β 10101013 mat transpose :

Tol 103 -^:

once solved disable ... isolate 'α' α x y β Tol LMSA :

1 time t0 - min 3.3

Max ...
location of Max
    ◄   width  ► α 30.98147 22.986893 15.631117 31 mat

α 30.98147 22.986893 15.631117 31 mat :

recast the fit x f α 1 el x α 2 el - (2^α 3 el / - exp * :

datum α 3 el :

zac u f datum - u 1232 solve : 19.7168 26.257 21 mat

x plot support x f x f x diff α 2 el 50 α 2 el 10 - 22 mat zac 1 el α 3 el zac 2 el α 3 el 22 mat 71 sys :

x plot

24 f α 2 el f 21 sys 29.0125 30.9815 21 sys

Smath ∫ ... acc 100
~ exact Gaussians x f x 344 int 217.1

General Gaussian ... shifted ►

β 2 el

Bar & Stem

Bar plot unicolor style, bar width user. data dt BarGraph v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

Spectrum stem plot for QuickPlot unicolor data Stem stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

Minimize Conjugate-Gradient

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

init iter Minimize XY f φ β LS X XY 1 col : Y XY 2 col : n XY rows : k x β φ length 1 - : vx-> raw initial vector fit i 1 X length range vx i el X i el β f eval : for M-> Cholesky GradientMatrix j 1 k range i 1 n range M i j el X i el β φ j 1 + el eval : for for CholeskySolver LeastSquares Φ M transpose M * : d M transpose Y vx - (* : Δβ Φ 1 -^eval d * : b β Δβ + : 131 line : β 1 el init : i 1 : i iter ≤ β i 1 + el XY f φ β i el LS : i i 1 + : 21 line while β β iter el : 41 line β 31 line :

Vector Gradient ... Optimiz Symbolic
seemingly a scalar kernel  not 
recognized in program loop x β f n φ ct 1 : i 1 n range PD ct el x β f (β ct el diff : stack f(x,β) & PD V 1 el x β f : V i 1 + el PD i el : 21 line ct ct 1 + : 41 line for V 31 line :

Error stem plot for QuickPlot unicolor data Error stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

L H N xd U 0 : dx H L - N / : i 1 N 1 + range U i el L dx - (dx i * + : for U 41 line :

Bar plot unicolor style, bar width user. data dt BarSize v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

corr Pearson correlation χ X Mean : μ Y Mean : n XY rows : 13 mat X i el χ - (Y i el μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum Y i el μ - (2^(i 1 n sum * sqrt / X i el χ - (X i el β f eval μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum X i el β f eval μ - (2^(i 1 n sum * sqrt / / 31 line :

correlation between [X,Y data]

correlation ... [X,fitted model]

Absolute  residuals observer
fetches β from calculated δ j 1 XY rows range Δ j el Y j el X j el β f - eval : for X Δ augment 21 line :

SSD = "Sum Square Differences" SSD Y i el X i el β f - (2^(i 1 XY rows sum :

Minimize Conjugate-Gradient

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

init iter Genfit XY f φ β LS X XY 1 col : Y XY 2 col : n XY rows : k x β φ length 1 - : vx-> raw initial vector fit i 1 X length range vx i el X i el β f eval : for M-> Cholesky GradientMatrix j 1 k range i 1 n range M i j el X i el β φ j 1 + el eval : for for CholeskySolver LeastSquares Φ M transpose M * : d M transpose Y vx - (* : Δβ Φ 1 -^eval d * : b β Δβ + : 131 line : β 1 el init : i 1 : i iter ≤ β i 1 + el XY f φ β i el LS : i i 1 + : 21 line while β β iter el : 41 line β 31 line :

Vector Gradient ... Optimiz Symbolic
seemingly a scalar kernel  not 
recognized in program loop x β f n φ ct 1 : i 1 n range PD ct el x β f (β ct el diff : stack f(x,β) & PD V 1 el x β f : V i 1 + el PD i el : 21 line ct ct 1 + : 41 line for V 31 line :

Error stem plot for QuickPlot unicolor data Error stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

L H N xd U 0 : dx H L - N / : i 1 N 1 + range U i el L dx - (dx i * + : for U 41 line :

Bar plot unicolor style, bar width user. data dt BarSize v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

corr Pearson correlation χ X Mean : μ Y Mean : n XY rows : 13 mat X i el χ - (Y i el μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum Y i el μ - (2^(i 1 n sum * sqrt / X i el χ - (X i el β f eval μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum X i el β f eval μ - (2^(i 1 n sum * sqrt / / 31 line :

correlation between [X,Y data]

correlation ... [X,fitted model]

Absolute  residuals observer
fetches β from calculated δ j 1 XY rows range Δ j el Y j el X j el β f - eval : for X Δ augment 21 line :

SSD = "Sum Square Differences" SSD Y i el X i el β f - (2^(i 1 XY rows sum :

t0 1 time :

Javier data XY 30507011014116319132122242925262618278293330360400440172 mat :

Symmetric Gaussian ...
General parametric user x β f β 1 el x β 2 el - (2^β 3 el / - exp * :

ConjugateGradient Vector Fnct x β φ x β f 3 φ :

X XY 1 col : Y XY 2 col : n XY rows : 13 mat support XY . 10 black plot : 21 sys

'init' is of same legth as β init 30201013 mat transpose :

iter 9 :

algo solving for β β init iter Genfit :

assign/isolate .. β β :

1 time t0 - 1899.5 s *

β 30.982767 22.987071 15.628292 31 mat

SSD 25.027

corr 0.834

β 30.982767 22.987071 15.628292 31 mat :

30.98147 22.986893 15.631117 31 mat

iter is adjusted on f(24,β) Max is @ β2 ... see 'plot' ........................ ►

24 β f 29.0140

β 2 el β f 30.983

x plot β 2 el 50 β 2 el 25 - 22 mat x β f support x β f x diff 61 sys :

x plot

sanity symmetry β 2 el 4 - β f β 2 el 4 + β f 21 mat 11.13 11.13 21 mat

Σ U 344 3.5 range : i 1 U rows range Cumul i el x β f x U 1 el U i el int eval : for U Cumul augment 31 line :

datum 12 :

Σ Σ eval :

xx β f datum - xx 344 solve 19.1369 26.8372 21 mat

Σ 217.1 2 1 sys

Closely exact ...
native Simpson  ∫ x β f x 344 int 217.1

Σ Σ rows row 44 217.1 12 mat

Carlos/Jean: 20230526