Dear community,
I've prepared two new commands to solve non-traditional types of linear programming via least squares (minimum norm general solution).
Many thanks to Kit Baum for uploading their latest versions to the SSC!
lppinv – a module to solve an under-, over- and identified linear problem without an objective function with the Moore-Penrose pseudoinverse and singular value decomposition (SVD).
The algorithm solves "hybrid" least squares linear programming (LS-LP) problems with the help of the Moore-Penrose inverse (pseudoinverse), calculated using singular value decomposition (SVD), with emphasis on estimation of non-typical constrained OLS (cOLS), Transaction Matrix (TM), and custom (user-defined) cases. Eventual regression analysis and a Monte-Carlo-based t-test of mean NRMSE is performed, the sample being drawn from a uniform or a user-specified distribution (Mata function).
tmpinv – a module to solve an under-, over- and identified Transaction Matrix (TM) problem with the help of the Moore-Penrose pseudoinverse, singular value decomposition (SVD), an F-test from linear regression/t-test of mean normalized RMSE, and results adjustment for extreme values.
The algorithm solves "hybrid" linear programming-least squares (LP-LS) Transaction Matrix (TM) problems with the help of the Moore-Penrose inverse (pseudoinverse), calculated using singular value decomposition (SVD). Estimation using 2x2 to 50x50 contiguous submatrices, repeated with compensatory slack variables until NRMSE is minimized in a given number of iterations, is followed by an F-test from linear regression/t-test of mean NRMSE from a pre-simulated distribution (Monte-Carlo, 50,000 iterations with matrices consisting of normal random variates). The result is adjusted for extreme values to match the RHS via shares of estimated row/column sums if the corresponding option is specified.
They can help solve a variety of interesting problems like calculating input-output tables, transaction matrices (like bilateral trade and investment tables), supply-use relations, a GDP (depvar) model with regards to a business cycle definition, etc.
If you are interested, you can find more information and examples in their help files.
I'm working on papers on the topic.
Have a nice day,
IB
Related Posts with New lppinv and tmpinv commands in Stata 16+
Percentile Calculation Panel DataHi I have a panel data set from 1987Q1 to 2016Q4. I want to assign each observation using deciles. B…
Generating quarterly date columnHi guys, A quick simple problem which is driving me nuts: I am trying to generate a quarterly time …
Finite Population specification for multi-level cross-classified models with repeated measuresHi all, I'm working with a three-level cross-classified model with repeated measures (i.e. multiple…
Stata - reproducible research (in 2021)Stata's update to 17 is impressive and tempting (not least because of increased speed and improved i…
How to generate a local list with numbers in a simple wayHi all, I have a simple question. Say I want to process filenames or variable names with consecutive…
Subscribe to:
Post Comments (Atom)
0 Response to New lppinv and tmpinv commands in Stata 16+
Post a Comment