在ggplot中绘制直方图

Plot histogram in ggplot

我有这个数据框。

dput(EF_Lat_Am)
structure(list(V1 = structure(c(4L, 3L, 5L, 6L, 1L, 2L, 7L, 8L, 
4L, 3L, 5L, 6L, 1L, 2L, 7L, 8L), .Label = c("Crop Agriculture", 
"Mining", "Mixed Agriculture", "Other land use", "Pasture", "Tree crops", 
"Urban", "Water"), class = "factor"), V2 = structure(c(1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = "Emission Factor", class = "factor"), 
    V3 = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L), .Label = c("2000", "2005"), class = "factor"), 
    V4 = c(77.0109486116396, 69.2454348145657, 73.684103657833, 
    71.0430911289891, 43.136201172115, 117.358146800995, 77.4653952935238, 
    89.0966064142874, 71.8286578413912, 67.9099357961953, 76.7438444998728, 
    67.4818461466729, 50.6468079101972, 117.799797611894, 78.7347377710757, 
    81.3020943196897)), .Names = c("V1", "V2", "V3", "V4"), row.names = c(NA, 
16L), class = "data.frame")

如您所见,对于 2000 年和 2005 年,我有每种土地利用类型的排放因子值。我想绘制一个直方图,x 轴为土地利用类型,y 轴为排放因子。此外,对于每种土地用途,我希望这两年的条形图相邻。我还想要一个图例,显示哪些年份对应于条形图(2000 年或 2005 年)。谢谢你的帮助。

答案在这里。

ggplot(EF_Lat_Am, aes(x=V1, y = V4, fill=V3, width=.85)) + geom_bar(position="dodge", stat="identity") +
  labs(x = "", y = "EF (T/ha)") +
  theme(axis.text=element_text(size=16),axis.title=element_text(size=20), 
        legend.title=element_text(size=20, face='bold'),legend.text=element_text(size=20), axis.line = element_line(colour = "black")) +
  scale_fill_grey("Period") + scale_y_continuous(limits=c(0,120)) + theme_classic(base_size = 20, base_family = "")