Compute Nagelkerke R2 ££ perhaps add: r2nag?
rnag(yobs = NULL, ypred = NULL)
yobs
is a vector of observed phenotypes
ypred
is a vector of predicted phenotypes