Здравствуйте, гость ( Вход | Регистрация )
30.10.2008 - 10:00
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 20 Регистрация: 29.10.2008 Из: Челябинск Пользователь №: 5459 |
Допустим имеем 2 группы по 6 крыс в каждой. У каждой крысы взято по 3 биоптата. Как анализировать?
Вариант 1 расчитать среднее по 3 биоптатам для каждой крысы и сравнить 6 и 6. (скажем, сделать Манна-Уитни) Вариант 2 сравнить 18 биоптатов одной группы с 18 другой. Правильно конечно делать первый вариант, поскольку при втором увеличивается вклад индивидуальной вариабельности и мы недооценим риск получить различия обусловленные случайной комбинацией крыс. Но чувствительность в первом варианте достаточно мала а хочется все-таки больше. Есть какие-нибудь ещё варианты? |
|
|
![]() |
![]() |
![]() |
6.11.2008 - 07:47
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 1219 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
Фактор может быть в анализе фиксированным или случайным исходя из целей эксперимента и контекста экспериментального плана. Хотя в данном случае выбор полей зрения, возможно, и не является абсолютно случайным в строгом понимании этого термина, но это - неподконтрольный эспериментатору фактор. Если бы анализировались поля зрения в центре и на периферии среза - фактор был бы фиксированным. Также проход препарата под микроскопом может осуществляться по разным схемам и т.о. случайная изменчивость опять-таки оказывается ограничена (фиксирована) какими-то рамками. В вашей же экспериментальной схеме фактор "поле зрения" однозначно является случайным. По поводу срезов. Грамотный анализ, безусловно, должен учитывать все особенности эксперимента и т.о. включать все ступени иерархии. Но давайте посмотрим приведет ли включение срезов к уменьшению P для эффекта "Группы". Для этого нужно знать структуру математических ожидаемых средних квадратов. Исходя из этой структуры проводится оценка значимости эффектов. В рассчитанном вами примере фактор «Группы» оценивался на фактор «Крысы внутри группы». Это можно увидеть по тому, на какой MS делился MS оцениваемого фактора. Фактор «Крысы внутри группы» оценивался на фактор «Поле зрения внутри Крысы», который выступал здесь Ошибкой всего комплекса. Если сейчас мы введем (а это нужно сделать) случайный фактор «Срез внутри Крысы», то MS «Крысы внутри группы» будет делиться на него, а сам он на MS «Поле зрения внутри среза». Таким образом, изменчивость, которая раньше уходила в Ошибку окажется разбитой на 2 части: изменчивость межу срезами и новую Ошибку, которая будет меньше прежней. Это приведет к тому, что эффект «Крысы внутри группы» будет оценен корректнее и точнее, P для него еще уменьшится, но это не значит, что уменьшится сам эффект - SS для этого фактора останется прежней. Введение фактора «Срез внутри Крысы» не затронет отношения средних квадратов MS для фактора «Группы» - он будет прежним. Т.о. иерархия никак не позволяет нам снизить индивидуальную изменчивость крыс, она по-прежнему остается высокой. Что имеем - то имеем. Выход - увеличение числа степеней свободы для фактора «Крыса». Тогда его MS уменьшится и величина F-критерия для «Групп» увеличится. Поскольку эксперимент проведен и количества животных не увеличить - давайте увеличивать степени свободы введением фактора «Время». Возвращаемся на старые позиции: попробуйте написать модель для такого анализа, а я проверю и напишу как это сделать в Statistica, хотя сильно подозреваю что она сдастся - предложить анализировать в другом модуле, что некорректно. (В SPSS принципиальным образом ничего не изменится). Думаю, с иерархическими комплексами мы разобрались
Сообщение отредактировал 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
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![]() ![]() |