#leastSquares

Marcin Paprzyckimarcinpaprzycki@masto.ai
2023-03-21

Using a #leastsquares method to estimate parameter values in the Lotka-Volterra model: β€œ#ParametersEstimation of a Lotka-Volterra Model in an Application for Market Graphics Processing Units” by D. Normatov, P. Mercorelli. ACSIS Vol. 30 p. 935–938; tinyurl.com/2p8s6kbu

2022-11-22

Last lecture πŸ”– of this semester is done! Interestingly the last topic was linear regression πŸ“ˆ and correlation. A fascinating topic with lots to say about #Leastsquares and scientific #fitting of data in general. But with little time βŒ› and little #mathematics it's hard to outline how useful this approach can be.

Looking forward to the mini projects results the students have to hand in soon. Hopefully they learned a little about #statistics πŸ˜…

2020-03-18

The #VectorAutoRegression stuff I've been doing can be summarized as
$$ y_t = \sum_{i=1}^p A_i y_{t - i} $$
where each $y_t$ is a $D$-vector and each $A_i$ is a $D \times D$ matrix.

Given an input series of $y$ values, the $A_i$ can be calculated by #LeastSquares minimization:

$$
Y = [ y_p ; y_1; \ldots ; y_{T - 1} ] \\
X = [ y_{p-1}, y_{p-2}, \ldots, y_0 ; y_p, y_{p-1}, \ldots, y_1 ; \ldots ; y_{T-2}, y_{T-3}, \ldots, y_{T-1 - p} ] \\
A = \left(X^T X\right)^{-1} X^T Y
$$
where the matrices have these initial dimensions:
Y : (T - p) Γ— D
X : (T - p) Γ— (p Γ— D)
A : (p Γ— D) Γ— D
then reshape $A$ to $p Γ— (D Γ— D)$

In my case all values are complex (and ^T is conjugate transpose) because each $y$ is an FFT of a block of input data - I'm using FFT size $256$ (making $D = 129 = 256/2+1$ unique bins for real input) overlapped 4x.

Client Info

Server: https://mastodon.social
Version: 2025.04
Repository: https://github.com/cyevgeniy/lmst