2015-08-03 2 views
3

Я сделал несколько линейных регрессий. Я попытался построить его с помощью этой команды.Регрессионный анализ ggplot2

layout(matrix(c(1,2,3,4),2,2)) 
plot(fit_ec_urban_franchise) 

После того, что у меня 'остатки vs. оборудованная' 4 'участков, масштаб против местоположение', 'нормальной С-С' и 'Разности vs.leverage'.

Возможно ли использовать ggplot2, чтобы поместить 4 участка в один?

+0

ищет пакеты с 'autoplot.lm' функций, таких как http://www.inside-r.org/packages/cran/eeptools/docs/autoplot.lm – scoa

+0

Спасибо за Помогите. Вы знаете другой способ сделать это? –

ответ

3

Решение должно использовать ggplot2::fortify. Вот код, который вы можете найти на его странице справки ?fortify. Я добавляю gridExtra, чтобы организовать 4 участка вместе.

library(ggplot2) 
library(gridExtra) 

mod <- lm(mpg ~ wt + cyl, data = mtcars) 

p1 <- qplot(.fitted, .resid, data = mod) + 
    geom_hline(yintercept = 0) + 
    geom_smooth(se = FALSE) 

p2 <- qplot(sample =.stdresid, data = mod, stat = "qq") + geom_abline() 

p3 <- qplot(.fitted, sqrt(abs(.stdresid)), data = mod) + geom_smooth(se = FALSE) 

p4 <- qplot(.hat, .stdresid, data = mod) + geom_smooth(se = FALSE) 


grid.arrange(p1,p2,p3,p4) 

enter image description here

+0

Как я могу добавить главный заголовок? Можно ли использовать 'top'. –

+0

+ ggtitle («ваше название здесь») – scoa

+0

Спасибо, что вам поможет. –

Смежные вопросы