Thank you for your help. I tried to use scale=TRUE but still not the correct answer. And I can't get the gamma on the result.