計算PRESS統計/預測Rsquared C＃

Question

我想創建一個方法，將使用C＃計算PRESS統計信息，但我顯然沒有做正確的事情，因爲即使我使用垃圾數據（它的rsquared值爲0.2 ）我仍然得到新聞統計的負面總數。有沒有人看到我在做什麼錯了？

參考如何計算按統計可以https://en.wikipedia.org/wiki/PRESS_statistic

public double CalculatePredictedRSquared(List<MultipleRegressionInfo> listMRInfo, Vector<double> vectorArray) 
    { 
     double predictedRSquared = 0, press = 0, tss = 0; 

     try 
     { 
      for (int i = 0; i < vectorArray.Count; i++) 
      { 
       var matrixArray = CreateMatrix.DenseOfColumnArrays(listMRInfo.ElementAt(0).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(1).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(2).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(3).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(4).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(5).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(6).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(7).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(8).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(9).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(10).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(11).ListValues.Where((v, j) => j != i).ToArray()); 
       var actualResult = vectorArray.ElementAt(i); 
       var newVectorArray = CreateVector.Dense(vectorArray.Where((v, j) => j != i).ToArray()); 
       var items = MultipleRegression.NormalEquations(matrixArray, newVectorArray); 

       var estimate = (items.ElementAt(0) * listMRInfo.ElementAt(0).ListValues.ElementAt(i)) + (items.ElementAt(1) * listMRInfo.ElementAt(1).ListValues.ElementAt(i)) + 
         (items.ElementAt(2) * listMRInfo.ElementAt(2).ListValues.ElementAt(i)) + (items.ElementAt(3) * listMRInfo.ElementAt(3).ListValues.ElementAt(i)) + 
         (items.ElementAt(4) * listMRInfo.ElementAt(4).ListValues.ElementAt(i)) + (items.ElementAt(5) * listMRInfo.ElementAt(5).ListValues.ElementAt(i)) + 
         (items.ElementAt(6) * listMRInfo.ElementAt(6).ListValues.ElementAt(i)) + (items.ElementAt(7) * listMRInfo.ElementAt(7).ListValues.ElementAt(i)) + 
         (items.ElementAt(8) * listMRInfo.ElementAt(8).ListValues.ElementAt(i)) + (items.ElementAt(9) * listMRInfo.ElementAt(9).ListValues.ElementAt(i)) + 
         (items.ElementAt(10) * listMRInfo.ElementAt(10).ListValues.ElementAt(i)) + (items.ElementAt(11) * listMRInfo.ElementAt(11).ListValues.ElementAt(i)); 

       press += actualResult - estimate; 
      } 

      tss += CalculateTotalSumOfSquares(vectorArray.ToList()); 
      predictedRSquared = 1 - (press/tss); 
     } 
     catch (Exception ex) 
     { 
      predictedRSquared = 0; 
      Console.WriteLine(ex.Message); 
      Console.WriteLine(ex.StackTrace); 
     } 

     return predictedRSquared; 
    } 

使用真實數據或垃圾數據按可變收益爲負數

Answer 1

答案很痛苦地簡單，我可以」可以看出我相信我沒有注意到它，但我並沒有將實際結果 - 估計和所以這裏是對任何有興趣的人的正確代碼

public double CalculatePredictedRSquared(List<MultipleRegressionInfo> listMRInfo, Vector<double> vectorArray) 
    { 
     double predictedRSquared = 0, press = 0, tss = 0; 

     try 
     { 
      for (int i = 0; i < vectorArray.Count; i++) 
      { 
       var matrixArray = CreateMatrix.DenseOfColumnArrays(listMRInfo.ElementAt(0).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(1).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(2).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(3).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(4).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(5).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(6).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(7).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(8).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(9).ListValues.Where((v, j) => j != i).ToArray(), listMRInfo.ElementAt(10).ListValues.Where((v, j) => j != i).ToArray(), 
        listMRInfo.ElementAt(11).ListValues.Where((v, j) => j != i).ToArray()); 
       var actualResult = vectorArray.ElementAt(i); 
       var newVectorArray = CreateVector.Dense(vectorArray.Where((v, j) => j != i).ToArray()); 
       var items = MultipleRegression.NormalEquations(matrixArray, newVectorArray); 
       var actualList = newVectorArray.ToList(); 
       var y = CalculateYIntercept(matrixArray, actualList, items); 

       var estimate = (items.ElementAt(0) * listMRInfo.ElementAt(0).ListValues.ElementAt(i)) + (items.ElementAt(1) * listMRInfo.ElementAt(1).ListValues.ElementAt(i)) + 
         (items.ElementAt(2) * listMRInfo.ElementAt(2).ListValues.ElementAt(i)) + (items.ElementAt(3) * listMRInfo.ElementAt(3).ListValues.ElementAt(i)) + 
         (items.ElementAt(4) * listMRInfo.ElementAt(4).ListValues.ElementAt(i)) + (items.ElementAt(5) * listMRInfo.ElementAt(5).ListValues.ElementAt(i)) + 
         (items.ElementAt(6) * listMRInfo.ElementAt(6).ListValues.ElementAt(i)) + (items.ElementAt(7) * listMRInfo.ElementAt(7).ListValues.ElementAt(i)) + 
         (items.ElementAt(8) * listMRInfo.ElementAt(8).ListValues.ElementAt(i)) + (items.ElementAt(9) * listMRInfo.ElementAt(9).ListValues.ElementAt(i)) + 
         (items.ElementAt(10) * listMRInfo.ElementAt(10).ListValues.ElementAt(i)) + (items.ElementAt(11) * listMRInfo.ElementAt(11).ListValues.ElementAt(i)) + y; 

       press += Math.Pow(actualResult - estimate, 2); 
      } 

      tss += CalculateTotalSumOfSquares(vectorArray.ToList()); 
      predictedRSquared = 1 - (press/tss); 
     } 
     catch (Exception ex) 
     { 
      predictedRSquared = 0; 
      Console.WriteLine(ex.Message); 
      Console.WriteLine(ex.StackTrace); 
     } 

     return predictedRSquared; 
    }