| Title: | Alternative Estimators to Adjusted R-Squared |
|---|---|
| Description: | Provides alternatives to the normal adjusted R-squared estimator for the estimation of the multiple squared correlation in regression models, as fitted by the lm() function. The alternative estimators are described in Karch (2020) <DOI:10.1525/collabra.343>. |
| Authors: | Julian Karch [aut, cre] |
| Maintainer: | Julian Karch <[email protected]> |
| License: | GPL-2 |
| Version: | 1.1.0 |
| Built: | 2024-09-11 05:52:35 UTC |
| Source: | https://github.com/karchjd/altr2 |
Returns different estimates of the multiple squared correlation.
altR2(lmOut)altR2(lmOut)
lmOut |
object of class "lm" as returned by the function |
A named vector with the different estimates
## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data. ctl <- c(4.17, 5.58, 5.18, 6.11, 4.50, 4.61, 5.17, 4.53, 5.33, 5.14) trt <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69) group <- gl(2, 10, 20, labels = c("Ctl", "Trt")) weight <- c(ctl, trt) lm.D9 <- lm(weight ~ group) estimates <- altR2(lm.D9)## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data. ctl <- c(4.17, 5.58, 5.18, 6.11, 4.50, 4.61, 5.17, 4.53, 5.33, 5.14) trt <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69) group <- gl(2, 10, 20, labels = c("Ctl", "Trt")) weight <- c(ctl, trt) lm.D9 <- lm(weight ~ group) estimates <- altR2(lm.D9)
Returns different estimates of the multiple squared correlation.
estimate_adj_R2(Rsquared, N, p)estimate_adj_R2(Rsquared, N, p)
Rsquared |
R-squared value |
N |
Number of observations |
p |
Number of predictors |
A named vector with the different estimates
## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data. ctl <- c(4.17, 5.58, 5.18, 6.11, 4.50, 4.61, 5.17, 4.53, 5.33, 5.14) trt <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69) group <- gl(2, 10, 20, labels = c("Ctl", "Trt")) weight <- c(ctl, trt) lm.D9 <- lm(weight ~ group) estimates <- estimate_adj_R2(summary(lm.D9)$r.squared, length(weight), 1)## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data. ctl <- c(4.17, 5.58, 5.18, 6.11, 4.50, 4.61, 5.17, 4.53, 5.33, 5.14) trt <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69) group <- gl(2, 10, 20, labels = c("Ctl", "Trt")) weight <- c(ctl, trt) lm.D9 <- lm(weight ~ group) estimates <- estimate_adj_R2(summary(lm.D9)$r.squared, length(weight), 1)