Я создаю линейную смешанную модель, которая включает термин взаимодействия и случайный эффект на основе сайта. Результатом является непрерывный балл, а предикторные переменные - все двоичные категориальные переменные - статус воздействия (да / нет) и группа (а / б). Я хочу знать разницу в средствах для лиц, подвергшихся воздействию, и не подвергавшихся воздействию отдельно для групп А и В, а также доверительный интервал для каждого из этих различий в средствах. Моя модель:
PROC MIXED DATA = work.data1;
CLASS exposure group (ref = "b");
MODEL score = exposure group exposure*group;
RANDOM site /subject = id type = un SOLUTION;
RUN;
Я пытался использовать операторы ESTIMATE, LSMEANS и LSMESTIMATE, но, похоже, я не могу получить стратифицированную разницу в средствах.