data R squared; *** Below, specify the sample size "n", the number of predictor variables "p", and the squared multiple correlation coefficient "Rsq" ***; n = 60; p = 5; Rsq = .31; ***********************************************************************; ERsq = ((p / (n-1))); a = ((ERsq / (1-ERsq))); *******************Darlington's F below***********************************; Dar_F = ((Rsq / (1-Rsq)) * ((n-p-1) / p)) / (1 + (a * ((n-1) / p))); ***********************************************************************; v = (((ERsq * (n-1)) + p))**2 / (((ERsq * (n-1) * (ERsq + 2)) + p)); Dar_p = ((1 - probf(Dar_F, v, n-p-1))); *************Huberty's adjusted Rsq value below****************************; HubadRsq = ((Rsq - ERsq) / (1 - ERsq)); *************Huberty's effect size index below*****************************; Hubefsiz = (HubadRsq - ERsq); ***********************************************************************; proc print noobs; var Dar_F Dar_p ERsq HubadRsq Hubefsiz; options nodate; options nocenter; title1 'NOTE: Dar_F=F test reported in Darlington (1990) for testing the'; title2 'null hypothesis that R-squared equals the expected value of R-squared.'; title3 'Dar_p=Probability value for Dar_F.'; title4 'ERsq=The expected value, or long-run mean, of R-squared under the null'; title5 'hypothesis that rho=0.'; title6 'HubadRsq=Adjusted R-squared value presented in Huberty (1994). Expressed'; title7 'as a proportional-reduction-in-error or improvement-over-chance statistic.'; title8 'Hubefsiz=Effect size index presented in Huberty (1994). Calculated by'; title9 "subtracting the expected value of R-squared from Huberty's adjusted"; title10 'R-squared value.'; run;