rsquaredstatistic<-function(valuesofdepvar, estimations){ nominator<-0 denominator<-0 for(i in 1:length(valuesofdepvar)){ nominator <- nominator + ((valuesofdepvar[i]-estimations[i])^2) denominator <- denominator + ((valuesofdepvar[i]-mean(valuesofdepvar))^2) } return(1-(nominator/denominator)) }