For 1) I don't know if we have better data respect Taiwan DSI i have to try all of them, for 2) I think that the negative values are numerical error due to the lack of a positive constrain in the pdf when the model fit the coefficients. We can clip the values of rtop and msd to solve the problem