/// <summary>
/// 1 - EMA of Close prices [EMA(Close)]
/// 2 - Double smooth [EMA(EMA(Close))]
/// 3 - Triple smooth [EMA(EMA(EMA(Close)))]
/// 4 - a) Calculation with percentage: [ROC(EMA(EMA(EMA(Close))))]
/// 4 - b) Calculation with percentage: [Momentum(EMA(EMA(EMA(Close))))]
/// </summary>
/// <see cref="http://stockcharts.com/school/doku.php?id=chart_school:technical_indicators:trix"/>
/// <see cref="http://www.fmlabs.com/reference/default.htm?url=TRIX.htm"/>
/// <returns></returns>
public override SingleDoubleSerie Calculate()
{
// EMA calculation
EMA ema = new EMA(Period, false);
ema.Load(OhlcList);
List <double?> emaValues = ema.Calculate().Values;
for (int i = 0; i < OhlcList.Count; i++)
{
OhlcList[i].Close = emaValues[i].HasValue ? emaValues[i].Value : 0.0;
}
// Double smooth
ema.Load(OhlcList.Skip(Period - 1).ToList());
List <double?> doubleSmoothValues = ema.Calculate().Values;
for (int i = 0; i < Period - 1; i++)
{
doubleSmoothValues.Insert(0, null);
}
for (int i = 0; i < OhlcList.Count; i++)
{
OhlcList[i].Close = doubleSmoothValues[i].HasValue ? doubleSmoothValues[i].Value : 0.0;
}
// Triple smooth
ema.Load(OhlcList.Skip(2 * (Period - 1)).ToList());
List <double?> tripleSmoothValues = ema.Calculate().Values;
for (int i = 0; i < (2 * (Period - 1)); i++)
{
tripleSmoothValues.Insert(0, null);
}
for (int i = 0; i < OhlcList.Count; i++)
{
OhlcList[i].Close = tripleSmoothValues[i].HasValue ? tripleSmoothValues[i].Value : 0.0;
}
// Last step
SingleDoubleSerie trixSerie = new SingleDoubleSerie();
if (CalculatePercentage)
{
ROC roc = new ROC(1);
roc.Load(OhlcList.Skip(3 * (Period - 1)).ToList());
trixSerie = roc.Calculate();
}
else
{
Momentum momentum = new Momentum();
momentum.Load(OhlcList.Skip(3 * (Period - 1)).ToList());
trixSerie = momentum.Calculate();
}
for (int i = 0; i < (3 * (Period - 1)); i++)
{
trixSerie.Values.Insert(0, null);
}
return(trixSerie);
}
public override ADXSerie Calculate()
{
ADXSerie adxSerie = new ADXSerie();
List <Ohlc> tempOhlcList = new List <Ohlc>();
for (int i = 0; i < OhlcList.Count; i++)
{
Ohlc tempOhlc = new Ohlc()
{
Close = OhlcList[i].High
};
tempOhlcList.Add(tempOhlc);
}
Momentum momentum = new Momentum();
momentum.Load(tempOhlcList);
List <double?> highMomentums = momentum.Calculate().Values;
tempOhlcList = new List <Ohlc>();
for (int i = 0; i < OhlcList.Count; i++)
{
Ohlc tempOhlc = new Ohlc()
{
Close = OhlcList[i].Low
};
tempOhlcList.Add(tempOhlc);
}
momentum = new Momentum();
momentum.Load(tempOhlcList);
List <double?> lowMomentums = momentum.Calculate().Values;
for (int i = 0; i < lowMomentums.Count; i++)
{
if (lowMomentums[i].HasValue)
{
lowMomentums[i] *= -1;
}
}
//DMIp <- ifelse( dH==dL | (dH< 0 & dL< 0), 0, ifelse( dH >dL, dH, 0 ) )
List <double?> DMIPositives = new List <double?>()
{
null
};
// DMIn <- ifelse( dH==dL | (dH< 0 & dL< 0), 0, ifelse( dH <dL, dL, 0 ) )
List <double?> DMINegatives = new List <double?>()
{
null
};
for (int i = 1; i < OhlcList.Count; i++)
{
if (highMomentums[i] == lowMomentums[i] || (highMomentums[i] < 0 & lowMomentums[i] < 0))
{
DMIPositives.Add(0);
}
else
{
if (highMomentums[i] > lowMomentums[i])
{
DMIPositives.Add(highMomentums[i]);
}
else
{
DMIPositives.Add(0);
}
}
if (highMomentums[i] == lowMomentums[i] || (highMomentums[i] < 0 & lowMomentums[i] < 0))
{
DMINegatives.Add(0);
}
else
{
if (highMomentums[i] < lowMomentums[i])
{
DMINegatives.Add(lowMomentums[i]);
}
else
{
DMINegatives.Add(0);
}
}
}
ATR atr = new ATR();
atr.Load(OhlcList);
List <double?> trueRanges = atr.Calculate().TrueRange;
adxSerie.TrueRange = trueRanges;
List <double?> trSum = wilderSum(trueRanges);
// DIp <- 100 * wilderSum(DMIp, n=n) / TRsum
List <double?> DIPositives = new List <double?>();
List <double?> wilderSumOfDMIp = wilderSum(DMIPositives);
for (int i = 0; i < wilderSumOfDMIp.Count; i++)
{
if (wilderSumOfDMIp[i].HasValue)
{
DIPositives.Add(wilderSumOfDMIp[i].Value * 100 / trSum[i].Value);
}
else
{
DIPositives.Add(null);
}
}
adxSerie.DIPositive = DIPositives;
// DIn <- 100 * wilderSum(DMIn, n=n) / TRsum
List <double?> DINegatives = new List <double?>();
List <double?> wilderSumOfDMIn = wilderSum(DMINegatives);
for (int i = 0; i < wilderSumOfDMIn.Count; i++)
{
if (wilderSumOfDMIn[i].HasValue)
{
DINegatives.Add(wilderSumOfDMIn[i].Value * 100 / trSum[i].Value);
}
else
{
DINegatives.Add(null);
}
}
adxSerie.DINegative = DINegatives;
// DX <- 100 * ( abs(DIp - DIn) / (DIp + DIn) )
List <double?> DX = new List <double?>();
for (int i = 0; i < OhlcList.Count; i++)
{
if (DIPositives[i].HasValue)
{
double?dx = 100 * (Math.Abs(DIPositives[i].Value - DINegatives[i].Value) / (DIPositives[i].Value + DINegatives[i].Value));
DX.Add(dx);
}
else
{
DX.Add(null);
}
}
adxSerie.DX = DX;
for (int i = 0; i < OhlcList.Count; i++)
{
if (DX[i].HasValue)
{
OhlcList[i].Close = DX[i].Value;
}
else
{
OhlcList[i].Close = 0.0;
}
}
EMA ema = new EMA(Period, true);
ema.Load(OhlcList.Skip(Period).ToList());
List <double?> emaValues = ema.Calculate().Values;
for (int i = 0; i < Period; i++)
{
emaValues.Insert(0, null);
}
adxSerie.ADX = emaValues;
return(adxSerie);
}