### Aptech Systems, Inc. Worldwide Headquarters

Address:

Aptech Systems, Inc.

2350 East Germann Road, Suite #21

Chandler, AZ 85286Phone: 360.886.7100

FAX: 360.886.8922**Ready to Get Started?**### Request Quote & Product Information

### Industry Solutions

### Products

### Resources

### Support

### Training & Events

Want more guidance while learning about the full functionality of GAUSS and its capabilities? Get in touch for in-person training or browse additional references below.

### Tutorials

Step-by-step, informative lessons for those who want to dive into GAUSS and achieve their goals, fast.

### What’s New in GAUSS 17

### Want to find out more?

### Have a Specific Question?

### Q&A: Register and Login

### Support Plans

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.

### User Forums

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 Tags

applications character vectors CML CMLMT Constrained Optimization datasets dates dlibrary dllcall econometrics Editor error error codes error handling errors Excel file i/o floating network GAUSS Engine GAUSS Light graphics GUI hotkeys installation license licensing linux loading data loop loops matrix manipulation Maximum Likelihood Maxlik MaxLikMT Memory multidimensional array optimization Optmum output PQG graphics procs random numbers strings structures threading### Recent Questions

- codes in 2016 and I cannot open them in 2017 version
- rows with selif
- How can I find the index of the largest element of a matrix
- Is there any function similar to rep in R?
- G0152 Variable not Initialized error message
- Does Gauss have similar function as " %in%" in R or "find" in matlab
- Shall I change from I7 6700K to Ryzen X1800
- inther procedure
- About precision of cdfMvn and multi-thread
- Axis fonts, LaTeX font in legend

### Features

### Time Series 2.0 MT

### Industry Solutions

### Find out more now

### Time Series MT 2.1

### Find out more now

### Find out more now

# Resources

# Can I define global level

Hi

I am doing a Monte Carlo simulation involving CML (constrained maximum likeliohood) procedure. In most cases, the codes are running well, however in certain DGP the inverse of a matrix, when calculating the log-likelihood, does not exist (non positive definiteness). I have set up a flag in the main body of the codes such that if there is any error in any particular Monte Carlo sample , that sample is dropped and new sample is drawn again. Now I set up trapchk and the trap in the following way

oldval = trapchk(1); trap 1,1; rhoi=invpd(rho); trap oldval,1; if scalerr(rhoi); print "wrong rho" rho; go to fap; endif;

Where Fap is the flag I set earlier. But it does not work. I tried with _fap, no luck. Can anybody help please?

## 3 Answers

Here is a simple working example that tries to convert the inverse of a non-positive definite matrix with the GAUSS command `invpd`. (NOTE that you can use the GAUSS command `inv` to compute the inverse of a non-positive definite matrix. Though, I assume in the context of your program this matrix must be positive definite).

//Create label mylabel: //non-positive definite matrix x = reshape(seqa(1, 1, 16), 4, 4); //store previous trap state old_trap = trapchk(1); trap 1; x_inv = invpd(x); trap old_trap; if scalinfnanmiss(x_inv); print "matrix not positive definite"; goto mylabel; endif;

It might be more clear to use a `do while` loop in this case, however. For example:

//Flag to track if matrix is positive definite matrix_ok = 1; do while (matrix_ok == 1); //non-positive definite matrix x = reshape(seqa(1, 1, 16), 4, 4); //store previous trap state old_trap = trapchk(1); trap 1; x_inv = invpd(x); trap old_trap; if scalinfnanmiss(x_inv); print "matrix not positive definite"; //Set flag to indicate matrix not pos. def. matrix_ok = 0; endif; endo;

Many thanks for this. But, suppose your non-positive definiteness problem appears within the user defined procedure (such as log-likelihood function) and you want to stop the program when the problem is identified by "trap" and resume again where it stopped.

To be more specific, suppose you are doing 2000 Monte Carlo simulations. You have identified the problem at 1401 replication (in your procedure). Then you want to drop this particular 1401-th sample, but keeping the previous 1400 replications results and drawing a new sample for 1401-th replication. To do this, I think I need to set a flag at the beginning of the simulation.

How can I use this flag within my procedure? or can I do that?

Many thanks in advance!

Are you using a loop for your simulation? Then use break. Doing it they way I have presented below you will get 2000 good simulations. The ones with singular rhos won't get counted.

simul = 1; do until simul > 2000; oldval = trapchk(1); trap 1,1; rhoi=invpd(rho); trap oldval,1; if scalerr(rhoi); print "wrong rho" rho; break; endif; . . . // cmlmt run . . simul = simul + 1; endo;

## Your Answer

## 3 Answers

Here is a simple working example that tries to convert the inverse of a non-positive definite matrix with the GAUSS command `invpd`. (NOTE that you can use the GAUSS command `inv` to compute the inverse of a non-positive definite matrix. Though, I assume in the context of your program this matrix must be positive definite).

//Create label mylabel: //non-positive definite matrix x = reshape(seqa(1, 1, 16), 4, 4); //store previous trap state old_trap = trapchk(1); trap 1; x_inv = invpd(x); trap old_trap; if scalinfnanmiss(x_inv); print "matrix not positive definite"; goto mylabel; endif;

It might be more clear to use a `do while` loop in this case, however. For example:

//Flag to track if matrix is positive definite matrix_ok = 1; do while (matrix_ok == 1); //non-positive definite matrix x = reshape(seqa(1, 1, 16), 4, 4); //store previous trap state old_trap = trapchk(1); trap 1; x_inv = invpd(x); trap old_trap; if scalinfnanmiss(x_inv); print "matrix not positive definite"; //Set flag to indicate matrix not pos. def. matrix_ok = 0; endif; endo;

Many thanks for this. But, suppose your non-positive definiteness problem appears within the user defined procedure (such as log-likelihood function) and you want to stop the program when the problem is identified by "trap" and resume again where it stopped.

To be more specific, suppose you are doing 2000 Monte Carlo simulations. You have identified the problem at 1401 replication (in your procedure). Then you want to drop this particular 1401-th sample, but keeping the previous 1400 replications results and drawing a new sample for 1401-th replication. To do this, I think I need to set a flag at the beginning of the simulation.

How can I use this flag within my procedure? or can I do that?

Many thanks in advance!

Are you using a loop for your simulation? Then use break. Doing it they way I have presented below you will get 2000 good simulations. The ones with singular rhos won't get counted.

simul = 1; do until simul > 2000; oldval = trapchk(1); trap 1,1; rhoi=invpd(rho); trap oldval,1; if scalerr(rhoi); print "wrong rho" rho; break; endif; . . . // cmlmt run . . simul = simul + 1; endo;