Many thanks for your reply! It definitely has to do something with specific model specification and indeed seems to be breaking some gradients. By the way, for this specific application my ‘bootstrapping approach’ shows we don’t need to worry that much about estimation error from the first stage. But sure would be interesting from a more technical perspective to look into some time!
1 Like