diff --git a/resolve/likelihood.py b/resolve/likelihood.py
index 269f622b17ef8205b045524dab8df05d236451e3..1dec1910f3e1a89094e95107bd8e75b485d6399c 100644
--- a/resolve/likelihood.py
+++ b/resolve/likelihood.py
@@ -120,7 +120,7 @@ def ImagingLikelihood(
             R = FullPolResponse(observation, sky_operator.target)
     else:
         if len(sdom.shape) == 3:
-            R = MfResponse(observation, sdom)
+            R = MfResponse(observation, sdom[0], sdom[1])
         else:
             R = StokesIResponse(observation, sdom)
     model_data = R @ sky_operator