Error log(log(ac05)/log(ac01)) : NaNs produced in GEE

[pepbioalerts]

Hi,
I get this error but I am not sure what it means or how to fix it. All coordinates are cartesian and integer and not NAs in the data.

Error in gee::gee(formula = formula, family = family, data = data, id = id, :
NA/NaN/Inf in foreign function call (arg 14)
In addition: Warning message:
In log(log(ac05)/log(ac01)) : NaNs produced

Any hint on what the issue could be ? I can give you a repex ifneeded

[levisc8]

Hi @pepbioalerts - sorry for the slow response! A reprex would be helpful, as I can't tell just from looking at this what the issue might be.

[yw3236]

@levisc8 @pepbioalerts - Hi Sam and pepbioalerts, I also encountered the same error when I ran spind GEE function. I specified "corstr" as "fixed". After checking the corresponding source code shown as the following:

if(corstr == "fixed"){
ac01 <- acfft(coord, res0, 1, 1.1, dmax = 1)
ac05 <- acfft(coord, res0, 5, 5.1, dmax = 1)
if(ac05 <= 0) v <- 1
if(ac05 > 0) v <- log(log(ac05) / log(ac01)) / log(5)
alpha <- ac01
para0 <- paste("n", nn,sep = "=")
para1 <- paste(", alpha", round(alpha, 3), sep = "=")
para2 <- paste(", v", round(v, 3), sep = "=")
A0 <- paste(para0, para1, para2)
id <- rep(1, nn)
coord <- cbind(x, y)
D <- as.matrix(stats::dist(coord))
R <- alpha ^ (D^v)
data <- data.frame(data, id)
suppressMessages(suppressWarnings(utils::capture.output({
mgee <- gee::gee(formula = formula, family = family,
data = data, id = id,R = R,corstr = "fixed",
scale.fix = scale.fix)
})))
...
...
...
}

I found ac01 is negative (-0.6596838), but ac05 is positive (0.3565027) in my case. Therefore, log(ac01) will produce NaNs. Could you please any hint to fix this issue? Thank you.