An implementation of the method published by Near et al MRM 73:44-50 (2015).
tdsr(mrs_data, ref = NULL, xlim = c(4, 0.5), max_t = 0.2)MRS data to be corrected.
optional MRS data to use as a reference, the mean of all dynamics is used if this argument is not supplied.
optional frequency range to perform optimisation, set to NULL to use the full range.
truncate the FID when longer than max_t to reduce time taken.
a list containing the corrected data; phase and shift values in units of degrees and Hz respectively.