Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

> Несколько биоптатов с крысы - как анализировать?
Varta
сообщение 30.10.2008 - 10:00
Сообщение #1





Группа: Пользователи
Сообщений: 20
Регистрация: 29.10.2008
Из: Челябинск
Пользователь №: 5459



Допустим имеем 2 группы по 6 крыс в каждой. У каждой крысы взято по 3 биоптата. Как анализировать?
Вариант 1 расчитать среднее по 3 биоптатам для каждой крысы и сравнить 6 и 6. (скажем, сделать Манна-Уитни)
Вариант 2 сравнить 18 биоптатов одной группы с 18 другой.
Правильно конечно делать первый вариант, поскольку при втором увеличивается вклад индивидуальной вариабельности и мы недооценим риск получить различия обусловленные случайной комбинацией крыс. Но чувствительность в первом варианте достаточно мала а хочется все-таки больше.
Есть какие-нибудь ещё варианты?

Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
 
Открыть тему
Ответов
nokh
сообщение 6.11.2008 - 07:47
Сообщение #2





Группа: Пользователи
Сообщений: 1219
Регистрация: 13.01.2008
Из: Челябинск
Пользователь №: 4704



Фактор может быть в анализе фиксированным или случайным исходя из целей эксперимента и контекста экспериментального плана. Хотя в данном случае выбор полей зрения, возможно, и не является абсолютно случайным в строгом понимании этого термина, но это - неподконтрольный эспериментатору фактор. Если бы анализировались поля зрения в центре и на периферии среза - фактор был бы фиксированным. Также проход препарата под микроскопом может осуществляться по разным схемам и т.о. случайная изменчивость опять-таки оказывается ограничена (фиксирована) какими-то рамками. В вашей же экспериментальной схеме фактор "поле зрения" однозначно является случайным. По поводу срезов. Грамотный анализ, безусловно, должен учитывать все особенности эксперимента и т.о. включать все ступени иерархии. Но давайте посмотрим приведет ли включение срезов к уменьшению P для эффекта "Группы". Для этого нужно знать структуру математических ожидаемых средних квадратов. Исходя из этой структуры проводится оценка значимости эффектов. В рассчитанном вами примере фактор «Группы» оценивался на фактор «Крысы внутри группы». Это можно увидеть по тому, на какой MS делился MS оцениваемого фактора. Фактор «Крысы внутри группы» оценивался на фактор «Поле зрения внутри Крысы», который выступал здесь Ошибкой всего комплекса. Если сейчас мы введем (а это нужно сделать) случайный фактор «Срез внутри Крысы», то MS «Крысы внутри группы» будет делиться на него, а сам он на MS «Поле зрения внутри среза». Таким образом, изменчивость, которая раньше уходила в Ошибку окажется разбитой на 2 части: изменчивость межу срезами и новую Ошибку, которая будет меньше прежней. Это приведет к тому, что эффект «Крысы внутри группы» будет оценен корректнее и точнее, P для него еще уменьшится, но это не значит, что уменьшится сам эффект - SS для этого фактора останется прежней. Введение фактора «Срез внутри Крысы» не затронет отношения средних квадратов MS для фактора «Группы» - он будет прежним. Т.о. иерархия никак не позволяет нам снизить индивидуальную изменчивость крыс, она по-прежнему остается высокой. Что имеем - то имеем. Выход - увеличение числа степеней свободы для фактора «Крыса». Тогда его MS уменьшится и величина F-критерия для «Групп» увеличится. Поскольку эксперимент проведен и количества животных не увеличить - давайте увеличивать степени свободы введением фактора «Время». Возвращаемся на старые позиции: попробуйте написать модель для такого анализа, а я проверю и напишу как это сделать в Statistica, хотя сильно подозреваю что она сдастся - предложить анализировать в другом модуле, что некорректно. (В SPSS принципиальным образом ничего не изменится). Думаю, с иерархическими комплексами мы разобрались smile.gif, давайте двинемся к перекрестно-иерархическим. (по терминологии Монтгомери это планы с группированными и пересекающимися факторами - стр. 257).

Сообщение отредактировал nokh - 6.11.2008 - 08:25
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Сообщений в этой теме
- 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
- - nokh   Дисперсионный анализ с повторными измерениями разр...   30.10.2008 - 20:56
|- - 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


Добавить ответ в эту темуОткрыть тему