 great work ,Can you share the model fitting code? I want to test whether my data is consistent, thank you