Лучший способ вывести эффекты взаимодействия на печать из линейной модели

Cmd +. , Profiles, Keys и выбирают Left option Key acts as: +Esc.

Работы для меня в emacs, хотя я не уверен, будет ли он иметь какие-либо другие нежелательные эффекты.

15
задан Ben Bolker 8 February 2015 в 17:00
поделиться

2 ответа

В пакете эффектов есть хорошие методы построения графиков для визуализации предсказанных значений регрессий.

thedata<-data.frame(x=rnorm(20),f=rep(c("level1","level2"),10))
thedata$y<-rnorm(20,,3)+thedata$x*(as.numeric(thedata$f)-1)

library(effects)
model.lm <- lm(formula=y ~ x*f,data=thedata)
plot(effect(term="x:f",mod=model.lm,default.levels=20),multiline=TRUE)
17
ответ дан 1 December 2019 в 03:24
поделиться

Ха - все еще пытаюсь осознать expand.grid () . Просто для сравнения, вот как я бы это сделал (используя ggplot2):

thedata <- data.frame(predict(thelm), thelm$model$x, thelm$model$f)

ggplot(thedata, aes(x = x, y = yhat, group = f, color = f)) + geom_line()

Я думаю, логика ggplot () довольно интуитивно понятна - группируйте и раскрашивайте линии по f. С увеличением числа групп необходимость указывать уровень для каждой становится все более полезной.

3
ответ дан 1 December 2019 в 03:24
поделиться
Другие вопросы по тегам:

Похожие вопросы: