我正在尝试使用循环生成图。
l1<-factor(rep(letters,4))
n1<-abs(rnorm(104))*10000
b1<-rep(c("1","2","3","4","5","6","7","8"),c(2,2,11,24,11,20,33,1))
k1<-rep((rep(c("A","B","C","D"),c(2,3,4,4))),8)
my.df<-data.frame(l1,b1,k1,n1) #make a dataframe
names(my.df)<-c("letter","branch","ltrtype","number") #factor names
library(ggplot2)
branch.list<-unique(my.df$branch)
sayi<-length(branch.list) # list of factor:letters
for (i in 1:sayi) {
branch.iter<-branch.list[i]
my.df.plot<-subset(my.df,my.df$branch==branch.iter,drop=T)
my.df.plot$branch<-factor(my.df.plot$branch) #So that unused levels don't show up
my.df.plot$letter<-factor(my.df.plot$letter) #So that unused levels don't show up
my.df.plot$ltrtype<-factor(my.df.plot$ltrtype) #So that unused levels don't show up
my.df.plot$number<-as.numeric(as.character(my.df.plot$number))
my.df.plot<-data.frame(my.df.plot)
myfilename<-paste(branch.iter,".jpeg",sep="")
jpeg(file=myfilename)
cizim<-ggplot(my.df.plot,aes(letter,number,fill=ltrtype))
cizim<-cizim + geom_bar(width = 1, position = "fill", binwidth = 1) + facet_grid(ltrtype~.)
cizim<-cizim + opts(title=branch.iter)
print(cizim)
dev.off()
}
(Q1):当x轴的级别数改变条的宽度时 如何防止这种情况并使每个图中的条宽相同?
alt text http://img411.imageshack.us/i/95325388.jpg/
alt text http://img411.imageshack.us/i/91510133.jpg/
(Q2):当i=7
R发出以下警告时:
(数据$ ymin == 0))警告(“当ymin!= 0时填充不明确”): 缺少需要TRUE / FALSE的值
我该怎么办呢?
(Q3):在这种情况下是否有更简单的方法来删除未使用的级别,因此我不必使用
my.df.plot$branch<-factor(my.df.plot$branch)
每次?
答案 0 :(得分:2)
你正在制作一些非常奇怪的情节。通过使用position="fill"
,您将每个条形伸展到高度为1(因为对应于该字母的一个观察值是对应于面板内字母的所有观察值的100%),完全失去了您试图绘制的任何信息。我的猜测是你的一些问题源于这个错误,但我不确定。
(Q1)您是否希望不同分支的条形图中的条形宽度相同?由于您要更改x变量的级别数,因此条形图必须更宽以填充图形。一些解决方案选项:
expand
的{{1}}选项让您的酒吧更窄并以图表为中心。因此,如果您有N个总x值(此处N = 26个字母),但某个特定情节仅使用其中的k个,则将scale_x_discrete
添加到您的绘图中。第一项是乘法扩展因子,这是默认值,第二项是加法因子。 (Q2)i = 7是唯一具有多个+ scale_x_discrete(expand=c(0.05, (N-k)/2))
值的组,对应于相同的字母/ ltrgroup组合。酒吧geom不知道该怎么做。我同意错误信息真的很神秘。
(Q3)一种选择是避免使用因子 - 在组合字符向量时使用number
,然后子集不会保留未使用的级别。
答案 1 :(得分:1)
(Q1)我认为不可能修复条宽。 Aniko建议保持所有级别对我来说最有意义。
(Q2)将binwidth = 1
替换为stat="identity"
,因为我认为您不需要stat="bin"
。
(Q3)其他选项包括drop.levels
中的gdata
- 包,以及dropUnusedLevels
中的Hmisc
- 包。