Implemented some functions to do the fitting, seems to work right. Also implemented log_inv_fit for Rmerge. 
Implemented some functions to do the fitting, seems to work right. Also implemented log_inv_fit for Rmerge. 
xia2/xia2,xia2/xia2