Aptech Systems, Inc. Worldwide Headquarters
Aptech Systems, Inc.
2350 East Germann Road, Suite #21
Chandler, AZ 85286
Ready to Get Started?
For Pricing and Distribution
Training & Events
Step-by-step, informative lessons for those who want to dive into GAUSS and achieve their goals, fast.
Have a Specific Question?
Q&A: Register and Login
Premier Support and Platinum Premier Support are annually renewable membership programs that provide you with important benefits including technical support, product maintenance, and substantial cost-saving features for your GAUSS System or the GAUSS Engine.
Join our community to see why our users are considered some of the most active and helpful in the industry!
Where to Buy
Available across the globe, you can have access to GAUSS no matter where you are.
Recent Tagsapplications character vectors CML CMLMT Constrained Optimization datasets dates dlibrary dllcall error error handling errors Excel FANPACMT file i/o floating network GAUSS Engine graphics GUI hotkeys installation Java API license licensing linux loading data matrices matrix matrix manipulation Maxlik MaxLikMT Memory Optmum output PQG graphics proc procs RAM random numbers string functions strings structures threading Time Series writing data
Time Series 2.0 MT
Find out more now
Time Series MT 2.1
Same results with equal optimization procedure
I am running a minimization procedure using optmum with a ML function in the context of a Kalman Filter model. When I get the outputs I use them to get the covariance matrix by inverting the hessian. Then I take square roots and get the standard errors.
The problem is that if I run the code several times I get different values for the standard errors. This is a bit confusing
I copy parts of my code, may be they help
1. the optimization call:
2. the inversion of the hessian:
3. standard error:
Thanks in advance!
What is the condition number of the Hessian? d = log(cond(_opfhess)) is an approximate measure of the number of digits lost in computing an inverse. If the Hessian in _opfhess is computed numerically, 8 digits are lost computing the Hessian. If d is greater than 8, you have lost all accuracy in the standard errors. It also indicates a catastrophic loss of precision in the calculations and under that condition the Hessian calculation can vary with very small differences in rounding errors.
I checked the condition of the Hessian and is 4.4.
But in any case, I don’t think the rounding problem will be the cause for the standard errors of all the coefficients to change and even be all equal in some cases.
May be in my previous post I was not very clear, but the thing is I am running several times the same code for the same data and I get different values for the standard errors of the estimated coefficients.
I believe it is a data problem. When I omit the first 16 iterations in my Kalman Filter loop this problem disappears.
However, the problem I have now is that for different data vectors (different countries but the same variable) the code does not converge.
Could you may be give me an idea which underlying problems may be there when convergence is not achieve?
Sorry if my question is too general.
Thanks in advance.
What is happening is that when the main optimization methods in Optmum get stuck, a random search is used. The random search is what is causing differences in the standard errors. If you set the random number seed with the rndseed command, that will give the same answer from run to run.