Комбинирование leveragePlots()

Question

Доброе утро,

новичок в R. Я пытаюсь объединить 9 участков, устраивая его 3 х 3.

У меня есть а) 7 leveragePlots() синтаксис, создает 1 график каждый, и b) 8-е плечоPlots() генерирует 2 графика, так как это две переменные линейные модели.

При использовании приведенных ниже первых 7 графиков, где расположены хорошо, но последний не вписывается в сетку 3x3, а открывает новое окно с двумя графиками.

Минимальный набор данных (sample.txt):

X Y 
276.67 5.20 
285.00 5.17 
296.67 4.96 
307.50 5.26 
307.50 5.11 
262.50 5.27 
278.33 5.20 
282.50 5.16 
267.50 5.17 
250.00 5.20 
275.00 5.18 
274.17 5.09 
311.67 5.03 
294.17 5.21 
279.17 5.29 
300.00 5.14 
292.50 5.09 
297.50 5.16 
290.00 4.99 
272.50 5.23 
289.17 5.32 
279.17 5.21 
325.83 5.14 
292.50 5.22 
255.83 5.15 
290.83 5.09 
283.33 5.09 
284.17 5.07 
298.33 5.27 
269.17 5.19 

пробовали использовать пар (mfrow = С (3,3)), а такжемакета (матрицу (C (1,2,3,4 , 5,6,7,8,8), 3, 3, byrow = TRUE)). Это не работает. Открывает новое окно для последнего синтаксиса leveragePlots(). Для того, чтобы повторить эти коды ниже,

library(car) 
sampledata <- read.table("D:/sample.txt", header=TRUE) 
y.1 <- sampledata$Y 
x.1 <- sampledata$X 
x.2 <- (sampledata$X)^2 

par(mfrow=c(3,3)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1 + x.2)) 

вы знаете, как я могу заставить заключительные 2 участка ли быть размещены в одной и той же компоновке 3x3? Если да, пожалуйста, посоветуйте мне.

В настоящее время я читаю на макете() в http://rfunction.com/archives/1538, не уверен, как использовать это для достижения вышеуказанного.

Спасибо за ваше руководство.

Answer 1

Исправленное решение, теперь с использованием предоставленных данных.

par(mfrow=c(3,3)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1)) 
leveragePlots(lm(y.1 ~ x.1 + x.2), layout = NA)