I am running a Regional I-O model. Although I get printed the Regional Transactions Matrix, I cannot get the Leontief inverse of it printed. Instead, I get only dots.

```
/* Step 3.4: Computation of Regional Matrices */
DRMRg=CILQ.*DRMNt;
/* Regional Direct Requirements Matrix mxm */
ImpRgCf1=(sumc(drmnt-drmrg));
/*Regional Imports Coefficients vector adjusted for CILQ mx1*/
Print; print "Direct Requirements Regional for THESS"; print DRMRg;
/* Step 3.5: Sectoral Aggregation – Definition of the Regional Classification Scheme*/
C3=eye(m);
C4=C3-DRMRg;
LEONTIEFRg=inv(C4);
/* Regional Leontief Inverse Matrix mxm */
PrdnRg=(PrdnNt).*(EmpRg./EmpNt);
/* Regional Sectoral Output mx1 */
TransRg=DRMRg*diagrv(zeros(m,m),PrdnRg);
/* Regional Transactions Matrix */
ImpRgf=ImpRgCf1.*PrdnRg;
/* Regional Imports Vector mx1, from coefficients to values */
Print; print "Leontief Regional for THESS"; print LEONTIEFRg;
Print; print "Regional Sectoral Output"; print PrdnRg;
format /rds 14,3;
Print; print "Regional Transactions Matrix for THESS"; print TransRg;
```

## 1 Answer

0

In GAUSS, missing values and NaN's are represented by a dots. So if you print a GAUSS matrix and you see one or more dots, this is the problem. Here is a quick example:

```
// Create 2x2 matrix with one missing value
a = { 1 .,
3 4 };
b = 2.5;
// Create a 1x1 missing value
c = { . };
print "a.*b = " a.*b;
print "a.*c = " a.*c;
```

The printed output will be:

a.b = 2.5000000 . 7.5000000 10.000000 a.c = . . . .

You can use the GAUSS function `isinfnanmiss`

to check to see if a matrix contains ANY missing values, NaN's or infinities. For example after the code above, we could check `a`

like this:

```
if isinfnanmiss(a);
print "'a' has infinities, NaN's or missing values";
else;
print "'a' is a finite matrix";
endif;
```

To find the cause of your problem, you need to go back and see where the matrix with the missing values, `TransRg`

, was created. From your code, we can see that it was created on this line:

```
/* Regional Sectoral Output mx1 */
TransRg=DRMRg*diagrv(zeros(m,m),PrdnRg);
```

This means that either `DRMRg`

or `PrdnRg`

has NaN's. You can use the `isinfnanmiss`

function and/or some print statements to find out which is the problem and then trace it back to the source of the NaN's to solve your problem.

## Your Answer

## 1 Answer

In GAUSS, missing values and NaN's are represented by a dots. So if you print a GAUSS matrix and you see one or more dots, this is the problem. Here is a quick example:

```
// Create 2x2 matrix with one missing value
a = { 1 .,
3 4 };
b = 2.5;
// Create a 1x1 missing value
c = { . };
print "a.*b = " a.*b;
print "a.*c = " a.*c;
```

The printed output will be:

a.b = 2.5000000 . 7.5000000 10.000000 a.c = . . . .

You can use the GAUSS function `isinfnanmiss`

to check to see if a matrix contains ANY missing values, NaN's or infinities. For example after the code above, we could check `a`

like this:

```
if isinfnanmiss(a);
print "'a' has infinities, NaN's or missing values";
else;
print "'a' is a finite matrix";
endif;
```

To find the cause of your problem, you need to go back and see where the matrix with the missing values, `TransRg`

, was created. From your code, we can see that it was created on this line:

```
/* Regional Sectoral Output mx1 */
TransRg=DRMRg*diagrv(zeros(m,m),PrdnRg);
```

This means that either `DRMRg`

or `PrdnRg`

has NaN's. You can use the `isinfnanmiss`

function and/or some print statements to find out which is the problem and then trace it back to the source of the NaN's to solve your problem.