Доверительный интервал долей в Excel Опции

[Naum]

Пожалуйста, научите считать доверительный интервал в Excel! понимаю, что очень глупая просьба, но....
Я изучала распространенность бронхиальной астмы в городе Н - получилось следующее: опрошенных 567, больных из них 37, распространенность 6,53%. По области, в которой находится город Н, распространенность от 5,7% до 5,9%. Рецензент, прочитав мою статью, прислала следующее: "При сопоставлении распространенности БА необходимо сравнивать известные популяционные величины и доверительный интервал данного исследования, а не точечную оценку. Так, 95% ДИ для распространенности БА 37/567=6.53% составляет 4.75-8.89. Данный интервал включает величины распространенности БА в городах области (5.7-5.9%), следовательно, различий в распространенности БА в сравниваемых регионах нет". Как получился этот ДИ 4,75-8,89? Недавно нашла на этом форуме как считать критерий Стьюдента в Excel - описано было очень подробно, шаг за шагом, доступно для таких чайников как я - супер! если можно, то в таком же ключе напишите

[Nikita]

Я перед Плав повинился, а Плав струсил. Хотя ничто не ценится так дорого и не стоит так дешево...
Сюда я больше не ходок.
На прощанье замечу, что смертельно обидевшись на безобидный полушутливый эпитет «небрежный», Плав сам здесь признал, что Clopper.xls «не программа, а иллюстрация к обсуждению». А народ ее уже скачивает (и наверное, судя по восторгам, активно пользует). Без четких разъяснений обязательно кто-то из пользователей воспользуется неадекватным ДИ, полученным с «помощью» «иллюстрации» Clopper.xls

Что касается того, что для получения надежных статистических значимостей на уровне 0,001 надо много работать, а это слишком дорого, то это из той оперы, когда скупой платит дважды. И в нашем случае платит пациент (и деньгами и здоровьем). А потом выясняется, что эффект, объявленный статистически значимым на уровне 0,05, не воспроизводится. Наука не даётся малой кровью. В научной литературе накоплены обширные сведения о том, что не менее 50% (а то и все 95%) «открытий», значимых на уровне 0,05, не воспроизводятся. См.: Ioannidis JPA (2005) Why Most Published Research Findings Are False. PLoS Med 2(8): e124. doi:10.1371/journal.pmed.0020124

Цифра низкой воспроизводимости таких «открытий» подтверждается и вычислениями по калибровке Р-значений и вероятностей воспроизведения.

Код

Pval      P(Ho|D)  P(H1|D)   Prep
0,05      > 30%     < 70%    < 50%
0,01      > 10%     < 90%    < 75%
0,001     > 2%      < 98%    < 90%

Pval - P-значение; Ho - нулевая гипотеза; H1 - альтернативная гипотеза; D - данные; Prep - вероятность воспроизведения.

Только при Р-значении 0,001 вероятность нулевой гипотезы H0 может приближаться к 2%, а вероятность воспроизведения наблюдаемого эффекта моежт приближаться к 90% !

Источники:
Thomas Sellke, M. J. Bayarri and James O. Berger. Calibration of p Values for Testing Precise Null Hypotheses. Source: The American Statistician, Vol. 55, No. 1 (Feb., 2001), pp. 62-71. Stable URL: http://www.jstor.org/stable/2685531 URL for draft: http://www.stat.duke.edu/~berger/papers/99-13b.ps

Формулы:
Нижняя граница Бейзова фактора (множителя, коэффициента) B для H0 по отношению к H1:
B=-epln(p) для p<1/e, где e - основание натурального логарифма, p - P-значение.
Для частотнической вероятности ошибки I рода (в отклонении H0):
Alpha(p)=B/(1+B)

Подробности и в частности мощное влияние на Р-значение объема выборки см.:
Adrian E. Raftery. Bayesian Model Selection in Social Research. Sociological Methodology, Vol. 25, (1995), pp. 111-163. Stable URL: http://www.jstor.org/stable/271063
Raftery, A.E. (1999). Bayes factors and BIC - Comment on "A critique of the Bayesian information criterion for model selection". Sociological Methods and Research, 27, 411-427.

Prep и другие можно вычислить с помощью программы LePAC или отдельного модуля LePREP:
http://lmrs.math.cnrs.fr/Persopage/Lecoutre/ErisA.html
См. также Сергиенко В.И., Бондарева И .Б. Математическая статистика в клинических исследованиях. ? М.: ГЭОТАР-Медиа, 2006. ? 304 с.

Сообщение отредактировал Nikita - 23.10.2008 - 19:36