Ivan Jacob Agaloos Pesigan 2024-10-01
Description
Generates nonparametric bootstrap confidence intervals (Efron & Tibshirani, 1993: https://doi.org/10.1201/9780429246593) for standardized regression coefficients (beta) and other effect sizes, including multiple correlation, semipartial correlations, improvement in R-squared, squared partial correlations, and differences in standardized regression coefficients, for models fitted by lm().
Installation
You can install the CRAN release of betaNB with:
install.packages("betaNB")You can install the development version of betaNB from GitHub with:
if (!require("remotes")) install.packages("remotes")
remotes::install_github("jeksterslab/betaNB")Example
In this example, a multiple regression model is fitted using program quality ratings (QUALITY) as the regressand/outcome variable and number of published articles attributed to the program faculty members (NARTIC), percent of faculty members holding research grants (PCTGRT), and percentage of program graduates who received support (PCTSUPP) as regressor/predictor variables using a data set from 1982 ratings of 46 doctoral programs in psychology in the USA (National Research Council, 1982). Confidence intervals for the standardized regression coefficients are generated using the BetaNB() function from the betaNB package.
df <- betaNB::nas1982Regression
Fit the regression model using the lm() function.
object <- lm(QUALITY ~ NARTIC + PCTGRT + PCTSUPP, data = df)Nonparametric Bootstrap
nb <- NB(object)Standardized Regression Slopes
BetaNB(nb, alpha = 0.05)
#> Call:
#> BetaNB(object = nb, alpha = 0.05)
#>
#> Standardized regression slopes
#> type = "pc"
#> est se R 2.5% 97.5%
#> NARTIC 0.4951 0.0724 5000 0.3544 0.6385
#> PCTGRT 0.3915 0.0765 5000 0.2311 0.5314
#> PCTSUPP 0.2632 0.0793 5000 0.1071 0.4171Other Effect Sizes
The betaNB package also has functions to generate nonparametric bootstrap confidence intervals for other effect sizes such as RSqNB() for multiple correlation coefficients (R-squared and adjusted R-squared), DeltaRSqNB() for improvement in R-squared, SCorNB() for semipartial correlation coefficients, PCorNB() for squared partial correlation coefficients, and DiffBetaNB() for differences of standardized regression coefficients.
Multiple Correlation Coefficients (R-squared and adjusted R-squared)
RSqNB(nb, alpha = 0.05)
#> Call:
#> RSqNB(object = nb, alpha = 0.05)
#>
#> R-squared and adjusted R-squared
#> type = "pc"
#> est se R 2.5% 97.5%
#> rsq 0.8045 0.0527 5000 0.6936 0.8978
#> adj 0.7906 0.0564 5000 0.6717 0.8905Improvement in R-squared
DeltaRSqNB(nb, alpha = 0.05)
#> Call:
#> DeltaRSqNB(object = nb, alpha = 0.05)
#>
#> Improvement in R-squared
#> type = "pc"
#> est se R 2.5% 97.5%
#> NARTIC 0.1859 0.0597 5000 0.0829 0.3165
#> PCTGRT 0.1177 0.0483 5000 0.0334 0.2203
#> PCTSUPP 0.0569 0.0341 5000 0.0087 0.1392Semipartial Correlation Coefficients
SCorNB(nb, alpha = 0.05)
#> Call:
#> SCorNB(object = nb, alpha = 0.05)
#>
#> Semipartial correlations
#> type = "pc"
#> est se R 2.5% 97.5%
#> NARTIC 0.4312 0.0697 5000 0.2880 0.5626
#> PCTGRT 0.3430 0.0727 5000 0.1828 0.4693
#> PCTSUPP 0.2385 0.0713 5000 0.0932 0.3731Squared Partial Correlation Coefficients
PCorNB(nb, alpha = 0.05)
#> Call:
#> PCorNB(object = nb, alpha = 0.05)
#>
#> Squared partial correlations
#> type = "pc"
#> est se R 2.5% 97.5%
#> NARTIC 0.4874 0.0984 5000 0.2860 0.6721
#> PCTGRT 0.3757 0.1090 5000 0.1561 0.5786
#> PCTSUPP 0.2254 0.1143 5000 0.0440 0.4817Differences of Standardized Regression Coefficients
DiffBetaNB(nb, alpha = 0.05)
#> Call:
#> DiffBetaNB(object = nb, alpha = 0.05)
#>
#> Differences of standardized regression slopes
#> type = "pc"
#> est se R 2.5% 97.5%
#> NARTIC-PCTGRT 0.1037 0.1318 5000 -0.1443 0.3721
#> NARTIC-PCTSUPP 0.2319 0.1246 5000 -0.0054 0.4825
#> PCTGRT-PCTSUPP 0.1282 0.1263 5000 -0.1196 0.3711Documentation
See GitHub Pages for package documentation.