Untitled

df<-read.csv("https://raw.githubusercontent.com/tuyenhavan/Statistics/Dataset/Area_Comparison.csv")

head(df)

##         Cover                   Product Percent
## 1 Agriculture        Landsat-based area   35.24
## 2       Urban        Landsat-based area    4.56
## 3      Mining        Landsat-based area    5.60
## 4      Forest        Landsat-based area   52.50
## 5       Water        Landsat-based area    1.69
## 6 Agriculture Local land use statistics   31.90

Produce graphics showing the difference and similarity between local land use reported and Landsat 8 derived area

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 3.3.3

ggplot(data=df,aes(x=Cover,y=Percent,fill=factor(Product,labels=c("Landsat Derived Area 2016","Local Statistics")))) + 
  geom_bar(stat="identity",position = "dodge") + scale_fill_discrete(name="Legend") + 
  xlab("Land Covers") + ylab("Area (percent)") + theme_bw()