# How can I calculate Nagelkerke pseudo r2 with glmmkin object from GMMAT package in R?

I have used glmmkin function from GMMAT package to fit a logistic mixed model with the binary phenotype 'disease', one fixed variable 'PRS' and one kinship matrix 'GRM' to model the covariance structure of the random effects exactly like in this example:

    glmmkin(disease ~ PRS, data = pheno, kins = GRM, id = "id", family =