DiagGelman.R finished