 How I can to obtain the variancia of the method sandwich to apply it in a model of Poisson with overdispersion in R ?

I need the algorithm 

	[[alternate HTML version deleted]]

