2013-11-29 1 views
2

У меня есть две группы лиц с повторными мерами (порядок мер не имеет значения [1,2] - это то же, что и [2,1]). Данные могут выглядеть (3 человека в группе, 6 Меры каждая):2 группы, 1 путь, повторная мера anova в MaLlab

groupA = [1 3 6 5 2 9; 2 5 3 4 5 8; 8 7 3 6 2 4]; 
groupB = [3 4 5 4 4 1; 2 8 4 2 1 2; 3 2 5 5 1 2]; 

Простым способом было бы сравнить 2 группы с помощью теста ranksum средних значений каждого человека:

meansA = mean(groupA, 2); % => [4.3 4.5 5.0] 
meansB = mean(groupB, 2); % => [3.5 3.2 3.0] 
[p, h] = ranksum(meansA, meansB) 

Однако этот тип анализа пренебрегает тем, что каждое из средних значений состоит из нескольких мер (и поэтому недооценивает значение).

Статистик сказал мне использовать «повторное измерение ANOVA», но ни одна из функций ANOVA MatLab, похоже, не делает то, что я хочу. Ближе всего, что я смог найти:

>> [p, atab] = anovan([1 3 6 5 2 9 2 5 3 4 5 8 8 7 3 6 2 4 3 4 5 4 4 1 2 8 4 2 1 2 3 2 5 5 1 2], {[zeros(1,18) ones(1,18)],[1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 6 6 6 6 6 6]}, 'varnames', {'individual', 'groupAorB'}, 'display', 'off') 
p = 
     NaN 
    0.9774 

Но это, кажется, не работает так, как я хочу его (значение NaN и нереалистичным р-значение). Я был бы рад за любые предложения о том, как выполнить соответствующий статистический тест по этим данным в MatLab.

ответ

0

Вы должны взглянуть на эту запись FileExchange, что имеет дело с односторонним движением повторяется меры ANOVA:

http://www.mathworks.com/matlabcentral/fileexchange/5576-rmaov1

Автор (Antonio Trujillo-Ortiz) сделал некоторые другие интересные записи для различных конструкций (2 и 3 ановас с повторными мерами).

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

+0

Хм этот код скорее проверяет, действует ли терапия с течением времени, и есть ли различия между людьми. В помощи говорится: «Девять испытуемых участвовали в релаксационной терапии с частотой недельной и недельной регистрации (w1 и w2) в течение нескольких недель и с несколькими неделями терапии (от w3 до w5). Интересно проверить, существуют ли различия в релаксационной терапии и внутри субъектов с уровнем значимости = 0,05. ". Поэтому он не проверяется на контрольную группу ... – user2952361

+0

Возможно, вас заинтересует другой дизайн - двухсторонняя anova, где Measure и Group могут быть вашими зависимыми переменными. – marsei

+0

Не могли бы вы привести пример, который должен выполнить это на описанных выше простых тестовых данных? – user2952361

0

nan означает, что модель, учитывающая INDIVIDUAL, учитывает всю дисперсию GROUP. В других словах, если вы поместили перехват для каждого INDIVIDUAL, попробуйте найти изменчивость из-за GROUP, у вас нет никакого различия. Модель «превышена».

Это потому, что вам действительно нужна модель смешанных эффектов - вы ищете эффект между субъектами. Для этого вам нужно сообщить MATLAB, что INDIVIDUAL вложен внутри GROUP. Таким образом, вы используете следующие параметры для anovan:

'random', INDIVIDUAL 
'nested', INDIVIDUAL_within_GROUP 

Сказав это, я не знаю, что ошибка ковариации предположения это делает - т.е. не принимает на себя диагональную матрицу ковариации?

Если вы хотите больше контролировать сделанные предположения, я предлагаю вам изучить NLMEFIT из панели инструментов статистики. Это позволяет использовать модели смешанных эффектов, определяющие ковариацию.