Доброго Вам времени суток!

Помогите, пожалуйста, решить животрепещущий вопрос!!!

Задача заключается в следующем: есть группа пациентов (n=500), у которых изучается ряд лабораторных показателей (концентрации веществ в крови, показатели не связаны между собой), эту группу нужно сравнить с контрольной (n=1209) но в контрольной группе известны лишь 95-е персентили (по сути верхняя граница нормы каждого лабораторного показателя) - исходных данных по контрольной группе нет, эти персентильные значения рекомендованы производителем реактивов. Если Вы знаете, как можно сравнить группы по значениям персентилей (или с использованием каких-либо других статистических параметров, которые можно вычислить у первой выборки), буду безмерно благодарна!!!

Единственное, что приходит в голову - это выстроить верхнюю доверительную границу для каждого конкретного параметра (лабораторного показателя) в изучаемой выборке и далее смотреть "на глаз", либо попробовать с помощью критерия Хотеллинга проверить гипотезу о многомерной средней.

Жду вашей безмерной благодарности Где можно получить?

Построить бутстреп оценку для Ваших перцентилей и сравнить с перцентилями производителя.

Прям здесь и сейчас )) Спасибо огромное за идею!
поскольку задача срочная - воспользуюсь для начала первым советом и представлю графиком для наглядности. Как правильно применить критерий Хотеллинга для гипотезы о многомерной средней мне пока непонятно...

Спасибо Вам за ответ!
А собственно как сделать такую оценку? и по какому критерию сравнивать в перцентилями производителя?

Вы просто делаете (например) 100000 перевыборок с возвращением из имеющейся выборки по конкретному показателю. Для каждой перевыборки рассчитываете нужный для сравнения перцентиль.

Затем сортируете полученный массив 100000 значений перцентиля, находите какое значение в этой отсортированной последовательности соответствует 2,5% (это по 2500 значений с края надо отсчитать).

Это будет граница 5% двустороннего критерия. Ну и если значение данное производителем вышло за эти границы то есть отклонение с 5% вероятностью ошибки.

Все эти расчеты легко выполнить в R, можете выложить пример данных.

Ремарка. Если вы пишите на русском, то и термины должны быть на русском. В русском языке нет ни слова "персент" (англ. - реrcent), ни "персентиль" (реrcentile), а есть процент (это то, что обозначают символом %) и процентиль. Дид ю андерстенд?

или перцентиль

http://dic.academic.ru/dic.nsf/enc_mathema...%98%D0%9B%D0%AC

PS но вариант с "про" конечно звучит профессионально

Of course, I've understood! Thank you for the comment!

Уже понятнее а какое количество перевыборок оптимально? от чего это зависит? опыт работы у меня в среде R отсутствует, к сожалению. Это можно сделать в Excel или в Statistica? если можно, то какими функциями, кроме сортировки и перцентилей, пользоваться. Не могу загрузить свой файл - вот пример по одному из показателей

Ala (исходные данные):
252,8
261,4
229,2
438,1
411,4
411,9
327,1
235,8
332,8
378,7
280,1
434,0
279,7
262,7
340,2
371,6
387,8
323,6
513,3
317,9
219,2
457,1
198,3
249,8
633,3
260,9
261,8
246,6
482,1
256,0
375,6
210,4
179,5
226,4 и т. д. - Контрольное значение = 643,0

данные на форум как я понял можно загрузить только в виде архива

помещаем в файл данные и импортируем их в R

> data<-read.table("data.txt", dec=",")

квантили считаются кучей алгоритмов, оптимально узнать что делал все таки изготовитель



по умолчанию метод ?7 построим бутстреп



Контрольное значение 643 достоверно больше (с вероятность ошибки не большей 0.05) значения которое дает данная выборка.

Общий смысл понятен! Спасибо! Испробую Ваш код в качестве первой задачи в R - как я понимаю, это достаточно гибкий инструмент для программируемой обработки,