Стандартное отклонение разницы средних |
Здравствуйте, гость ( Вход | Регистрация )
Стандартное отклонение разницы средних |
15.05.2008 - 05:55
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 36 Регистрация: 19.10.2007 Из: Москва Пользователь №: 4460 |
Есть среднее +/- стандартное отклонение показателя до лечения, есть среднее +\- после лечения. Как расчитать разницу средних и стандартное отклонение? Прочитал. что можно по формуле (корень из суммы квадратов станд отклонений средних значений до и после), но когда использовал формулу имея реальные данные оказалось, что величины отклонений (реального и рассчитаного) сильно различаются. Например: уровень креатинина до 92+/-24, после 88+/-23, разница по реальным данным вышла -3,6+/-12,8. По формуле, соответственно, +/-33, 2. Подскажите, пожалуйста, что не так ...
|
|
15.05.2008 - 06:45
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 1114 Регистрация: 10.04.2007 Пользователь №: 4040 |
А может, взять исходные выборки до и после лечения (очевидно, данные парные?)? Из них составить одну выборку разностей вариант. По полученной выборке рассчитать искомые статистические показатели.
Ebsignasnan prei wissant Deiws ainat! As gijwans! Sta ast stas arwis!
|
|
15.05.2008 - 10:29
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 36 Регистрация: 19.10.2007 Из: Москва Пользователь №: 4460 |
Еще раз!
Анализирую материал из статьи. Есть 2 группы (активное лечение и плацебо). Известны исходные и итоговые средние значения некоего показателя +/- стандартное отклонение в обеих группах. Необходимо оценить среднее значение изменения показателя (отдельно для каждой группы) и определить стандартное отклонение этого изменения (опять же отдельно для каждой группы). Я думал все просто, считаешь стандартное отклонение разности средних и все ... Но когда проверил результат на реальных данных, различия реального стандартного отклонения разности средних (итоговое минус исходное) от расчитываемого, были существенными. Видимо не понимаю чего-то ну очень простого ... |
|
15.05.2008 - 12:53
Сообщение
#4
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Вы используете не ту формулу SD для оценки средней разности между выборками с попарно связанными вариантами. Вы взяли просто формулу для оценки ошибки разности средних для равночисленных, но не связанных выборок, очевидно из Лакина стр 115. Но перелистнув 2 страницы, есть подробное описание с примером как оценить среднюю разность при связанных выборках и дана формула расчета ошибки средней разности. В формулу входит средняя разность и сумма квадратов разностей. Чтобы получить SD умножте результат на корень иp n. Пример из Лакина и решение с формулами - в пришпиленном файле, но формулы могут слететь, файлы с 2007 и 2003 wordom.
Прикрепленные файлы
|
|
15.05.2008 - 19:09
Сообщение
#5
|
|
Группа: Пользователи Сообщений: 36 Регистрация: 19.10.2007 Из: Москва Пользователь №: 4460 |
Правильно говорят, один дурак задаст вопрос и 100 мудрецов не ответят
Книгу уважаемого Лакина не встречал, к сожалению. Но все равно спасибо. Еще раз, но уже совсем просто. У больных с заболеванием А уровень креатинина до начала лечения 92,4 +/- 24,5 мкмоль/л, после лечения 88,9 +/- 23,3 мкмоль/л, объем выборки - 378 пациентов. Единственное чего нет, так это матрицы данных, чтобы посчитать изменение +/- стандартное отклонение. так вот, можно ли посчитать это самое отклонение без данных, имея только результат (исходно/после лечения)? |
|
15.05.2008 - 21:16
Сообщение
#6
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
На самом деле, вы задали очень важный вопрос. Очень часто, читая литературу, хочется сравнить свои результаты с результатами других методик, или сравнить несколько методов. Разность до и после лечения дает возможность сравнить несколько методик лечения посредством сравнения доверительных интервалов (ДИ) , но для расчета ДИ нужно знать SD или m. Все программы при сравнении средних значений связанных выборок предоставляют не только SD для каждого ряда, но и SD разницы средних, которые почему то в работах не приводят. К сожалению не имея рядов нельзя посчитать сумму квадратов разностей, которая входит в формулу расчета m.
|
|
16.05.2008 - 09:35
Сообщение
#7
|
|
Группа: Пользователи Сообщений: 1013 Регистрация: 4.10.2006 Пользователь №: 1933 |
Рассчитать по приведенным данным стандартное отклонение для разностей нельзя. Ибо формула расчета
SD^2=SD1^2+SD2^2-2*r*SD1*SD2, где r - корреляция между значениями первого и второго определений. Соответственно, не зная этого r SD суммарное не определяется (так что мудрецы-то могут ответить). Да, и кстати запись +/- это запись не для стандартного отклонения, а для стандартной ошибки, которая значительно меньше стандартного отклонения (в корень из числа наблюдений раз). Так что надо аккуратнее с данными. Для тех, кто еще будет читать эту ветку: обратите внимание, что из суммы квадратов вычитается произведение стандартных отклонений. Поэтому квадрат стандартного отклонения для парных случаев всегда меньше суммы квадратов стандартных отклонений и равен ему лишь в случае независимости наблюдений (что и наблюдается в случае непарных наблюдений). |
|
16.05.2008 - 11:09
Сообщение
#8
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
На всякого мудреца - довольно простоты
Автор постера ясно указал, что приводится среднее +-стандартное отклонение, что принято в качестве стандарта представления данных в зарубежных публикациях. Обратите внимание, что в отчетах, которые выдают обе программы, в пришпиленном файле, средняя разница приводится именно со значением SD, а не ошибки. Соглашусь с Плавом только в том, что действительно, отсутствуют данные для расчет SD разницы. Но не соглашусь с тем, что по приведенной формуле можно рассчитать SD. Оно получается значительно выше, чем есть на самом деле. Продолжение расчета по примеру из Лакина - в пришпиленном файле. Интересно, а что SAS выдаст по рассматриваемому примеру?
Прикрепленные файлы
|
|
16.05.2008 - 12:51
Сообщение
#9
|
|
Группа: Пользователи Сообщений: 1013 Регистрация: 4.10.2006 Пользователь №: 1933 |
Всегда приятно, когда есть люди, которые могут легко найти ошибку, правда не могут обнаружить где она скрыта. Для любителей все пересчитывать (формула не предназначалась для пересчетов, но это так) - пропущен множитель 2 перед произведениями SD. Добавьте и все получите с точностью до округления. Это метод не мой, это - вообще-то - основы анализа связанных переменных, так что вместо предложений по поводу считать в SAS (все программы дают одинаковые результаты) могли бы (если бы знали) просто указать на пропущенный множитель.
Что же касается +/- стандартное отклонение в зарубежных публикациях, попробуйте найти такое представление в зарубежных публикациях. Не надо смешивать французский с нижегородским. Стандартное представление M(SD), а не +/-, которое означает интервал, причем 68%. |
|
16.05.2008 - 15:50
Сообщение
#10
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
И мне приятно, что модератор легко ее может исправить. А насчет французского с нижегородским см. ссылку.
Прикрепленные файлы
|
|
16.05.2008 - 17:13
Сообщение
#11
|
|
Группа: Пользователи Сообщений: 36 Регистрация: 19.10.2007 Из: Москва Пользователь №: 4460 |
Коллеги, большое спасибо!
Открыли глаза еще на одну деталь. А по поводу +/- ... В зарубежных публикациях указывают по разному, и +/- и (SD). Основное правило - это подписать, что означает +/- или что находится в скобочках. Хотя соглашусь, в последнее время доминирует (SD). |
|
18.05.2008 - 11:17
Сообщение
#12
|
|
Группа: Пользователи Сообщений: 1013 Регистрация: 4.10.2006 Пользователь №: 1933 |
И мне приятно, что модератор легко ее может исправить. А насчет французского с нижегородским см. ссылку. Не хотел ввязываться в спор, но уж... Вы бы еще статью из континетального Китая прикрепили, а еще лучше из International Medical Journal, издаваемого в Москве. Не все "зарубежные" журналы созданы равными, при ссылках на то как надо делать следует смотреть на журналы с высоким воздействием impact, а не просто на то, что напечатано на английском. |
|
18.05.2008 - 21:11
Сообщение
#13
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Dear Sir Плав, Вы не ввязались в спор, а Вы его затеяли.
А спорить не о чем, в нашей специальности IOVS - ведущий международный журнал и он не виноват, что китайцы, и не только из Тайваня, там чаще печатаются, чем нижегородцы. Клинические испытания, результаты которых публикуются в этом журнале должны быть зарегистрированы по следующим ссылкам: http://www.actr.org.au (Australia) http://www.clinicaltrials.gov (US) http://isrctn.org (UK) http://www.umin.ac.jp/ctr/index/htm (Japan) http://www.trialregister.nl/trialreg/index.asp (Netherlands) Investigative Ophthalmology & Visual Science (IOVS), published monthly in print and online, is an official journal of the Association for Research in Vision and Ophthalmology (ARVO), an international organization whose purposes are to encourage and assist research, training, publication, and dissemination of knowledge in vision and ophthalmology. Included are original contributions that emphasize clinical and laboratory hypothesis-based research with statistically good results that clearly advance the fields of ophthalmic and vision research. IOVS de-emphasizes purely descriptive research. Contact Information for IOVS Editorial Office Address: 12300 Twinbrook Parkway Suite 250 Email: iovs@arvo.org Если у Вас есть международные правила оформления работ, то поделитесь первоисточником, пожалуйста. |
|
19.05.2008 - 11:11
Сообщение
#14
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Издательство ELSEVIER не нуждается в представлении. У них во всех работах, которые мне встречались среднее +- SD. В постсоветстком пространстве, видимо, можно и то и другое, потому и возникают вопросы.
Прикрепленные файлы
|
|
19.05.2008 - 12:41
Сообщение
#15
|
|
Группа: Пользователи Сообщений: 1114 Регистрация: 10.04.2007 Пользователь №: 4040 |
В постсоветстком пространстве, видимо, можно и то и другое, потому и возникают вопросы. Причем тут постсоветское пространство? Наука интернациональна. А национальная самоидентификация не должна идти далее культуры и языка, которые каждому следует беречь и хранить. Теперь по делу. Уважаемый доктор, посмотрите монографию "Тейлор Дж. Введение в теорию ошибок. - М.: Мир, 1985", с. 24 и далее. Там поясняется, почему после +- должна стоять величина, характеризующая именно ошибку, а не рассеяние. В цитируемой же работе далее даны подробные рассуждения о манипуляциях с величинами, рассматриваемыми в данной теме. Как раз ошибки (а не стандартные отклонения) позволяют данные манипуляции провести, не привлекая дополнительной информации об исходных данных, которая часто недоступна. Если возникнет затруднение с поиском печатного варианта книги, поищите в Интернете электронный вариант. Где-то я его видел Ebsignasnan prei wissant Deiws ainat! As gijwans! Sta ast stas arwis!
|
|