R: В anova.lm (г): F-тесты ANOVA по существу идеально подходят ненадежно - PullRequest
4 голосов
/ 18 декабря 2011

Я соединяю интерактивные руководства со старым текстом, чтобы выучить R (стр. 182 - http://cran.r -project.org / doc / contrib / Faraway-PRA.pdf ).Когда я использую данные из пакета из R (как в примерах учебника), это не проблема.Однако, когда я использую данные из своего текста, я всегда заканчиваю без F-значения и предупреждения.

Взгляните:

данные в data.frame:

car.noise <- data.frame( speed = c("idle", "0-60mph", "over 60"), chrysler = c(41,65,76), 
bmw = c(45,67,72), ford = c(44,66,76), chevy = c(45,66,77), subaru = c(46,76,64))

проверьте data.frame:

car.noise
    speed chrysler bmw ford chevy subaru
1    idle       41  45   44    45     46
2 0-60mph       65  67   66    66     76
3 over 60       76  72   76    77     64

данные расплава.frame:

mcar.noise<- melt(car.noise, id.var="speed")

проверка расплавленных данных.

1: я неправильно кодирую 2: текстовые примеры слишком «идеальны» для подбора, так как они пытаются показать четкий пример

1 Ответ

7 голосов
/ 18 декабря 2011

Вы пытаетесь подобрать модель, которая дает отдельное среднее значение для каждой комбинации переменной * скорости. С имеющимися у вас данными это означает, что у вас вообще нет репликации. Это все равно что пытаться сравнить две группы, когда у вас есть только одно значение из каждой группы.

Если вы посмотрите на строку «Остатки» в своей таблице анова, вы должны заметить, что у вас там нет степеней свободы, и ваши суммы квадратов также равны 0. Вы можете попытаться приспособить модель без взаимодействия, если считаете, что это уместно, но у вас недостаточно данных для согласования модели с взаимодействием.

...