Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Continuously add lines to ggplot with for loop
1
gravatar for frcamacho
3.2 years ago by
frcamacho • 180
United States
frcamacho • 180 wrote:

Hi there,

I have an object called gPlot and I want to continuously add geom_lines to the gPlot object. 

df is a dataframe 
sstart ssend
167    2637
552    8273

My following code is : 

gPlot<-ggplot() + geom_line(aes(x=1:6549, y=1, colour="red")) + geom_line(aes(x=295:1218, y=.99, colour="blue")) + geom_line(aes(x=3084:4340, y=.99, colour="green"))+coord_cartesian(xlim = c(0, 6549), ylim = c(0,1)) +  xlab("coordinates") + ylab("regions")

geomLine <- NULL
for(i in 1:nrow(df)){
  rowDF<-df[i,]
  dfstart <-rowDF$sstart
  dfend<-rowDF$send
  geomLine<-gPlot+geom_line(aes(x=dfstart:dfend, y=.95, colour="black"))

}

However, when I plot the new ranges from the df, the plot does not add all the ranges but, just puts the last number range. Any help would be appreciated.

Thanks!
ggplot2 R • 6.1k views
ADD COMMENT • link •
modified 3.2 years ago • written 3.2 years ago by frcamacho • 180

did you find the reason why this happen in the for loop?

ADD REPLY • link written 10 months ago by lffu_0032 • 40
4
gravatar for frcamacho
3.2 years ago by
frcamacho • 180
United States
frcamacho • 180 wrote:

I finally solved the above question after staring at the computer. Using aes_string() made it work. 

for(i in 1:nrow(stool_polyketide_analysis_df_vfam)){
  rowDF<-stool_polyketide_analysis_df_vfam[i,]
  dfstart <-rowDF$sstart
  dfend<-rowDF$send
  gPlot<-gPlot+geom_line(aes_string(x=dfstart:dfend, y=.95-.005*i))

}
ADD COMMENT • link modified 2.5 years ago • written 3.2 years ago by frcamacho • 180
1

Believe it or not, this is the best answer to solving the looping problem with ggplot. Most of the answers have a lot of elaborate schemes that I was not able to get to work. This simple one works!! Thanks

ADD REPLY • link written 7 months ago by godcent70 • 10

Cool, aes_string is exactly what I was looking for :)

ADD REPLY • link written 5 weeks ago by ATpoint ♦ 32k
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1283 users visited in the last hour