INAR with negative binomial thinning

[MUHAMMED_AHAMMED_E_S]

How can I simulate an INAR model with the negative binomial thinning operator?