Looking at more_replacement, this :
sample_shifts = approximation.sample_node(shifts, more_replacements={batches: fulldata[:128,:]}, size=100)
posteriorshifts = sample_shifts.eval()
…does return something. I’ll check if the data makes any sense and update.
Thanks!