Are there any Julia packages that allow one to specify spatial correlation (or any correlation structure) when implementing generalizes least squares regression? I’m looking for something akin to the gls() function in the nlme package in R, which allows one to specify variables that dictate correlation, along with correlation structure (e.g., AR1, symmetric, etc).

I’ve looked around a bunch but haven’t found anything even remotely close. Is anyone aware of something that might work?

I have a work-in-progress: MeasurementErrorModels.jl. It allows lag-covariance matrices for both the predictor noise and response noise. I’ve added a github gist which illustrates a spatial model example. I’m currently working on adding models that allow both noisy and noiseless predictors, and Fisher information derived confidence intervals. Currently I assume that the noise *covariances* are known, but I’m interested in implementing models where e.g. only the *correlation* structure is known too.

1 Like