I have a dataset consisting of mobility data at the county-level for the US for about one year. So the number of observations is >1m. Apart from the county code, the date, and the mobility index, the dataset also has a number of socio-economic variables at the county level, which do not vary temporally (i.e. there is one single value for all day-county combinations) and also do not have any missing values. The mobility data contain a lot of missing values clustered both temporally and spatially, meaning that some counties have more missing values than others and some days have more missing values than others.
I'm looking for an imputation method that takes into account temporal, spatial correlation as well as the association with some sociodemographic county-level variables such as income, education, and so forth. I'm thinking that multiple imputation would be the way to go (like Hmisc or mice) but how would I incorporate the spatial and temporal correlation? Is there an R package that does what I'm looking for?