Performs matrix skew-t parameter estimation, Gallaugher and McNicholas (2017) <doi: 10.1002/sta4.143>.