Code used in the analysis is here