Как получить четыре графика, сгенерированные из регрессионного анализа в R? - PullRequest
Новая
       21

Как получить четыре графика, сгенерированные из регрессионного анализа в R?

0 голосов
/ 22 мая 2018

Как получить четыре сгенерированных графика (один отделенный от другого / один за раз) из регрессионного анализа, не используя консоль / приглашение (Hit <Return> to see next plot:)? 

if(!require("truncnorm")) install.packages("truncnorm") ; library(truncnorm)


FACTORA <- c(0,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0)
FACTORB <- c(0,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0)
FACTORC <- c(0,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0)
FACTORD <- c(0,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0)
FACTORE <- c(0,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0)

planDOE <- data.frame(FACTORA, 
                  FACTORB, 
                  FACTORC, 
                  FACTORD, 
                  FACTORE, 
                  rtruncnorm(n=32, a=0, mean=61.16, sd=31.32))

rg <- lm(planDOE)
plot(rg)

Я пыталсяследующие способы и не увенчались успехом (с использованием true или false в командах аргумента ask ask). 

plot(rg)
par(ask=TRUE)

# 

par(ask=FALSE)
plot(rg)

# 

plot(rg, ask = FALSE)

# 

par(mfrow=c(2,2))
plot(rg)
# In one image it can put together the graphic information.

# 

png('test2.png', units="in", width=13, height=7, res=300)
plot(rg)
dev.off()

В этом последнем способе я понимаю, что можно получить график отдельно, но только сохранить график остатков.

Ответ: 

plot(rg, which = 1)
plot(rg, which = 2)
#...

Ответы [ 2 ]

0 голосов
/ 22 мая 2018

Вы действительно можете увидеть это, если вы сделали ?plot.lm

mod <- lm(1:10 ~ rnorm(10))

Теперь вы можете построить каждый из графиков отдельно, указав аргумент which в plot.lm.Например, 

png('plot_1.png', units="in", width=13, height=7, res=300)
 plot(mod, which = 1)
dev.off()

png('plot_2.png', units="in", width=13, height=7, res=300)
 plot(mod, which = 2)
dev.off()

plot(mod, which = 3)
plot(mod, which = 5)

Они будут отображать фигуру, которую вы обычно получаете, если вы используете plot(mod).

Естественно, вы можете сохранять после каждого графика или только те участки, которые предпочитаете.

ПРИМЕЧАНИЕ : я не использовал which = 4.

0 голосов
/ 22 мая 2018

установить номинал (mfrow = params).используя mtcars, вот пример.После этого вы можете сохранить график в виде файла .png 

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

par(mfrow=c(2,2))
plot(rg)
...