Я должен найти 95% CI на медиане и других процентилях. Я не знаю, как подойти к этому. Я в основном использую R как инструмент программирования.
r
confidence-interval
median
Доминик Комтуа
источник
источник
library(boot)
появляется для подтверждения этого:> boot.ci (boot (x, function (x, i) медиана (x [i]), R = 1000)) Интервалы: Уровень Нормальный Базовый 95% (74.42, 78.22) (75.00 , 78,49) Процентный уровень BCa 95% (73,51, 77,00) (73,00, 77,00)Другой подход основан на квантилях биномиального распределения.
например:
источник
Проверьте бутстрап ресэмплинга. Поиск R справка для функции загрузки. В зависимости от ваших данных с передискретизацией вы можете оценить доверительные интервалы для чего угодно.
источник
wilcox.test(..., conf.int=TRUE)
функцией R.Есть и другие подходы: один из них основан на тесте суммы рангов Уилкоксона, примененном к одному образцу с коррекцией непрерывности. В R это может быть поставлено как:
И здесь есть CI Дэвида Олив для медианы:
CI для медианы
источник
Результат, основанный на подходе qbinom, не подходит для небольших выборок. Предположим, что x имеет 10 компонентов. Тогда qbinom (c (.025, .975), 10, .5) дает 2 и 8. Результирующий интервал не обрабатывает статистику порядка в нижнем хвосте симметрично статистике из верхнего хвоста; Вы должны получить либо 2 и 9, либо 3 и 8. Правильный ответ - 2 и 9. Вы можете проверить однофакторность в SAS. Поймайте здесь, вам нужно не более 0,025 вероятности ниже и выше; нижний квантиль не делает этого, так как он дает по крайней мере 0,025 в или ниже. Вы спасены снизу, потому что счетчик, который должен быть 1, должен быть сопоставлен со статистикой второго порядка, считая 0, и поэтому «выключено одним» отменяется. Эта случайная отмена не происходит сверху, и поэтому вы получаете неправильный ответ здесь. Сортировка кода (x) [qbinom (c (.025, .975), длина (x) ,. 5) + c (0,1)] почти работает, и .5 можно заменить другими значениями квантиля, чтобы получить доверительные интервалы для других квантилей, но это не будет правильно, когда существует такое, что P [X <= a ] =. 025. См., Например, Хиггинс, Непараметрическая статистика.
источник