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))
}