diff --git a/+nla/+qualityControl/checkHeadMotion.m b/+nla/+qualityControl/checkHeadMotion.m
index c937609..505a5b4 100755
--- a/+nla/+qualityControl/checkHeadMotion.m
+++ b/+nla/+qualityControl/checkHeadMotion.m
@@ -15,6 +15,10 @@ function checkHeadMotion(fig, input_struct, motion, remove_index)
     distances = nla.helpers.euclidianDistanceROIs(network_atlas);
     prog.Value = 0.75;
     [r_vec, p_vec] = corr(motion, functional_connectivity.v', 'type', 'Pearson');
+    if any(isnan(r_vec))
+        msgbox("NaNs appearing in correlation between motion and functional connectivity. Please choose another column for motion")
+        return
+    end
     
     prob = nla.TriMatrix(network_atlas.numROIs());
     r = nla.TriMatrix(network_atlas.numROIs());