Биостатистика, фармакокинетика и клинические исследования.
Ранее уже писал, что есть некоторое несоответствие SPSS и SAS по расчетам степеней свободы. Не могу точно сказать – какие причины приводят к различиям, но они есть. Для иллюстрации можно рассмотреть доступный датасет, который входит в состав Phoenix WinNonlin RAEBE Template – Example Data EMA 618604 2008 Rev 7.xls Выполним код SAS для этого датасета:
Ниже представлен код для IBM SPSS для расчета 90% доверительных интервалов исследования биоэквивалентности. Код соответствует референсному коду SAS при условии положительной матрицы Гессе. Т.к. FDA допускает использование COVTYPE – CSH (Heterogeneous Compound Symmetry), а FA0() в SPSS отсутствует, то используется именно этот тип. Надо отметить, что в SPSS при использовании GENLINMIXED структура CSH – недоступна,…
IBM SPSS дает возможность разработки собственных команд синтаксиса. Это средство совместно с применением Python и R Project дает очень широкие возможности применения и расширения функционала SPSS (и это действительно важно, ведь SPSS не дает, к примеру, человеческим образом получить доверительный интервал для разности пропорций). А если вы захотите найти руководство, то задав поисковику множество запросов…
В IBM SPSS для расчтета доверительных интервалов разницы лог-преобразованных параметров используется процедура GLM. За основу был взять код предложенный пользователем Oiinkie на форуме bebac.at Небольшая модификация связана с явным указанием фактора для контраста: /CONTRAST(Trt)=Simple(1) – требуется с вниманием указать какой уровень фактора будет референсным. Дополнительный участок кода использует результаты таблицы контрастов для вычисления искомых доверительных…