Синтаксис BayesFactor novaBF

Я хочу получить коэффициенты Байеса для ANOVA, которые аналогичны классическим F-тестам, и я просто хочу убедиться, что правильно понимаю, как писать синтаксис, особенно в отношении идентификаторов субъектов.

Например, у меня есть независимые переменные между субъектами a_between и b_between и внутрисубъектные переменные c_within и d_within, зависимая переменная values, где subject_id идентифицирует каждого субъекта; в наборе данных my_data.

Если я правильно понимаю, для полного дисперсионного анализа я должен использовать:

anovaBF(values~a_between*b_between*c_within*d_within+subject_id, data = my_data, whichModels = "bottom", whichRandom = "subject_id") # and I assume the order of variables does not matter, e.g. it could also be d_within*a_between*c_within*b_between+subject_id

Только для дисперсионного анализа внутри субъектов я должен использовать:

anovaBF(values~c_within*d_within+subject_id, data = my_data, whichModels = "bottom", whichRandom = "subject_id")

Только для дисперсионного анализа между субъектами я должен использовать:

anovaBF(values~a_between*b_between, data = my_data, whichModels = "bottom", whichRandom = "subject_id")

Так что в последнем случае у меня нет +subject_id - иначе я получаю Error in base::try(expression, silent = silent) : not enough observations. (Возможно, потому что для каждого subject_id есть только одна строка?)

Два основных вопроса:

1. Независимо от причин, являются ли приведенные выше решения правильными?
2. Если решения верны, почему я должен дважды указывать идентификатор субъекта (один раз как whichRandom и один раз в начале как +subject_id) для внутритематических переменных, и почему нет, когда есть только межтематические переменные?

(К вашему сведению, есть связанный вопрос с ответом, но не совсем то, что я хочу знать: https://stats.stackexchange.com/questions/230224/mixed-bayesian-anova-using-bayesfactor-package-in-r)