воскресенье, 17 октября 2010 г.

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

<< Главная страница