Check for a y value being Inf and print error

[moritzpschwarz]

This happened to us for -Inf (logging a 0 variable):

library(gets)
set.seed(124)
x <- rnorm(100)
y <- rnorm(100) + 10
y[20] <- Inf

isat(y, mxreg = x, iis = TRUE)
IIS block 1 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

IIS block 2 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

IIS block 3 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

IIS block 4 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

GETS of union of retained IIS variables... 

SIS block 1 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

SIS block 2 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

SIS block 3 of 4:
25 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

SIS block 4 of 4:
24 path(s) to search
Searching: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 

GETS of union of retained SIS variables... 

GETS of union of ALL retained variables...
Error in if (ar.LjungBox$p.value <= ar.LjungB[2]) { : 
  missing value where TRUE/FALSE needed

For x variables this is fine:

library(gets)
set.seed(124)
x <- rnorm(100)
y <- rnorm(100) + 10
x[20] <- Inf

isat(y, mxreg = x, iis = TRUE)
Error in qr.default(x, tol = tol, LAPACK = LAPACK) : 
  NA/NaN/Inf in foreign function call (arg 1)