Здравствуйте, гость ( Вход | Регистрация )
30.10.2008 - 10:00
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 20 Регистрация: 29.10.2008 Из: Челябинск Пользователь №: 5459 |
Допустим имеем 2 группы по 6 крыс в каждой. У каждой крысы взято по 3 биоптата. Как анализировать?
Вариант 1 расчитать среднее по 3 биоптатам для каждой крысы и сравнить 6 и 6. (скажем, сделать Манна-Уитни) Вариант 2 сравнить 18 биоптатов одной группы с 18 другой. Правильно конечно делать первый вариант, поскольку при втором увеличивается вклад индивидуальной вариабельности и мы недооценим риск получить различия обусловленные случайной комбинацией крыс. Но чувствительность в первом варианте достаточно мала а хочется все-таки больше. Есть какие-нибудь ещё варианты? |
|
|
![]() |
![]() |
![]() |
30.10.2008 - 20:56
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 1219 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
Дисперсионный анализ с повторными измерениями разрабатывался и традиционно используется именно для вашего случая. Во-первых он позволяет учесть зависимый характер между измерениями, а во-вторых - разделить ошибку эксперимента на 2 части: в нашем случае - на биологическую и техническую. Изменчивость крыс внутри группы - биологическая изменчивость, а изменчивость биоптатов внутри крыс -методическая (техническая). Такой анализ будет работать и при изучении динамики (если бы 2 ваши группы характеризовали 2 временные точки), но отнюдь им не ограничивается. Дисперсионный анализ это - единственное что вы можете сделать, чтобы максимально использовать имеющуюся информацию. Работа с усредненными значениями - грубый подход снижающий мощность исследования и делающий во-многом бессмысленным вашу затею с уменьшением ошибки.
В иерархическом дисперсионном комплексе будут следующие факторы. (1) Группы. Фиксированный фактор, отражающий различия двух групп по числу нейтрофилов при усреднении данных по крысам и биоптатам внутри групп. (2) Крыса внутри группы. Случайный фактор, отражающий изменчивость крыс по числу нейтрофилов при усреднении данных по биоптатам от одной крысы - т.е. «внутри» крысы. Взаимодействия факторов (Группы х Крыса внутри группы) быть не может, т.к. в каждой группе крысы свои. (Если бы эксперимент был организован так, что и в первой и во второй группах были использованы одни и те же животные - мы имели бы другой экспериментальный план (experimental design) в котором бы было такое взаимодействие). В ошибку дисперсиооного комплекса и всего эксперимента уйдет техническая изменчивость между биоптатами внутри крысы. Таким образом, в результате единственного анализа общая изменчивость оказывается разложенной на все возможные составляющие компоненты. Это хорошо и грамотно, поскольку если бы вы использовали только один биоптат от животного, то в ошибку эксперимента ушла бы техническая изменчивость (ошибка анализа) и биологическая изменчивость между крысами. Ошибка была бы больше и обнаружить различия межу группами при фиксированном числе экспериментальных животных было бы труднее. А непараметрика здесь вообще грубо говоря все сведет на нет. Распределение сосчитанных значений близко к логарифмически нормальному, поэтому перед анализом исходные данные - количество нейтрофилов в поле зрения - можно преобразовать логарифмированием. Напишите какие статистические программы у вас есть. Я могу подсказать как делать такой анализ в Statistica, SYSTAT и SPSS. P.S. Думаю, что методическая изменчивость будет очень велика если смотреть только 1 поле зрения. Она может перекрыть (и скорее всего перекроет) именчивость между кусочками ткани из разных мест раны. В этом случае изменчивость между биоптатами будет почти целиком отражать неравномерность содержания нейтрофилов в поле зрения и смысл во взятии 3х биоптатов теряется. Сколько полей зрения рекомендуется просматривать по методике которой вы работаете? Сообщение отредактировал nokh - 31.10.2008 - 06:35 |
|
|
![]() |
![]() |
Varta Несколько биоптатов с крысы - как анализировать? 30.10.2008 - 10:00
плав Цитата(Varta @ 30.10.2008 - 10:00) Д... 30.10.2008 - 11:56
DrgLena А что и как измерили в биоптатах? Если например, ... 30.10.2008 - 12:16
Varta Плав, я поняла вашу точку зрения. Я неполно сформу... 30.10.2008 - 18:17
DoctorStat Наблюдения для одной крысы зависимы, поэтому нельз... 30.10.2008 - 20:02
DoctorStat Цитата(nokh @ 30.10.2008 - 20:56) Ди... 31.10.2008 - 13:11
nokh Цитата(DoctorStat @ 31.10.2008 - 16... 31.10.2008 - 14:27
DoctorStat Цитата(nokh @ 31.10.2008 - 14:27) По... 31.10.2008 - 15:44
плав Цитата(DoctorStat @ 31.10.2008 - 15... 31.10.2008 - 17:52
Varta Спасибо за ценную информацию.
При пристальном изуч... 3.11.2008 - 08:28
плав Цитата(Varta @ 3.11.2008 - 08:28) 1)... 3.11.2008 - 12:36
Varta repeated measures ANOVA rat это все-таки анализ по... 3.11.2008 - 14:09
плав Цитата(Varta @ 3.11.2008 - 14:09) re... 3.11.2008 - 16:51
nokh Выше я неправильно назвал комплекс перекрестно-иер... 3.11.2008 - 18:18
плав А почему не хотите использовать процедуру MIXED?
M... 3.11.2008 - 20:36
nokh С mixed я еще не разобрался. Ваш код выдает резуль... 3.11.2008 - 22:06
плав Цитата(nokh @ 3.11.2008 - 22:06) При... 3.11.2008 - 22:11
Varta Таблица из Статистики 6.0
взяты группа - крыса -... 3.11.2008 - 22:55
nokh >плав. Ясно, значит это терминологические тонко... 4.11.2008 - 00:16
Varta Преобразование ситуацию не спасет, я для того, что... 4.11.2008 - 07:10
Varta плав
ЦитатаЕсли бы Вы потрудились посмотреть ссылк... 4.11.2008 - 08:39
плав Цитата(Varta @ 4.11.2008 - 08:39) пл... 4.11.2008 - 10:27
nokh Вы вступаете в неконструктивную полемику с человек... 4.11.2008 - 09:44
Varta ничего биологически значимого не произошло.
В данн... 4.11.2008 - 11:23
плав Мы вообще-то об одном говорим?
Цитата из первого п... 4.11.2008 - 18:00
Varta Когда я лет пять назад делала эксперимент в аспира... 6.11.2008 - 05:53
nokh Фактор может быть в анализе фиксированным или случ... 6.11.2008 - 07:47
Varta Фактически я не могу проанализировать данные по по... 7.11.2008 - 17:46
nokh С модулем Repeated Measures не работаю по уже указ... 7.11.2008 - 21:32
Varta Спасибо за объяснения, правда вникаю с трудом пока... 9.11.2008 - 10:23
nokh Цитата(Varta @ 9.11.2008 - 13:23) ti... 10.11.2008 - 19:26
Varta Цитата(nokh @ 10.11.2008 - 19:26) Но... 11.11.2008 - 13:31
DrgLena Позволю себе реплику, не по поводу статистических ... 11.11.2008 - 15:34
Varta Цитата(DrgLena @ 11.11.2008 - 15:34)... 12.11.2008 - 14:07![]() ![]() |