Re: [R] Nagelkerke Pseudo R-squared
Already answered many thanks Frank. Best Regards SV Envoyé de mon iPhone Le 18 juil. 2015 à 20:39, varin sacha varinsa...@yahoo.fr a écrit : 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. __ 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.
Re: [R] Nagelkerke Pseudo R-squared
Thanks David for your response. Best regards, SV - Mail original - De : David Winsemius dwinsem...@comcast.net À : varin sacha varinsa...@yahoo.fr Cc : R-help Mailing List r-help@r-project.org Envoyé le : Samedi 18 juillet 2015 3h33 Objet : Re: [R] Nagelkerke Pseudo R-squared On Jul 17, 2015, at 4:33 PM, varin sacha wrote: Dear R-Experts, I have fitted an ordinal logistic regression with just 1 explanatory variable for the reproducible example here below. Everything is working, now I try to calculate the Nagelkerke Pseudo R-squared. I have found a package BaylorEdPsych providing many Pseudo R-squared, but the example shown in the package is for GLM (binary logistic regression) not for ordinal logistic regression. How can I calculate the Nagelkerke Pseudo R-squared for ordinal logistic regression ? polr-objects have a deviance node. If this has statistical value (which I have some doubts regarding) then just apply the usual formula to compare nested models. -- David. Many thanks as usual for your precious help. Reproducible example : install.packages(MASS) library(MASS) a=factor(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) m - polr(a ~ b, Hess=TRUE) summary(m) __ 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. David Winsemius Alameda, CA, USA __ 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.
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.
Re: [R] Nagelkerke Pseudo R-squared
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.
Re: [R] Nagelkerke Pseudo R-squared
On Jul 17, 2015, at 4:33 PM, varin sacha wrote: Dear R-Experts, I have fitted an ordinal logistic regression with just 1 explanatory variable for the reproducible example here below. Everything is working, now I try to calculate the Nagelkerke Pseudo R-squared. I have found a package BaylorEdPsych providing many Pseudo R-squared, but the example shown in the package is for GLM (binary logistic regression) not for ordinal logistic regression. How can I calculate the Nagelkerke Pseudo R-squared for ordinal logistic regression ? polr-objects have a deviance node. If this has statistical value (which I have some doubts regarding) then just apply the usual formula to compare nested models. -- David. Many thanks as usual for your precious help. Reproducible example : install.packages(MASS) library(MASS) a=factor(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) m - polr(a ~ b, Hess=TRUE) summary(m) __ 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. David Winsemius Alameda, CA, USA __ 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.