chisq2S | R Documentation |
This function computes the robust effect size index from Vandekar, Tao, & Blume (2020). Vector arguments are accepted. If different length arguments are passed they are dealt with in the usual way of R. For mixed effects models, RESI is conditional on the average correlation structure within subjects.
chisq2S(chisq, df, n)
chisq |
The chi-square statistic for the parameter of interest. |
df |
Number of degrees of freedom of the chi-square statistic. |
n |
Number of independent samples. |
The formula for converting a Chi-square statistic to RESI is:
S = \sqrt(max( 0, (chisq - df)/n))
Returns a scalar or vector argument of the robust effect size index estimate.
# obtain Chi-sq value by fitting an lm and running a Wald test
mod = lm(charges ~ region * age + bmi + sex, data = RESI::insurance)
# run a Wald test with robust variance
wt = lmtest::waldtest(mod, vcov = sandwich::vcovHC, test = "Chisq")
# get Chi-sq value and degrees of freedom
chisq = wt$Chisq[2]
df = abs(wt$Df[2])
# run chisq2S to convert to RESI
chisq2S(chisq, df = df, n = nrow(mod$model))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.