# fit a simple linear regression with a binary predictor
mod = lm(charges ~ sex, data = RESI::insurance)
# calculate t-value
t = summary(mod)$coefficients[2, "t value"]
# calculate RESI (S)
S = t2S(t, n = 1338, rdf = 1336)
# convert S to R^2
S2Rsq(S)
Run the code above in your browser using DataLab