В настоящее время я использую пакет statnet для ERGMS в R, и я установил модель и проверил ее вырожденность (используя mcmc diagnostics()).Теперь я хочу оценить качество соответствия:
fit1 <- ergm(ExampleData ~ edges + mutual + ctriple, set.seed(123))
gof_fit1 <- gof(fit1)
summary(gof_fit1)
jpeg("GOF model1.jpg", width = 600, height = 600)
par(mfrow=c(2,3))
plot(gof_fit1)
dev.off()
Проблема в том, что plot(gof_fit1) выдает пять графиков (один для степени, степени, общих партнеров по краю, минимального геодезического расстояния и моделистатистика каждого).Поэтому, если я открою «GOF model1.jpg», он покажет мне только последний график (статистику модели).
Обычно я использую следующие команды для сохранения нескольких графиков в одном файле jpg, и это работает, как каждый графикзаписано так:
jpeg("Example1.jpg")
par(mfrow=c(2,2))
plot(w)
plot(x)
plot(y)
plot(z)
dev.off()
Однако, в моем примере выше, одна команда (plot(gof_fit1)) создает несколько графиков.Как сохранить все пять графиков, сгенерированных командой plot(gof_fit1), в одном файле jpg?
Заранее большое спасибо!