wagedata = read.table("C:\\Users\\DELL\\Desktop\\wagedata.csv")

regression_36 <- lm(log(wage) ~ educ, data = wagedata)
summary(regression_36) #beta1(intercept) 0.583773 beta2 0.082744 (slope)
summary(regression_36)$r.squared

resid <- regression_36$residuals
y <- log(wagedata[1])
r_squared_bar <- 1 - (sum(resid^2)/sum(y^2))
print (r_squared_bar) # 0.9212892

#0.1858065 <=  0.9212892