使用gganimate
和ggplot
制作动画饼图时遇到问题。
我想每年都有正常的馅饼,但我的产量完全不同。
您可以使用mtcars
:
library(ggplot2)
library(gganimate)
#Some Data
df<-aggregate(mtcars$mpg, list(mtcars$cyl,mtcars$carb), sum)
colnames(df)<-c("X","Y","Z")
bp<- ggplot(df, aes(x="", y=Z, fill=X, frame=Y))+
geom_bar(width = 1, stat = "identity") + coord_polar("y", start=0)
gganimate(pie, "output.gif")
这是输出:
frame
只有一个级别时效果很好:
答案 0 :(得分:3)
ggplot代码创建一个单独的堆积条形图,其中包含df
中每一行的一个部分。使用coord_polar
,这将成为单个饼图,其中数据框中的每一行都有一个楔形。然后,当您使用gg_animate
时,每个框架仅包含与给定级别Y
对应的楔形。这就是为什么每次只获得完整饼图的一部分的原因。
如果你需要为Y
的每个级别设置一个完整的饼图,那么一个选项是为Y
的每个级别创建一个单独的饼图,然后将这些饼图组合成一个GIF。这里有一些假数据的例子(我希望)与您的真实数据相似:
library(animation)
# Fake data
set.seed(40)
df = data.frame(Year = rep(2010:2015, 3),
disease = rep(c("Cardiovascular","Neoplasms","Others"), each=6),
count=c(sapply(c(1,1.5,2), function(i) cumsum(c(1000*i, sample((-200*i):(200*i),5))))))
saveGIF({
for (i in unique(df$Year)) {
p = ggplot(df[df$Year==i,], aes(x="", y=count, fill=disease, frame=Year))+
geom_bar(width = 1, stat = "identity") +
facet_grid(~Year) +
coord_polar("y", start=0)
print(p)
}
}, movie.name="test1.gif")
以上GIF中的馅饼大小都相同。但您也可以根据count
的每个级别Year
的总和来更改馅饼的大小(代码改编自this SO answer):
library(dplyr)
df = df %>% group_by(Year) %>%
mutate(cp1 = c(0, head(cumsum(count), -1)),
cp2 = cumsum(count))
saveGIF({
for (i in unique(df$Year)) {
p = ggplot(df %>% filter(Year==i), aes(fill=disease)) +
geom_rect(aes(xmin=0, xmax=max(cp2), ymin=cp1, ymax=cp2)) +
facet_grid(~Year) +
coord_polar("y", start=0) +
scale_x_continuous(limits=c(0,max(df$cp2)))
print(p)
}
}, movie.name="test2.gif")
如果我可以暂时编辑,虽然动画很酷(但饼图很酷,所以也许动画一堆饼图只会增加侮辱伤害),数据可能会更容易理解一个普通的静态线图。例如:
ggplot(df, aes(x=Year, y=count, colour=disease)) +
geom_line() + geom_point() +
scale_y_continuous(limits=c(0, max(df$count)))
或许这个:
ggplot(df, aes(x=Year, y=count, colour=disease)) +
geom_line() + geom_point(show.legend=FALSE) +
geom_line(data=df %>% group_by(Year) %>% mutate(count=sum(count)),
aes(x=Year, y=count, colour="All"), lwd=1) +
scale_y_continuous(limits=c(0, df %>% group_by(Year) %>%
summarise(count=sum(count)) %>% max(.$count))) +
scale_colour_manual(values=c("black", hcl(seq(15,275,length=4)[1:3],100,65)))