Print ggplot2 Plot within for-Loop in R (Example)
In this article you’ll learn how to draw ggplot2 plots within a for-loop in the R programming language.
The post looks as follows:
- Introducing Example Data
- Example: Plotting ggplot2 Plots within Loop
- Video, Further Resources & Summary
You’re here for the answer, so let’s get straight to the programming part…
Introducing Example Data
We will use the following data as basement for the example of this R programming tutorial:
set.seed(159159) # Create example data data <- data.frame(x = 1:100, y1 = rnorm(100), y2 = rnorm(100), y3 = rnorm(100)) |
Our example data contains 100 rows and four columns. The variable x is ranging from 1 to 100 and the columns y1, y2, and y2 contain 100 random normally distributed numbers each.
We also need to install and load the ggplot2 add-on package to RStudio:
install.packages("ggplot2") # Install and load ggplot2 library("ggplot2") |
Now, we can draw a basic ggplot2 graph as follows:
ggplot(data, aes(x = x, y = y1)) + # Basic ggplot2 plot of x & y1 geom_point() |
Figure 1: Basic Scatterplot Created by ggplot2 Package.
In Figure 1, you can see the result of the previous R code: A scatterplot of x and y1.
Now, let’s assume we want to create a ggplot2 plot of each combination of x and y1, y2, and y3 respectively. In such a scenario, we may want to use a for-loop:
for(i in 2:ncol(data)) { # ggplot within for-loop ggplot(data, aes(x = x, y = data[ , i])) + geom_point() Sys.sleep(2) } |
Damn! No graphic was returned!
So why does ggplot2 not work if it is inside a for-loop although it works outside of it?!
That’s what I’m going to show you next…
Example: Plotting ggplot2 Plots within Loop
If we want to draw a plot within a for-loop, we need to wrap the print function around the R code creating the plot. Have a look at the following R syntax:
for(i in 2:ncol(data)) { # Printing ggplot within for-loop print(ggplot(data, aes(x = x, y = data[ , i])) + geom_point()) Sys.sleep(2) } |
Figure 2: Showing ggplot2 Plots within for-Loop using print() Function.
After running the previous R code, you will see three ggplot2 graphs popping up at the bottom right of RStudio with a delay of 2 seconds.
Looks good!
Video, Further Resources & Summary
Would you like to know more about the ggplot2 package in R? Then you could watch the following video of my YouTube channel. In the video, I explain the R codes of this article in RStudio:
Please accept YouTube cookies to play this video. By accepting you will be accessing content from YouTube, a service provided by an external third party.
If you accept this notice, your choice will be saved and the page will refresh.
Furthermore, you might have a look at the related articles of this website. A selection of articles can be found here.
This tutorial illustrated how to use the ggplot2 package within while- or for-loops in R. Don’t hesitate to let me know in the comments, in case you have additional questions. In addition, please subscribe to my email newsletter in order to receive updates on the newest articles.
I’m Joachim Schork. On this website, I provide statistics tutorials as well as code in Python and R programming.
43 Comments. Leave new
Hi , Thanks for the tutorial on ggplot loop! Is there a way to dynamically assign the variable names as Y-Axis labels?
Hi Anoop,
Thank you for the comment.
Please have a look at the following tutorials:
https://statisticsglobe.com/name-variables-in-for-loop-dynamically-in-r
https://statisticsglobe.com/assign-function-in-r
They explain how to assign names dynamically. You can also use this logic to assign Y-axis labels dynamically.
Regards,
Joachim
Hi Joachim
This was an excellent description of how to create the for loop for ggplotting.
I have a follow up question related to Anoop’s question that I hope you can help with, because I don’t seem to be able to work it out from the links you have shared.
I would like to assign my column names to be the x labels (my y variable is always the same), but when I try with the code that I put together from the links you shared I only ever get the column number and not the column name. I have tried passing a character list of column names, however it only ever uses the first name in the list.
Can you help clear this up a little for me please, as I am not sure where I am going wrong.
Hey,
Thanks a lot for the kind comment, it’s great to hear that you find this tutorial helpful!
Regarding your question, could you share your code and illustrate how your data set looks like?
Regards,
Joachim
Thank you.
I have a dataframe of 36 observations of 17 variables, but in this iteration of what I am doing I only need to plot column 2 (Total_Erosion) against columns 8-17.
Am using the code below at the moment:
for (i in 8:ncol(Bank1Variables)) {
print(ggplot(data = Bank1Variables, aes_string(x = Bank1Variables[ , i], y = Total_Erosion*-1)) +
geom_point() +
labs(x = assign(paste0(” “, i), value = i), y = names(Bank1Variables[2])) +
geom_smooth(method = “lm”, se = FALSE))
}
Which produces my plots beautifully, and names my Y axis correctly but my X axes are only ever named based on their row number and not their name.
I have also created a character vector of the column names that I am using
w <- c("Mean_Stage", "Max_Stage", "Peaks_AboveQ10", "Peaks_AboveQ50", "Hours_AboveQ10", "Hours_AboveQ50")
But when I include the vector in the code as follows:
for (i in 8:ncol(Bank1Variables)) {
print(ggplot(data = Bank1Variables, aes_string(x = Bank1Variables[ , i], y = Total_Erosion*-1)) +
geom_point() +
labs(x = assign(paste0(" ", w), value = w), y = names(Bank1Variables[2])) +
geom_smooth(method = "lm", se = FALSE))
}
It only ever uses the first name of Mean_Stage and labels all plots the same and then if i replace the w with w[i] it labels them all as NA
It's likely to be something extremely simply and silly, but I have been working on it for a while now and can't get it to work.
Thank you for your time. I really appreciate it!
Thanks for the clarifications! You may simply add the xlab and ylab functions to your ggplot2 code.
Have a look at the following example. It shows how to modify the yaxis label dynamically:
Let me know in case you need further help! 🙂
Regards,
Joachim
Hi,
I am getting following error.
Error in data[, i] : object of type ‘closure’ is not subsettable.
Any solution please? Thanks
Hi Moh,
It seems like you are trying to subset something that is not a data object.
Have a look here for more details: https://statisticsglobe.com/error-object-of-type-closure-is-not-subsettable-in-r
I hope that helps!
Joachim
Hi and thanks for the tutorial! Is there a way to show all graphs at the same time? In R, without using a pdf (it’s only 9 graphs for me)
Hi Ale,
Thank you for the comment!
You could use the patchwork package to draw your plots in RStudio side-by-side: https://statisticsglobe.com/patchwork-r-package
Regards
Joachim
Hi Joachim
This is great. Been having this issue a few times. But why does it work? What does print() do that is needed here?
Hi Christer,
Thank you for the nice comment, glad it helped!
You can read more about the reasons in this article, specifically: “An interesting side-effect of this is that ggplots are only rendered when explicitly print()ed/plot()ed within a loop, as only the last return value in a sequence of calls gets its print method invoked.”
Regards
Joachim
THis website was very helpful for me. I have always tried several times to learn creating loops in R and the books have not been good. But, your website made me like to go through all your examples and even try some other new options just to see how the outputs changes.
Thank you very much!!!!!!!!!!!!!!!!!!!!!!!!!
Thanks a lot for this wonderful feedback Richard! I’m very glad to hear that my tutorials helped to improve your R programming skills 🙂
Regards
Joachim
I spent days on this problem, but with your help it will be a big time-saver in the future! I especially like how you showed the initial lack of result (“Damn! No graphic…”), just as I experienced it (many times). I can’t thank you enough.
Thanks a lot for the kind comment Littleag! Glad it helped! 🙂
Regards
Joachim
HI ！
i fellowed your code ,and it work fine, when i try to save each plot of them,i meet so error.
———————————————————————————————————————————
for (i in 2:ncol(data)) {
p=( ggplot(data=data) + geom_bar(aes(x=tissue,y=data[,i]),stat=”identity”)+theme_classic())
plot_list[[i]] = p}
———————————————————————————————————————————
Rstudio did not repot error,but the plot was not the plot when i use code as you descripted
————————————————————————————————————————
for (i in 2:ncol(data)) {
print( ggplot(data=data) + geom_bar(aes(x=tissue,y=data[,i]),stat=”identity”)+theme_classic())
Sys.sleep(2)}
——————————–
thanks u
Hey Joe,
You mean you want to export these plots as PDF or PNG? I’m not certain if I get your question correctly.
Regards
Joachim
I think that’s what the poster was asking, and it’s what I really want to know as well.
Hey Mike,
You can export all plots as PDFs as shown below:
I hope that helps!
Joachim
Hi Joachim, very useful script. How ever I have a question could you please answer this.
for(i in 2:ncol(data))
Can you please explain why you have put 2. Thanks
Hey Sameer,
Thank you, glad it is useful!
The 2:ncol is used because the first column in our data frame is the x-variable, and we do only want to loop over the y-variables (i.e. positions 2:ncol(data)).
Regards,
Joachim
Hi Joachim,
I am using for-loops to produce multiple graphs in RMarkdown. I assign them to a list so that they are able to print or to be exported to html or word document. My problem is that i would like to use for-loops then export results of each single code chunk to a separate document i.e. 5 RMarkdown code chunks = 5 separate word documents. I can do this with officer package when the plots are not produced by for-loops or assigned to a list but if i do this for plots on a list produced by for-loop, i only get 1 plot only.
Is there a way to do this in a simple way?
Hey Joshua,
Unfortunately, I’m not an expert on RMarkdown. However, could you share your code? I could have a look and try to fx it.
Regards,
Joachim
Hi Joachim
For some reason I couldn’t reply to your last message on our thread but I just wanted to say that your solution worked perfectly.
I am so grateful for your help! Thank you for your time and effort.
It is much appreciated!
BB
This is really great to hear, glad it was useful! 🙂
Regards,
Joachim
Thanks for the information – really helpful to learn so far!
I’d really like to be able to do the same steps, but have each plot combine into the same plot – is there a way to adapt this to work?
Specifically, I am trying to plot several normal distributions from a data.table onto the same plot. Any advice appreciated!
Hey Sam,
Thanks for the kind comment, glad it’s useful.
Regarding your question, did you already have a look at the patchwork package? It’s a great package for combining multiple plots.
Regards,
Joachim
Hi,
Thank you very much for your video. I tried to apply your code to my data with the following script :
for (i in 2:ncol(df)) {
print(ggplot(df, aes(x = Date_Time, y = df[ , i]))+
geom_point())
Sys.sleep(2)
}
But it’s sending an error message :
Error in `[.data.table`(df, , i) :
j (the 2nd argument inside […]) is a single symbol but column name ‘i’ is not found. Perhaps you intended DT[, ..i]. This difference to data.frame is deliberate and explained in FAQ 1.1.
I don’t understand why it’s not working, if it’s working on your session.
Hey Alice,
It seems like you are working with a data.table instead of a data.frame. Could you please try if it works when you convert your data.table to the data.frame class first? See here for more details.
Regards,
Joachim
Hi Joachim,
Thanks a lot for there instructions. I’m still running into a small problem with my own data.
I would like to create multiple plots, one for each of my participants.
So my dataframe looks like this:
x= subjectID
y=time point
z=concentration of a biomarker.
I want to create a plot of y vs. z, for every x. So i need to loop over x, but I cant figure out how to incorporate a third column in this code..
Hi Niels,
Thank you for the kind comment.
Please have a look at the following example code. It reproduces your question:
Regards,
Joachim
Hi Joachim,
thanks so much for your answer. I was able to adjust a little to fit my data and everything worked fine!
However, after this I tried to incorporate an earlier answer from you in order to save my plots as PDF. But this resulted in the following error message:
Error in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
invalid font type
I’ve tried several thing to fix this, all not working sadly. Strange thing now is when I remove all the PDF lines, my original code is also returning the same error…
Do you have any suggestions?
My complete code looks like this:
for(i in 1:length(unique(TropT$Record_ID))) {
print(ggplot(TropT[TropT$Record_ID == unique(TropT$Record_ID)[i], ], aes(x = phase, y = Lab_TropT, group=1)) +
ggtitle(paste(TropT$Record_ID == unique(TropT$Record_ID)[i])) +
theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 40,
hjust = 0))+
ylab(“Troponine T (ng/L)”) + theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 20))+
xlab(“Meetmoment”) + theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 20))+
geom_line(colour = “#333333″)+
geom_point(shape=21, color=”black”, fill=”#69b3a2″, size=6)+
scale_y_continuous(
limits = c(0, 50),
breaks = seq(from = 0, to = 50, by = 5))+
theme_ipsum())
Sys.sleep(2)
}
Previous question includes the wrong code..
This is the correct version:
for(i in 1:length(unique(TropT$Record_ID))) {
pdf(paste0(“TropT plots.pdf”))
print(ggplot(TropT[TropT$Record_ID == unique(TropT$Record_ID)[i], ], aes(x = phase, y = Lab_TropT, group=1)) +
ggtitle(paste(TropT$Record_ID[i])) +
theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 40,
hjust = 0))+
ylab(“Troponine T (ng/L)”) + theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 20))+
xlab(“Meetmoment”) + theme(plot.title = element_text(
face = “bold”,
colour = “#333333”,
size = 20))+
geom_line(colour = “#333333″)+
geom_point(shape=21, color=”black”, fill=”#69b3a2″, size=6)+
theme_ipsum())
dev.off()
}
Hi Nils,
I’m sorry for the late reply, unfortunately I couldn’t respond earlier. Do you still need help with this error message?
Regards,
Joachim
Hi Joachim,
I would like to follow up on my previous questions with an new problem.
For data exploring, I would like to create multiple graphs, one for each Markers, showing the concentration of that marker on the Y axis, with time on the X-axis with multiple ID’s in the same graph.
somewhat example would be: https://i2.wp.com/michaeltoth.me/figures/20190416_ggplot_geom_line/color_aes-1.png?zoom=1.5&w=578&ssl=1
I can’t figure out how to create te loop. Dataset example would be:
set.seed(8374636)
data <- data.frame(ID = 1:5,
time = 1:5,
markerA = rnorm(100)),
markerB = rnorm(100),
markerC = rnorm(100),
head(data)
ggplot(data, aes(x=time,
y=marker (A, B, C etc),
col=ID,
group=ID,
text = ID))+
geom_line(
size=1)+
geom_point(
size=2)
I won't bother you with the design part of the code, so I left that out.
Hope you can help me once again. Thanks in advance!
Hi Niels,
I’m sorry for the delayed response. I was on a long vacation, so unfortunately I wasn’t able to get back to you earlier. Do you still need help with your syntax?
Regards,
Joachim
Hi Joachim,
No problem. But yes, I would still like your help please!
Hi Niels,
I’m uncertain if I got it correctly which component of your data should be shown at which position. However, I have created an updated example, which hopefully helps to find a solution for your situation.
I think the crucial point of the following code is that I’m reshaping the data from wide to long format using the pivot_longer function.
Consider the example code and its resulting graph below:
You may read more on how to draw multiple time series lines in the same graph in this tutorial.
I hope this helps!
Joachim
Hi Joachim,
Unfortunatelly, this is not what I intended..
Maybe an easier way to explain my question.
I want to create the same plot you’ve made. But then with multiple subjects in one graph in stead of markers.
Now I’m looking for an easy way to create the same plot, but then for marker A, B, C etc.
It would save me a lot of time if I don’t have to copy paste a few lines of code and edit the Marker value.
Hope you wan’t to help me once more!
Hi Niels,
Thanks for the further explanations. I’m still uncertain if I got your question correctly. However, is the following code doing what you are looking for? It creates three separate images for the three different markers.
Regards,
Joachim
Hi Joachim,
Yes! I was able to adjust some things, but pivoting to a long format apparently did the trick !
Thanks again!
Hi Niels,
This is great to hear, glad you found a solution!
Regards,
Joachim