using Dew.Math;
using Dew.Math.Units;
using Dew.Stats.Units;
namespace Dew.Examples
{
  private void Example()
  {
    Matrix data = new Matrix(0,0);
    data.LoadFromFile("process_data.mtx");

    Vector sh = new Vector(0);
    Vector sl = new Vector(0);
    Vector s = new Vector(0);
    Vector outofcontrol = new Vector(0);
    // estimate k=0.32, h = 4.7
    double k = 0.32;
    double h = 4.7;
    // estimate process mean and sigma
    double m = 230.3;
    double s = 5.2;
    StatControlCharts.QCCumSumChart(data,s,k,h,m,s,sh,sl);
    // find point indexes which exceed h (out-of-control points)
    outofcontrol.FindIndexes(shigh,">",h);
  }
}