static void Main(string[] args)
{
Mtb.Application mtbApp = new Mtb.Application();
Mtb.Project proj = mtbApp.ActiveProject;
Mtb.Worksheet ws = proj.ActiveWorksheet;
mtbApp.UserInterface.Visible = true;
mtbApp.UserInterface.DisplayAlerts = false;
proj.ExecuteCommand("rand 10 c1");
for (int i = 0; i < 1100; i++)
{
Console.WriteLine(i + 1);
proj.ExecuteCommand(string.Format("Note line {0}", i + 1));
proj.Commands.Delete();
}
Console.WriteLine("Done");
Console.ReadKey();
}
private void buttonTest_Click(object sender, EventArgs e)
{
Mtb.Application MtbApp = new Mtb.Application();
Mtb.Project MtbProject = MtbApp.ActiveProject;
Mtb.Worksheet MtbWorksheet = MtbProject.ActiveWorksheet;
MtbApp.UserInterface.Visible = false;
MtbApp.UserInterface.DisplayAlerts = false;
try
{
//need to exe like %'c:\...\...\xx.mac' "99000R99" "1"
string macro = "c:\\vm\\test.mac";
string robotname = "24010R01";
string appl = "Hand";
string subgroup = "Clamp";
MtbProject.ExecuteCommand("%'" + macro + "' \"" + robotname + "%\" \"" + appl + "%\" \"" + subgroup + "%\"");
//MtbProject.ExecuteCommand("%'" + textBoxMacro.Text + "' \"" + textBoxRobot.Text + "\"");
}
catch (COMException ex)
{
MessageBox.Show(ex.Message, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
textBoxOutput.Text = "Commands count: " + MtbProject.Commands.Count + Environment.NewLine;
for (int i = 1; i <= MtbProject.Commands.Count; i++)
{
textBoxOutput.Text += "*******************************************************************************" + Environment.NewLine;
textBoxOutput.Text += "*******************************************************************************" + Environment.NewLine;
textBoxOutput.Text += "Command " + i + Environment.NewLine;
textBoxOutput.Text += "Command name: " + MtbProject.Commands.Item(i).Name + Environment.NewLine;
textBoxOutput.Text += "Outputs count: " + MtbProject.Commands.Item(i).Outputs.Count + Environment.NewLine;
for (int j = 1; j <= MtbProject.Commands.Item(i).Outputs.Count; j++)
{
textBoxOutput.Text += "-------------------------------------------------------------------------------" + Environment.NewLine;
textBoxOutput.Text += "Output " + j + Environment.NewLine;
textBoxOutput.Text += "Output type: " + MtbProject.Commands.Item(i).Outputs.Item(j).OutputType + Environment.NewLine;
textBoxOutput.Text += "Output text: " + MtbProject.Commands.Item(i).Outputs.Item(j).Text + Environment.NewLine;
if (MtbProject.Commands.Item(i).Outputs.Item(j).OutputType == Mtb.MtbOutputTypes.OTGraph)
{
MtbProject.Commands.Item(i).Outputs.Item(j).Graph.CopyToClipboard();
pictureBoxChart.Image = Clipboard.GetImage();
}
}
}
while (MtbProject.Commands.Count > 0)
{
MtbProject.Commands.Remove(1);
}
}
public override void Execute(Mtb.Project proj)
{
// 將資料匯入 Minitab
if (RawData == null || RawData.Rows.Count == 0)
{
return;
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉
proj.Commands.Delete();
proj.Worksheets.Delete();
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
Mtb.Column bkCol = ws.Columns.Item("BK"); bkCol.Name = "斷絲率";
Mtb.Column grouCol = ws.Columns.Item("GROUP_DESC"); grouCol.Name = "紡位群";
grouCol.SetValueOrder(Mtb.MtbValueOrderTypes.WorksheetOrder);
//如果沒有資料或欄位長度不一致就離開
if ((int)bkCol.DataType == 3 || (int)grouCol.DataType == 3 ||
(bkCol.RowCount != grouCol.RowCount))
{
return;
}
StringBuilder cmnd = new StringBuilder();
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
//繪製Boxplot
cmnd.AppendFormat("boxplot ({0})*{1};\r\n", bkCol.SynthesizedName, grouCol.SynthesizedName);
string gpath = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("posidiff_{0}.jpg", _rawdata.TableName));
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendLine("iqrbox;");
cmnd.AppendLine("outl.");
//計算統計量
cmnd.AppendFormat("stat {0};\r\n", bkCol.SynthesizedName);
cmnd.AppendFormat("by {0};\r\n", grouCol.SynthesizedName);
string[] colstr = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 5, Mtblib.Tools.MtbVarType.Column);
Mtb.Column
gvalCol = ws.Columns.Item(colstr[0]),
meanCol = ws.Columns.Item(colstr[1]),
stdevCol = ws.Columns.Item(colstr[2]),
minCol = ws.Columns.Item(colstr[3]),
maxCol = ws.Columns.Item(colstr[4]);
cmnd.AppendFormat("mean {0};\r\n", meanCol.SynthesizedName);
cmnd.AppendFormat("stdev {0};\r\n", stdevCol.SynthesizedName);
cmnd.AppendFormat("mini {0};\r\n", minCol.SynthesizedName);
cmnd.AppendFormat("maxi {0};\r\n", maxCol.SynthesizedName);
cmnd.AppendFormat("gval {0}.\r\n", gvalCol.SynthesizedName);
string mPath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("Exec \"{0}\" 1", mPath), ws);
if (File.Exists(gpath)) //如果檔案存在就加入報表內容
{
_rptLst.Add(
new RptOutput
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath)
}
);
}
if (gvalCol.RowCount > 0 && (int)gvalCol.DataType != 3) // 有產出報表就加入報表
{
DataTable descStatTable = Mtblib.Tools.MtbTools.GetDataTableFromMtbCols(
new Mtb.Column[] { gvalCol, meanCol, stdevCol, minCol, maxCol });
if (descStatTable != null)
{
descStatTable.Columns[0].ColumnName = "紡位群組";
descStatTable.Columns[1].ColumnName = "平均";
descStatTable.Columns[2].ColumnName = "標準差";
descStatTable.Columns[3].ColumnName = "最小值";
descStatTable.Columns[4].ColumnName = "最大值";
_rptLst.Add(
new RptOutput
{
OType = MtbOType.TABLE,
OutputInByteArr = Tool.ConvertDataSetToByteArray(descStatTable)
}
);
}
}
}
public void GeneratePicturesInstance(Mtb.Application MtbApp, Models.Instance modelInstance, Models.RowData modelRowData)
{
Mtb.Project MtbProj = MtbApp.ActiveProject;
parseColumn(modelInstance, modelRowData);
//计算标准差
CalculateSTDVE(column3, modelInstance);
Mtb.Columns MtbColumns = MtbProj.ActiveWorksheet.Columns;
Mtb.Column MtbColumn1 = MtbColumns.Add(null, null, 1);
MtbColumn1.SetData(column1.ToArray());
Mtb.Column MtbColumn2 = MtbColumns.Add(null, null, 1);
MtbColumn2.SetData(column2.ToArray());
Mtb.Column MtbColumn3 = MtbColumns.Add(null, null, 1);
MtbColumn3.SetData(column3.ToArray());
try
{
string imgPath = Path.Combine(modelRowData.filePath, parseFileName(modelInstance.title) + " Process Capability");
string lspec = Double.IsNaN(modelInstance.LCL) ? "" : "Lspec " + modelInstance.LCL.ToString("f3") + ";";
string uspec = Double.IsNaN(modelInstance.UCL) ? "" : "Uspec " + modelInstance.UCL.ToString("f3") + ";";
MtbProj.ExecuteCommand(" Capa C3 " + column1.Count + "; " + lspec + uspec + " Pooled; AMR; UnBiased; OBiased; Toler 6; Within; Percent; Title \"" + getPictureTitle(0, modelInstance) + "\"; CStat.");
Mtb.Graph MtbGraph = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1).Graph;
MtbGraph.SaveAs(imgPath, true, Mtb.MtbGraphFileTypes.GFPNGHighColor, 768, 531);
modelInstance.pictures.Add(imgPath + ".png");
}
catch
{
modelInstance.pictures.Add(null);
}
try
{
string imgPath = Path.Combine(modelRowData.filePath, parseFileName(modelInstance.title) + " Individual Polt");
MtbProj.ExecuteCommand(" Indplot ( C3 ) * C1; Title \"" + getPictureTitle(1, modelInstance) + "\"; Individual.");
Mtb.Graph MtbGraph2 = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1).Graph;
MtbGraph2.SaveAs(imgPath, true, Mtb.MtbGraphFileTypes.GFPNGHighColor, 768, 531);
modelInstance.pictures.Add(imgPath + ".png");
}
catch
{
modelInstance.pictures.Add(null);
}
try
{
string imgPath = Path.Combine(modelRowData.filePath, parseFileName(modelInstance.title) + " Scatter Plot");
MtbProj.ExecuteCommand(" Plot C3*C2; Symbol C1; Title \"" + getPictureTitle(2, modelInstance) + "\"; JITTER.");
Mtb.Graph MtbGraph3 = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1).Graph;
MtbGraph3.SaveAs(imgPath, true, Mtb.MtbGraphFileTypes.GFPNGHighColor, 768, 531);
modelInstance.pictures.Add(imgPath + ".png");
}
catch
{
modelInstance.pictures.Add(null);
}
try
{
string imgPath = Path.Combine(modelRowData.filePath, parseFileName(modelInstance.title) + " Probability Plot");
MtbProj.ExecuteCommand(" PPlot C3; Normal; Symbol; FitD; NoCI; Grid 2; Grid 1; MGrid 1; Title \"" + getPictureTitle(3, modelInstance) + "\".");
Mtb.Graph MtbGraph4 = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1).Graph;
MtbGraph4.SaveAs(imgPath, true, Mtb.MtbGraphFileTypes.GFPNGHighColor, 768, 531);
modelInstance.pictures.Add(imgPath + ".png");
}
catch
{
modelInstance.pictures.Add(null);
}
MtbProj.Delete();
MtbApp.New();
}
public override void Execute(Mtb.Project proj)
{
// 將資料匯入 Minitab
if (_rawdata == null || _rawdata.Rows.Count == 0)
{
throw new ArgumentNullException("查無對應資料");
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉
proj.Commands.Delete();
proj.Worksheets.Delete();
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
StringBuilder cmnd = new StringBuilder();
/*
* Create TSPlot
* 該圖會將數據於多個平面中顯示(k*1的方式)，所以要先設定 Panel 的安排方式
*
*/
cmnd.AppendLine("Unstack ('DateTime' 'Temperature');");
cmnd.AppendLine("Subscripts 'Area';");
cmnd.AppendLine("After;");
cmnd.AppendLine("VarNames.");
for (int i = 0; i < CountPlotArea; i++)
{
cmnd.AppendLine(string.Format("NAME C{0} '{1}'.", 4 + 2 * (i + 1), AreaArray[i]));
}
string tmpCmnd = "Tsplot ";
for (int i = 0; i < CountPlotArea; i++)
{
tmpCmnd = tmpCmnd + string.Format("C{0} ", (i + 3) * 2);
}
tmpCmnd = tmpCmnd + ";";
cmnd.AppendLine(tmpCmnd);
string gpath_tsplot = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("trend_{0}.jpg", _rawdata.TableName));
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath_tsplot);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendLine("Scale 2;");
cmnd.AppendLine("Min 100;");
cmnd.AppendLine("Max 600;");
cmnd.AppendLine("Min 100;");
cmnd.AppendLine("Max 600;");
cmnd.AppendLine("Same 2;");
cmnd.AppendLine("Panel;");
cmnd.AppendLine("Index;");
cmnd.AppendLine("Stamp C3;");
cmnd.AppendLine("Connect;");
cmnd.AppendLine("Symbol;");
cmnd.AppendLine("Type 6;");
cmnd.AppendLine("Color 4;");
cmnd.AppendLine("Size 0.7;");
cmnd.AppendLine("Title;");
cmnd.AppendLine("Footnote;");
cmnd.AppendLine("FPer;");
cmnd.AppendLine("Footnote;");
cmnd.AppendLine("FPanel;");
cmnd.AppendLine("NoDTitle;");
cmnd.AppendLine("NoPerFootnote.");
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
string fpath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", fpath), ws);
//將檔案轉為二進位陣列
if (File.Exists(gpath_tsplot))
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath_tsplot),
Tag = "TSPlot"
});
}
//Tsplot 'Temperature_1' 'Temperature_2';
//Scale 2;
// Min 100;
// Max 600;
// Min 100;
// Max 600;
//Same 2;
//Panel;
//Index;
//Stamp 'Plant_2';
//Connect;
//Symbol;
//Title;
//Footnote;
// FPer;
//Footnote;
// FPanel;
//NoDTitle;
//NoPerFootnote.
this.Tag = "TSPlot";
}
public override void Execute(Mtb.Project proj)
{
if (_rawdata == null || _rawdata.Rows.Count == 0)
{
return;
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
List <string> varnames = new List <string>();
List <Mtb.Column> varCols = new List <Mtb.Column>(); //變數的欄位集合
for (int i = 0; i < _rawdata.Columns.Count; i++)
{
DataColumn col = _rawdata.Columns[i];
switch (col.ColumnName)
{
case "TIMESTAMP":
case "CHART_PARA_INDEX":
break;
default:
varnames.Add(col.ColumnName);
varCols.Add(ws.Columns.Item(col.ColumnName));
break;
}
}
Mtb.Column timeCol = ws.Columns.Item("TIMESTAMP");
foreach (var col in varCols)
{
if ((int)col.DataType == 3 || col.MissingCount == col.RowCount)
{
throw new ArgumentNullException(string.Format("[{0}]查無資料-多變量管制圖", col.Name));
}
}
StringBuilder cmnd = new StringBuilder();
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
cmnd.AppendFormat("fdate {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("format(dtyyyy-MM-dd hh:mm).");
cmnd.AppendFormat("tsplot {0};\r\n", string.Join(" ", varCols.Select(x => x.SynthesizedName)));
double gHeight = 4;
if (varCols.Count > 4)
{
gHeight = Math.Min(30, 4 + ((double)varCols.Count - 4) * 0.8);
}
cmnd.AppendFormat("graph 8 {0};\r\n", gHeight);
string gpath = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("overlytrnd_{0}.jpg", _rawdata.TableName));
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendFormat("stamp {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("axlab 1;");
cmnd.AppendLine("adis 0;");
cmnd.AppendLine("scale 1;");
cmnd.AppendLine(" psize 7;");
cmnd.AppendLine(" angle 90;");
cmnd.AppendLine("scale 2;");
cmnd.AppendLine(" psize 8;");
cmnd.AppendLine("symbol;");
cmnd.AppendLine("conn;");
cmnd.AppendLine("pane;");
cmnd.AppendFormat("rc {0} 1;\r\n", varCols.Count);
cmnd.AppendLine("label;");
cmnd.AppendLine(" psize 8;");
cmnd.AppendLine("noal;");
cmnd.AppendLine("nodt.");
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
string macroPath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", macroPath), ws);
//將檔案轉為二進位陣列
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath)
});
}
/// <summary>
/// 計算穩定狀態下的參數組(單一觀測值)，使用一般的 Variance-Covariance Matrix 去計算，而非 Sullivan-Woodall
/// </summary>
/// <param name="data">欲使用的資料表</param>
/// <returns></returns>
public TsquareParameters Execute(Algebra.Matrix <double> data, Mtb.Project proj)
{
//確認基本資訊能否計算
int p = data.ColumnCount;
int m = data.RowCount;
if (m - p <= 0)
{
throw new Exception(string.Format("樣本數不可低於{0}，請重新確認。", p));
}
//組合管制界線指令並計算 Phase1 的 limit
double ucl;
StringBuilder cmnd = new StringBuilder();
cmnd.AppendLine("invCdf 0.9973 k1;");
cmnd.AppendFormat("Beta {0} {1}.\r\n", (double)p / 2, (double)(m - p - 1) / 2);
cmnd.AppendFormat("let k1=(({0}-1)**2)/{0}*k1\r\n", m);
try
{
Mtb.Worksheet ws = proj.Worksheets.Add(1);
string fpath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("Exec \"{0}\" 1", fpath), ws);
ucl = ws.Constants.Item("k1").GetData();
ws.Delete();
proj.Commands.Delete();
}
catch (Exception e)
{
throw new Exception(e.Message);
}
double miss = Mtblib.Tools.MtbTools.MISSINGVALUE;
var M = Algebra.Matrix <double> .Build;
//計算 Mean vector & Covariance matrix
IEnumerable <Algebra.Vector <double> > cols = data.EnumerateColumns();
var mean = M.DenseOfColumnMajor(1, data.ColumnCount, cols.Select(x => x.Where(o => o < miss).Average()));
var subData = M.DenseOfRows((data.EnumerateRows().Where(x => x.All(y => y < miss)).ToArray()));
var diff = subData - M.Dense(subData.RowCount, 1, (i, j) => 1).Multiply(mean);
var cov = diff.Transpose().Multiply(diff) / (subData.RowCount - 1);
var invS = cov.Inverse();
TsquareParameters tsquarePara = new TsquareParameters()
{
Mean = mean,
Covariance = cov,
SampleSize = m,
SubgroupSize = 1
};
//計算出資料的Tsquare value
List <double> t2 = new List <double>();
for (int i = 0; i < data.RowCount; i++)
{
if (data.Row(i).ToArray().All(x => x < miss))
{
t2.Add(Tool.CalculateTSquare(data.Row(i).ToRowMatrix(), mean, invS));
}
else
{
t2.Add(miss);
}
}
if (t2.Any(x => x >= ucl && x < miss))
{
int[] oocRow = t2.Select((x, i) => new { Tsquare = x, Index = i }).Where(x => x.Tsquare >= ucl && x.Tsquare < miss).Select(x => x.Index).OrderByDescending(x => x).ToArray();
subData = data.Clone();
//逐列移除OOC的資料列
for (int i = 0; i < oocRow.Length; i++)
{
subData = subData.RemoveRow(oocRow[i]);
}
tsquarePara = Execute(subData, proj);
}
return(tsquarePara);
}
public override void Execute(Mtb.Project proj)
{
// 將資料匯入 Minitab
if (_rawdata == null || _rawdata.Rows.Count == 0)
{
throw new ArgumentNullException("查無對應資料");
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉，因為 Cotent 唯讀
proj.Commands.Delete();
proj.Worksheets.Delete();
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
List <Mtb.Column> varCols = new List <Mtb.Column>(); //變數的欄位集合
varCols.Add(ws.Columns.Item("VALUE"));
Mtb.Column timeCol = ws.Columns.Item("TEST_DATE");
#region Create Minitab command
StringBuilder cmnd = new StringBuilder();
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
cmnd.AppendFormat("fdate {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("format(dtyyyy-MM-dd).");
cmnd.AppendFormat("tsplot {0};\r\n", string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
cmnd.AppendFormat(" stamp {0};\r\n", timeCol.SynthesizedName);
string propName = ws.Columns.Item("PROD_NAME").GetData(1, 1);
string compName = ws.Columns.Item("ITEM_NAME").GetData(1, 1);
string unitName = ws.Columns.Item("ITEM_UNIT").GetData(1, 1);
string gpath = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("trend_{0}_{1}.jpg", propName, compName));
cmnd.AppendFormat(" gsave \"{0}\";\r\n", gpath);
cmnd.AppendLine(" repl;");
cmnd.AppendLine(" jpeg;");
cmnd.AppendFormat(" axlab 2 \"{0}({1})\";\r\n", compName, unitName);
cmnd.AppendLine(" scale 1;");
cmnd.AppendLine(" angle 90;");
cmnd.AppendLine(" axlab 1;");
cmnd.AppendLine(" adis 0;");
cmnd.AppendLine(" conn;");
cmnd.AppendLine(" symb;");
cmnd.AppendLine(" size 2;");
cmnd.AppendLine(" nodt.");
cmnd.AppendLine("brief 2");
cmnd.AppendFormat("Desc {0};\r\n", string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
//cmnd.AppendFormat(" by {0};\r\n");
cmnd.AppendLine(" n;");
cmnd.AppendLine(" mean;");
cmnd.AppendLine(" stdev;");
cmnd.AppendLine(" mini;");
cmnd.AppendLine(" maxi;");
cmnd.AppendLine(" ztag \"DESC\".");
cmnd.AppendLine("brief 0");
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
#endregion
string fpath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", fpath), ws);
//取得圖和敘述統計量
try
{
if (File.Exists(gpath))
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath),
Tag = "Trend"
});
}
}
catch (Exception ex)
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = null,
Tag = "Trend"
});
}
#region 取得敘述統計量表格
try
{
//將報表結果放到 DataTable
DataTable dt = new DataTable();
//取得結果的 html string
string htmlCmnd = proj.Commands.Cast <Mtb.Command>().Where(x => x.Tag == "DESC").Select(x => x.Outputs).FirstOrDefault().Item(2).HTMLText;
var htmlDoc = new HtmlDocument();
htmlDoc.LoadHtml(htmlCmnd);
//取得並設定表頭
htmlDoc.LoadHtml(htmlDoc.DocumentNode.SelectSingleNode("//tr").OuterHtml);
var htmlHeader = htmlDoc.DocumentNode.SelectSingleNode("//tr");
foreach (var cell in htmlHeader.SelectNodes("//td"))
{
if (cell.InnerText != "\n")
{
dt.Columns.Add(cell.InnerText);
}
}
//取得表格內容
htmlDoc.LoadHtml(htmlCmnd);
var htmlRows = htmlDoc.DocumentNode.SelectNodes("//table/tr");
DataRow dr;
HtmlDocument tmpHtmlDoc = new HtmlDocument();
for (int i = 1; i < htmlRows.Count; i++)
{
dr = dt.NewRow();
var row = htmlRows[i];
tmpHtmlDoc.LoadHtml(row.OuterHtml);
var xx = tmpHtmlDoc.DocumentNode.SelectNodes("//td").Select(x => x.InnerText).ToArray();
dr.ItemArray = xx;
dt.Rows.Add(dr);
}
//轉成 DataTable
this.Contents.Add(new RptOutput()
{
OType = MtbOType.TABLE,
OutputInByteArr = Tool.ConvertDataSetToByteArray(dt),
Tag = "Table"
});
}
catch (Exception ex)
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.TABLE,
OutputInByteArr = null,
Tag = "Table"
});
}
#endregion
}
public Dictionary <string, object> MtbAnalysis(string chartType, Dictionary <string, string> rules)
{
try
{
foreach (Process proc in Process.GetProcessesByName("Mtb"))
{
proc.Kill();
}
}
catch (Exception ex)
{}
System.IO.FileInfo imgFile = new System.IO.FileInfo(imgDir + this.controlId + "_" + chartType + ".PNG");
if (imgFile.Exists)
{
imgFile.Delete();
}
Mtb.Application MtbApp = new Mtb.Application();
Dictionary <string, object> resultObject = new Dictionary <string, object>();
try
{
MtbApp.UserInterface.Visible = false;
Mtb.Project MtbProj = MtbApp.ActiveProject;
Mtb.Columns MtbColumns = MtbProj.ActiveWorksheet.Columns;
//加载dt所有数据到mtb
for (int i = 0; i < endColumn; i++)
{
Mtb.Column MtbColumn1 = MtbColumns.Add(null, null, 1);
ArrayList columnData = new ArrayList();
for (int j = 0; j < dt.Rows.Count; j++)
{
if (dt.Rows[j][i].ToString().Length < 1)
{
columnData.Add(DBNull.Value);
}
else
{
columnData.Add(dt.Rows[j][i]);
}
}
MtbColumn1.SetData(columnData.ToArray());
}
int rawColumnCount = MtbColumns.Count;
//Minitab 控制图 or 能力分析命令。
//命令先生成，先调用checkRule()方法检查规则数量，以及对应生成的检验结果列数
string execAnalysis = this.GetExcuteCommands(chartType, rules);
//Minitab 分析结果存储列命令
string varCommand = this.GetVarCommands(chartType);
if (!string.IsNullOrEmpty(varCommand))
{
MtbProj.ExecuteCommand(varCommand);
}
if (!string.IsNullOrEmpty(execAnalysis))
{
MtbProj.ExecuteCommandAsync(execAnalysis);
}
Mtb.Output MtbOut = null;
try {
MtbOut = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1);
Console.WriteLine("Normal output!Normal output!Normal output!");
}
catch (Exception ex)
{
MtbOut = MtbProj.Commands.Item(MtbProj.Commands.Count).Outputs.Item(1);
Console.WriteLine("Unormal output!Unormal output!Unormal output!");
}
//Minitab 图形保存
Mtb.Graph MtbGraph = MtbOut.Graph;
MtbGraph.SaveAs(imgDir + this.controlId + "_" + chartType, true, Mtb.MtbGraphFileTypes.GFPNGHighColor, 790, 500);
//Minitab 工作区表保存
Mtb.Worksheet MtbSheet = MtbProj.ActiveWorksheet;
MtbSheet.SaveAs(imgDir + this.controlId + "_" + chartType, true);
for (int i = 1; i <= MtbSheet.Columns.Count; i++)
{
resultObject.Add(MtbSheet.Columns.Item(i).Name, MtbSheet.Columns.Item(i).GetData());
}
if (chartType.Contains("Chart"))
{
resultObject.Add("RuleOut", MtbProj.Commands.OutputDocument.Text);
}
MtbProj.Delete();
MtbApp.Quit();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
finally { MtbApp.Quit(); }
return(resultObject);
}
public override void Execute(Mtb.Project proj)
{
// 將資料匯入 Minitab
if (_rawdata == null || _rawdata.Rows.Count == 0)
{
throw new ArgumentNullException("查無對應資料");
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉
proj.Commands.Delete();
proj.Worksheets.Delete();
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
List <Mtb.Column> varCols = new List <Mtb.Column>(); //變數的欄位集合
List <Mtb.Column> bkCol = new List <Mtb.Column>(); //斷絲率的集合
Mtb.Column timeCol = null;
for (int i = 0; i < _rawdata.Columns.Count; i++)
{
DataColumn col = _rawdata.Columns[i];
switch (col.ColumnName)
{
case "TIMESTAMP":
timeCol = ws.Columns.Item(col.ColumnName);
break;
case "RPT_TIMEHOUR":
case "GROUP_ID":
break;
case "BK":
bkCol.Add(ws.Columns.Item(col.ColumnName));
break;
default:
varCols.Add(ws.Columns.Item(col.ColumnName));
break;
}
}
foreach (var col in bkCol)
{
if ((int)col.DataType == 3 || col.MissingCount == col.RowCount)
{
throw new ArgumentNullException(string.Format("[{0}]查無資料", col.Name));
}
}
foreach (var col in varCols)
{
if ((int)col.DataType == 3 || col.MissingCount == col.RowCount)
{
throw new ArgumentNullException(string.Format("[{0}]查無資料", col.Name));
}
}
StringBuilder cmnd = new StringBuilder();
/*
* Create TSPlot
* 該圖會將數據於多個平面中顯示(k*1的方式)，所以要先設定 Panel 的安排方式
*
*/
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
cmnd.AppendFormat("fdate {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("format(dtyyyy-MM-dd hh:mm).");
cmnd.AppendFormat("tsplot {0} {1};\r\n", string.Join(" &\r\n", bkCol.Select(x => x.SynthesizedName)),
string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
double gHeight = 4;
if (varCols.Count > 3)
{
gHeight = Math.Min(30, 4 + ((double)varCols.Count - 3) * 0.8);
}
cmnd.AppendFormat("graph 8 {0};\r\n", gHeight);
string gpath_trend = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("trend_{0}.jpg", _rawdata.TableName));
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath_trend);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendLine("panel;");
cmnd.AppendFormat("rc {0} 1;\r\n", varCols.Count + bkCol.Count);
cmnd.AppendLine("noal;");
cmnd.AppendLine("label;");
cmnd.AppendLine(" psize 8;");
cmnd.AppendLine("scale 1;");
cmnd.AppendLine(" psize 8;");
cmnd.AppendLine("axlab 1;");
cmnd.AppendLine("lshow;");
cmnd.AppendLine(" angle 90;");
cmnd.AppendLine("scale 2;");
cmnd.AppendLine(" psize 8;");
if (timeCol != null && (int)timeCol.DataType != 3)
{
cmnd.AppendFormat("stamp {0};\r\n", timeCol.SynthesizedName);
}
cmnd.AppendLine("symb;");
cmnd.AppendLine("conn;");
//cmnd.AppendFormat("foot \"建立時間: {0}\";\r\n", DateTime.Now);
cmnd.AppendLine("nodt.");
//Create Matrix Plot
cmnd.AppendFormat("matrixplot ({0})*({1});\r\n", string.Join(" &\r\n", bkCol.Select(x => x.SynthesizedName)),
string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
double gWidth = 8;
if (varCols.Count > 5)
{
gWidth = Math.Min(30, 8 + ((double)varCols.Count - 5) * 0.8);
}
cmnd.AppendFormat("graph {0} 5;\r\n", gWidth);
string gpath_corr = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("corr_{0}.jpg", _rawdata.TableName));
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath_corr);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendLine("scale 1;");
cmnd.AppendLine(" psize 8;");
cmnd.AppendLine("nojitter;");
cmnd.AppendLine("full;");
cmnd.AppendLine("bound;");
cmnd.AppendLine("noal;");
cmnd.AppendLine("regr;");
cmnd.AppendLine("symb;");
cmnd.AppendLine("nodt.");
//Calculate correlation coefficient
string[] matIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 1, Mtblib.Tools.MtbVarType.Matrix);
cmnd.AppendFormat("corr {0} {1} {2}\r\n", string.Join(" &\r\n", bkCol.Select(x => x.SynthesizedName)),
string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)),
string.Join(" &\r\n", matIds)
);
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
string fpath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", fpath), ws);
//將檔案轉為二進位陣列
if (File.Exists(gpath_trend))
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath_trend),
Tag = "Trend"
});
}
if (File.Exists(gpath_corr))
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath_corr),
Tag = "Scatter"
});
}
//Create Correlation table
DataTable corrTable = new DataTable();
corrTable.Columns.Add("熔爐項目", typeof(string));
corrTable.Columns.Add("與斷絲率相關係數", typeof(double));
var corrMat = ws.Matrices.Item(matIds[0]).GetData();
double[] subCorr = new double[varCols.Count];
double[] corr = corrMat as double[];
if (corr != null)
{
Array.Copy(corr, 1, subCorr, 0, varCols.Count);
string[] labNames = varCols.Select(x => x.Name).ToArray();
DataRow dr;
for (int i = 0; i < subCorr.Length; i++)
{
dr = corrTable.NewRow();
dr[0] = labNames[i];
dr[1] = (object)Math.Round(subCorr[i], 3);
corrTable.Rows.Add(dr);
}
}
this.Contents.Add(new RptOutput()
{
OType = MtbOType.TABLE,
OutputInByteArr = Tool.ConvertDataSetToByteArray(corrTable)
});
}
public override void Execute(Mtb.Project proj)
{
// 將資料匯入 Minitab
if (_rawdata == null || _rawdata.Rows.Count == 0)
{
return;
}
_rptLst = new List <IRptOutput>(); //重新建立一個分析結果列舉
Mtb.Worksheet ws = proj.Worksheets.Add(1); //新增工作表
Mtblib.Tools.MtbTools.InsertDataTableToMtbWs(_rawdata, ws); //匯入資料至 Minitab
List <string> varnames = new List <string>();
List <Mtb.Column> varCols = new List <Mtb.Column>(); //變數的欄位集合
for (int i = 0; i < _rawdata.Columns.Count; i++)
{
DataColumn col = _rawdata.Columns[i];
switch (col.ColumnName)
{
case "TIMESTAMP":
case "CHART_PARA_INDEX":
break;
default:
varnames.Add(col.ColumnName);
varCols.Add(ws.Columns.Item(col.ColumnName));
break;
}
}
foreach (var col in varCols)
{
if ((int)col.DataType == 3 || col.MissingCount == col.RowCount)
{
throw new ArgumentNullException(string.Format("[{0}]查無資料-多變量管制圖", col.Name));
}
}
StringBuilder cmnd = new StringBuilder();
List <TsquareParameters> tmpParaList = new List <TsquareParameters>(); //用於計算 decomposition
string[] colIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 4, Mtblib.Tools.MtbVarType.Column); //Plot point, CL, UCL and Test column
// 指定繪圖時需要的欄位變數
Mtb.Column
pplotCol = ws.Columns.Item(colIds[0]),
clCol = ws.Columns.Item(colIds[1]),
uclCol = ws.Columns.Item(colIds[2]),
oocCol = ws.Columns.Item(colIds[3]),
timeCol = ws.Columns.Item("TIMESTAMP");
// 計算 T2 plot points
if (_parameters == null || _parameters.Count == 0)
{
//自己算的時候需要的變數和子命令
#region Phase I
cmnd.AppendLine("macro");
cmnd.AppendLine("myt2calculator x.1-x.p;");
cmnd.AppendLine("variance cov ssiz;");
cmnd.AppendLine("meanvect location; ");
cmnd.AppendLine("ppoint t2;");
cmnd.AppendLine("climit cl ucl;");
cmnd.AppendLine("test ooc;");
cmnd.AppendLine("siglevel alpha.");
cmnd.AppendLine("mcolumn x.1-x.p t2");
cmnd.AppendLine("mcolumn loc.1-loc.p xx.1-xx.p tmp.1-tmp.3 ooc ucl cl");
cmnd.AppendLine("mconstant ssiz m alpha conf a1 a2 kk");
cmnd.AppendLine("mmatrix location diff cov tdiff invCov");
cmnd.AppendLine("default alpha = 0.0027");
cmnd.AppendLine("rnmiss x.1-x.p tmp.1");
cmnd.AppendLine("copy x.1-x.p xx.1-xx.p;");
cmnd.AppendLine("exclude;");
cmnd.AppendFormat("where \"tmp.1>0\".\r\n");
cmnd.AppendLine("cova xx.1-xx.p cov"); // Get Covariance matrix
cmnd.AppendLine("let ssiz = n(xx.1)"); // Get the sample size of covariance calculation, ALSO..THAT IS THE NUMBER OF NONMISSING OBSERVATIONS
cmnd.AppendLine("stat x.1-x.p;");
cmnd.AppendLine("mean loc.1-loc.p."); //Get mean vector.
/*
* Check if there is enough nonmissing observation to calculate control limit under
* Tsquare command (m>p+1, where m=#observation, include missing obs, p=#items)
* If no, you still need to draw plot points on graph..
* The trick is we given parameters (Mean & Covariance) and given a fake sample size, then
* you get t-sq value, becasue we don't need the CL & UCL from Minitab.
* (WE USE REGULAR COVARIANCE INSTEAD OF COVARIANCE BY Sullivan & Woodall)
*
*/
cmnd.AppendLine("if(ssiz < p+1)");
cmnd.AppendLine("let kk = p+1");
cmnd.AppendLine("else");
cmnd.AppendLine("copy ssiz kk\r\n");
cmnd.AppendLine("endif");
cmnd.AppendLine("tsquare x.1-x.p 1;");
cmnd.AppendLine("mu loc.1-loc.p;");
cmnd.AppendLine("sigma cov;");
cmnd.AppendLine("number kk;");
cmnd.AppendLine("sampsize tmp.1;"); //Get subgroup size.
cmnd.AppendLine("ppoint t2."); //Get Tsquare value.
cmnd.AppendLine("copy loc.1-loc.p location"); //Copy column to mean vector
cmnd.AppendLine("let m = sum(tmp.1)"); //Get the actual sample size (include missing observations)
cmnd.AppendLine("if (m <= p+1)"); // Check if there is enough obs to calc contril limit...
cmnd.AppendLine("let cl[m]=miss()"); //no center line
cmnd.AppendLine("let ucl[m]=miss()"); //no control limit
cmnd.AppendLine("set ooc"); //no ooc
cmnd.AppendLine("(0)m");
cmnd.AppendLine("end");
cmnd.AppendLine("else");
cmnd.AppendLine("let conf = 1-alpha");
cmnd.AppendLine("let a1 = p/2");
cmnd.AppendLine("let a2 = (m-p-1)/2");
// Calculate UCL
cmnd.AppendLine("set tmp.1");
cmnd.AppendLine("(conf)m");
cmnd.AppendLine("end");
cmnd.AppendLine("invcdf tmp.1 ucl;");
cmnd.AppendLine(" beta a1 a2.");
cmnd.AppendLine("let ucl = (m-1)**2/m*ucl"); //Get upper control limit
// Calculate CL
cmnd.AppendLine("set tmp.1");
cmnd.AppendLine("(0.5)m");
cmnd.AppendLine("end");
cmnd.AppendLine("invcdf tmp.1 cl;");
cmnd.AppendLine(" beta a1 a2.");
cmnd.AppendLine("let cl = (m-1)**2/m*cl"); //Get center line
cmnd.AppendLine("let ooc = if(t2>ucl and t2<>MISS(),1,0)"); //Get OOC info
cmnd.AppendLine("endif");
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
cmnd.AppendLine("endmacro");
string macPath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mytsquare.mac", cmnd.ToString());
string[] matIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 2, Mtblib.Tools.MtbVarType.Matrix); //紀錄Mean vecot & Covariance matrix
string[] constIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 1, Mtblib.Tools.MtbVarType.Constant); // Sample size of Covariance matrix
cmnd.Clear();
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
cmnd.AppendFormat("%\"{0}\" {1};\r\n", macPath, string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
cmnd.AppendFormat("variance {0} {1};\r\n", matIds[0], constIds[0]);
cmnd.AppendFormat("meanvect {0};\r\n", matIds[1]);
cmnd.AppendFormat("ppoint {0};\r\n", pplotCol.SynthesizedName);
cmnd.AppendFormat("climit {0} {1};\r\n", clCol.SynthesizedName, uclCol.SynthesizedName);
cmnd.AppendFormat("test {0}.\r\n", oocCol.SynthesizedName);
string path = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", path), ws);
#endregion
//取得參數組
TsquareParameters tmpPara = new TsquareParameters();
Mtb.Matrix mat;
mat = ws.Matrices.Item(matIds[1]);
tmpPara.Mean = LinearAlgebra.Matrix <double> .Build.DenseOfColumnMajor(1, mat.ColumnCount, mat.GetData());
mat = ws.Matrices.Item(matIds[0]);
tmpPara.Covariance = LinearAlgebra.Matrix <double> .Build.DenseOfColumnMajor(mat.RowCount, mat.ColumnCount, mat.GetData());
tmpPara.SampleSize = ws.Constants.Item(constIds[0]).GetData();
tmpPara.SubgroupSize = 1;
tmpPara.ParaID = null;
tmpParaList.Add(tmpPara);
}
else
{
//有指定參數的時候
#region Phase II
int nParas = ws.Columns.Item("CHART_PARA_INDEX").GetNumDistinctRows(); //宣告多組參數組
cmnd.AppendLine("macro");
cmnd.AppendLine("myt2calculator x.1-x.p;");
cmnd.AppendFormat("variance cov.1-cov.{0} ssiz.1-ssiz.{0};\r\n", nParas);
cmnd.AppendFormat("meanvect mvect.1-mvect.{0};\r\n", nParas);
cmnd.AppendLine("ppoint t2;");
cmnd.AppendLine("climit cl ucl;");
cmnd.AppendLine("test ooc;");
cmnd.AppendLine("siglevel alpha.");
cmnd.AppendLine("mcolumn x.1-x.p xx.1-xx.p xxx.1-xxx.p loc.1-loc.p");
cmnd.AppendLine("mcolumn tmp t2 ucl cl ooc tmpT2 tmpUCL tmpCL tmpOOC");
cmnd.AppendFormat("mmatrix mvect.1-mvect.{0} cov.1-cov.{0}\r\n", nParas);
cmnd.AppendFormat("mconstant m histm alpha conf a1 a2 ssiz.1-ssiz.{0} kk\r\n", nParas);
cmnd.AppendLine("default alpha = 0.0027");
var _chartparaindex = _rawdata.AsEnumerable()
.Select(x => x.Field <string>("CHART_PARA_INDEX")).ToArray();
var _distinctIds = _chartparaindex.Distinct().ToArray();
for (int i = 0; i < _distinctIds.Length; i++) //對每個參數與對應的資料計算 Tsquare 和管制界限
{
string id = _distinctIds[i];
cmnd.AppendLine("copy x.1-x.p xx.1-xx.p;");
cmnd.AppendLine("include;");
cmnd.AppendFormat("rows {0}.\r\n",
string.Join(" &\r\n",
_chartparaindex.Select((x, rowid) => new { Value = x, RowId = rowid + 1 })
.Where(x => x.Value == id).Select(x => x.RowId).ToArray()));
#region Set parameter information into Minitab macro
TsquareParameters para = new TsquareParameters();
para.Mean = null; para.Covariance = null; para.SampleSize = MISSINGVALUE; para.SubgroupSize = 1;
if (id != null)
{
para = _parameters.Where(x => x.ParaID == id).First();
}
if (para.Mean != null) //把已知的 Mean vector 寫入 Minitab
{
cmnd.AppendLine("read loc.1-loc.p");
cmnd.AppendFormat("{0}\r\n", string.Join(" &\r\n", para.Mean.Enumerate()));
cmnd.AppendLine("end");
}
else
{
cmnd.AppendLine("stat xx.1-xx.p;");
cmnd.AppendLine("mean loc.1-loc.p.");
}
if (para.Covariance != null) //把已知的 Covariance 寫入 Minitab
{
cmnd.AppendFormat("read {1} {2} cov.{0} \r\n", i + 1, para.Covariance.RowCount, para.Covariance.ColumnCount);
List <LinearAlgebra.Vector <double> > valuesByRow = para.Covariance.EnumerateRows().ToList();
for (int r = 0; r < para.Covariance.RowCount; r++)
{
cmnd.AppendFormat("{0}\r\n", string.Join(" &\r\n", valuesByRow[r]));
}
cmnd.AppendLine("end");
cmnd.AppendFormat("let ssiz.{0}={1}\r\n", i + 1, para.SampleSize);
}
else
{
cmnd.AppendLine("rnmiss xx.1-xx.p tmp");
cmnd.AppendLine("copy xx.1-xx.p xxx.1-xxx.p;");
cmnd.AppendLine("exclud;");
cmnd.AppendFormat("where \"tmp>0\".\r\n");
cmnd.AppendFormat("cova xxx.1-xxx.p cov.{0}\r\n", i + 1);
cmnd.AppendFormat("let ssiz.{0} = count(xxx.1)\r\n", i + 1);
}
#endregion
//Tsquare command
cmnd.AppendFormat("if(ssiz.{0} < p+1)\r\n", i + 1);
cmnd.AppendLine("let kk = p+1");
cmnd.AppendLine("else");
cmnd.AppendFormat("copy ssiz.{0} kk\r\n", i + 1);
cmnd.AppendLine("endif");
if (id == null && _chartparaindex.Where(x => x == id).Count() == 1) //如果只有一個觀測值且沒有參數設定
{
cmnd.AppendLine("let tmp=1");
cmnd.AppendLine("let tmpT2=miss()");
}
else
{
cmnd.AppendFormat("tsquare xx.1-xx.p {0};\r\n", para.SubgroupSize);
cmnd.AppendLine(" mu loc.1-loc.p;");
cmnd.AppendFormat(" sigma cov.{0};\r\n", i + 1);
cmnd.AppendLine(" number kk;");
cmnd.AppendLine("sampsize tmp;"); //Get subgroup size.
cmnd.AppendLine("ppoint tmpT2."); //Get Tsquare value.
}
cmnd.AppendFormat("copy loc.1-loc.p mvect.{0}\r\n", i + 1); //Copy column to mean vector
/*
* Get the value of m which use to calculate UCL, CL.
* In phase II, m is the sample size of the data used to calculate historical
* covariance matrix.
* In phase I, m is the number of observations.
*
*/
cmnd.AppendLine("let m = sum(tmp)"); //Get the actual sample size (include missing data)
if (para.Covariance != null)
{
cmnd.AppendFormat("let histm = ssiz.{0}\r\n", i + 1);
}
else
{
cmnd.AppendLine("copy m histm");
}
// Calculate control limits
if (id == null && _chartparaindex.Where(x => x == id).Count() == 1)
{
cmnd.AppendLine("let tmpUCL=miss()");
cmnd.AppendLine("let tmpCL=miss()");
}
else
{
cmnd.AppendLine("let conf = 1-alpha");
cmnd.AppendLine("let a1 = p");
cmnd.AppendLine("let a2 = histm");
#region Get UCL
cmnd.AppendLine("set tmp");
cmnd.AppendLine("(conf)m");
cmnd.AppendLine("end");
if (id == null) //phase I case
{
cmnd.AppendLine("invcdf tmp tmpUCL;");
cmnd.AppendLine(" beta a1 a2.");
cmnd.AppendLine("let tmpUCL = (m-1)**2/m*tmpUCL"); //Get upper control limit
}
else //phase II case
{
cmnd.AppendLine("invcdf tmp tmpUCL;");
cmnd.AppendLine(" f a1 a2.");
cmnd.AppendLine("let tmpUCL = p*(histm+1)*(histm-1)/histm/(histm-p)*tmpUCL"); //Get upper control limit
}
#endregion
#region Get CL
cmnd.AppendLine("set tmp");
cmnd.AppendLine("(0.5)m");
cmnd.AppendLine("end");
if (id == null)
{
cmnd.AppendLine("invcdf tmp tmpCL;");
cmnd.AppendLine(" beta a1 a2.");
cmnd.AppendLine("let tmpCL = (m-1)**2/m*tmpCL"); //Get upper control limit
}
else
{
cmnd.AppendLine("invcdf tmp tmpCL;");
cmnd.AppendLine(" f a1 a2.");
cmnd.AppendLine("let tmpCL = p*(histm+1)*(histm-1)/histm/(histm-p)*tmpCL"); //Get center line
}
#endregion
}
cmnd.AppendLine("let tmpOOC = if(tmpT2>tmpUCL AND tmpT2<>MISS(),1,0)"); //Get OOC info
if (i == 0)
{
cmnd.AppendLine("copy tmpT2 tmpUCL tmpCL tmpOOC t2 ucl cl ooc");
}
else
{
cmnd.AppendLine("stack (t2 ucl cl ooc) (tmpT2 tmpUCL tmpCL tmpOOC) (t2 ucl cl ooc)");
}
}
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
cmnd.AppendLine("endmacro");
string macPath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mytsquare.mac", cmnd.ToString());
string[] matMeanIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, nParas, Mtblib.Tools.MtbVarType.Matrix);
string[] matCovIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, nParas, Mtblib.Tools.MtbVarType.Matrix);
string[] constIds = Mtblib.Tools.MtbTools.CreateVariableStrArray(ws, 1 * nParas, Mtblib.Tools.MtbVarType.Constant);
cmnd.Clear();
cmnd.AppendLine("notitle");
cmnd.AppendLine("brief 0");
cmnd.AppendFormat("%\"{0}\" {1};\r\n", macPath, string.Join(" &\r\n", varCols.Select(x => x.SynthesizedName)));
cmnd.AppendFormat("variance {0} {1};\r\n", string.Join(" &\r\n", matCovIds), string.Join(" &\r\n", constIds));
cmnd.AppendFormat("meanvect {0};\r\n", string.Join(" &\r\n", matMeanIds));
cmnd.AppendFormat("ppoint {0};\r\n", pplotCol.SynthesizedName);
cmnd.AppendFormat("climit {0} {1};\r\n", clCol.SynthesizedName, uclCol.SynthesizedName);
cmnd.AppendFormat("test {0}.\r\n", oocCol.SynthesizedName);
string path = Mtblib.Tools.MtbTools.BuildTemporaryMacro("mymacro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", path), ws);
#endregion
//取得參數組
TsquareParameters tmpPara;
Mtb.Matrix mat;
for (int i = 0; i < _distinctIds.Length; i++)
{
tmpPara = new TsquareParameters();
tmpPara.ParaID = _distinctIds[i];
mat = ws.Matrices.Item(matMeanIds[i]);
tmpPara.Mean = LinearAlgebra.Matrix <double> .Build.DenseOfColumnMajor(mat.RowCount, mat.ColumnCount, mat.GetData());
mat = ws.Matrices.Item(matCovIds[i]);
tmpPara.Covariance = LinearAlgebra.Matrix <double> .Build.DenseOfColumnMajor(mat.RowCount, mat.ColumnCount, mat.GetData());
tmpPara.SampleSize = ws.Constants.Item(constIds[i]).GetData();
tmpParaList.Add(tmpPara);
}
}
//繪圖
//會使用 TIMESTAMP, PPOINT, CL, UCL, TEST
#region 繪製 T2 Control Chart (TSPLOT)
string gpath = System.IO.Path.Combine(Environment.GetEnvironmentVariable("tmp"), "Minitab", string.Format("tsquare_{0}.jpg", _rawdata.TableName));
cmnd.Clear();
cmnd.AppendFormat("fdate {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("format(dtyyyy-MM-dd hh:mm).");
cmnd.AppendFormat("tsplot {0} {1} {2};\r\n", pplotCol.SynthesizedName, uclCol.SynthesizedName, clCol.SynthesizedName);
cmnd.AppendFormat("gsave \"{0}\";\r\n", gpath);
cmnd.AppendLine("repl;");
cmnd.AppendLine("jpeg;");
cmnd.AppendLine("over;");
cmnd.AppendFormat("symb {0};\r\n", oocCol.SynthesizedName);
cmnd.AppendLine("type &");
double[] oocData = oocCol.GetData();
if (oocData.Any(x => x == 0))
{
cmnd.AppendLine("6 &");
}
if (oocData.Any(x => x == 1))
{
cmnd.AppendLine("12 &");
}
cmnd.AppendLine("0 0 0 0;");
cmnd.AppendLine("size 1;");
cmnd.AppendLine("color &");
if (oocData.Any(x => x == 0))
{
cmnd.AppendLine("1 &"); //r17 color 64
}
if (oocData.Any(x => x == 1))
{
cmnd.AppendLine("2 &");
}
cmnd.AppendLine(";");
cmnd.AppendLine("conn;");
cmnd.AppendLine("type 1 1 1;");
cmnd.AppendLine("color 1 2 120;"); //r17 conn:64 cl: 9, climit:8
//cmnd.AppendLine("graph;");
//cmnd.AppendLine("color 22;");
cmnd.AppendLine("nole;");
cmnd.AppendFormat("stamp {0};\r\n", timeCol.SynthesizedName);
cmnd.AppendLine("scale 1;");
cmnd.AppendFormat("tick 1:{0}/{1};\r\n", pplotCol.RowCount,
pplotCol.RowCount > 35 ? Math.Ceiling((double)pplotCol.RowCount / 35) : 1);
cmnd.AppendLine("axla 1;");
cmnd.AppendLine("adis 0;");
cmnd.AppendLine("axla 2;");
cmnd.AppendLine("adis 0;");
string ttlString = string.Join(",", varnames);
if (ttlString.Length > 23)
{
ttlString = ttlString.Substring(0, 20) + "...";
}
cmnd.AppendLine("graph 8 4;");
cmnd.AppendFormat("title \"T2管制圖 {0}\";\r\n", ttlString);
cmnd.AppendFormat("footn \"更新時間: {0}\";\r\n", DateTime.Now);
cmnd.AppendFormat("ZTag \"{0}\";\r\n", "_T2CHART");
cmnd.AppendLine(".");
//刪除所有圖形
cmnd.AppendLine("gmana;");
cmnd.AppendLine("all;");
cmnd.AppendLine("close;");
cmnd.AppendLine("nopr.");
string t2MacroPath = Mtblib.Tools.MtbTools.BuildTemporaryMacro("myt2macro.mtb", cmnd.ToString());
proj.ExecuteCommand(string.Format("exec \"{0}\" 1", t2MacroPath), ws);
#endregion
//將檔案轉為二進位陣列
this.Contents.Add(new RptOutput()
{
OType = MtbOType.GRAPH,
OutputInByteArr = File.ReadAllBytes(gpath)
});
//計算 Decomposition
if (oocData.Any(x => x == 1))
{
DataTable tmpDataTable = _rawdata.Copy();
tmpDataTable.Columns.Add("OOC", typeof(int));
for (int r = 0; r < tmpDataTable.Rows.Count; r++)
{
DataRow dr = tmpDataTable.Rows[r];
dr["OOC"] = oocData[r];
}
//建立 Decomposition 的表格
DataTable decoTable = new DataTable();
decoTable.Columns.Add("TIMESTAMP", typeof(DateTime));
foreach (var item in varnames)
{
decoTable.Columns.Add(item, typeof(double));
}
//將OOC的項目值取出
var subData = tmpDataTable.Select("OOC=1").CopyToDataTable();
var oocParaSet = subData.AsEnumerable().Select(dr => dr.Field <string>("CHART_PARA_INDEX")).Distinct().ToArray(); //OOC 的參數組有哪些
foreach (var item in tmpParaList)
{
if (oocParaSet.Contains(item.ParaID)) //該參數組的觀測值是OOC
{
//把對應的觀測值轉換成 List<double[]>，其中每個 double[] 是每個 row 的數據
var subsubData = subData.Select(string.Format("CHART_PARA_INDEX {0}", item.ParaID == null ? "IS NULL" : "='" + item.ParaID + "'"))
.CopyToDataTable();
var obsArray = subsubData.DefaultView.ToTable(false, varnames.ToArray())
.AsEnumerable().Select(x => x.ItemArray.Select(o => Convert.ToDouble(o)).ToArray()).ToList();
LinearAlgebra.Matrix <double> obs = LinearAlgebra.Matrix <double> .Build.DenseOfRowArrays(obsArray);
LinearAlgebra.Matrix <double> t2deco = Tool.T2Decomposition(obs, item.Mean, item.Covariance);
var t2decoByRow = t2deco.EnumerateRows().ToArray();
for (int r = 0; r < t2decoByRow.Count(); r++)
{
DataRow dr = decoTable.NewRow();
object[] o = new object[1 + t2deco.ColumnCount];
t2decoByRow[r].ToArray().CopyTo(o, 1);
o[0] = subsubData.Rows[r].Field <DateTime>("TIMESTAMP");
decoTable.Rows.Add(o);
}
}
}
cmnd.Clear();
//LinearAlgebra.Matrix<double> obs = LinearAlgebra.Matrix<double>.Build.DenseOfRowArrays(subData);
this.Contents.Add(new RptOutput()
{
OType = MtbOType.TABLE,
OutputInByteArr = Tool.ConvertDataSetToByteArray(decoTable)
});
}
else // 沒有 OOC 就不做
{
this.Contents.Add(new RptOutput()
{
OType = MtbOType.TABLE,
OutputInByteArr = null
});
}
//Console.ReadKey();
}