diff --git a/knoten/bundle.py b/knoten/bundle.py
index a5cfe5a05cd2f667f10a5b1f5c8c7d775ea83d39..d59250afe73e279038ce3ec60fb33562ef001632 100644
--- a/knoten/bundle.py
+++ b/knoten/bundle.py
@@ -475,8 +475,10 @@ def compute_sigma0(V, dX, W_parameters, W_observations):
     
     if (dof > 0):
         sigma0 = VTPV/dof
-    else:
+    elif (dof == 0):
         sigma0 = VTPV
+    else:
+        raise ValueError(f"Computed degrees of freedom [{dof}] is invalid.")
 
     return np.sqrt(sigma0)