public override void SetCalibrationModel(IMultivariateCalibrationModel calib)
{
if(calib is PLS2CalibrationModel)
_calib = (PLS2CalibrationModel) calib;
else
throw new ArgumentException("Expecting argument of type PLS2CalibrationModel, but actual type is " + calib.GetType().ToString());
}
public override void Import(
IMultivariateCalibrationModel calibrationSet,
DataTable table)
{
PLS2CalibrationModel calib = (PLS2CalibrationModel)calibrationSet;
// store the x-loads - careful - they are horizontal in the matrix
for (int i = 0; i < calib.XLoads.Rows; i++)
{
Altaxo.Data.DoubleColumn col = new Altaxo.Data.DoubleColumn();
for (int j = 0; j < calib.XLoads.Columns; j++)
{
col[j] = calib.XLoads[i, j];
}
table.DataColumns.Add(col, GetXLoad_ColumnName(i), Altaxo.Data.ColumnKind.V, 0);
}
// now store the y-loads - careful - they are horizontal in the matrix
for (int i = 0; i < calib.YLoads.Rows; i++)
{
Altaxo.Data.DoubleColumn col = new Altaxo.Data.DoubleColumn();
for (int j = 0; j < calib.YLoads.Columns; j++)
{
col[j] = calib.YLoads[i, j];
}
table.DataColumns.Add(col, GetYLoad_ColumnName(i), Altaxo.Data.ColumnKind.V, 1);
}
// now store the weights - careful - they are horizontal in the matrix
for (int i = 0; i < calib.XWeights.Rows; i++)
{
Altaxo.Data.DoubleColumn col = new Altaxo.Data.DoubleColumn();
for (int j = 0; j < calib.XWeights.Columns; j++)
{
col[j] = calib.XWeights[i, j];
}
table.DataColumns.Add(col, GetXWeight_ColumnName(i), Altaxo.Data.ColumnKind.V, 0);
}
// now store the cross product vector - it is a horizontal vector
{
Altaxo.Data.DoubleColumn col = new Altaxo.Data.DoubleColumn();
for (int j = 0; j < calib.CrossProduct.Columns; j++)
{
col[j] = calib.CrossProduct[0, j];
}
table.DataColumns.Add(col, GetCrossProduct_ColumnName(), Altaxo.Data.ColumnKind.V, 3);
}
}
public override void SetCalibrationModel(IMultivariateCalibrationModel calib)
{
if (calib is PLS2CalibrationModel)
{
_calib = (PLS2CalibrationModel)calib;
}
else
{
throw new ArgumentException("Expecting argument of type PLS2CalibrationModel, but actual type is " + calib.GetType().ToString());
}
}
public override void Reset()
{
_calib = new PLS2CalibrationModel();
base.Reset();
}
public override void Reset()
{
_calib = new PLS2CalibrationModel();
base.Reset();
}
/// <summary>
/// Exports a table to a PLS2CalibrationSet
/// </summary>
/// <param name="table">The table where the calibration model is stored.</param>
/// <param name="calibrationSet"></param>
public static void Export(
DataTable table,
out PLS2CalibrationModel calibrationSet)
{
int numberOfX = GetNumberOfX(table);
int numberOfY = GetNumberOfY(table);
int numberOfFactors = GetNumberOfFactors(table);
calibrationSet = new PLS2CalibrationModel();
calibrationSet.NumberOfX = numberOfX;
calibrationSet.NumberOfY = numberOfY;
calibrationSet.NumberOfFactors = numberOfFactors;
MultivariatePreprocessingModel preprocessSet = new MultivariatePreprocessingModel();
MultivariateContentMemento plsMemo = table.GetTableProperty("Content") as MultivariateContentMemento;
if (plsMemo != null)
{
preprocessSet.PreprocessOptions = plsMemo.SpectralPreprocessing;
}
calibrationSet.SetPreprocessingModel(preprocessSet);
Altaxo.Collections.AscendingIntegerCollection sel = new Altaxo.Collections.AscendingIntegerCollection();
Altaxo.Data.DataColumn col;
col = table[GetXOfX_ColumnName()];
if (col == null || !(col is INumericColumn))
{
NotFound(GetXOfX_ColumnName());
}
preprocessSet.XOfX = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector((INumericColumn)col, numberOfX);
col = table[GetXMean_ColumnName()];
if (col == null)
{
NotFound(GetXMean_ColumnName());
}
preprocessSet.XMean = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector(col, numberOfX);
col = table[GetXScale_ColumnName()];
if (col == null)
{
NotFound(GetXScale_ColumnName());
}
preprocessSet.XScale = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector(col, numberOfX);
sel.Clear();
col = table[GetYMean_ColumnName()];
if (col == null)
{
NotFound(GetYMean_ColumnName());
}
sel.Add(table.DataColumns.GetColumnNumber(col));
preprocessSet.YMean = DataColumnWrapper.ToROVector(col, numberOfY);
sel.Clear();
col = table[GetYScale_ColumnName()];
if (col == null)
{
NotFound(GetYScale_ColumnName());
}
sel.Add(table.DataColumns.GetColumnNumber(col));
preprocessSet.YScale = DataColumnWrapper.ToROVector(col, numberOfY);
sel.Clear();
for (int i = 0; i < numberOfFactors; i++)
{
string colname = GetXWeight_ColumnName(i);
col = table[colname];
if (col == null)
{
NotFound(colname);
}
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.XWeights = DataTableWrapper.ToRORowMatrix(table.DataColumns, sel, numberOfX);
sel.Clear();
for (int i = 0; i < numberOfFactors; i++)
{
string colname = GetXLoad_ColumnName(i);
col = table[colname];
if (col == null)
{
NotFound(colname);
}
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.XLoads = DataTableWrapper.ToRORowMatrix(table.DataColumns, sel, numberOfX);
sel.Clear();
for (int i = 0; i < numberOfFactors; i++)
{
string colname = GetYLoad_ColumnName(i);
col = table[colname];
if (col == null)
{
NotFound(colname);
}
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.YLoads = DataTableWrapper.ToRORowMatrix(table.DataColumns, sel, numberOfY);
sel.Clear();
col = table[GetCrossProduct_ColumnName()];
if (col == null)
{
NotFound(GetCrossProduct_ColumnName());
}
sel.Add(table.DataColumns.GetColumnNumber(col));
calibrationSet.CrossProduct = DataTableWrapper.ToRORowMatrix(table.DataColumns, sel, numberOfFactors);
}
/// <summary>
/// Exports a table to a PLS2CalibrationSet
/// </summary>
/// <param name="table">The table where the calibration model is stored.</param>
/// <param name="calibrationSet"></param>
public static void Export(
DataTable table,
out PLS2CalibrationModel calibrationSet)
{
int numberOfX = GetNumberOfX(table);
int numberOfY = GetNumberOfY(table);
int numberOfFactors = GetNumberOfFactors(table);
calibrationSet = new PLS2CalibrationModel();
calibrationSet.NumberOfX = numberOfX;
calibrationSet.NumberOfY = numberOfY;
calibrationSet.NumberOfFactors = numberOfFactors;
MultivariatePreprocessingModel preprocessSet = new MultivariatePreprocessingModel();
MultivariateContentMemento plsMemo = table.GetTableProperty("Content") as MultivariateContentMemento;
if(plsMemo!=null)
preprocessSet.PreprocessOptions = plsMemo.SpectralPreprocessing;
calibrationSet.SetPreprocessingModel(preprocessSet);
Altaxo.Collections.AscendingIntegerCollection sel = new Altaxo.Collections.AscendingIntegerCollection();
Altaxo.Data.DataColumn col;
col = table[GetXOfX_ColumnName()];
if(col==null || !(col is INumericColumn)) NotFound(GetXOfX_ColumnName());
preprocessSet.XOfX = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector((INumericColumn)col,numberOfX);
col = table[GetXMean_ColumnName()];
if(col==null) NotFound(GetXMean_ColumnName());
preprocessSet.XMean = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector(col,numberOfX);
col = table[GetXScale_ColumnName()];
if(col==null) NotFound(GetXScale_ColumnName());
preprocessSet.XScale = Altaxo.Calc.LinearAlgebra.DataColumnWrapper.ToROVector(col,numberOfX);
sel.Clear();
col = table[GetYMean_ColumnName()];
if(col==null) NotFound(GetYMean_ColumnName());
sel.Add(table.DataColumns.GetColumnNumber(col));
preprocessSet.YMean = DataColumnWrapper.ToROVector(col,numberOfY);
sel.Clear();
col = table[GetYScale_ColumnName()];
if(col==null) NotFound(GetYScale_ColumnName());
sel.Add(table.DataColumns.GetColumnNumber(col));
preprocessSet.YScale = DataColumnWrapper.ToROVector(col,numberOfY);
sel.Clear();
for(int i=0;i<numberOfFactors;i++)
{
string colname = GetXWeight_ColumnName(i);
col = table[colname];
if(col==null) NotFound(colname);
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.XWeights = DataTableWrapper.ToRORowMatrix(table.DataColumns,sel,numberOfX);
sel.Clear();
for(int i=0;i<numberOfFactors;i++)
{
string colname = GetXLoad_ColumnName(i);
col = table[colname];
if(col==null) NotFound(colname);
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.XLoads = DataTableWrapper.ToRORowMatrix(table.DataColumns,sel,numberOfX);
sel.Clear();
for(int i=0;i<numberOfFactors;i++)
{
string colname = GetYLoad_ColumnName(i);
col = table[colname];
if(col==null) NotFound(colname);
sel.Add(table.DataColumns.GetColumnNumber(col));
}
calibrationSet.YLoads = DataTableWrapper.ToRORowMatrix(table.DataColumns,sel,numberOfY);
sel.Clear();
col = table[GetCrossProduct_ColumnName()];
if(col==null) NotFound(GetCrossProduct_ColumnName());
sel.Add(table.DataColumns.GetColumnNumber(col));
calibrationSet.CrossProduct = DataTableWrapper.ToRORowMatrix(table.DataColumns,sel,numberOfFactors);
}
