Discussion about twinsb_fit