beta.div function

Beta diversity computed as Var(Y)