private void mahalanobisDistanceToolStripMenuItem_Click(object sender, EventArgs e)
{
FormForMahalanobisDistance MainWindow = new FormForMahalanobisDistance();
PanelForClassSelection ClassSelectionPanel = new PanelForClassSelection(true, eClassType.WELL);
ClassSelectionPanel.Height = MainWindow.panelForSourceCloud.Height;
MainWindow.panelForSourceCloud.Controls.Add(ClassSelectionPanel);
PanelForClassSelection HitClassSelectionPanel = new PanelForClassSelection(false, eClassType.WELL);
HitClassSelectionPanel.Height = MainWindow.panelHitClass.Height;
HitClassSelectionPanel.ListRadioButtons[1].Checked = true;
MainWindow.panelHitClass.Controls.Add(HitClassSelectionPanel);
if (MainWindow.ShowDialog() != System.Windows.Forms.DialogResult.OK) return;
cExtendedList ListClassSelected = new cExtendedList();
for (int i = 0; i < ClassSelectionPanel.ListCheckBoxes.Count; i++)
if (ClassSelectionPanel.ListCheckBoxes[i].Checked) ListClassSelected.Add(1);
else
ListClassSelected.Add(0);
int IdxClassForOutliers = 0;
for (int i = 0; i < HitClassSelectionPanel.ListRadioButtons.Count; i++)
if (HitClassSelectionPanel.ListRadioButtons[i].Checked)
{
IdxClassForOutliers = i;
break;
}
if (ListClassSelected.Sum() < 1)
{
MessageBox.Show("At least one classe has to be selected.", "Error !", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
// start by computing the inversed covariance matrix
//if (checkedListBoxActiveDescriptors.CheckedItems.Count <= 1)
//{
// MessageBox.Show("At least two descriptors have to be selected", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
// return;
//}
cCovarianceMatrix CM = new cCovarianceMatrix();
cExtendedTable NewTable = null;
if (this.ProcessModeplateByPlateToolStripMenuItem.Checked)
{
//foreach (cPlate TmpPlate in CompleteScreening.ListPlatesActive)
//{
// cListWells ListWellsToProcess = new cListWells(null);
// foreach (cWell item in TmpPlate.ListActiveWells)
// if ((item.GetClassIdx() != -1) && (ListClassSelected[item.GetClassIdx()] == 1)) ListWellsToProcess.Add(item);
// cExtendedTable NewTable = new cExtendedTable(ListWellsToProcess, true);
// CM.SetInputData(NewTable);
//}
}
else if (this.ProcessModeEntireScreeningToolStripMenuItem.Checked)
{
cListWells ListWellsToProcess = new cListWells(null);
foreach (cPlate TmpPlate in cGlobalInfo.CurrentScreening.ListPlatesActive)
foreach (cWell item in TmpPlate.ListActiveWells)
if ((item.GetCurrentClassIdx() != -1) && (ListClassSelected[item.GetCurrentClassIdx()] == 1)) ListWellsToProcess.Add(item);
NewTable = new cExtendedTable(ListWellsToProcess, true);
// cTwoSampleFTest CM = new cTwoSampleFTest();
CM.SetInputData(NewTable);
}
else
{
cListWells ListWellsToProcess = new cListWells(null);
foreach (cWell item in cGlobalInfo.CurrentScreening.GetCurrentDisplayPlate().ListActiveWells)
if ((item.GetCurrentClassIdx() != -1) && (ListClassSelected[item.GetCurrentClassIdx()] == 1)) ListWellsToProcess.Add(item);
NewTable = new cExtendedTable(ListWellsToProcess, true);
// cTwoSampleFTest CM = new cTwoSampleFTest();
CM.SetInputData(NewTable);
}
CM.Run();
cInverse cI = new cInverse();
cI.SetInputData(CM.GetOutPut());
cI.Run();
// get the cloud center
cStatistics cstat = new cStatistics();
cstat.UnselectAll();
cstat.IsMean = true;
cstat.SetInputData(NewTable);
cstat.Run();
if (cstat.GetOutPut() == null) return;
cExtendedList ListMeans = cstat.GetOutPut().GetRow(0);
cDescriptorType MahalanobisType = new cDescriptorType("Mahalanobis Distance", true, 1);
#region Compute the Threshold
int DegreeOfFreedom = cGlobalInfo.CurrentScreening.ListDescriptors.GetActiveDescriptors().Count;
cExtendedTable T = new cExtendedTable();
T.Add(new cExtendedList());
T[0].Add((double)MainWindow.numericUpDownPValue.Value);
cFunctions F = new cFunctions();
F.SetInputData(T);
F.IsInverse = true;
F.DegreeOfFreedom = DegreeOfFreedom;
F.Run();
double ThresholdForMahalanobis = Math.Sqrt(F.GetOutPut()[1][0]);
#endregion
int IdxClassForNonOutliers = 1;
foreach (cPlate TmpPlate in cGlobalInfo.CurrentScreening.ListPlatesAvailable)
{
for (int Col = 0; Col < cGlobalInfo.CurrentScreening.Columns; Col++)
for (int Row = 0; Row < cGlobalInfo.CurrentScreening.Rows; Row++)
{
cWell TmpWell = TmpPlate.GetWell(Col, Row, false);
if (TmpWell == null) continue;
double ValueDistance = TmpWell.GetAverageValuesList(true)[0].Dist_Mahalanobis(ListMeans, cI.GetOutPut());
if (ValueDistance > ThresholdForMahalanobis || double.IsNaN(ValueDistance) ) // that's an outlier
{
TmpWell.SetClass(IdxClassForOutliers);
}
else
{
TmpWell.SetClass(IdxClassForNonOutliers);
}
if (MainWindow.checkBoxDistAsDesc.Checked)
{
cListSignature LDesc = new cListSignature();
cSignature NewDesc = new cSignature(ValueDistance, MahalanobisType, cGlobalInfo.CurrentScreening);
LDesc.Add(NewDesc);
TmpWell.AddSignatures(LDesc);
}
}
}
if (MainWindow.checkBoxDistAsDesc.Checked)
{
cGlobalInfo.CurrentScreening.ListDescriptors.AddNew(MahalanobisType);
}
cGlobalInfo.CurrentScreening.ListDescriptors.UpDateDisplay();
cGlobalInfo.CurrentScreening.UpDatePlateListWithFullAvailablePlate();
for (int idxP = 0; idxP < cGlobalInfo.CurrentScreening.ListPlatesActive.Count; idxP++)
cGlobalInfo.CurrentScreening.ListPlatesActive[idxP].UpDataMinMax();
}
private void ToolStripMenuItem_OperationsStatistics(object sender, EventArgs e)
{
cExtendedTable ET = new cExtendedTable();
ET.Name = " Statistics(" + this.Input.Name + ")";
cLinearize L = new cLinearize();
L.SetInputData(this.Input);
L.Run();
cStatistics S = new cStatistics();
S.SetInputData(L.GetOutPut());
S.Run();
cDisplayExtendedTable CDT = new cDisplayExtendedTable();
CDT.SetInputData(S.GetOutPut());
CDT.Run();
}
private void statisticsToolStripMenuItem1_Click_1(object sender, EventArgs e)
{
string NameFunction = "";
if (_StatCVItem.Checked)
NameFunction = "Coeff. of Variation";
else if (_StatMeanItem.Checked)
NameFunction = "Mean";
else if (_StatSumItem.Checked)
NameFunction = "Sum";
else if (_StatJarqueBeraItem.Checked)
NameFunction = "Jarque-Bera";
cGUI_ListClasses GUI_ListClasses = new cGUI_ListClasses();
GUI_ListClasses.IsCheckBoxes = true;
GUI_ListClasses.IsSelectAll = true;
if (GUI_ListClasses.Run().IsSucceed == false) return;
cExtendedTable ListClassSelected = new cExtendedTable(GUI_ListClasses.GetOutPut());// GetOutPut();
int IdxClass = -1;
for (int IdxC = 0; IdxC < ListClassSelected[0].Count; IdxC++)
{
if (ListClassSelected[0][IdxC] == 1) IdxClass = IdxC;
}
#region single plate and plate by plate
cDesignerTab DT = new cDesignerTab();
if ((ProcessModeCurrentPlateOnlyToolStripMenuItem.Checked) || (ProcessModeplateByPlateToolStripMenuItem.Checked)/*||(ProcessModeEntireScreeningToolStripMenuItem.Checked)*/)
{
List<cPlate> ListPlatesToProcess = new List<cPlate>();
if ((ProcessModeplateByPlateToolStripMenuItem.Checked)/*||(ProcessModeEntireScreeningToolStripMenuItem.Checked)*/)
{
foreach (cPlate TmpPlate in cGlobalInfo.CurrentScreening.ListPlatesActive)
ListPlatesToProcess.Add(TmpPlate);
}
else
ListPlatesToProcess.Add(cGlobalInfo.CurrentScreening.GetCurrentDisplayPlate());
foreach (cPlate TmpPlate in ListPlatesToProcess)
{
cListWells ListWellsToProcess1 = new cListWells(null);
foreach (cWell item in TmpPlate.ListActiveWells)
{
if (item.GetCurrentClassIdx() != -1)
{
if (ListClassSelected[0][item.GetCurrentClassIdx()] == 1)
ListWellsToProcess1.Add(item);
}
}
cExtendedTable NewTable1 = new cExtendedTable(ListWellsToProcess1, true);
if ((NewTable1.Count == 0) || (NewTable1[0].Count < 3))
{
if (ProcessModeCurrentPlateOnlyToolStripMenuItem.Checked)
{
MessageBox.Show("Insufficient number of control wells", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
else
continue;
}
cExtendedList ListValues = new cExtendedList();
List<cDescriptorType> ListDescs = new List<cDescriptorType>();
List<string> ListNames = new List<string>();
int RealIdx = 0;
for (int IDxDesc = 0; IDxDesc < cGlobalInfo.CurrentScreening.ListDescriptors.Count; IDxDesc++)
{
if (!cGlobalInfo.CurrentScreening.ListDescriptors[IDxDesc].IsActive()) continue;
cExtendedTable TableForValues = new cExtendedTable();
TableForValues.Add(NewTable1[RealIdx]);
RealIdx++;
if (_StatJarqueBeraItem.Checked)
{
cNormalityJarqueBera JB = new cNormalityJarqueBera();
JB.SetInputData(TableForValues);
JB.Run();
ListValues.Add(JB.GetOutPut()[0][0]);
}
else
{
cStatistics CS = new cStatistics();
CS.UnselectAll();
if (_StatCVItem.Checked)
CS.IsCV = true;
else if (_StatMeanItem.Checked)
CS.IsMean = true;
else if (_StatSumItem.Checked)
CS.IsSum = true;
CS.SetInputData(TableForValues);
CS.Run();
ListValues.Add(CS.GetOutPut()[0][0]);
}
ListDescs.Add(cGlobalInfo.CurrentScreening.ListDescriptors[IDxDesc]);
}
cExtendedTable ET = new cExtendedTable(new cExtendedTable(ListValues));
ET[0].ListTags = new List<object>();
ET[0].ListTags.AddRange(ListDescs);
ET.Name = TmpPlate.GetName() + "\n" + NameFunction + " - " + cGlobalInfo.ListWellClasses[IdxClass].Name + " (" + NewTable1[0].Count + " wells)";
ET[0].Name = ET.Name;
cSort S = new cSort();
S.SetInputData(ET);
S.ColumnIndexForSorting = 0;
S.Run();
//ZFactorList.Sort(delegate(cSimpleSignature p1, cSimpleSignature p2) { return p1.AverageValue.CompareTo(p2.AverageValue); });
cViewerGraph1D VG1 = new cViewerGraph1D();
VG1.SetInputData(S.GetOutPut());
VG1.Chart.LabelAxisY = NameFunction;
VG1.Chart.LabelAxisX = "Descriptor";
VG1.Chart.IsZoomableX = true;
VG1.Chart.IsBar = true;
VG1.Chart.IsBorder = true;
VG1.Chart.IsDisplayValues = true;
VG1.Chart.IsShadow = true;
VG1.Chart.MarkerSize = 4;
VG1.Title = TmpPlate.GetName();
VG1.Run();
DT.SetInputData(VG1.GetOutPut());
}
DT.Run();
cDisplayToWindow CDW = new cDisplayToWindow();
CDW.SetInputData(DT.GetOutPut());//VG1.GetOutPut());
if (ProcessModeCurrentPlateOnlyToolStripMenuItem.Checked)
CDW.Title = NameFunction + " - " + ListPlatesToProcess[0].GetName();
else
CDW.Title = NameFunction + " - " + ListPlatesToProcess.Count + " plates";
CDW.Run();
CDW.Display();
}
#endregion
#region entire screening
else if (ProcessModeEntireScreeningToolStripMenuItem.Checked)
{
List<cPlate> ListPlatesToProcess = new List<cPlate>();
foreach (cPlate TmpPlate in cGlobalInfo.CurrentScreening.ListPlatesActive)
ListPlatesToProcess.Add(TmpPlate);
cExtendedList ListZ = new cExtendedList();
List<cPlate> ListPlatesForZFactor = new List<cPlate>();
foreach (cPlate TmpPlate in ListPlatesToProcess)
{
cListWells ListWellsToProcess1 = new cListWells(null);
foreach (cWell item in TmpPlate.ListActiveWells)
{
if (item.GetCurrentClassIdx() != -1)
{
if (ListClassSelected[0][item.GetCurrentClassIdx()] == 1)
ListWellsToProcess1.Add(item);
}
}
cExtendedTable NewTable1 = new cExtendedTable(ListWellsToProcess1, cGlobalInfo.CurrentScreening.ListDescriptors.GetDescriptorIndex(cGlobalInfo.CurrentScreening.ListDescriptors.GetActiveDescriptor()));
if ((NewTable1.Count == 0) || (NewTable1[0].Count < 3))
{
if (ProcessModeCurrentPlateOnlyToolStripMenuItem.Checked)
{
MessageBox.Show("Insufficient number of control wells", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
else
continue;
}
if (_StatJarqueBeraItem.Checked)
{
cNormalityJarqueBera JB = new cNormalityJarqueBera();
JB.SetInputData(NewTable1);
JB.Run();
ListZ.Add(JB.GetOutPut()[0][0]);
}
else
{
cStatistics CS = new cStatistics();
CS.UnselectAll();
if (_StatCVItem.Checked)
CS.IsCV = true;
else if (_StatMeanItem.Checked)
CS.IsMean = true;
else if (_StatSumItem.Checked)
CS.IsSum = true;
CS.SetInputData(NewTable1);
CS.Run();
ListZ.Add(CS.GetOutPut()[0][0]);
}
ListPlatesForZFactor.Add(TmpPlate);
}
#endregion
cExtendedTable ET = new cExtendedTable(new cExtendedTable(ListZ));
ET[0].ListTags = new List<object>();
ET[0].ListTags.AddRange(ListPlatesForZFactor);
ET.Name = NameFunction + " - " + cGlobalInfo.CurrentScreening.ListDescriptors[cGlobalInfo.CurrentScreening.ListDescriptors.CurrentSelectedDescriptorIdx].GetName();// +" - " +cGlobalInfo.ListWellClasses[IdxClassNeg].Name + " (" + NewTable1[0].Count + " wells) vs. " +cGlobalInfo.ListWellClasses[IdxClassPos].Name + " (" + NewTable2[0].Count + " wells)";
ET[0].Name = ET.Name;
cViewerGraph1D VG1 = new cViewerGraph1D();
VG1.SetInputData(ET);
VG1.Chart.LabelAxisY = NameFunction;
VG1.Chart.LabelAxisX = "Plate";
VG1.Chart.IsZoomableX = true;
VG1.Chart.IsBar = true;
VG1.Chart.IsBorder = true;
VG1.Chart.IsDisplayValues = true;
VG1.Chart.IsShadow = true;
VG1.Chart.MarkerSize = 4;
VG1.Title = cGlobalInfo.CurrentScreening.ListDescriptors[cGlobalInfo.CurrentScreening.ListDescriptors.CurrentSelectedDescriptorIdx].GetName();
VG1.Run();
cDisplayToWindow CDW = new cDisplayToWindow();
CDW.SetInputData(VG1.GetOutPut());
CDW.Title = NameFunction + " - " + ListPlatesToProcess.Count + " plates";
CDW.Run();
CDW.Display();
}
}
private void ToolStripMenuItem_Statistics(object sender, EventArgs e)
{
cExtendedTable CET = new cExtendedTable();
foreach (DataGridViewColumn item in GridView.SelectedColumns)
CET.Add(Input[item.Index]);
cStatistics S = new cStatistics();
S.SetInputData(CET);
S.Run();
cDisplayExtendedTable CDT = new cDisplayExtendedTable();
CDT.SetInputData(S.GetOutPut());
CDT.Run();
}