Продолжаю ковыряться в 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]
<< Главная страница