Uses MtxExpr, Regress, Math387;
procedure Example;
var Data: Matrix;
  ANOVARes : TANOVA2Result;
  pRows, pCols, pInteract : double;
begin
  Data.SetIt(4,4,false,
            [2.5, 3.2, 4.2, 3.9,
             1.9, 3.5, 3.6, 3.7,
             3.4, 3.5, 3.3, 3.4,
             4.2, 5.0, 3.1, 2.4]);
  ANOVARes := ANOVA2(Data,pRows,pCols,pInteract,1); // no interaction !!
  // ANOVARes =(SS1:0.505; SS2:1.37; SS3:NAN; SS4:6.555; SSTotal:8.43 ;
  // Deg1: 3; Deg2:3 ; Deg3:+NAN; Deg4:9 ; DegTotal:15 ;
  // MS1:0.16833333333 ; MS2:0.45666666667 ; MS3:+NAN ; MS4:0.72833333333 ;
  //FDist1:0.12260277778 ; FDist2:0.33260555556 ; FDist3:+NAN ;
  // FCrit1:3.8625483576 ; FCrit2:3.8625483576 ; FCrit3:+NAN ;
  // pRows: 0.94442876898
  // pCols: 0.80219818603
  // pInteract: 2.0820886898e-307
end;

#include "MtxExpr.hpp"
#include "Regress.hpp"
#include "Math.hpp"
void __fastcall Example()
{
  sMatrix Data;
  Data.SetIt(4,4,false,OPENARRAY(double,(2.5, 3.2, 4.2, 3.9,
         1.9, 3.5, 3.6, 3.7,
         3.4, 3.5, 3.3, 3.4,
         4.2, 5.0, 3.1, 2.4)));
  double prows,pcols,pinteract;
  TANOVA2Result ar = ANOVA2(Data,prows,pcols,pinteract,1,0.05);
  // ar =(SS1:0.505; SS2:1.37; SS3:NAN; SS4:6.555; SSTotal:8.43 ;
  // Deg1: 3; Deg2:3 ; Deg3:+NAN; Deg4:9 ; DegTotal:15 ;
  // MS1:0.16833333333 ; MS2:0.45666666667 ; MS3:+NAN ; MS4:0.72833333333 ;
  //FDist1:0.12260277778 ; FDist2:0.33260555556 ; FDist3:+NAN ;
  // FCrit1:3.8625483576 ; FCrit2:3.8625483576 ; FCrit3:+NAN ;
  // pRows: 0.94442876898
  // pCols: 0.80219818603
  // pInteract: 2.0820886898e-307
}