Another way to visualize the distribution, is using histograms.
par(mfrow=c(1,2)) hist(BEPAS$MVPA_acc [BEPAS$Walkability=="Low walkability"]) hist(BEPAS$MVPA_acc [BEPAS$Walkability=="High walkability"])