Spatial prediction code & paper with gIVH, prediction variance using delta method.  Currently spatial gIVH programmed with RSR-ICAR model which doesn't work well because you can't limit K to just those cells that are sampled since Q is still full dimensional. Need to switch to process convolution for gIVH computation.