R 의사-R²

1 개요[ | ]

R 의사-R²
df = data.frame(
  Hours = c(0.50,0.75,1.00,1.25,1.50,1.75,1.75,2.00,2.25,2.50,2.75,3.00,3.25,3.50,4.00,4.25,4.50,4.75,5.00,5.50),
  Pass = c(0,0,0,0,0,0,1,0,1,0,1,0,1,0,1,1,1,1,1,1)
)
model = glm(Pass ~ Hours, df, family=binomial)

library(DescTools)
PseudoR2(model, c("CoxSnell","Nagelkerke","McFadden"))
df = data.frame(
  Hours = c(0.50,0.75,1.00,1.25,1.50,1.75,1.75,2.00,2.25,2.50,2.75,3.00,3.25,3.50,4.00,4.25,4.50,4.75,5.00,5.50),
  Pass = c(0,0,0,0,0,0,1,0,1,0,1,0,1,0,1,1,1,1,1,1)
)
model = glm(Pass ~ Hours, df, family=binomial)

# R-square
nullmodel = glm(Pass ~ 1, df, family=binomial)
cat("R² =", 1-logLik(model)/logLik(nullmodel), "\n")

# -2 log-likelihood
cat("-2 log-likelihood =", -2*logLik(model), "\n")

2 같이 보기[ | ]

문서 댓글 ({{ doc_comments.length }})
{{ comment.name }} {{ comment.created | snstime }}