Корреляция Опции

[riba4kasonja]

Здравствуйте, сразу расскажу пример. ДЕти со сниженной пиковой скоростью выдоха Было обследовано 64, из них снижена пиковая скорость выдоха у 22. из этих 22 детей, 11 инфицированы хламидией и микоплазмой. Соответственное, есть ли между этим связь и каким способом лучше посчитать (или используя какой критерий доказать достоверность связи)?

[DoctorStat]

Здравствуйте, сразу расскажу пример. ДЕти со сниженной пиковой скоростью выдоха Было обследовано 64, из них снижена пиковая скорость выдоха у 22. из этих 22 детей, 11 инфицированы хламидией и микоплазмой. Соответственное, есть ли между этим связь и каким способом лучше посчитать (или используя какой критерий доказать достоверность связи)?

Для исследования зависимости между снижением пиковой скорости выдоха и инфицированностью в группе детей с пониженным слухом у Вас не хватает данных по инфицированности детей, у которых отсутствует снижение пиковой скорости выдоха.

Сообщение отредактировал DoctorStat - 4.01.2013 - 14:04

[riba4kasonja]

Для исследования зависимости между снижением пиковой скорости выдоха и инфицированностью в группе детей с пониженным слухом у Вас не хватает данных по инфицированности детей, у которых отсутствует снижение пиковой скорости выдоха.

Посчитала, у детей с нормальной скоростью выдоха маркеры хламидийной и микоплазменной инфекции обнаружились у 4 из 42 детей

[DoctorStat]

Посчитала, у детей с нормальной скоростью выдоха маркеры хламидийной и микоплазменной инфекции обнаружились у 4 из 42 детей

Делим всех детей с пониженным слухом на 4 группы в зависимости от пониженной скорости выдоха и инфицированности. Заносим для удобства полученный данные в таблицу Excel (см. левый рисунок). Анализируем полученную 2х2 таблицу сопряженности в любой программе, умеющей вычислять точный 2-х сторонний критерий Фишера. В качестве такой программы я выбрал электронную таблицу DoctorStat (см. правый рисунок). Из рисунка видно, что значимость P-VALUE=0.000544289, что меньше альфа=0.001. Отсюда делаем вывод, что инфицированность и скорость выдоха в группе детей с пониженном слухом связана с высокой степенью значимости альфа=0.001. Наличие инфекций (хламидии и микоплазмы) способствует пониженной пиковой скорости выдоха, т.к. отношение шансов иметь пониженную скорость выдоха в группе инфицированных по отношению к неинфицированным равно Odds_Ratio = (11*38)/(4*11)=9.5 > 1 .

Сообщение отредактировал DoctorStat - 4.01.2013 - 16:44

Эскизы прикрепленных изображений

 

[riba4kasonja]

Спасибо, а когда применяют корреляционные коэффициенты?

[Larina Tatjana]

Спасибо, а когда применяют корреляционные коэффициенты?

Коллега riba4kasonja!

Существуют более 10 видов коэффициентов корреляции. Поэтому и здесь следует уточнять о каких конкретно коэффициента Вы спрашиваете. Почитайте прекрасную книгу Фёрстера и Рёнца "Методы корреляционного и регрессионного анализа", и там найдёте многие ответы. Применительно же к Вашим данным можно говорить о корреляции (взаимосвязи) между двумя качественными признаками. Ниже привожу полученные по моей просьбе результаты детального анализа Вашей таблицы:


Statistics for Table of A by B

Statistic DF Value Prob
------------------------------------------------------
Chi-Square 1 13.1815 0.0003
Likelihood Ratio Chi-Square 1 12.7813 0.0004
Continuity Adj. Chi-Square 1 11.0223 0.0009
Mantel-Haenszel Chi-Square 1 12.9755 0.0003
Phi Coefficient 0.4538
Contingency Coefficient 0.4133
Cramer's V 0.4538

Здесь 3 последних показателя как раз и отражают интенсивность связи двух качественных признаков. Это Фи-коэффициент, коэффициент контингенции, и коэффициент Крамера. Рекомендую использовать именно его.
Кроме этих 3-х показателей интенсивности связи есть немало и других. Ниже привожу ещё одну таблицу с такими показателями:


Statistic Value ASE
------------------------------------------------------
Gamma 0.8095 0.1166
Kendall's Tau-b 0.4538 0.1184
Stuart's Tau-c 0.3652 0.1087

Somers' D C|R 0.4048 0.1158
Somers' D R|C 0.5088 0.1288

Pearson Correlation 0.4538 0.1184
Spearman Correlation 0.4538 0.1184

Lambda Asymmetric C|R 0.0000 0.3127
Lambda Asymmetric R|C 0.3182 0.1454
Lambda Symmetric 0.1892 0.1914

Uncertainty Coefficient C|R 0.1834 0.0957
Uncertainty Coefficient R|C 0.1552 0.0835
Uncertainty Coefficient Symmetric 0.1681 0.0888

Эти показатели хорошо описаны в книге Афифи и Эйзена.

Теперь относительно результатов проведения точного анализа по Фишеру. Ниже привожу результаты с использованием метода Монте-Карло с объёмом выборок в 10 млн.


Monte Carlo Estimate for the Exact Test

Pr >= ChiSq 5.544E-04
99% Lower Conf Limit 5.352E-04
99% Upper Conf Limit 5.736E-04

Number of Samples 10000000

Вывод: существует СТАТИСТИЧЕСКИ ЗНАЧИМАЯ взаимосвязь двух изучаемых качественных признаков.

Более глубокий анализ этих результатов требует вникания и в другие показатели. Вот, к примеру, оценки относительных рисков:

Estimates of the Relative Risk (Row1/Row2)

Type of Study Value 95% Confidence Limits
-----------------------------------------------------------------
Case-Control (Odds Ratio) 9.5000 2.5211 35.7985
Cohort (Col1 Risk) 5.2500 1.8902 14.5815
Cohort (Col2 Risk) 0.5526 0.3598 0.8489

Видим, что значения относительных рисков для разных колонок разные. В частности, максимальный вклад в эту зависимость вносит первая колонка. Тогда как вторая колонка вносит гораздо меньший вклад. Чтобы понять смысл этих различий, необходимо понимать ВСЮ ПРОЦЕДУРУ вычисления этих статистик. Вот тогда и станет ясно, в каких клетках положительная связь, а в каких - отрицательная.

Вывод: необходимо учиться, если действительно желаете извлечь ценную и полезную информацию. В противном случае всё то, что напишите далее с помощью участников форума, будет лишь создавать видимость знания.
Кстати, один из аспектов Вашего исследования вообще остался "за бортом". Ведь у всех пациентов помимо фиксации 2-х признаков есть и показатели по анамнезу. И тогда встаёт вопрос об использовании многомерных методов. Которые и помогут выяснить детали того, почему такова реакция организма.

Желаю успешного овладения знаниями в этой области.

[riba4kasonja]

Здравствуйте, Larina Tatjana. Спасибо вам огромное, что разъяснили мне хоть немного, заданный мной вопрос и потратили свое время. Как показывает опыт даж. многие, написавшие докторскую, только создают видимость знания. У меня конечно имеются книги по статистике, но про которую вы писали нет. Просто у всех коэффициентов корреляции есть много критериев когда не следует расчитывать тот или иной коэффициент и я не всегда могу понять (в силу неглубоких знаний) есть они у меня или нет. Про анамнез вы правы, просто он не у всех был настолько подробно собран в плане явлений бронхообструкции, поэтомы этот вопрос я немного опустила

Сообщение отредактировал riba4kasonja - 5.01.2013 - 19:49