library(ggplot2)

a <- data.frame(
landward_part=c('4001740','4001740','4001740','7001503','7001503','7001503','7003775','7003775','7003775'),
seaward_part=c('4001740','7001503','7003775','4001740','7001503','7003775','4001740','7001503','7003775'),
erosion_volume=c(73.97,74.36,80.91,210.70,215.41,226.88,103.63,106.35,115.02))

setwd('E:/Sources/OpenEarthModels/tudelft/dune_erosion_sensitivity_bathy/R/')

win.metafile("desb_seaward_landward_erosion.emf",
width=5, height=4)
ggplot(data=a,aes(landward_part,seaward_part, fill=erosion_volume))  + geom_tile() +
xlab('Landward part') +
ylab('Seaward part')
dev.off()

win.metafile("desb_seaward_erosion.emf",
width=5, height=4)
ggplot(a, aes(seaward_part,erosion_volume)) +
geom_jitter(data=a,aes(seaward_part,color=landward_part,shape=landward_part),position=position_jitter(width=.1)) +
xlab('Seaward part') +
ylab(expression(paste('Erosion volume [',m^3,'/m]')))
dev.off()

win.metafile("desb_landward_erosion.emf",
width=5, height=4)
ggplot(a, aes(landward_part,erosion_volume)) +
geom_jitter(data=a,aes(landward_part,color=seaward_part,shape=seaward_part),position=position_jitter(width=.1)) +
xlab('Landward part') +
ylab(expression(paste('Erosion volume [',m^3,'/m]')))
dev.off()