IBM SPSS Bioequivalence Full Replicate (Код IBM SPSS для репликативного дизайна биоэквивалентности)

Ниже представлен код для 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).

Добавить комментарий