在for循环中用ggplot绘图不起作用 [英] plot with ggplot in for-loop doesn't work

查看:182
本文介绍了在for循环中用ggplot绘图不起作用的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!

问题描述

我做了一个ggplot,里面介绍了一个物种的种群数量。这个脚本运作良好。所以我做了for-loop连锁产生25个物种这个情节,但它不起作用。

下面是一个物种的例子:

pre $ A < data.frame(sp =A,x = c(rep(0.1,3),rep(2,3),rep(4.5:8.5,4),6.5),width = c(rep(NA,3) ,rep(4,3),rep(1,20),5),y = c(rep(1:3,2),rep(1:4,each = 5),5),col = c (white,21),rep(grey80,6)),name1 = c(CI(%),CN(%),rep(,25)),name1i = c (n,7),1.33,0.53,0.52,0.42,NA,27.12,10.90,10.55,8.56,rep(n,1) (NA,11)),NB2 = C(REP(NA,16),24754,13558,5449,5273,4281,代表(NA,6)),NB3 = C(REP(NA,6),2.43,代表(NA,4),49.51,rep(NA,15)),name2 = c(rep(,21)J,F,M,A,M ))

library(ggplot2)
library(gridExtra)
table < - A
windows()
ggplot(table,aes(x = x,y = y,width = width,fill = col))+
geom_tile(color =black)+
scale_fill_manual(values = c(white=white,grey80 =grey80),
breaks = levels(table $ col),na.value =white)+
geom_text(aes(label = name1),size = 6,hjust = 0) +
geom_t分机(AES(标记= name1i),大小= 6,hjust = 0,fontface = 3)+
geom_text(AES(标记= NB1),大小= 6)+
geom_text(AES(标记= (aes(label = nb3),size = 6,fontface =bold)+
geom_text(aes(label = name2),size = 6,fontface =bold)+
coord_cartesian(xlim = c(0,9),ylim = c(0.5,5.5))+
labs(x =,y = )+
scale_x_discrete(expand = c(0,0))+
scale_y_discrete(expand = c(0,0))+
theme(axis.ticks = element_blank(),
axis.text = element_blank(),
legend.position =none,
panel.background = element_rect(fill ='white',color ='white'),
plot .margin = unit(c(0.5,0.5,-1,-1),lines))



有for循环的类:

  A<  -  data.frame (sp =A,x = c(rep(0.1,3),rep(2,3),rep(4.5:8.5,4),6.5),width = c(rep(NA,3),rep 4),rep(1,20),5),y = c(rep(1:3,2),rep(1:4,each = 5),5),col = c(rep ,21),rep(grey80,6)),name1 = c(CI (%)  C.N。 (%),rep(,25)),name1i = c(rep(,2)Total,rep(,24)),nb1 = c(rep(NA,7), 1.33,0.53,0.52,0.42,NA,27.12,10.90,10.55,8.56,代表(NA,11)),NB2 = C(REP(NA,16),24754,13558,5449,5273,4281,代表(NA ,6)),nb3 = c(rep(NA,6),2.43,rep(NA,4),49.51,rep(NA,15)),name2 = c(rep(,21)J ,F,M,A,M,Month))
