diff --git a/docs/pyerrors/fits.html b/docs/pyerrors/fits.html
index 98dd2d7b..a90800ae 100644
--- a/docs/pyerrors/fits.html
+++ b/docs/pyerrors/fits.html
@@ -367,7 +367,7 @@
256 if sorted(list(funcd.keys())) != key_ls:
257 raise ValueError('x and func dictionaries do not contain the same keys.')
258
-259 x_all = np.concatenate([np.array(xd[key]) for key in key_ls])
+259 x_all = np.concatenate([np.array(xd[key]).transpose() for key in key_ls]).transpose()
260 y_all = np.concatenate([np.array(yd[key]) for key in key_ls])
261
262 y_f = [o.value for o in y_all]
@@ -1351,7 +1351,7 @@ Hotelling t-squared p-value for correlated fits.
257 if sorted(list(funcd.keys())) != key_ls:
258 raise ValueError('x and func dictionaries do not contain the same keys.')
259
-260 x_all = np.concatenate([np.array(xd[key]) for key in key_ls])
+260 x_all = np.concatenate([np.array(xd[key]).transpose() for key in key_ls]).transpose()
261 y_all = np.concatenate([np.array(yd[key]) for key in key_ls])
262
263 y_f = [o.value for o in y_all]