standardize vcv

[daijiang]

@arives should we standardize the vcv matrix here too? Without doing so, sometimes, the diag values are too large, and then rr2::R2_resid will have value of 1 because the scal <- prod(diag(s2 * phyV))^(1/n) will be Inf.

phyr/R/pglmm_compare.R

Line 292 in d0f91d9

re.1 <- list(covar = vcv(phy))

[arives]

This sounds like a good idea, and I can't see any reason not to.