B < - data.frame(sp =B,x = c(rep(0.1,3) ,rep(2,3),rep(4.5:8.5,4),6.5),width = c(rep(NA,3),rep(4,3),rep(1,20),5) c(rep(1:3,2),rep(1:4,each = 5),5),col = c(rep(white,21),rep(grey80,6) c(CI(%),CN(%),rep(,25)),name1i = c(rep(,2)Total,rep(,24) NB1 = C(REP(NA,6),4.5,6.74,3.43,3.56,NA,19.89,29.82,15.18,15.75,代表(​​NA,12)),NB2 = C(REP(NA,16),3282, (NA,10),7.30,rep(NA,4),32.29,rep(NA,11)),name2 = c( (A,B)$ b $($,$)
data $ (b)b







$ table< - data [data $ sp == i,]
windows()
ggplot(table,aes(x = x,y = y,width = width,fill = col))+
geom_tile(color =black)+
scale_fill_manual(values = c(white=white,grey80 =grey80),
breaks = levels(table $ col),na.value =white)+
geom_text(aes(label = name1),size = 6,hjust = 0) +
geom_text(AES(标记= name1i),大小= 6,hjust = 0,fontface = 3)+
geom_text(AES(标记= NB1),大小= 6)+
geom_text (aes(label = nb2),size = 6,fontface = 3)+
geom_text(aes(label = nb3),size = 6,fontface =bold)+
geom_text (x = c(0,9),ylim = c(0.5,5.5))+
labs(x =),size = 6,fontface =bold)+
coord_cartesian ,y =)+
scale_x_discrete(expand = c(0,0))+
scale_y_discrete(expand = c(0,0))+
theme(axis.ticks = element_blank (),
axis.text = element_blank(),
legend.position =none,
panel.background = element_rect(fill ='white',color ='w (c(0.5,0.5,-1,-1),lines))
}
pre>

情节是空的,没有错误信息,所以我不明白我的错误在哪里。这只是一个ggplot问题?事实上,我不这么认为,因为其他的ggplot可以和for-loop一起工作...
你有什么想法吗?
感谢您的帮助!

解决方案

问题在于for循环。您需要在循环中使用打印。

  for(i in sp){
table< - data [ data $ sp ==A,]
windows()
ggp< - ggplot(...)+ ...
print(ggp)
}

考虑这个简单的例子:

 库(GGPLOT2)
DF = data.frame(X = 1:10,Y = RNORM(10))#样品数据
ggplot(DF)+ geom_point(AES(X ,Y))#渲染ggplot
为(I在1:2)ggplot(DF)+ geom_point(AES(X,Y))#没什么
为(I在1:2)照片(ggplot (df)+ geom_point(aes(x,y)))#渲染

另外,@ user229552说,你同时使用同一张表。


I made ggplot which presents population numbers of one species. This script works well. So I made for-loop to chain produce this plot for 25 species, but it doesn't works.

Here's an example for one species :

A <- data.frame(sp="A",x=c(rep(0.1,3),rep(2,3),rep(4.5:8.5,4),6.5), width=c(rep(NA,3), rep(4,3), rep(1,20),5),y=c(rep(1:3,2), rep(1:4, each=5),5),col=c(rep("white",21), rep("grey80",6)), name1=c("C.I. (%)","C.N. (%)", rep("",25)), name1i=c(rep("",2), "Total",rep("",24)), nb1=c(rep(NA,7),1.33,0.53,0.52,0.42,NA,27.12,10.90,10.55,8.56,rep(NA,11)),nb2=c(rep(NA,16),24754,13558,5449,5273,4281,rep(NA,6)),nb3=c(rep(NA,6),2.43,rep(NA,4),49.51, rep(NA,15)), name2=c(rep("",21),"J","F","M","A","M","Month"))

library(ggplot2)
library(gridExtra)
  table <- A
  windows()
  ggplot(table, aes(x=x, y=y, width=width, fill=col))+
    geom_tile(colour="black")+
    scale_fill_manual(values = c("white"="white", "grey80"="grey80"),    
                      breaks = levels(table$col), na.value = "white")+
    geom_text(aes(label=name1), size=6, hjust=0)+
    geom_text(aes(label=name1i), size=6, hjust=0, fontface=3)+
    geom_text(aes(label=nb1), size=6)+
    geom_text(aes(label=nb2), size=6, fontface=3)+
    geom_text(aes(label=nb3), size=6, fontface="bold")+
    geom_text(aes(label=name2), size=6, fontface="bold")+
    coord_cartesian(xlim=c(0,9), ylim=c(0.5,5.5))+
    labs(x = "",y = "") +                   
    scale_x_discrete(expand = c(0, 0)) +   
    scale_y_discrete(expand = c(0, 0))+     
    theme(axis.ticks = element_blank(),     
          axis.text = element_blank(),      
          legend.position = "none",        
          panel.background = element_rect(fill='white', colour='white'),   
          plot.margin = unit(c(0.5,0.5,-1,-1), "lines"))  

And for two species with for-loop :

A <- data.frame(sp="A",x=c(rep(0.1,3),rep(2,3),rep(4.5:8.5,4),6.5), width=c(rep(NA,3), rep(4,3), rep(1,20),5),y=c(rep(1:3,2), rep(1:4, each=5),5),col=c(rep("white",21), rep("grey80",6)), name1=c("C.I. (%)","C.N. (%)", rep("",25)), name1i=c(rep("",2), "Total",rep("",24)), nb1=c(rep(NA,7),1.33,0.53,0.52,0.42,NA,27.12,10.90,10.55,8.56,rep(NA,11)),nb2=c(rep(NA,16),24754,13558,5449,5273,4281,rep(NA,6)),nb3=c(rep(NA,6),2.43,rep(NA,4),49.51, rep(NA,15)), name2=c(rep("",21),"J","F","M","A","M","Month"))
B <- data.frame(sp="B",x=c(rep(0.1,3),rep(2,3),rep(4.5:8.5,4),6.5), width=c(rep(NA,3), rep(4,3), rep(1,20),5),y=c(rep(1:3,2), rep(1:4, each=5),5),col=c(rep("white",21), rep("grey80",6)), name1=c("C.I. (%)","C.N. (%)", rep("",25)), name1i=c(rep("",2), "Total",rep("",24)), nb1=c(rep(NA,6),4.5,6.74,3.43,3.56,NA,19.89,29.82,15.18,15.75,rep(NA,12)),nb2=c(rep(NA,16),3282,4920,2504,2599,5328,rep(NA,6)),nb3=c(rep(NA,10),7.30,rep(NA,4),32.29, rep(NA,11)), name2=c(rep("",21),"J","F","M","A","M","Month"))
data <- rbind(A,B)
sp <- c("A","B")

library(ggplot2)
library(gridExtra)
for(i in sp){
  table <- data[data$sp==i,]
  windows()
  ggplot(table, aes(x=x, y=y, width=width, fill=col))+
    geom_tile(colour="black")+
    scale_fill_manual(values = c("white"="white", "grey80"="grey80"),    
                      breaks = levels(table$col), na.value = "white")+
    geom_text(aes(label=name1), size=6, hjust=0)+
    geom_text(aes(label=name1i), size=6, hjust=0, fontface=3)+
    geom_text(aes(label=nb1), size=6)+
    geom_text(aes(label=nb2), size=6, fontface=3)+
    geom_text(aes(label=nb3), size=6, fontface="bold")+
    geom_text(aes(label=name2), size=6, fontface="bold")+
    coord_cartesian(xlim=c(0,9), ylim=c(0.5,5.5))+
    labs(x = "",y = "") +                   
    scale_x_discrete(expand = c(0, 0)) +   
    scale_y_discrete(expand = c(0, 0))+    
    theme(axis.ticks = element_blank(),    
          axis.text = element_blank(),     
          legend.position = "none",         
          panel.background = element_rect(fill='white', colour='white'),   
          plot.margin = unit(c(0.5,0.5,-1,-1), "lines"))  
}

Plot is empty without error message, so I do not understand where is my error. Is it just a ggplot problem? In fact I don't think so, because other ggplot works well with for-loop... Do you have any idea? Thanks for your help!

解决方案

The problem is the for loop. You need to use print in a loop.

for (i in sp) {
  table <- data[data$sp=="A",]
  windows()
  ggp <- ggplot(...) + ...
  print(ggp)
}

Consider this simple example:

library(ggplot2)
df=data.frame(x=1:10,y=rnorm(10))                       # sample data
ggplot(df)+geom_point(aes(x,y))                         # render ggplot
for (i in 1:2) ggplot(df)+geom_point(aes(x,y))          # nothing
for (i in 1:2) print(ggplot(df) + geom_point(aes(x,y))) # renders

Also, as @user229552 says, you are using the same table both times.

这篇关于在for循环中用ggplot绘图不起作用的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!

查看全文
登录 关闭
扫码关注1秒登录
发送“验证码”获取 | 15天全站免登陆