calculate.var.np function

calculates closed form variance estimate for R