hotfix: plot_rm_corr crash with specific column names (#351)

* hotfix: plot_rm_corr crash with specific column names

* missed a formula

pingouin/parametric.py

@@ -1164,6 +1164,8 @@ def anovan(data=None, dv=None, between=None, ss_type=2, effsize="np2"):
     # C marks the data as categorical
     # Q allows to quote variable that do not meet Python variable name rule
     # e.g. if variable is "weight.in.kg" or "2A"
+    assert dv not in ["C", "Q"], "`dv` must not be 'C' or 'Q'."
+    assert all(fac not in ["C", "Q"] for fac in between), "`between` must not contain 'C' or 'Q'."
     formula = "Q('%s') ~ " % dv
     for fac in between:
         formula += "C(Q('%s'), Sum) * " % fac
@@ -1709,6 +1711,9 @@ def ancova(data=None, dv=None, between=None, covar=None, effsize="np2"):
 
     # Fit ANCOVA model
     # formula = dv ~ 1 + between + covar1 + covar2 + ...
+    assert dv not in ["C", "Q"], "`dv` must not be 'C' or 'Q'."
+    assert between not in ["C", "Q"], "`between` must not be 'C' or 'Q'."
+    assert all(c not in ["C", "Q"] for c in covar), "`covar` must not contain 'C' or 'Q'."
     formula = "Q('%s') ~ C(Q('%s'))" % (dv, between)
     for c in covar:
         formula += " + Q('%s')" % (c)

pingouin/plotting.py

@@ -1001,6 +1001,9 @@ def plot_rm_corr(
     # C marks the data as categorical
     # Q allows to quote variable that do not meet Python variable name rule
     # e.g. if variable is "weight.in.kg" or "2A"
+    assert x not in ["C", "Q"], "`x` must not be 'C' or 'Q'."
+    assert y not in ["C", "Q"], "`y` must not be 'C' or 'Q'."
+    assert subject not in ["C", "Q"], "`subject` must not be 'C' or 'Q'."
     formula = "Q('%s') ~ C(Q('%s')) + Q('%s')" % (y, subject, x)
     model = ols(formula, data=data).fit()