Индикатор на основе индикатора tradingview.com "Кривая линейной регрессии"

Нужно написать такой же для MQ4.

Вот есть текущие наработки по нему на языке MQ4, но неправильно рисует:

#property strict

//--- input parameters

input int      len=100;

input bool     useUpperDev=true;

input bool     useLowerDev=true;

input bool     showPearson=true;

input bool     extendLines=false;

//--- indicator buffers

double         Label1Buffer[];

int OnInit() {

//--- indicator buffers mapping

   SetIndexBuffer(0,Label1Buffer);

   return(INIT_SUCCEEDED);

}

void calcDev(double &src[], int len, double slope, double average, double intercept, double &result[]) {

   double upDev = 0.0;

   double dnDev = 0.0;

   double stdDevAcc = 0.0;

   double dsxx = 0.0;

   double dsyy = 0.0;

   double dsxy = 0.0;

   

   int periods = len - 1;

   

   double daY = intercept + (slope * periods) / 2;

   double val = intercept;

   

   for (int i = 0; i <= periods; i++) {

      double price = iHigh(NULL, 0, i) - val;

      if (price > upDev) {

         upDev = price;

      }

   

      price = val - iLow(NULL, 0, i);

      if (price > dnDev) {

         dnDev = price;

      }

   

      price = iClose(NULL, 0, i);

      double dxt = price - average;

      double dyt = val - daY;

     

      price = price - val;

      stdDevAcc = stdDevAcc + price * price;

      dsxx = dsxx + dxt * dxt;

      dsyy = dsyy + dyt * dyt;

      dsxy = dsxy + dxt * dyt;

      val = val + slope;

   }

   

   double stdDev = sqrt(stdDevAcc / (periods == 0 ? 1 : periods));

   double pearsonR = (dsxx == 0 || dsyy == 0 ? 0 : dsxy / sqrt(dsxx * dsyy));

   

   result[0] = stdDev;

   result[1] = pearsonR;

   result[2] = upDev;

   result[3] = dnDev;

}

void calcSlope(int len, double &result[]) {

   double sumX = 0.0;

   double sumY = 0.0;

   double sumXSqr = 0.0;

   double sumXY = 0.0;

   //for (int i = 0; i < len-1; i++) {

   for (int i = len-1; i > 0; i--) {

      double val = iClose(NULL, 0, i);

      double per = i + 1.0;

      sumX = sumX + per;

      sumY = sumY + val;

      sumXSqr = sumXSqr + per * per;

      sumXY = sumXY + val * per;

   }

   double slope = (len * sumXY - sumX * sumY) / (len * sumXSqr - sumX * sumX);

   double average = sumY / (len);

   double intercept = average - slope * sumX / len + slope;

   

   result[0] = slope;

   result[1] = average;

   result[2] = intercept;

}

int OnCalculate(const int rates_total,

                const int prev_calculated,

                const datetime &time[],

                const double &open[],

                const double &high[],

                const double &low[],

                const double &close[],

                const long &tick_volume[],

                const long &volume[],

                const int &spread[]) {

   double calcSlopeResults[3];

   double calcDevResults[4];

   

   calcSlope(len, calcSlopeResults);

   double s = calcSlopeResults[0];

   double a = calcSlopeResults[1];

   double i = calcSlopeResults[2];

   

   double startPrice = i + s * (len - 1);

   double endPrice = i;

   

   //baseLine := line.new(bar_index - len + 1, startPrice, bar_index, endPrice, width=1, extend=extend, color=color.red)

   ObjectDelete(0, "LR Base Line");

   datetime time1 = iTime(NULL, NULL, 1);

   datetime time2 = iTime(NULL, NULL, len);

   // Рисуем base line (должна быть посередине канала):

   if (!ObjectCreate(0, "LR Base Line", OBJ_TREND, 0, time1, startPrice, time2, endPrice)) { 

      Print(__FUNCTION__, ": не удалось создать линию. Код ошибки = ", GetLastError());

      return(false);

   }

   ObjectSetInteger(0, "LR Base Line", OBJPROP_RAY_RIGHT, false); 

   Print("time1: " + time1 + "; time2: " + time2 + "; startPrice: " + startPrice + "; endPrice: " + endPrice);

   return(rates_total);

}

Нужно доработать/исправить этот скрипт, чтобы он рисовал канал линейной регрессии как на tradingview.com.

Чтобы проверить правильность работы индикатора, нужно в tradingview.com открыть график любой валютной пары, нажать кнопку «индикаторы» и добавить индикатор «Кривая линейной регрессии».