Bar graph and box-and-whisker plot sorted over binary variable

#1
Hi,

I have a dataset that looks like that below.


I want to create graphs like the example below:


The variable names var1-var4 in my dataset would replace the variable names on the x axis in the graph. My variable "status" is what is labeled Unstable/Stable in the graph. Pairid is a variable that describes which pairid each row belongs to, as this is a matched case-control study with two controls per case.

Have you got tips on how to create such a graph?
Thank you in advance.