zvar <- function(yi, vi) {
  vi/(1 - yi^2)^2
}