In line 152 of the RMSD_DIF() function both mean and standard deviation are assigned to variance estimate from the model, COV_11. This seems like a typo.
if(dentype %in% c('norm', 'snorm')){
mu <- sv$value[sv$name == "MEAN_1"]
sigma2 <- mu <- sv$value[sv$name == "COV_11"]
In line 152 of the
RMSD_DIF()
function both mean and standard deviation are assigned to variance estimate from the model, COV_11. This seems like a typo.