source
# Analyze production measured in tonnes
prod <- production[production$measure == "Q_tlw" & production$value > 0,]
prod <- merge(prod, source)
# Production by source in 2023, in million tonnes
x <- aggregate(value~source_name, prod, sum, subset=year==2023)
transform(x, value=round(value/1e6))
Run the code above in your browser using DataLab