Implement the bias reduced pim equations for the vcov matrix (sandwich) estimator

[JoFAM]

The option bias.reduced = TRUE should give correct vcov estimates. For now the sandwich estimator gives the unbiased estimates.

Either the default should be the sandwich estimator, or the one based on the experimental likelihood should be integrated in the pim function.