Variance estimation in complex survey sampling for generalized linear models