Generate samples using model function itself?

[pookie]

Thank you very much @martinjankowiak, this works!
I see it will ignore the data values in obs=data, and sample from the distribution.