*(Вопрос) Как мне построить 95% доверительный интервал для медианы в SPSS? * (Ответ) Автор: Marta Garcia-Granero, 11.3.2002. *Существует 2 способа. * Пример данных. INPUT PROGRAM. - VECTOR V(1). - LOOP #I = 1 TO 100. - COMPUTE V(1) = uniform(10). - END CASE. - END LOOP. - END FILE. END INPUT PROGRAM. execute. * Решение #1 (начиная с SPSS 11). TEMPORARY. COMPUTE constant=1. RATIO STATISTICS v1 WITH constant /PRINT = CIN(95) MEDIAN . * Решение #2 (любая версия SPSS). DEFINE CIMED95(!POSITIONAL !TOKENS(1)). * Сохраним исходный файл на будущее. SAVE OUTFILE='c:\\temp\\cimed__.sav'. SORT CASES BY !1 (A). COMPUTE ID=$CASENUM. RANK VARIABLES=!1 (A) /N INTO n /PRINT=NO . * Ранг нижней границы***. COMPUTE l=TRUNC(((n+1)/2-(n<=284)*(0.9789*sqrt(n)))-(n>284)*(0.98*sqrt(n))). * Обработка исключений (2 случая). IF n=17 l=5. IF n=67 l=26. * Ранг верхней границы. COMPUTE u=n+1-l. * Оставляем случаи, соответствующие нижней и верхней границам. SELECT IF((id=l) or (id=u)). EXECUTE . * Вывод результатов. STRING limits (A12). IF ($CASENUM=1) limits = 'НГ 95% ДИ:' . IF ($CASENUM=2) limits = 'ВГ 95% ДИ:' . REPORT FORMAT=LIST AUTOMATIC ALIGN(CENTER) /VARIABLES=limits !1 /TITLE 'Приблизительный 95% ДИ для медианы'. * Восстанавливаем исходный файл. GET FILE='c:\\temp\\cimed__.sav'. !ENDDEFINE. * Вызов макроса. cimed95 v1. *** С практической точки зрения имеет смысл рассмотреть 2 особых случая: когда размер выборки не превосходит 284 и когда он превосходит это число. Формула работает при размерах выборок, начиная с 6. Существуют 2 исключения, когда n=17 и n=67. Для них требуется корректировка расчётных показателей. Подробнее см. на http://www.math.unb.ca/~knight/utility/MedInt95.htm