Changes in version 2.0 Major Changes There are a few major changes in output and input that will likely break existing dependencies. Here are the highlights - the main output (lfsr, lfdr, etc) is rearranged into a dataframe, called result. - so, for example, the lfsr is now a$result$lfsr instead of a$lfsr - or, better, use the accessor function get_lfsr(a) to extract the lfsr etc - I added accessor functions get_lfsr, get_lfdr, get_pm, get_psd etc to access the lfsr, lfdr, posterior mean and posterior standard deviation. Using these functions to access results will help ensure your code remains valid if I happen to change the internal structure of the results again (although not anticipated...) - output fitted.g is renamed fitted_g, and flash.data becomes flash_data to make the whole package convention more consistent. Also fit becomes fit_details. - function prefixes comppost and compdens replaced with comp_post and comp_dens, again for consistency. - nonzeromode option is replaced with the option mode to specify mode. Or use mode="estimate" to estimate the mode. - more flexible control of output. For example, you can say you want only the logLR output by specifying outputlevel = c("lfsr","logLR"), or only posterior mean by outputlevel = c("PosteriorMean").