Ivan Jacob Agaloos Pesigan 2024-05-05
Description
Generates Monte Carlo confidence intervals 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()
. betaMC
combines ideas from Monte Carlo confidence intervals for the indirect effect (Pesigan and Cheung, 2023: http://doi.org/10.3758/s13428-023-02114-4) and the sampling covariance matrix of regression coefficients (Dudgeon, 2017: http://doi.org/10.1007/s11336-017-9563-z) to generate confidence intervals effect sizes in regression.
Installation
You can install the CRAN release of betaMC
with:
install.packages("betaMC")
You can install the development version of betaMC
from GitHub with:
if (!require("remotes")) install.packages("remotes")
remotes::install_github("jeksterslab/betaMC")
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 BetaMC()
function from the betaMC
package.
df <- betaMC::nas1982
Regression
Fit the regression model using the lm()
function.
object <- lm(QUALITY ~ NARTIC + PCTGRT + PCTSUPP, data = df)
Monte Carlo Sampling Distribution of Parameters
Normal-Theory Approach
mvn <- MC(object, type = "mvn")
Asymptotic distribution-free Approach
adf <- MC(object, type = "adf")
Heteroskedasticity Consistent Approach (HC3)
hc3 <- MC(object, type = "hc3")
Standardized Regression Slopes
Normal-Theory Approach
BetaMC(mvn, alpha = 0.05)
#> Call:
#> BetaMC(object = mvn, alpha = 0.05)
#>
#> Standardized regression slopes
#> type = "mvn"
#> est se R 2.5% 97.5%
#> NARTIC 0.4951 0.0750 20000 0.3414 0.6333
#> PCTGRT 0.3915 0.0767 20000 0.2370 0.5388
#> PCTSUPP 0.2632 0.0748 20000 0.1185 0.4093
Asymptotic distribution-free Approach
BetaMC(adf, alpha = 0.05)
#> Call:
#> BetaMC(object = adf, alpha = 0.05)
#>
#> Standardized regression slopes
#> type = "adf"
#> est se R 2.5% 97.5%
#> NARTIC 0.4951 0.0675 20000 0.3522 0.6148
#> PCTGRT 0.3915 0.0708 20000 0.2431 0.5201
#> PCTSUPP 0.2632 0.0769 20000 0.1072 0.4096
Heteroskedasticity Consistent Approach (HC3)
BetaMC(hc3, alpha = 0.05)
#> Call:
#> BetaMC(object = hc3, alpha = 0.05)
#>
#> Standardized regression slopes
#> type = "hc3"
#> est se R 2.5% 97.5%
#> NARTIC 0.4951 0.0795 20000 0.3230 0.6334
#> PCTGRT 0.3915 0.0820 20000 0.2187 0.5379
#> PCTSUPP 0.2632 0.0849 20000 0.0936 0.4254
Other Effect Sizes
The betaMC
package also has functions to generate Monte Carlo confidence intervals for other effect sizes such as RSqMC()
for multiple correlation coefficients (R-squared and adjusted R-squared), DeltaRSqMC()
for improvement in R-squared, SCorMC()
for semipartial correlation coefficients, PCorMC()
for squared partial correlation coefficients, and DiffBetaMC()
for differences of standardized regression coefficients.
Multiple Correlation Coefficients (R-squared and adjusted R-squared)
RSqMC(hc3, alpha = 0.05)
#> Call:
#> RSqMC(object = hc3, alpha = 0.05)
#>
#> R-squared and adjusted R-squared
#> type = "hc3"
#> est se R 2.5% 97.5%
#> rsq 0.8045 0.0614 20000 0.6487 0.8866
#> adj 0.7906 0.0658 20000 0.6236 0.8785
Improvement in R-squared
DeltaRSqMC(hc3, alpha = 0.05)
#> Call:
#> DeltaRSqMC(object = hc3, alpha = 0.05)
#>
#> Improvement in R-squared
#> type = "hc3"
#> est se R 2.5% 97.5%
#> NARTIC 0.1859 0.0688 20000 0.0496 0.3197
#> PCTGRT 0.1177 0.0543 20000 0.0251 0.2346
#> PCTSUPP 0.0569 0.0374 20000 0.0068 0.1484
Semipartial Correlation Coefficients
SCorMC(hc3, alpha = 0.05)
#> Call:
#> SCorMC(object = hc3, alpha = 0.05)
#>
#> Semipartial correlations
#> type = "hc3"
#> est se R 2.5% 97.5%
#> NARTIC 0.4312 0.0867 20000 0.2227 0.5654
#> PCTGRT 0.3430 0.0829 20000 0.1583 0.4844
#> PCTSUPP 0.2385 0.0775 20000 0.0823 0.3852
Squared Partial Correlation Coefficients
PCorMC(hc3, alpha = 0.05)
#> Call:
#> PCorMC(object = hc3, alpha = 0.05)
#>
#> Squared partial correlations
#> type = "hc3"
#> est se R 2.5% 97.5%
#> NARTIC 0.4874 0.1193 20000 0.1753 0.6497
#> PCTGRT 0.3757 0.1154 20000 0.1074 0.5528
#> PCTSUPP 0.2254 0.1122 20000 0.0288 0.4550
Differences of Standardized Regression Coefficients
DiffBetaMC(hc3, alpha = 0.05)
#> Call:
#> DiffBetaMC(object = hc3, alpha = 0.05)
#>
#> Differences of standardized regression slopes
#> type = "hc3"
#> est se R 2.5% 97.5%
#> NARTIC-PCTGRT 0.1037 0.1416 20000 -0.1752 0.3763
#> NARTIC-PCTSUPP 0.2319 0.1326 20000 -0.0420 0.4789
#> PCTGRT-PCTSUPP 0.1282 0.1366 20000 -0.1532 0.3813
Documentation
See GitHub Pages for package documentation.