Thanks Frank for your response. I have used the rms package. Reproducible example here below :
install.packages("rms") library(rms) a=c("tres grand", "grand", "petit","petit","tres grand","grand","petit","petit","tres grand","grand") b=c("homme", "homme", "femme", "femme", "femme", "homme", "homme", "homme", "femme", "femme") h <- orm(a ~ b) h The Nagelkerke R2 = 0.074. Is it right ? Best Regards, SV ----- Mail original ----- De : Frank Harrell <f.harr...@vanderbilt.edu> À : r-help@r-project.org Cc : Envoyé le : Samedi 18 juillet 2015 19h51 Objet : Re: [R] Nagelkerke Pseudo R-squared It's implemented in the R rms package's lrm and orm functions. ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Nagelkerke-Pseudo-R-squared-tp4710014p4710031.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.