Вычисление 95% доверительного интервала для медианы

*(Вопрос) Как мне построить 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