R：McFadden＆Nagelkerke之間的關係

Question

我想知道是否有一種方法可以根據產生的輸出計算Nagelkerke R-square。我知道我可以直接計算McFadden R-square。但是Nagelkerke產生了我們認爲是模型更精確的力量。 如果你有這樣的想法，我沒有把運氣添加到我的設置上。

謝謝。

Answer 1

這個問題還沒有定義，所以我會盡我所能，假設「產生的輸出」是一個glm對象。此函數應用於glm對象時應產生適當的僞R方形。

Nagelkerke <- function(mod) { 
    l_full <- exp(logLik(mod)) 
    l_intercept <- exp(logLik(update(mod, . ~ 1))) 
    N <- length(mod$y) 
    r_2 <- (1 - (l_intercept/l_full)^(2/N))/(1 - l_intercept^(2/N)) 
    return(as.numeric(r_2)) 
} 

實施例：

model <- glm(formula = vs ~ mpg + disp, family = binomial("logit"), data = mtcars); 
Nagelkerke(model); 
#[1] 0.6574295