Продолжаю ковыряться в R
Я, наверное, мазохист, но я продолжаю разбираться с R. Вот, даже осилил провести Repeated Measures Anova:
1. Запустим Rcmdr
 library (Rcmdr)
 2. Загрузим файл для обработки:
 e1 <- read.table («http://volokhonsky.ru/data/ex01.csv», header=TRUE)
 3. Сформируем массив данных из двух зависимых переменных:
 dvm <- with (e1, cbind (o1,o2))
 4. Опишем многомерную линейную модель:
 mlm1 <- lm (dvm ~ 1)
 5. Сформируем внутригрупповой фактор:
 rfactor <- factor (c («o1»,"o2"))
 6. Загрузим библиотеку car
 library (car)
 7. Выполним дисперсионный анализ:
mlm1.aov <- Anova (mlm1, idata = data.frame (rfactor), idesign = ~rfactor, type="III")
 summary (mlm1.aov, multivariate=FALSE)
Теперь осталось понять, как в этот синтаксис запихнуть межгрупповой фактор. Остальные варианты синтаксиса, которые я находил, позволяли запихнуть межгрупповой фактор, но они скорее выглядели как раз наоборот — как попытка протащить повторные измерения через одномерный ANOVA...
То ли его на четвёртом шаге запихивать, то ли в конце... Ща буду копать мануалы, почувствую себя манулом.
Оригинальная запись в блоге Волохонского: 
Продолжаю ковыряться в R.

Комментарии: 0:
Отправить комментарий
Подпишитесь на каналы Комментарии к сообщению [Atom]
<< Главная страница