Ниже представлен код для IBM SPSS для расчета 90% доверительных интервалов исследования биоэквивалентности. Код соответствует референсному коду SAS при условии положительной матрицы Гессе. Т.к. FDA допускает использование COVTYPE – CSH (Heterogeneous Compound Symmetry), а FA0() в SPSS отсутствует, то используется именно этот тип. Надо отметить, что в SPSS при использовании GENLINMIXED структура CSH – недоступна, поэтому остается только MIXED. В SAS использование FA0() и CSH дает идентичные результаты. При отрицательном гессиане меняется расчет степеней свободы и результаты уже не совпадают, но разница и ошибка остаются прежними. Об этом в другой записи. По умолчанию расчет степеней свободы в MIXED выполняется по Satterthwaite.
1 2 3 4 5 6 7 8 | MIXED logAUC BY Period Formulation Sequence Subject /CRITERIA=CIN(90) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE) /FIXED=Period Formulation Sequence | SSTYPE(3) /METHOD=REML /RANDOM=Formulation | SUBJECT(Subject) COVTYPE(CSH) /REPEATED=Formulation | SUBJECT(Subject*Period) COVTYPE(DIAG) /EMMEANS=TABLES(Formulation) COMPARE REFCAT(LAST) ADJ(LSD). |