Двойной цикл для врезанных графиков с использованием ggplot в R - PullRequest
Новая
       28

Двойной цикл для врезанных графиков с использованием ggplot в R

1 голос
/ 09 января 2020

Я пытаюсь запустить al oop для следующего кода: 

library(tidyverse)
library(cowplot)

#Setup of Variables
N <- 199
K <- N+1
x <- rep(0,N)
x[1] <- 0.5
g <- 3.84
time <- c(1:K)

for (t in 1:N){
  x[t+1] = g*x[t]*(1-x[t])
}

A <- data.frame(time,x)

#create separate plots 1 & 2 --> combine into plot3 using ggdraw()
plot1 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="blue") + xlim(0,150) + ylab("Fish Stock") + theme_classic()

plot2 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="forestgreen") + xlim(8,15) + theme(legend.position="none", axis.title.x=element_blank(), axis.title.y=element_blank(), panel.border=element_blank(), panel.grid.major=element_blank(), panel.grid.minor=element_blank(), axis.line=element_line(colour="black"), panel.background=element_rect(fill="white",colour="white"))

plot3 <- 
  ggdraw() +
  draw_plot(plot1) + 
  draw_plot(plot2, .68, .68, .3, .3) #insetplot, x, y, width, height

plot3

, который создает следующие графики c:

enter image description here

То, что я хочу сделать с этим, это вложить этот l oop в другой, который проходит через вектор значений для g. Моя попытка этого кода выглядит следующим образом: 

G <- c(2.7, 2.9, 3.0, 3.5, 3.82, 3.83, 3.84, 3.85)

#Loop --> creation of x, dataframe, insetplots
for (g in G) {
  for (t in 1:N){
    x[t+1] = g*x[t]*(1-x[t])
  }

  A <- data.frame(time,x)

  plot1 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="blue") + xlim(0,150) + ylab("Fish Stock") + theme_classic()
  plot2 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="forestgreen") + xlim(8,15) + theme(legend.position="none", axis.title.x=element_blank(), axis.title.y=element_blank(), panel.border=element_blank(), panel.grid.major=element_blank(), panel.grid.minor=element_blank(), axis.line=element_line(colour="black"), panel.background=element_rect(fill="white",colour="white"))
  plot3 <- 
  ggdraw() +
  draw_plot(plot1) + 
  draw_plot(plot2, .68, .68, .3, .3) #insetplot, x, y, width, height

  png(file = paste0("Fish-Stock-Inset_", g, ".png"))
  plot3
  dev.off()
}

Однако я получаю серию пустых изображений. Сохранение графиков происходит правильно, так как я получаю правильно названные файлы, но сами графики отсутствуют. Я почти уверен, что проблема связана с моим размещением A <- data.frame(time,x) после начального for (t in 1:N) l oop.

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

1 Ответ

1 голос
/ 09 января 2020

Ваша проблема возникает из-за того, что plot3 оценивается только в конце для l oop (так называемая отложенная оценка в R).

Попробуйте это ниже, где я использую ggsave для рендеринга графика и примечание для объявления вашего x внутри l oop: 

G <- c(2.7, 2.9, 3.0, 3.5, 3.82, 3.83, 3.84, 3.85)

#Loop --> creation of x, dataframe, insetplots
for (g in G) {
  x <- rep(0,N)
  x[1] <- 0.5
  for (t in 1:N){
    x[t+1] = g*x[t]*(1-x[t])
  }

  A <- data.frame(time,x)

  plot1 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="blue") + xlim(0,150) + ylab("Fish Stock") + theme_classic()
  plot2 <- ggplot(data=A, aes(x=time, y=x)) + geom_line(color="forestgreen") + xlim(8,15) + 
  theme(legend.position="none", axis.title.x=element_blank(), 
  axis.title.y=element_blank(), panel.border=element_blank(), 
  panel.grid.major=element_blank(), panel.grid.minor=element_blank(), 
  axis.line=element_line(colour="black"), 
  panel.background=element_rect(fill="white",colour="white"))
  plot3 <- 
    ggdraw() +
    draw_plot(plot1) + 
    draw_plot(plot2, .68, .68, .3, .3) #insetplot, x, y, width, height

  ggsave(plot3,file = paste0("Fish-Stock-Inset_", g, ".png"))
}
...